Dear Friends,
I need to 3D render the AAL parcellation volume, with labels (in annotation file) for each ROI.
My AAL volume is on MNI152 standard space, with 2mm 91x109x91, as well my results volumes.
I used the following command to generate my surface files from my AAL volume.
mri_vol2surf --mov aal2std.nii.gz --reg cvs_avg35_inMNI152/mri/transforms/reg.mni152.2mm.dat --projdist-max 0 1 0.1 --interp nearest --hemi lh --out lh.aal2std-1mm_surf.mgz —reshape
mri_vol2surf --mov aal2std.nii.gz --reg cvs_avg35_inMNI152/mri/transforms/reg.mni152.2mm.dat --projdist-max 0 1 0.1 --interp nearest --hemi rh --out rh.aal2std-1mm_surf.mgz --reshape
Then I tried to create my labels using the following command:
mri_cor2label --i lh.aal2std-1mm_surf.mgz --surf cvs_avg35_inMNI152 lh --id 75 --l lh.Pallidum_L
I run this command for each one of the 90 ROIs on AAL volume.
My problem is that when I execute the last command specifically for Pallidum, both Left and Right (IDs 75 and 76), I found this error:
$Id: mri_cor2label.c,v 1.12 2011/03/02 00:04:14 nicks Exp $ Loading mri lh.aal2std-1mm_surf.mgz Loading cvs_avg35_inMNI152/surf/lh.white Scanning the volume Found 0 label voxels ERROR: found no voxels matching id 75
It seems that I loose some voxels when I convert from volume voxels to surface vertexes.
Do you know what could be happening?
I appreciate any help you could give to us.
Thanks very much,
best regards,
-- João V. Dornas Phd student in Neuroscience Cognitive Biology Group, Institute of Biology Otto von Guericke University, Magdeburg Germany