Dear Surfers,
I (successfully) resampled a functional data on the surface map using
mri_vol2surf:
> mri_vol2surf --mov func.nii --reg register.dat --projfrac 0.5 --interp nearest --hemi lh --o lh.func.mgh
Then, I normalized this image to a standard space (fsaverage) using
mri_surf2surf:
> mri_surf2surf --srcsubject sub1 --trgsubject fsaverage --sval sub1/lh.func.mgh --tval sub1/lh.func.norm.mgh --fwhm 5 --cortex --hemi lh
The dimension of the output file, lh.func.norm.mgh, was 163842 x 1 x 1
(identified from mri_info). This file can be opened from Matlab, and I
can successfully get access to 163842 vertices. What I want now is to
take the aparc+aseg information from the average brain (fsaverage) and
be able identify which vertex falls under which ROI label. I decided
to resample aparc+aseg.mgz in fsaverage on the surface map, hoping to
get label information at the dimension of 163842 x 1 x 1. So I ran:
> mri_vol2surf --mov fsaverage/mri/aparc+aseg.mgz --regheader fsaverage --projfrac 0.5 --interp nearest --hemi lh --o fsaverage/surf/lh.roi.mgh
> mri_vol2surf --mov fsaverage/mri/aparc+aseg.mgz --regheader fsaverage --projfrac 0.5 --interp nearest --hemi rh --o fsaverage/surf/rh.roi.mgh
The resulting mgh files, however, had a dimension of 23406 x 1 x 7.
1. As I said, I wanted to get label information (roi number from the
FreeSurferColorLUT) for each vertex of the surface image in the
fsaverage space. Is this the right approach? And did I do this
correctly (I'm concerned about many options above)?
2. If so, why is the resulting dimension different from the other one?
I assume that it's probably a reshape issue, but how I can reshape the
dimension to match 163842 x 1 x 1?
Thanks in advance.
Joon
--
Joonkoo Park
Ph.D. Candidate, Cognition & Cognitive Neuroscience
Department of Psychology, University of Michigan