|
Size: 3503
Comment:
|
← Revision 41 as of 2018-01-04 15:07:58 ⇥
Size: 5133
Comment:
|
| Deletions are marked like this. | Additions are marked like this. |
| Line 1: | Line 1: |
| [wiki:Self:FreeSurferWorkFlows top] | [[FreeSurferWorkFlows|top]] |
| Line 11: | Line 11: |
| FreeSurfer now runs automated labeling of the brain volume and this is included in the stable v3.0 release, [wiki:Self:ReconAllDevTable during the -autorecon2 stage]. However, if you processed your anatomical data using previous versions and you wish to obtain the automated labels, you can just run the subcortical segmentation separately. To obtain automatically segmented volumes for the first time (assuming the [wiki:Self:ReconAllDevTable -autorecon1 stage] has completed), run: {{{ recon-all -subcortseg -subjid <subject name> }}} |
FreeSurfer runs automated labeling of the brain volume [[ReconAllDevTable|during the -autorecon2 stage]]. |
| Line 21: | Line 15: |
| It is important to note that the stable v3.0 release of FreeSurfer by default uses the automatically segmented brain volume (ASEG - in mri/aseg.mgz) to segment the white matter volume (WM - mri/wm). In doing this the edits are different than the manual edits you may have already done in older versions of freesurfer (i.e., the ventricles are filled in completely) this creates a different folding pattern and uses different atlases. You should be sure that all of your subjects are edited in the same manner, using the same atlases (i.e., don't mix freesurfer versions). If you wish to process one subject further (i.e., to get final surfaces) and wish to maintain your edits you can use the ''' -noedit_wm_with_aseg''' flag along with the '''-old-atlas''' flag. This will preserve your manual edits and utilize the old atlases. This is not ideal, as you should process all your subjects fully with the same version of FreeSurfer. | |
| Line 23: | Line 16: |
| To view the segmentation, use this command: | To view just the segmentation, use this command: |
| Line 26: | Line 19: |
| tkmedit <subject name> norm.mgz -segmentation mri/aseg.mgz $FREESURFER_HOME/FreeSurferColorLUT.txt | tkmedit <subject name> norm.mgz -aseg |
| Line 31: | Line 24: |
| If the voxels are incorrectly labeled (this is more likely to happen along the medial aspect of the brain), then you can re-label them yourself. Refer to the TkMeditGuide/TkMeditWorkingWithData/TkMeditSelectionsLabels page for detailed information on this procedure. | If the voxels are incorrectly labeled, then you can re-label them yourself although we suggest contacting the freesurfer mailing list first to see if there is an automatic procedure that can be used to re-label them.. Refer to the TkMeditGuide/TkMeditWorkingWithData/TkMeditSelectionsLabels page for detailed information on manually editing the aseg. |
| Line 37: | Line 30: |
| * [https://surfer.nmr.mgh.harvard.edu/ftp/articles/fischl02-labeling.pdf Whole Brain Segmentation: Automated Labeling of Neuroanatomical Structures in the Human Brain], Fischl et al., (2002). Neuron, 33:341-355. | * [[https://surfer.nmr.mgh.harvard.edu/ftp/articles/fischl02-labeling.pdf|Whole Brain Segmentation: Automated Labeling of Neuroanatomical Structures in the Human Brain]], Fischl et al., (2002). Neuron, 33:341-355. |
| Line 39: | Line 32: |
| The atlas included with the Freesurfer distribution is found in the 'average' directory, and is called 'RB_all_2006-02-15.gca'. It is possible to construct your own atlas. This is described next. | The atlas included with the Freesurfer distribution is found in the 'average' directory, and is called 'RB_all_2008-03-26.gca'. It is possible to construct your own atlas. This is described next. |
| Line 46: | Line 39: |
| * Run ["mri_ca_train"] to create the atlas. | * Run [[mri_ca_train]] to create the atlas. |
| Line 49: | Line 42: |
See also AtlasSubjects For accuracy evaluations see also SubcorticalSegmentationAccuracy = Frequently Asked Questions = == Are there any files which contain the left and right hemisphere segmentation of aseg.mgz separately? == No, but you could accomplish this in MATLAB by setting all the left (or right) segments to zero. For example, you can use the following code as a template and set each of the labels equal to zero. {{{ a = MRIread('aseg.mgz'); ind = find(a.vol == 17); % Left hippo a.vol(ind) = 0; % Set all left hippo to 0 }}} and so on for the other segments = References = The original CMA segmentation scheme used for subcortical segmentation training is defined by (Filipek, et al, Cerebral Cortex, 1994) ([[https://www.ncbi.nlm.nih.gov/pubmed/7950308|paper]]). However, with better imaging resolution, anatomists Verne Caviness and Nikos Mkris developed the separation of the "thalamus proper" from a "ventral diencephalon" region that subtends many of the smaller nuclei and structures in the area inferior to the thalamus, such as hypothalamus, red nuclei, later and medial geniculate, ect. The first definition of the ventral diencephalon method is provided in Am J Med Genet. 1997 Sep 19;74(5):507-14. ([[https://www.ncbi.nlm.nih.gov/pubmed/9342202|paper]]). As defined by Makris in Makris, et al. Biol Psychiatry. 2008 Aug 1;64(3):192-202. ([[https://www.ncbi.nlm.nih.gov/pmc/articles/PMC2572710/|paper]] "...ventral diencephalon (49), which according to our morphometric definition contains the hypothalamus, basal forebrain, and sublenticular extended amygdala (SLEA), as well as a large portion of ventral tegmentum (which is included in our ventral diencephalon region by convention although part of midbrain)". Also, from BIOL PSYCHIATRY L.J. Seidman et al 1999;46:941–954 ([[https://www.ncbi.nlm.nih.gov/pubmed/10509177|paper]]) "We used the same segmentation procedure of the thalamus developed in our preliminary MRI report (Seidman et al 1997). The central gray nuclei were divided into superior and inferior segments, the thalamus proper and the “ventral diencephalon,” respectively. The thalamus proper (including all thalamic nuclei except the lateral and medial geniculate bodies) was defined on all coronal slices on which it was present, from the first slice containing the anterior nucleus through the last slice containing the pulvinar nucleus. The structure was bounded medially by the third ventricle, and laterally by the internal capsule. The superior border was the body of the lateral ventricle, and the inferior border, the hypothalamic fissure. This structure definition differed from the central gray nuclei of Filipek and co-workers (1994) only in its inferior boundary." The 5th ventricle, referred to as the "cavum septum pallucidum," is a bit of an enigma that tends to only occur on occasion. It's not reported in any manuscripts, but it was confirmed by neuroradiologist, P. Ellen Grant. |
Subcortical Segmentation
In automatic subcortical segmentation, each voxel in the normalized brain volume is assigned one of about 40 labels, including:
- Cerebral White Matter, Cerebral Cortex, Lateral Ventricle, Inferior Lateral Ventricle, Cerebellum White Matter, Cerebellum Cortex, Thalamus, Caudate, Putamen, Pallidum, Hippocampus, Amygdala, Lesion, Accumbens area, Vessel, Third Ventricle, Fourth Ventricle, Brain Stem, Cerebrospinal Fluid
FreeSurfer runs automated labeling of the brain volume during the -autorecon2 stage.
The automatic subcortical segmentation can take many (11+) hours to complete.
To view just the segmentation, use this command:
tkmedit <subject name> norm.mgz -aseg
Float your cursor over any voxel and the label assigned to it will be displayed in the TkMeditTools window.
If the voxels are incorrectly labeled, then you can re-label them yourself although we suggest contacting the freesurfer mailing list first to see if there is an automatic procedure that can be used to re-label them.. Refer to the TkMeditGuide/TkMeditWorkingWithData/TkMeditSelectionsLabels page for detailed information on manually editing the aseg.
Aseg Atlas
Automatic subcortical segmentation of a brain volume is based upon the existence of an atlas containing probablistic information on the location of structures. This is decribed here:
Whole Brain Segmentation: Automated Labeling of Neuroanatomical Structures in the Human Brain, Fischl et al., (2002). Neuron, 33:341-355.
The atlas included with the Freesurfer distribution is found in the 'average' directory, and is called 'RB_all_2008-03-26.gca'. It is possible to construct your own atlas. This is described next.
Constructing an Aseg Atlas
Basic steps:
- Using tkmedit, label each volume in the brain. Repeat for all subjects to be included in the atlas.
Run mri_ca_train to create the atlas.
See also rebuild_gca_atlas.csh script in $FREESURFER_HOME/bin.
See also AtlasSubjects
For accuracy evaluations see also SubcorticalSegmentationAccuracy
Frequently Asked Questions
Are there any files which contain the left and right hemisphere segmentation of aseg.mgz separately?
No, but you could accomplish this in MATLAB by setting all the left (or right) segments to zero. For example, you can use the following code as a template and set each of the labels equal to zero.
a = MRIread('aseg.mgz');
ind = find(a.vol == 17); % Left hippo
a.vol(ind) = 0; % Set all left hippo to 0and so on for the other segments
References
The original CMA segmentation scheme used for subcortical segmentation training is defined by (Filipek, et al, Cerebral Cortex, 1994) (paper). However, with better imaging resolution, anatomists Verne Caviness and Nikos Mkris developed the separation of the "thalamus proper" from a "ventral diencephalon" region that subtends many of the smaller nuclei and structures in the area inferior to the thalamus, such as hypothalamus, red nuclei, later and medial geniculate, ect. The first definition of the ventral diencephalon method is provided in Am J Med Genet. 1997 Sep 19;74(5):507-14. (paper). As defined by Makris in Makris, et al. Biol Psychiatry. 2008 Aug 1;64(3):192-202. (paper "...ventral diencephalon (49), which according to our morphometric definition contains the hypothalamus, basal forebrain, and sublenticular extended amygdala (SLEA), as well as a large portion of ventral tegmentum (which is included in our ventral diencephalon region by convention although part of midbrain)".
Also, from BIOL PSYCHIATRY L.J. Seidman et al 1999;46:941–954 (paper) "We used the same segmentation procedure of the thalamus developed in our preliminary MRI report (Seidman et al 1997). The central gray nuclei were divided into superior and inferior segments, the thalamus proper and the “ventral diencephalon,” respectively. The thalamus proper (including all thalamic nuclei except the lateral and medial geniculate bodies) was defined on all coronal slices on which it was present, from the first slice containing the anterior nucleus through the last slice containing the pulvinar nucleus. The structure was bounded medially by the third ventricle, and laterally by the internal capsule. The superior border was the body of the lateral ventricle, and the inferior border, the hypothalamic fissure. This structure definition differed from the central gray nuclei of Filipek and co-workers (1994) only in its inferior boundary."
The 5th ventricle, referred to as the "cavum septum pallucidum," is a bit of an enigma that tends to only occur on occasion. It's not reported in any manuscripts, but it was confirmed by neuroradiologist, P. Ellen Grant.
