There's a lot of stuff going on there, so it's hard to say withhout more
detail. But it sounds like you might not be taking into account that the
voxel coords in matlab are 1-based.
On 6/26/2020 12:04 PM, Franziska Klein wrote:
External Email - Use Caution
Hello,
I am using the isosurface() function in order to create a mesh based
on the T1.mgz and on the ribbon.mgz. I need those meshes for
the coregistration of some channels.
After coregistration I do project the channels onto the cortex and
extract the corresponding voxel coordinates. Until here, everything
works fine.
Next, I am using MRIwrite in order to save the T1.mgz as *.nii.mgz
which I am using to generate a deformation field with SPM12 to
transform the voxel coordinates to normalized MNI coordinates.
However, when I put the voxel coordinates (not normalized yet)
into SPM12 while displaying my generated *.nii image it happens that
the coordinates on the right hemisphere (MNI space) are slightly
outside of the cortex which is of course not the case on my cortex
mesh after projection.
Now I am wondering whether there is a mismatch between the mesh and
the *.nii image (one voxel maybe)? And whether that is introduced by
MRIwrite? If so, do I have to shift all points of the meshes by +1 in
order to get correct results the way I want to? If not, any ideas on
where that mismatch could come from?
Thank you and kind regards,
Franziska
_______________________________________________
Freesurfer mailing list
Freesurfer@nmr.mgh.harvard.edu
https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer
_______________________________________________
Freesurfer mailing list
Freesurfer@nmr.mgh.harvard.edu
https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer
