|
Size: 2416
Comment: renamed from GEMS
|
Size: 3919
Comment:
|
| Deletions are marked like this. | Additions are marked like this. |
| Line 1: | Line 1: |
| Segmentation of hippocampal subfields | = Segmentation of hippocampal subfields = ''Author: Koen Van Leemput'' == DEPRECATED - This functionality is found only in v5.2 and v5.3. == '''See [[HippocampalSubfields]] for the functionality found in v6.0 and downloadable.''' See also: [[GEMS]] |
| Line 5: | Line 13: |
| * Automated Segmentation of Hippocampal Subfields from Ultra-High Resolution In Vivo MRI. K. Van Leemput, A. Bakkour, T. Benner, G. Wiggins, L.L. Wald, J. Augustinack, B.C. Dickerson, P. Golland, B. Fischl. Hippocampus, vol. 19, no. 6, pp. 549-557, June 2009 | * [[https://surfer.nmr.mgh.harvard.edu/pub/articles/VanLeemputHippocampus2009.pdf|Automated Segmentation of Hippocampal Subfields from Ultra-High Resolution In Vivo MRI]]. K. Van Leemput, A. Bakkour, T. Benner, G. Wiggins, L.L. Wald, J. Augustinack, B.C. Dickerson, P. Golland, B. Fischl. Hippocampus, vol. 19, no. 6, pp. 549-557, June 2009 |
| Line 7: | Line 15: |
| To generate the hippocampal subfield segmentations as part of the normal FreeSurfer processing stream, simply add the "-hippo-subfields" flag to the recon-all script: | === recon-all stream === |
| Line 9: | Line 17: |
| To generate the hippocampal subfield segmentations as part of the normal !FreeSurfer processing stream, simply add the "-hippo-subfields" flag to the recon-all script: {{{ |
|
| Line 10: | Line 20: |
| }}} It is also possible to only run the hippocampal subfield segmentation for a subject that has already undergone the standard volumetric !FreeSurfer pipeline (in particular, the program needs the !FreeSurfer-generated files "aseg.mgz", "nu.mgz" and "talairach.xfm"): {{{ recon-all -s bert -hippo-subfields }}} The output consists of a collection of images $SUBJECTS_DIR/bert/mri/posterior_* that indicate each voxel's posterior probability of belonging to different subregions. The estimated volumes of these subregions can be obtained using the "kvlQuantifyHippocampalSubfieldSegmentations.sh" script. Doing {{{ kvlQuantifyHippocampalSubfieldSegmentations.sh }}} will create the text files "nonPartialVolumeStatsLeft.txt" and "nonPartialVolumeStatsRight.txt" in the directory $SUBJECTS_DIR, containing tables that summarize the volumes of the hippocampal subfield segmentations of the left and right side, respectively, in all of the subjects in $SUBJECTS_DIR. The reported volumes are given in number of voxels, where each voxel is of the size "0.5x0.5x0.5mm^3". It is expected that the "-hippo-subfields" flag has been used for each subject in $SUBJECTS_DIR; if this is not the case the "kvlQuantifyHippocampalSubfieldSegmentations.sh" script will fail and no (or incomplete) output files will be generated. |
|
| Line 11: | Line 31: |
| It is also possible to only run the hippocampal subfield segmentation for a subject that has already undergone the standard volumetric FreeSurfer pipeline (in particular, the program needs the FreeSurfer-generated files "aseg.mgz", "nu.mgz" and "talairach.xfm"): | The estimated hippocampal subfield volumes can also be computed for one specific subject only: For the right hippocampus, the command is |
| Line 13: | Line 33: |
| recon-all -s bert -hippo-subfields The output consists of a collection of images $SUBJECTS_DIR/bert/mri/posterior_* that indicate each voxel's posterior probability of belonging to different subregions. The estimated volumes of these subregions can be obtained using the "kvlQuantifyHippocampalSubfieldSegmentations.sh" script. Doing cd $SUBJECTS_DIR kvlQuantifyHippocampalSubfieldSegmentations.sh will create the text files "nonPartialVolumeStatsLeft.txt" and "nonPartialVolumeStatsRight.txt" in the directory $SUBJECTS_DIR, containing tables that summarize the volumes of the hippocampal subfield segmentations of the left and right side, respectively, in all of the subjects in $SUBJECTS_DIR. The reported volumes are given in number of voxels, where each voxel is of the size "0.5x0.5x0.5mm^3". It is expected that the "-hippo-subfields" flag has been used for each subject in $SUBJECTS_DIR. |
{{{ cd $SUBJECTS_DIR/bert/mri kvlQuantifyPosteriorProbabilityImages $FREESURFER_HOME/data/GEMS/compressionLookupTable.txt \ posterior_right_* posterior_Right-Hippocampus.mgz }}} and for the left {{{ cd $SUBJECTS_DIR/bert/mri kvlQuantifyPosteriorProbabilityImages $FREESURFER_HOME/data/GEMS/compressionLookupTable_left.txt \ posterior_left_* posterior_Left-Hippocampus.mgz }}} |
| Line 23: | Line 46: |
| {{{ cd $SUBJECTS_DIR/bert/mri freeview nu.mgz \ -p-labels posterior_left_* posterior_Left-Hippocampus.mgz \ -p-labels posterior_right_* posterior_Right-Hippocampus.mgz \ -p-prefix posterior_ -p-lut $FREESURFER_HOME/FreeSurferColorLUT.txt }}} |
|
| Line 24: | Line 54: |
| cd $SUBJECTS_DIR/bert/mri freeview nu.mgz -p-labels posterior_left_* \ posterior_right_* \ posterior_Left-Hippocampus.mgz \ posterior_Right-Hippocampus.mgz \ -p-prefix posterior_ \ -p-lut $FREESURFER_HOME/data/GEMS/FreeSurferColorLUT_edited.txt |
The segmentation results are given as probability maps: a (unsigned 8-bit) value of 0 means a probability is 0, 225 means a probability of 1, and everthing in between scales linearly. |
| Line 32: | Line 56: |
| Report bugs to <freesurfer@nmr.mgh.harvard.edu> | === Log files === Since the console output generated during subfields segmentation is long, the log is written in an alternate file - `$SUBJECTS_DIR/bert/scripts/hippo-subfields.log`. The file `recon-all.log` simply shows the invocation of various commands during the subfields generation. Report bugs to <freesurfer@nmr.mgh.harvard.edu> |
Segmentation of hippocampal subfields
Author: Koen Van Leemput
DEPRECATED - This functionality is found only in v5.2 and v5.3.
See HippocampalSubfields for the functionality found in v6.0 and downloadable.
See also: GEMS
Automated segmentation of the subfields of the hippocampus, using models built from manual segmentations of the right hippocampus in 0.38x0.38x0.8 mm^3 in vivo MRI scans in 10 subjects. For the left hippocampus, the model for the right side is simply mirrored. The exact algorithm is described in the paper:
Automated Segmentation of Hippocampal Subfields from Ultra-High Resolution In Vivo MRI. K. Van Leemput, A. Bakkour, T. Benner, G. Wiggins, L.L. Wald, J. Augustinack, B.C. Dickerson, P. Golland, B. Fischl. Hippocampus, vol. 19, no. 6, pp. 549-557, June 2009
recon-all stream
To generate the hippocampal subfield segmentations as part of the normal FreeSurfer processing stream, simply add the "-hippo-subfields" flag to the recon-all script:
recon-all -s bert -all -hippo-subfields
It is also possible to only run the hippocampal subfield segmentation for a subject that has already undergone the standard volumetric FreeSurfer pipeline (in particular, the program needs the FreeSurfer-generated files "aseg.mgz", "nu.mgz" and "talairach.xfm"):
recon-all -s bert -hippo-subfields
The output consists of a collection of images $SUBJECTS_DIR/bert/mri/posterior_* that indicate each voxel's posterior probability of belonging to different subregions. The estimated volumes of these subregions can be obtained using the "kvlQuantifyHippocampalSubfieldSegmentations.sh" script. Doing
kvlQuantifyHippocampalSubfieldSegmentations.sh
will create the text files "nonPartialVolumeStatsLeft.txt" and "nonPartialVolumeStatsRight.txt" in the directory $SUBJECTS_DIR, containing tables that summarize the volumes of the hippocampal subfield segmentations of the left and right side, respectively, in all of the subjects in $SUBJECTS_DIR. The reported volumes are given in number of voxels, where each voxel is of the size "0.5x0.5x0.5mm^3". It is expected that the "-hippo-subfields" flag has been used for each subject in $SUBJECTS_DIR; if this is not the case the "kvlQuantifyHippocampalSubfieldSegmentations.sh" script will fail and no (or incomplete) output files will be generated.
The estimated hippocampal subfield volumes can also be computed for one specific subject only: For the right hippocampus, the command is
cd $SUBJECTS_DIR/bert/mri kvlQuantifyPosteriorProbabilityImages $FREESURFER_HOME/data/GEMS/compressionLookupTable.txt \ posterior_right_* posterior_Right-Hippocampus.mgz
and for the left
cd $SUBJECTS_DIR/bert/mri kvlQuantifyPosteriorProbabilityImages $FREESURFER_HOME/data/GEMS/compressionLookupTable_left.txt \ posterior_left_* posterior_Left-Hippocampus.mgz
Visualizing the subregion segmentations on top of the MRI data can be achieved by doing
cd $SUBJECTS_DIR/bert/mri freeview nu.mgz \ -p-labels posterior_left_* posterior_Left-Hippocampus.mgz \ -p-labels posterior_right_* posterior_Right-Hippocampus.mgz \ -p-prefix posterior_ -p-lut $FREESURFER_HOME/FreeSurferColorLUT.txt
The segmentation results are given as probability maps: a (unsigned 8-bit) value of 0 means a probability is 0, 225 means a probability of 1, and everthing in between scales linearly.
Log files
Since the console output generated during subfields segmentation is long, the log is written in an alternate file - $SUBJECTS_DIR/bert/scripts/hippo-subfields.log. The file recon-all.log simply shows the invocation of various commands during the subfields generation.
Report bugs to <freesurfer@nmr.mgh.harvard.edu>