How does it work if you just use vanilla recon-all v6 beta with the T2 refinement flag and the hires option? We haven¹t yet updated the HCP scripts for FreeSurfer v6, but this is on our to do list.
Peace,
Matt.
From: freesurfer-bounces@nmr.mgh.harvard.edu on behalf of Antonin Skoch ansk@ikem.cz Reply-To: Freesurfer support list freesurfer@nmr.mgh.harvard.edu Date: Sunday, November 27, 2016 at 7:03 PM To: freesurfer@nmr.mgh.harvard.edu Subject: [Freesurfer] freesufer v6 beta - T2 pial refinement performance and nsigma_below setting
Dear FreeSurfer experts,
I have processed a group of subjects by the freesurfer v6 beta version (with recent upgrade of mris_make_surfaces from 2016/11/21).
It seems to me that in contrast to v5.3 version (with the same nsigma_above 2, nsigma_below 3) the T2 pial refinement cuts out too much gray matter in regions where the gray matter has rather low intensity in T2 images. I have observed this behavior quite often in my subjects.
I have uploaded the example subject cut_gm_T2pial.tar.gz to the server. The problematic area has RAS location 3.9, -70, 0. Also see the screenshot.
My command line was:
mris_make_surfaces -nsigma_above 2 -nsigma_below 3 -aseg aseg.hires -filled filled.hires -wm wm.hires -mgz -sdir $my_subject_dir -orig white.deformed -nowhite -orig_white white.deformed -orig_pial pial.preT2.pass2 -T2 my_subject/mri/T2w_hires.norm -T1 T1w_hires.norm -output .postT2.pass2 my_subject rh
and for FS 5.3
mris_make_surfaces -nsigma_above 2 -nsigma_below 3 -aseg aseg.hires -filled filled.hires -wm wm.hires -mgz -sdir $my_subject_dir -orig white.deformed -nowhite -orig_white white.deformed -orig_pial pial.preT2.pass2 -T2 my_subject/mri/T2w_hires.norm -T1 T1w_hires.norm -output .postT2.pass2.FS_5.3 my_subject rh
Could you please comment on how to proceed?
It seems that in contrast to v5.3 it is necessary to increase nsigma_below value since the nsigma values seem to be more "sensitive" in v6beta version. I tried to increase nsigma_below (nsigma_below for T2 in recon-all is 5) but it seemed to me that it leads to "leak" of the pial surface out of the gray matter in other areas, mainly in regions adjacent to cerebellum. See the example leaking the pial surface to the cerebellum in left hemisphere with nsigma_below 5 in the second screenshot in attachment ( I have uploaded this subject as file cut_gm_T2pial2.tar.gz, his cut-out gray matter RAS location is -33.7, -18, 45, command line is similar as above, but for left hemisphere). Yellow surface is with nsigma_below 5, blue is nsigma_below 3. Btw, the -aseg option in mris_make_surfaces does not prevent to leak surfaces into the cerebellum?
Could also be the reason of these problems, that I used HCP scripts which do not do the T2 normalization with anatomy-informed mri_normalize (using aseg and white surface), but only global normalization and bias field correction?
Regards,
Antonin Skoch _______________________________________________ Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer The information in this e-mail is intended only for the person to whom it is addressed. If you believe this e-mail was sent to you in error and the e-mail contains patient information, please contact the Partners Compliance HelpLine at http://www.partners.org/complianceline . If the e-mail was sent to you in error but does not contain patient information, please contact the sender and properly dispose of the e-mail.