Dear Freesurfer experts,
I'm running into some problems when running the long_mris_slopes command (longitudinal 2 stage model) to get the rate of change for thickness and area for each hemisphere.
This is the command I ran:
for i in thickness area ; do for j in rh ; do
long_mris_slopes --qdec $SUBJECTS_DIR/qdec/39kids-covariates-MT2-1.dat --meas ${i} --hemi ${j} --do-rate --do-stack --do-label --time years --qcache 81kids_template --sd $SUBJECTS_DIR
done done
For the left hemisphere everything works fine but for the right hemisphere it stops working at my second subject and I get the following error:
ERROR -6 : mris_label_calc intersect did not work?
This is the terminal output for said subject: Subject-Template: BJ4K.base
INFO: 2 TPs in BJ4K.base , mean age: 0.041095
=============================================================================== SUBJECT BJ4K.base Intersecting Within-Subject Cortex Label
cp /data/p_npsy001_neurotrain/kids_data/39kids-sem-con-3goodscans/semcon-3goodscans2-1/BJ4K_1.long.BJ4K.base/label/rh.cortex.label /data/p_npsy001_neurotrain/kids_data/39kids-sem-con-3goodscans/semcon-3goodscans2-1/BJ4K.base/label/rh.long.cortex.label
mris_label_calc intersect /data/p_npsy001_neurotrain/kids_data/39kids-sem-con-3goodscans/semcon-3goodscans2-1/BJ4K_2.long.BJ4K.base/label/rh.cortex.label /data/p_npsy001_neurotrain/kids_data/39kids-sem-con-3goodscans/semcon-3goodscans2-1/BJ4K.base/label/rh.long.cortex.label /data/p_npsy001_neurotrain/kids_data/39kids-sem-con-3goodscans/semcon-3goodscans2-1/BJ4K.base/label/rh.long.cortex.label
mris_label_calc: malloc.c:2392: sysmalloc: Assertion `(old_top == initial_top (av) && old_size == 0) || ((unsigned long) (old_size) >= MINSIZE && prev_inuse (old_top) && ((unsigned long) old_end & (pagesize - 1)) == 0)' failed. ERROR -6 : mris_label_calc intersect did not work?
I've googled a bit and suggestions were a memory problem but I don't think that's the case if it is hemisphere specific and the first subject in the list works fine...
I've also attached the input file.
What could be the problem and how can I fix it? Cheers Clara