Overview:
The human cortical surface is a highly complex,
folded structure. Sulci, the spaces between the folds,
define location on the cortex and provide a parcellation
into anatomically distinct areas. A topic that has
recently received increased attention is the segmentation
of these sulci from magnetic resonance images,
with most work focusing on extracting either the sulcal
spaces between the folds or curve representations
of sulci. Unlike these methods, we propose a technique
that extracts actual regions of the cortical surface that
surround sulci, which we call "sulcal regions." The
method is based on a watershed algorithm applied to a
geodesic depth measure on the cortical surface. A
well-known problem with the watershed algorithm is
a tendency toward over segmentation, meaning that a
single region is segmented as several pieces. To address
this problem, we propose a post processing algorithm
that merges appropriate segments from the watershed
algorithm. The sulcal regions are then
manually labeled by simply selecting the appropriate
regions with a mouse click and a preliminary study of
sulcal depth is reported. Finally, a scheme is presented
for computing a complete parcellation of the
cortical surface.
Introduction:
Quantitative anatomic studies of the human cortex
are challenging due to its highly complex, convoluted
folding pattern. The various folds, called gyri, and the
spaces between the folds, called sulci, define location
on the cortex and provide a parcellation of the cortex
into anatomically distinct areas. With the advancement
of magnetic resonance (MR) imaging techniques,
high-resolution, high-contrast three-dimensional images
of the brain can now be routinely acquired in vivo.
As a result, methods for modeling the cortical surface from
these images have emerged, providing a means
for furthering the understanding of morphometric variability
in human populations. A topic that has recently
received increased attention is the study of sulci, in
particular, their segmentation from MR images.
Previous work in the segmentation of sulci has focused
on fitting a surface, finding a set of
points, or extracting volumetric regions within sulcal spaces.
Other work has focused on
extracting curve representations of the sulci. Unlike these
methods, we propose a technique that segments the
actual cortical regions surrounding sulci. The advantage
of segmenting actual cortical regions is that it
allows for a direct geometric study of the cortical surface
and provides a means for mapping functional activation
sites. For ease of terminology, we refer to our
segmented regions as sulcal regions, which define the
buried regions of cortex surrounding the sulcal spaces.
Another advantage of the proposed method is that it
segments sulcal regions on the medial surface as well
as the lateral and inferior surfaces. In addition, this
segmentation method is completely automated and can
be used as a visualization tool for manual labeling of
individual sulcal regions.
Method:
A cortical surface is extracted from volumetric MR data using
CRUISE
and the surface is mapped to the sphere for computational and
visualization purposes. A typical surface and its spherical map
(displaying mean curvature) is shown below.
Sulcal regions are "buried cortices" as shown below
Watershed segmentation proceeds by computationally immersing the underside of the cortex into
a water bath as illustrated below
To accomplish this, the hemispheres are first separated and then a deformable surface "shrinkwrap"
algorithm is run on each hemisphere, thereby defining the outer surface of each hemisphere.
All nodes of the triangle mesh representing a cortex that are touching the shrinkwrap
surface are deemed to be "gyral regions." Two distance maps are computed: a 3D Euclidean map and a
geodesic distance map. The geodesic map is computed using a fast marching algorithm run on the triangle
mesh. All nodes greater (Euclidean distance) from the gyral regions are considered to be sulcal regions.
A Euclidean distance map and a sulcal segmentation are shown below:
A geodesic map is shown below:
A 2D watershed algorithm is then run on the sulcal regions that were previously
segmented, yielding finer regions, as shown below:
This result is an "over segmentation" - very typical of watershed algorithms - and therefore
a merging of catchment basins operation is done, yielding the result below:
|
Results:
Over 1,000 brains have been processed using this algorithm, which is fully automatic.
A typical result is shown in cross sections below:
A measure of average geodesic depth in the Central, Sylvian, and Superior Temporal Sulci
was computed for 15 subjects, and is given below:
Geodesic Depth in Milimeters (N = 15)
Sulcus | Mean | SD | Max | Min |
---|
Ce (L) | 20.5 | 2.27 | 24.5 | 16.1 |
Ce (R) | 20.1 | 1.20 | 22.0 | 18.1 |
Syl (L) | 40.8 | 2.71 | 45.0 | 36.7 |
Syl (R) | 39.3 | 2.89 | 46.2 | 35.7 |
ST (L) | 19.1 | 1.57 | 22.1 | 16.4 |
ST (R) | 21.2 | 2.64 | 25.6 | 16.0 |
Conclusion:
Watershed segmentation of sulcal regions can be used for automatic identification
of sulcal regions in the analysis of changes due to aging and disease or for the
automatic generation of cortical landmarks for applications such as intersubject
registration. These methods have been used to analyze changes in cortical
thickness and shape in (Rettmann et al., 2006).
Publications
- M. E. Rettmann, C. Xu, D. L. Pham, D. N. Yu, and J. L. Prince,
"On Automated Segmentation of Buried Gyri," 5th
International Conference on Functional Mapping of the Human Brain,
June 1999. NeuroImage, vol. 9, no. 6, pg. 143, 1999.
- M. E. Rettmann, C. Xu, D. L. Pham, and J. L. Prince, "Automated
Segmentation of Sulcal Regions", Second International Conference
on Medical Image Computing and Computer Assisted Interventions
(MICCAI'99), pp. 158-167, September 1999.
- X. Han, M. E. Rettmann, C. Xu, and J. L. Prince, "Morphology on
Triangle Meshes Using Geodesic Distance," Conf. on
Info. Sci. Sys., Princeton University March 15-17, 2000.
- M. E. Rettmann, X. Han, D. L. Pham, J. L. Prince, "Geodesics for
Sulcal Segmentation and Depth Measurements," NeuroImage Human
Brain Mapping 2000 Meeting, Poster No. 667, NeuroImage Vol. 11, No. 5,
May 2000.
- M. E. Rettmann, X. Han, and J. L. Prince, "Watersheds on the
cortical surface for automated sulcal segmentation," Proceedings
IEEE Workshop on Mathematical Methods in Biomedical Image Analysis,
Hilton Head Island, SC, Jun 11-12, p.20-27, June 2000.
- X. Han, M. E. Rettmann, C. Xu, and J. L. Prince, "Automatic
segmentation editing for cortical surface reconstruction," Proc. SPIE
Medical Imaging, Conf. 4322, Paper 22, San Diego, February 2001.
- M. E. Rettmann, X. Han, J. L. Prince, "Automated Parcellation of
the Cortical Surface for Computation of Regional Gyrification
Indices," Abstract 230, Human Brain Mapping, Brighton UK, 2001.
- X. Tao, X. Han, M.E. Rettmann, J. L. Prince, and C. A. Davatzikos,
"Statistical Study on Cortical Sulci of Human Brains," Information
Processing in Medical Imaging, 17th Int'l Conference, Davis CA,
pp.475-487, June 2001.
- M. E. Rettman, X. Han, and J. L. Prince, "Automated Sulcal
Segmentation Using Watersheds on the Cortical Surface," NeuroImage,
vol. 15, no. 2, pp. 329-344, February 1, 2002
- M.E. Rettmann, X. Tao, J. L. Prince, "Assisted labeling
techniques for the human brain cortex," Proc. SPIE Medical Imaging
2002, vol. 4684, pp. 179-190, San Diego, CA, Feb. 23-28, 2002.
- K. Behnke, M. Rettmann, D. Pham, D. Shen, S. Resnick, C. Davatzikos,
and J. L. Prince, "Automatic classification of sulcal regions of
the human brain cortex using pattern recognition," in Proc. SPIE's
Medical Imaging, San Diego, CA, Feb. 15-20, 2003.
- M. E. Rettmann, J. L. Prince, and S. M. Resnick. "Analysis of
Sulcal Shape Changes Associated with Aging," in Proc. Human Brain
Mapping, New York, NY, June 18-22, 2003.
- D. Tosun, M. E. Rettmann, and J. L. Prince, "Mapping Techniques
for Aligning Sulci Across Multiple Brains," Conference on Medical
Image Computing and Computer Assisted Intervention, Montreal, 14-18
November 2003.
- D. Tosun, M.E. Rettmann, and J. L. Prince, "Mapping Techniques
for Aligning Sulci Across Multiple Brains," Medical Image Analysis,
vol.8, pp. 295-309, 2004.
- D. Tosun, M. E. Rettmann, X. Han, X. Tao, C. Xu, S. M. Resnick,
D. Pham, and J. L. Prince, "Cortical Surface Segmentation and
Mapping," NeuroImage, vol. 23, pp. S108-S118, 2004.
- X. Han, D.L. Pham, D. Tosun, M.E. Rettmann, C. Xu, and J. L. Prince,
"CRUISE: Cortical Reconstruction Using Implicit Surface
Evolution," NeuroImage, vol. 23, pp. 997-1012, 2004.
- M.E. Rettmann, D. Tosun, X. Tao, S.M. Resnick, and J.L. Prince,
"Program for the Assisted Labeling of Sulcal Regions (PALS):
Description and Reliability," NeuroImage, volume 24, issue 2,
pp. 398-416, 15 January 2005.
- M.E. Rettmann, D. Tosun , X. Tao, S.M. Resnick, and J.L. Prince,
"Mapping Cross-Sectional Differences in Cortical Thickness During
Aging," 11th Annual Meeting of the Organization for Human Brain
Mapping (HBM), Toronto, Ontario, Canada, June 12-16, 2005.
- M. E. Rettmann, M. A. Kraut, J. L. Prince, and S. M. Resnick,
"Cross-Sectional and Longitudinal Analyses of Anatomical Sulcal
Changes Associated with Aging," Cerebral Cortex, 16:1584-1594,
November 2006.