Hello all Freesurfer users!
I am planning to do some cortical thickness analysis in children affected by psychiatric disorders, and after a few months of reflection, I still have some trouble to obtain a correct segmentation in the parietal region with my images.
Actually, I am especially interested in the parietal lobe where we already observed a decrease in the gray matter volume in this sample of patients, but I have some trouble with the segmentation of the gray-white matter in this region. In fact, my orig and smoothwm surfaces look pretty good, except in about 10 slices in the parietal lobe where parietal cortex is so bright that it is classified as white matter. I first thought that it was either due to the quality of the scan (1.5 Tesla, slices 1.5mm) or related to a bias in the magnetic field. So first I tried to apply different bias correctors, but it has no effect on the segmentation. Finally, I found that forcing the high gray and white low values during the segmentation (in the Wm Filter Expert Preferences) was really interesting to obtain a better orig surface on my images. What I did is I applied a Gaussian Fitting algorithm of segmentation on the normalized and skull stripped image (brain), and took the limit values of the partial volume intensity for the high gray and low white value, so that the segmentation included the tissue classified as partial volume (about 97 for white low and 104 for high gray, depending on the subjects). This had a quite benefic effect on the parietal cortex, but as you can probably imagine, it excluded the subcortical structure of my wm images, and increase drastically the time spent for manual editing. Further, most of the time I have to delineate the insular cortex manually, so I just displaced the segmentation problem from the upper to the medial structures.
I am quite perturbed by this failure to obtain a correct cortical ribbon in the parietal region with my images, any advice or idea will be greatly appreciated…
Are my images not enough good for what I want to do? Did someone already notice that parietal cortical thickness was sometimes difficult to assess? (I checked the bert’s one, which also show a really thin parietal cortex?)
How can I solve this problem? Is there any way to apply differentially the rules for the segmentation according to the region in Freesurfer? For example, based on a reference atlas that specifies that partial volume have to be included in the parietal region, and excluded in the subcortical region? Or referring to the talairach boxes, that means specify the high gray and low white value for all the regions that correspond to the coordinates of the parietal talairach boxes?
Thank you very much for your help,
Best regards,
Marie Schaer Neuroimaging Laboratory Child Psychiatry Department Geneva Faculty of Medicine Switzerland
PS: I attached the images with my first message, but it was refused because too big for the mailing list, but I can easily provide you the images again.
freesurfer@nmr.mgh.harvard.edu