External Email - Use Caution
Hi all
Regarding hippocampal subfield segmentation in FS7- how can I move the segmentation result to native T2w space? I plan to compare the result with a manual segmentation which was done in this space, using scripts which require nifti inputs.
I have run FS7 specifying with the -cm flag that this takes place in my T1 (MP2RAGE) 0.65mm isotropic resolution: recon-all -all -s <subject> -i T1w.nii -T2 T2w.nii -T2pial -gcut -cm
Then I ran the Hc subfield segmentation, specifying that my T2w volume (3D-SPACE, 0.55mm isotropic) should be used for this: segmentHA_T2.sh <subject> T2w.nii 3DT2 0
The output that I want to move back to native 3DT2 space is the labelled volume lh.hippoAmygLabels-3DT2.v21.mgz (or the nifti of this which of course I can create using mri_convert). However this volume is 0.333mm isotropic and has a limited field of view - from the documentation: "Note that [lr]h.hippoAmygLabels-T1.v21.mgz and [lr]h.hippoAmygLabels-T1.v21.[hierarchy].mgz cover only a patch around the hippocampus, at a higher resolution than the input image."
The most relevant guidance I can find in the documentation is "[lr]h.hippoAmygLabels-T1.v21.FSvoxelSpace.mgz live in the same voxel space as the other Freesurfer volumes (e.g., orig.mgz, nu.mgz, aseg.mgz), so you can use it directly to produce masks for further analyses, but its resolution is lower (1 mm vs 0.333 mm)." In my case FSvoxelSpace is at 0.65mm resolution, but this is still lower than my T2w volume (0.55mm) and so would involve some downsampling (0.333mm --> 0.65mm --> 0.55mm) which I would like to avoid.
I have tried various linear registration approaches but I always end up with a result which, although in 0.55mm space, is clearly not aligned with the hippocampus when viewed together with my T2w.nii in ITK-SNAP.
Any suggestions would be much appreciated!
Many thanks
Mike