External Email - Use Caution
Dear experts,
I am running mri_aparc2aseg using two annotations as input (one for each hemisphere), each with 159 labels, and wish to end up with a single output segmentation of 318 labels. However, when I run the below, this is not the case:
mri_aparc2aseg \ --s subject1 \ --annot 500_sym.aparc \ --wmparc-dmax 2 \ --labelwm \ --hypo-as-wm \ --o subject1/mri/parcellation.nii.gz
In fact, the output has 680 labels rather than 318 labels. A few questions.
- How can I see what exactly the 680 label names are (how to extract from the .nii.gz output file)? Currently, I have just quantified it to be 680 by doing this in matlab: gunzip('parcellation.nii.gz') a = MRIread('parcellation.nii'); u = unique(a.vol(:)); length(u) % total number of labels
- I believe the extra labels are from the Freesurfer LUT. How can I make sure they are not added when running mri_aparc2aseg? (such that only my annotation labels are included in the output)
Thank you very much,
Ellen
***** Ellen Ji, PhD Postdoctoral Research Fellow Psychiatric University Hospital University of Zürich ellen.ji@bli.uzh.ch homanlab.github.io/ellen/