Extracting line representations of sulcal and gyral patterns in MR images of the human brain

Cortical complexity alterations in motor subtypes of Parkinson's disease: A surface-based morphometry analysis of fractal dimension

Based on motor symptoms, Parkinson's disease (PD) can be classified into tremor dominant (TD) and postural instability gait difficulty (PIGD) subtypes. Few studies have examined cortical complexity differences in PD motor subtypes. This study aimed to investigate differences in cortical complexity and grey matter volume (GMV) between TD and PIGD. We enrolled 36 TD patients, 27 PIGD patients and 66 healthy controls (HC) from the PPMI (Parkinson's Progression Markers Initiative) database. Voxel-based morphometry (VBM) and surface-based morphometry (SBM) were utilized to assess differences in GMV, cortical thickness and cortical complexity. The structural MRI data of participants was analysed using CAT12/SPM12 (p < 0.05, FDR corrected). Additionally, correlations between clinical data and structural changes were examined (p < 0.05, Holm-Bonferroni corrected). In comparison to both HC and TD groups, PIGD patients exhibited a significant fractal dimension (FD) decrease in many cortical regions. A significant negative correlation between age and FD was observed in the left insula for the PIGD patients and in the bilateral insula for the TD patients. However, no significant differences were found in GMV, cortical thickness or other complexity indices. Altered FD in the bilateral insula indicates that postural instability and gait disturbances may result from a failure to integrate information from various structures, whereas parkinsonian rest tremor is not associated with this integration. Also, widespread decreases in cortical FD demonstrate that FD is more sensitive than other complexity measures and can serve as a novel biomarker for identifying subtle changes in cortical morphology in the PIGD subtype.

Toward a better understanding of how a gyrified brain develops

The size and shape of the cerebral cortex have changed dramatically across evolution. For some species, the cortex remains smooth (lissencephalic) throughout their lifetime, while for other species, including humans and other primates, the cortex increases substantially in size and becomes folded (gyrencephalic). A folded cortex boasts substantially increased surface area, cortical thickness, and neuronal density, and it is therefore associated with higher-order cognitive abilities. The mechanisms that drive gyrification in some species, while others remain lissencephalic despite many shared neurodevelopmental features, have been a topic of investigation for many decades, giving rise to multiple perspectives of how the gyrified cerebral cortex acquires its unique shape. Recently, a structurally unique germinal layer, known as the outer subventricular zone, and the specialized cell type that populates it, called basal radial glial cells, were identified, and these have been shown to be indispensable for cortical expansion and folding. Transcriptional analyses and gene manipulation models have provided an invaluable insight into many of the key cellular and genetic drivers of gyrification. However, the degree to which certain biomechanical, genetic, and cellular processes drive gyrification remains under investigation. This review considers the key aspects of cerebral expansion and folding that have been identified to date and how theories of gyrification have evolved to incorporate this new knowledge.

Automatic Quantification of Normal Brain Gyrification Patterns and Changes in Fetuses with Polymicrogyria and Lissencephaly Based on MRI

BACKGROUND AND PURPOSE

The current imaging assessment of fetal brain gyrification is performed qualitatively and subjectively using sonography and MR imaging. A few previous studies have suggested methods for quantification of fetal gyrification based on 3D reconstructed MR imaging, which requires unique data and is time-consuming. In this study, we aimed to develop an automatic pipeline for gyrification assessment based on routinely acquired fetal 2D MR imaging data, to quantify normal changes with gestation, and to measure differences in fetuses with lissencephaly and polymicrogyria compared with controls.

MATERIALS AND METHODS

We included coronal T2-weighted MR imaging data of 162 fetuses retrospectively collected from 2 clinical sites: 134 controls, 12 with lissencephaly, 13 with polymicrogyria, and 3 with suspected lissencephaly based on sonography, yet with normal MR imaging diagnoses. Following brain segmentation, 5 gyrification parameters were calculated separately for each hemisphere on the basis of the area and ratio between the contours of the cerebrum and its convex hull. Seven machine learning classifiers were evaluated to differentiate control fetuses and fetuses with lissencephaly or polymicrogyria.

RESULTS

In control fetuses, all parameters changed significantly with gestational age (P < .05). Compared with controls, fetuses with lissencephaly showed significant reductions in all gyrification parameters (P ≤ .02). Similarly, significant reductions were detected for fetuses with polymicrogyria in several parameters (P ≤ .001). The 3 suspected fetuses showed normal gyrification values, supporting the MR imaging diagnosis. An XGBoost-linear algorithm achieved the best results for classification between fetuses with lissencephaly and control fetuses (n = 32), with an area under the curve of 0.90 and a recall of 0.83. Similarly, a random forest classifier showed the best performance for classification of fetuses with polymicrogyria and control fetuses (n = 33), with an area under the curve of 0.84 and a recall of 0.62.

CONCLUSIONS

This study presents a pipeline for automatic quantification of fetal brain gyrification and provides normal developmental curves from a large cohort. Our method significantly differentiated fetuses with lissencephaly and polymicrogyria, demonstrating lower gyrification values. The method can aid radiologic assessment, highlight fetuses at risk, and may improve early identification of fetuses with cortical malformations.

On the definition, construction, and presentation of the human cerebral sulci: A morphology-based approach

Although the term sulcus is known for almost four centuries, its formal, precise, consistent, constructive, and quantitative definition is practically lacking. As the cerebral sulci (and gyri) are vital in cortical anatomy which, in turn, is central in neuroeducation and neuroimage processing, a new sulcus definition is needed. The contribution of this work is threefold, namely to (1) propose a new, morphology-based definition of the term sulcus (and consequently that of gyrus), (2) formulate a constructive method for sulcus calculation, and (3) provide a novel way for the presentation of sulci. The sulcus is defined here as a volumetric region on the cortical mantle between adjacent gyri separated from them at the levels of their gyral white matter crest lines. Consequently, the sulcal inner surface is demarcated by the crest lines of the gyral white matter of its adjacent gyri. Correspondingly, the gyrus is defined as a volumetric region on the cortical mantle separated from its adjacent sulci at the level of its gyral white matter crest line. This volumetric sulcus definition is conceptually simple, anatomy-based, educationally friendly, quantitative, and constructive. Considering the sulcus as a volumetric object is a major differentiation from other works. Based on the introduced sulcus definition, a method for volumetric sulcus construction is proposed in two, conceptually straightforward, steps, namely, sulcal intersection formation followed by its propagation which steps are to be repeated for every sulcal segment. These sulcal and gyral constructions can be automated by applying existing methods and public tools. As a volumetric sulcus forms an imprint into the white matter, this enables prominent sulcus presentation. Since this type of presentation is novel yet unfamiliar to the reader, also a dual surface presentation was proposed here by employing the spatially co-registered white matter and cortical surfaces. The results were presented as dual surface labeled sulci on eight standard orthogonal views, anterior, left lateral, posterior, right lateral, superior, inferior, medial left, and medial right by using a 3D brain atlas. Moreover, additional 108 labeled images were created with sulcus-oriented views for 27 individual left and right sulci forming 54 dual white matter-cortical surface images strengthening in this way the educational value of the proposed approach. These images were included for public use in the NOWinBRAIN neuroimage repository with over 7700 3D images available at www.nowinbrain.org. The results demonstrated the superiority of white matter surface sulci presentation over the standard cortical surface and cross-sectional presentations in terms of sulcal course, continuity, size, shape, width, depth, side branches, and pattern. To my best knowledge, this is the first work ever presenting the labeling of sulci on all cerebral white matter surfaces as well as on dual white matter-cortical surfaces. Additionally to neuroeducation, three other applications of the proposed approach were discussed, sulcal reference maps, sulcus quantification in terms of new parameters introduced here (sulcal volume, wall skewness, and the number of white matter basins), and an atlas-assisted tool for exploration and studying of cerebral sulci and gyri .

Aberrant Changes in Cortical Complexity in Right-Onset Versus Left-Onset Parkinson's Disease in Early-Stage

There is increasing evidence to show that motor symptom lateralization in Parkinson’s disease (PD) is linked to non-motor features, progression, and prognosis of the disease. However, few studies have reported the difference in cortical complexity between patients with left-onset of PD (LPD) and right-onset of PD (RPD). This study aimed to investigate the differences in the cortical complexity between early-stage LPD and RPD. High-resolution T1-weighted magnetic resonance images of the brain were acquired in 24 patients with LPD, 34 patients with RPD, and 37 age- and sex-matched healthy controls (HCs). Cortical complexity including gyrification index, fractal dimension (FD), and sulcal depth was analyzed using surface-based morphometry via CAT12/SPM12. Familywise error (FWE) peak-level correction at p < 0.05 was performed for significance testing. In patients with RPD, we found decreased mean FD and mean sulcal depth in the banks of the left superior temporal sulcus (STS) compared with LPD and HCs. The mean FD in the left superior temporal gyrus (STG) was decreased in RPD compared with HCs. However, in patients with LPD, we did not identify significantly abnormal cortical complex change compared with HCs. Moreover, we observed that the mean FD in STG was negatively correlated with the 17-item Hamilton Depression Scale (HAMD) among the three groups. Our findings support the specific influence of asymmetrical motor symptoms in cortical complexity in early-stage PD and reveal that the banks of left STS and left STG might play a crucial role in RPD.

Fundamental functional differences between gyri and sulci: implications for brain function, cognition, and behavior

Folding of the cerebral cortex is a prominent characteristic of mammalian brains. Alterations or deficits in cortical folding are strongly correlated with abnormal brain function, cognition, and behavior. Therefore, a precise mapping between the anatomy and function of the brain is critical to our understanding of the mechanisms of brain structural architecture in both health and diseases. Gyri and sulci, the standard nomenclature for cortical anatomy, serve as building blocks to make up complex folding patterns, providing a window to decipher cortical anatomy and its relation with brain functions. Huge efforts have been devoted to this research topic from a variety of disciplines including genetics, cell biology, anatomy, neuroimaging, and neurology, as well as involving computational approaches based on machine learning and artificial intelligence algorithms. However, despite increasing progress, our understanding of the functional anatomy of gyro-sulcal patterns is still in its infancy. In this review, we present the current state of this field and provide our perspectives of the methodologies and conclusions concerning functional differentiation between gyri and sulci, as well as the supporting information from genetic, cell biology, and brain structure research. In particular, we will further present a proposed framework for attempting to interpret the dynamic mechanisms of the functional interplay between gyri and sulci. Hopefully, this review will provide a comprehensive summary of anatomo-functional relationships in the cortical gyro-sulcal system together with a consideration of how these contribute to brain function, cognition, and behavior, as well as to mental disorders.

Associations Between Individual Differences in Mathematical Competencies and Surface Anatomy of the Adult Brain

Previously conducted structural magnetic resonance imaging (MRI) studies on the neuroanatomical correlates of mathematical abilities and competencies have several methodological limitations. Besides small sample sizes, the majority of these studies have employed voxel-based morphometry (VBM)-a method that, although it is easy to implement, has some major drawbacks. Taking this into account, the current study is the first to investigate in a large sample of typically developed adults the associations between mathematical abilities and variations in brain surface structure by using surface-based morphometry (SBM). SBM is a method that also allows the investigation of brain morphometry by avoiding the pitfalls of VBM. Eighty-nine young adults were tested with a large battery of psychometric tests to measure mathematical competencies in four different areas: (1) simple arithmetic; (2) complex arithmetic; (3) higher-order mathematics; and (4) numerical intelligence. Also, we asked participants for their mathematics grades for their final school exams. Inside the MRI scanner, we collected high-resolution T1-weighted anatomical images from each subject. SBM analyses were performed with the computational anatomy toolbox (CAT12) and indices for cortical thickness, for cortical surface complexity, for gyrification, and sulcal depth were calculated. Further analyses revealed associations between: (1) the cortical surface complexity of the right superior temporal gyrus and numerical intelligence; (2) the depth of the right central sulcus and adults' ability to solve complex arithmetic problems; and (3) the depth of the left parieto-occipital sulcus and adults' higher-order mathematics competence. Interestingly, no relationships with previously reported brain regions were observed, thus, suggesting the importance of similar research to confirm the role of the brain regions found in this study.

Determinants of structural segregation and patterning in the human cortex

This study aimed at uncovering mechanisms that govern the spatio-temporal patterning of the human cortex and its structural variability, and drawing links between fetal brain development and variability in adult brains. A data-driven analytic approach based on structural MR images revealed the following findings: (1) The cortical surface can be subdivided into 13 independent regions ("communities") based on macroscopic features. (2) Thirty centers of low inter-subject variability were found in major sulci on the cortical surface. Their variability showed a strong positive correlation with the known time points at which they appear in fetal development. Centers forming early induce a higher inter-subject regularity in a larger local vicinity, while those forming later result in smaller regions of higher variability. (3) The layout of sulcal and gyral patterns within a community is governed typically by two centers. Depending on the relative variability of each center, communities can be classified into structural sub-types. (4) Sub-types across ipsi-lateral communities are independent, but associated with the sub-type of the same community on the contra-lateral side. Results shown here integrate well with current knowledge about macroscopic, microscopic, and genetic determinants of brain development.

Gyral net: A new representation of cortical folding organization

One distinct feature of the cerebral cortex is its convex (gyri) and concave (sulci) folding patterns. Due to the remarkable complexity and variability of gyral/sulcal shapes, it has been challenging to quantitatively model their organization patterns. Inspired by the observation that the lines of gyral crests can form a connected graph on each brain hemisphere, we propose a new representation of cortical gyri/sulci organization pattern - gyral net, which models cortical architecture from a graph perspective, starting with nodes and edges obtained from the reconstructed cortical surfaces. A novel computational framework is developed to efficiently and automatically construct gyral nets from surface meshes, and four measurements are devised to quantify the folding patterns. Using an MRI dataset for autism study as a test bed, we identified reduced local connectivity cost and increased curviness of gyral net bilaterally on the parietal lobe, occipital lobe, and temporal lobe in autistic patients. Additionally, we found that the cortical thickness and the gyral straightness of gyral joints are higher than the rest of gyral crest regions. The proposed representation offers a new tool for a comprehensive and reliable characterization of the cortical folding organization. This novel computational framework will enable large-scale analyses of cortical folding patterns in the future.

Simultaneous Cortical Surface Labeling and Sulcal Curve Extraction

Automatic labeling of the gyri and sulci on the cortical surface is important for studying cortical morphology and brain functions within populations. A method to simultaneously label gyral regions and extract sulcal curves is proposed. Assuming that the gyral regions parcellate the whole cortical surface into contiguous regions with certain fixed topology, the proposed method labels the subject cortical surface by deformably registering a network of curves that form the boundary of gyral regions to the subject cortical surface. In the registration process, the curves are encouraged to follow the fine details of the sulcal geometry and to observe the shape statistics learned from training data. Using the framework of probabilistic point set registration methods, the proposed algorithm finds the sulcal curve network that maximizes the posterior probability by Expectation-Maximization (EM). The automatic labeling method was evaluated on 15 cortical surfaces using a leave-one-out strategy. Quantitative error analysis is carried out on both labeled regions and major sulcal curves.

Automated extraction of nested sulcus features from human brain MRI data

Extracting objects related to a fold in the cerebral cortex ("sulcus features") from human brain magnetic resonance imaging data has applications in morphometry, landmark-based registration, and anatomical labeling. In prior work, sulcus features such as surfaces, fundi and pits have been extracted separately. Here we define and extract nested sulcus features in a hierarchical manner from a cortical surface mesh having curvature or depth values. Our experimental results show that the nested features are comparable to features extracted separately using other methods, and that they are consistent across subjects and with manual label boundaries. Our open source feature extraction software will be made freely available as part of the Mindboggle project (http://www.mindboggle.info).

Multimodality data integration in epilepsy

An important goal of software development in the medical field is the design of methods which are able to integrate information obtained from various imaging and nonimaging modalities into a cohesive framework in order to understand the results of qualitatively different measurements in a larger context. Moreover, it is essential to assess the various features of the data quantitatively so that relationships in anatomical and functional domains between complementing modalities can be expressed mathematically. This paper presents a clinically feasible software environment for the quantitative assessment of the relationship among biochemical functions as assessed by PET imaging and electrophysiological parameters derived from intracranial EEG. Based on the developed software tools, quantitative results obtained from individual modalities can be merged into a data structure allowing a consistent framework for advanced data mining techniques and 3D visualization. Moreover, an effort was made to derive quantitative variables (such as the spatial proximity index, SPI) characterizing the relationship between complementing modalities on a more generic level as a prerequisite for efficient data mining strategies. We describe the implementation of this software environment in twelve children (mean age 5.2 +/- 4.3 years) with medically intractable partial epilepsy who underwent both high-resolution structural MR and functional PET imaging. Our experiments demonstrate that our approach will lead to a better understanding of the mechanisms of epileptogenesis and might ultimately have an impact on treatment. Moreover, our software environment holds promise to be useful in many other neurological disorders, where integration of multimodality data is crucial for a better understanding of the underlying disease mechanisms.

Consistent sulcal parcellation of longitudinal cortical surfaces

Automated accurate and consistent sulcal parcellation of longitudinal cortical surfaces is of great importance in studying longitudinal morphological and functional changes of human brains, since longitudinal cortical changes are normally very subtle, especially in aging brains. However, applying the existing methods (which were typically developed for cortical sulcal parcellation of a single cortical surface) independently to longitudinal cortical surfaces might generate longitudinally-inconsistent results. To overcome this limitation, this paper presents a novel energy function based method for accurate and consistent sulcal parcellation of longitudinal cortical surfaces. Specifically, both spatial and temporal smoothness are imposed in the energy function to obtain consistent longitudinal sulcal parcellation results. The energy function is efficiently minimized by a graph cut method. The proposed method has been successfully applied to sulcal parcellation of both real and simulated longitudinal inner cortical surfaces of human brain MR images. Both qualitative and quantitative evaluation results demonstrate the validity of the proposed method.

An improved representation of regional boundaries on parcellated morphological surfaces

Establishing the correspondences of brain anatomy with function is important for understanding neuroimaging data. Regional delineations on morphological surfaces define anatomical landmarks and help to visualize and interpret both functional data and morphological measures mapped onto the cortical surface. We present an efficient algorithm that accurately delineates the morphological surface of the cerebral cortex in real time during generation of the surface using information from parcellated 3D data. With this accurate region delineation, we then develop methods for boundary-preserved simplification and smoothing, as well as procedures for the automated correction of small, misclassified regions to improve the quality of the delineated surface. We demonstrate that our delineation algorithm, together with a new method for double-snapshot visualization of cortical regions, can be used to establish a clear correspondence between brain anatomy and mapped quantities, such as morphological measures, across groups of subjects.

Cortical sulci recognition and spatial normalization

Brain mapping techniques pair similar anatomical information across individuals. In this context, spatial normalization is mainly used to reduce inter-subject differences to improve comparisons. These techniques may benefit from anatomically identified landmarks useful to drive the registration. Automatic labeling, classification or segmentation techniques provide such labels. Most of these approaches depend strongly on normalization, as much as normalization depends on landmark accuracy. We propose in this paper a coherent Bayesian framework to automatically identify approximately 60 sulcal labels per hemisphere based on a probabilistic atlas (a mixture of spam models: Statistical Probabilistic Anatomy Map) estimating simultaneously normalization parameters. This way, the labelization method provides also with no extra computational costs a new automatically constrained registration of sulcal structures. We have limited our study to global affine and piecewise affine registration. The suggested global affine approach outperforms significantly standard affine intensity-based normalization techniques in term of sulci alignments. Further, by combining global and local joint labeling, a final mean recognition rate of 86% has been obtained with much more reliable labeling posterior probabilities. The different methods described in this paper have been integrated since the release version 3.2.1 of the BrainVISA software platform (Riviére et al., 2009).

Spectral-based automatic labeling and refining of human cortical sulcal curves using expert-provided examples

We present a spectral-based method for automatically labeling and refining major sulcal curves of a human cerebral cortex. Given a set of input (unlabeled) sulcal curves automatically extracted from a cortical surface and a collection of expert-provided examples (labeled sulcal curves), our objective is to identify the input major sulcal curves and assign their neuroanatomical labels, and then refines these curves based on the expert-provided example data, without employing any atlas-based registration scheme as preprocessing. In order to construct the example data, neuroanatomists manually labeled a set of 24 major sulcal curves (12 each for the left and right hemispheres) for each individual subject according to a precise protocol. We collected 30 sets of such curves from 30 subjects. Given the raw input sulcal curve set of a subject, we choose the most similar example curve to each input curve in the set to label and refine the latter according to the former. We adapt a spectral matching algorithm to choose the example curve by exploiting the sulcal curve features and their relationship. The high dimensionality of sulcal curve data in spectral matching is addressed by using their multi-resolution representations, which greatly reduces time and space complexities. Our method provides consistent labeling and refining results even under high variability of cortical sulci across the subjects. Through experiments we show that the results are comparable in accuracy to those done manually. Most output curves exhibited accuracy values higher than 80%, and the mean accuracy values of the curves in the left and the right hemispheres were 84.69% and 84.58%, respectively.

An automated pipeline for cortical sulcal fundi extraction

In this paper, we propose a novel automated pipeline for extraction of sulcal fundi from triangulated cortical surfaces. This method consists of four consecutive steps. Firstly, we adopt a finite difference method to estimate principal curvatures, principal directions and curvature derivatives, along the principal directions, for each vertex. Then, we detect the sulcal fundi segment in each triangle of the cortical surface based on curvatures and curvature derivatives. Afterwards, we link the sulcal fundi segments into continuous curves. Finally, we connect breaking sulcal fundi and smooth bumping sulcal fundi by using the fast marching method on the cortical surface. The proposed method can find the accurate sulcal fundi using curvatures and curvature derivatives without any manual interaction. The method was applied to 10 normal brain MR images on inner cortical surfaces. We quantitatively evaluated the accuracy of the sulcal fundi extraction method using manually labeled sulcal fundi by experts. The average difference between automatically extracted major sulcal fundi and the expert labeled results is consistently around 1.0mm on 10 subject images, indicating the good performance of the proposed method.

The development of gyrification in childhood and adolescence

Gyrification is the process by which the brain undergoes changes in surface morphology to create sulcal and gyral regions. The period of greatest development of brain gyrification is during the third trimester of pregnancy, a period of time in which the brain undergoes considerable growth. Little is known about changes in gyrification during childhood and adolescence, although considering the changes in gray matter volume and thickness during this time period, it is conceivable that alterations in the brain surface morphology could also occur during this period of development. The formation of gyri and sulci in the brain allows for compact wiring that promotes and enhances efficient neural processing. If cerebral function and form are linked through the organization of neural connectivity, then alterations in neural connectivity, i.e., synaptic pruning, may also alter the gyral and sulcal patterns of the brain. This paper reviews developmental theories of gyrification, computational techniques for measuring gyrification, and the potential interaction between gyrification and neuronal connectivity. We also present recent findings involving alterations in gyrification during childhood and adolescence.

Automatic cortical sulcal parcellation based on surface principal direction flow field tracking

Automatic parcellation of cortical surfaces into sulcal based regions is of great importance in structural and functional mapping of human brain. In this paper, a novel method is proposed for automatic cortical sulcal parcellation based on the geometric characteristics of the cortical surface including its principal curvatures and principal directions. This method is composed of two major steps: 1) employing the hidden Markov random field model (HMRF) and the expectation maximization (EM) algorithm on the maximum principal curvatures of the cortical surface for sulcal region segmentation, and 2) using a principal direction flow field tracking method on the cortical surface for sulcal basin segmentation. The flow field is obtained by diffusing the principal direction field on the cortical surface. The method has been successfully applied to the inner cortical surfaces of twelve healthy human brain MR images. Both quantitative and qualitative evaluation results demonstrate the validity and efficiency of the proposed method.

Automatic extraction of sulcal lines on cortical surfaces based on anisotropic geodesic distance

Analyzing cortical sulci is important for studying cortical morphology and brain functions. Although sulcal lines on cortical surfaces can be defined in various ways, it is critical in a neuroimaging study to define a sulcal line along the valley of a cortical surface with a high curvature within a sulcus. To extract the sulcal lines automatically, we present a new geometric algorithm based on the computation of anisotropic skeletons of sulcal regions. Because anisotropic skeletons are highly adaptive to the anisotropic nature of the surface shape, the resulting sulcal lines lie accurately on the valleys of the sulcal areas. Our sulcal lines remain unchanged under local shape variabilities in different human brains. Through experiments, we show that the errors of the sulcal lines for both synthetic data and real cortical surfaces were nearly as constant as the function of random noise. By measuring the changes in sulcal shape in Alzheimer's disease (AD) patients, we further investigated the effectiveness of the accuracy of our sulcal lines using a large sample of MRI data. This study involved 70 normal controls (n [men/women]: 29/41, age [mean+/-SD]: 71.7+/-4.9 years), and 100 AD subjects (37/63, 72.3+/-5.5). We observe significantly lower absolute average mean curvature and shallower sulcal depth in AD subjects, where the group difference becomes more significant if we measure the quantities along the sulcal lines rather than over the entire sulcal area. The most remarkable difference in the AD patients was the average sulcal depth (control: 11.70 and AD: 11.34).

Automatic cortical sulcal parcellation based on surface principal direction flow field tracking

The human cerebral cortex is a highly convoluted structure composed of sulci and gyri, corresponding to the valleys and ridges of the cortical surface respectively. Automatic parcellation of the cortical surface into sulcal regions is of great importance in structural and functional mapping of the human brain. In this paper, a novel method is proposed for automatic cortical sulcal parcellation based on the geometric characteristics of cortical surface including its principal curvatures and principal directions. This method is composed of two major steps: 1) employing the hidden Markov random field model (HMRF) and the expectation maximization (EM) algorithm on the maximum principal curvatures of the cortical surface for sulcal region segmentation, and 2) using a principal direction flow field tracking method on the cortical surface for sulcal basin segmentation. The flow field is obtained by diffusing the principal direction field on the cortical surface mesh. A unique feature of this method is that the automatic sulcal parcellation process is quite robust and efficient, and is independent of any external guidance such as atlas-based warping. The method has been successfully applied to the inner cortical surfaces of twelve healthy human brain MR images. Both quantitative and qualitative evaluation results demonstrate the validity and efficiency of the proposed method.

Joint sulcal detection on cortical surfaces with graphical models and boosted priors

In this paper, we propose an automated approach for the joint detection of major sulci on cortical surfaces. By representing sulci as nodes in a graphical model, we incorporate Markovian relations between sulci and formulate their detection as a maximum a posteriori (MAP) estimation problem over the joint space of major sulci. To make the inference tractable, a sample space with a finite number of candidate curves is automatically generated at each node based on the Hamilton-Jacobi skeleton of sulcal regions. Using the AdaBoost algorithm, we learn both individual and pairwise shape priors of sulcal curves from training data, which are then used to define potential functions in the graphical model based on the connection between AdaBoost and logistic regression. Finally belief propagation is used to perform the MAP inference and select the joint detection results from the sample spaces of candidate curves. In our experiments, we quantitatively validate our algorithm with manually traced curves and demonstrate the automatically detected curves can capture the main body of sulci very accurately. A comparison with independently detected results is also conducted to illustrate the advantage of the joint detection approach.

A novel method for cortical sulcal fundi extraction

Sulcal fundi are 3D curves along the bottom of sulcal regions of the human cerebral cortex. In this paper, we propose a novel automatic method for extraction of sulcal fundi from triangulated cortical surface. Compared to existing methods, the proposed method can find accurate sulcal fundi using curvatures and curvature derivatives without manual interaction. Given a triangulated cortical surface, our method is composed of four steps: estimating curvatures and curvature derivatives for each vertex, detecting the sulcal fundi segments in each triangle, linking the sulcal fundi segments and combining of adjacent sulcal fundi, and connecting breaking sulcal fundi and smoothing using the fast marching method on the cortical surface. The proposed sulcal fundi extraction method is applied to ten normal brain inner cortical surfaces. We quantitatively validated the proposed method of sulcal fundi extraction using manually labeled sulcal fundi by experts as the ground truth.

Hamilton-Jacobi skeleton on cortical surfaces

In this paper, we propose a new method to construct graphical representations of cortical folding patterns by computing skeletons on triangulated cortical surfaces. In our approach, a cortical surface is first partitioned into sulcal and gyral regions via the solution of a variational problem using graph cuts, which can guarantee global optimality. After that, we extend the method of Hamilton-Jacobi skeleton [1] to subsets of triangulated surfaces, together with a geometrically intuitive pruning process that can trade off between skeleton complexity and the completeness of representing folding patterns. Compared with previous work that uses skeletons of 3-D volumes to represent sulcal patterns, the skeletons on cortical surfaces can be easily decomposed into branches and provide a simpler way to construct graphical representations of cortical morphometry. In our experiments, we demonstrate our method on two different cortical surface models, its ability of capturing major sulcal patterns and its application to compute skeletons of gyral regions.

Automatic trajectory planning for deep brain stimulation: a feasibility study

DBS for Parkinson's disease involves an extensive planning to find a suitable electrode implantation path to the selected target. We have investigated the feasibility of improving the conventional planning with an automatic calculation of possible paths in 3D. This requires the segmentation of anatomical structures. Subsequently, the paths are calculated and visualized. After selection of a suitable path, the settings for the stereotactic frame are determined. A qualitative evaluation has shown that automatic avoidance of critical structures is feasible. The participating neurosurgeons estimate the time gain to be around 30 minutes.

A geometric method for automatic extraction of sulcal fundi

Sulcal fundi are 3-D curves that lie in the depths of the cerebral cortex and, in addition to their intrinsic value in brain research, are often used as landmarks for downstream computations in brain imaging. In this paper, we present a geometric algorithm that automatically extracts the sulcal fundi from magnetic resonance images and represents them as spline curves lying on the extracted triangular mesh representing the cortical surface. The input to our algorithm is a triangular mesh representation of an extracted cortical surface as computed by one of several available software packages for performing automated and semi-automated cortical surface extraction. Given this input we first compute a geometric depth measure for each triangle on the cortical surface mesh, and based on this information we extract sulcal regions by checking for connected regions exceeding a depth threshold. We then identify endpoints of each region and delineate the fundus by thinning the connected region while keeping the endpoints fixed. The curves, thus, defined are regularized using weighted splines on the surface mesh to yield high-quality representations of the sulcal fundi. We present the geometric framework and validate it with real data from human brains. Comparisons with expert-labeled sulcal fundi are part of this validation process.

Joint Sulci Detection Using Graphical Models and Boosted Priors

In this paper we propose an automated approach for joint sulci detection on cortical surfaces by using graphical models and boosting techniques to incorporate shape priors of major sulci and their Markovian relations. For each sulcus, we represent it as a node in the graphical model and associate it with a sample space of candidate curves, which is generated automatically using the Hamilton-Jacobi skeleton of sulcal regions. To take into account individual as well as joint priors about the shape of major sulci, we learn the potential functions of the graphical model using AdaBoost algorithm to select and fuse information from a large set of features. This discriminative approach is especially powerful in capturing the neighboring relations between sulcal lines, which are otherwise hard to be captured by generative models. Using belief propagation, efficient inferencing is then performed on the graphical model to estimate each sulcus as the maximizer of its final belief. On a data set of 40 cortical surfaces, we demonstrate the advantage of joint detection on four major sulci: central, precentral, postcentral and the sylvian fissure.

An automated labeling system for subdividing the human cerebral cortex on MRI scans into gyral based regions of interest

In this study, we have assessed the validity and reliability of an automated labeling system that we have developed for subdividing the human cerebral cortex on magnetic resonance images into gyral based regions of interest (ROIs). Using a dataset of 40 MRI scans we manually identified 34 cortical ROIs in each of the individual hemispheres. This information was then encoded in the form of an atlas that was utilized to automatically label ROIs. To examine the validity, as well as the intra- and inter-rater reliability of the automated system, we used both intraclass correlation coefficients (ICC), and a new method known as mean distance maps, to assess the degree of mismatch between the manual and the automated sets of ROIs. When compared with the manual ROIs, the automated ROIs were highly accurate, with an average ICC of 0.835 across all of the ROIs, and a mean distance error of less than 1 mm. Intra- and inter-rater comparisons yielded little to no difference between the sets of ROIs. These findings suggest that the automated method we have developed for subdividing the human cerebral cortex into standard gyral-based neuroanatomical regions is both anatomically valid and reliable. This method may be useful for both morphometric and functional studies of the cerebral cortex as well as for clinical investigations aimed at tracking the evolution of disease-induced changes over time, including clinical trials in which MRI-based measures are used to examine response to treatment.

A learning based algorithm for automatic extraction of the cortical sulci

This paper presents a learning based method for automatic extraction of the major cortical sulci from MRI volumes or extracted surfaces. Instead of using a few pre-defined rules such as the mean curvature properties, to detect the major sulci, the algorithm learns a discriminative model by selecting and combining features from a large pool of candidates. We used the Probabilistic Boosting Tree algorithm to learn the model, which implicitly discovers and combines rules based on manually annotated sulci traced by neuroanatomists. The algorithm almost has no parameters to tune and is fast because of the adoption of integral volume and 3D Haar filters. For a given approximately registered MRI volume, the algorithm computes the probability of how likely it is that each voxel lies on a major sulcus curve. Dynamic programming is then applied to extract the curve based on the probability map and a shape prior. Because the algorithm can be applied to MRI volumes directly, there is no need to perform preprocessing such as tissue segmentation or mapping to a canonical space. The learning aspect makes the approach flexible and it also works on extracted cortical surfaces.

A construction of an averaged representation of human cortical gyri using non-linear principal component analysis

Because of the complex shape of human cortical gyri and great variation between individuals, development of effective representation schemes which allow establishment of correspondence between individuals, extraction of average structure of a population, and co-registration has proved very difficult. We introduce an approach which extracts line representations of gyri at different depths from high resolution MRI, labels main gyri semi-automatically, and extracts a template from a population using non-linear principal component analysis. The method has been tested on data from 96 healthy human volunteers. The model captures the most salient shape features of all major cortical gyri, and can be used for inter-subject registration, for investigating regionalized inter-subject variability, and for inter-hemispheric comparisons.

Mindboggle: a scatterbrained approach to automate brain labeling

Mindboggle (http://www.binarybottle.com/mindboggle.html) is a fully automated, feature matching approach to label cortical structures and activity anatomically in human brain MRI data. This approach does not assume that the existence of component structures and their relative spatial relationship is preserved from brain to brain, but instead disassembles a labeled atlas and reassembles its pieces to match corresponding pieces in an unlabeled subject brain before labeling. Mindboggle: (1) converts linearly coregistered subject and atlas MRI data into sulcus pieces, (2) matches each atlas piece with a combination of subject pieces by minimizing a cost function, (3) transforms atlas label boundaries to the matching subject pieces, (4) warps atlas labels to their transformed boundaries, and (5) propagates labels to fill remaining gaps in a mask derived from the subject brain. We compared Mindboggle with four registration methods: linear registration, and nonlinear registration using SPM2, AIR, and ANIMAL. Automated labeling by all of the nonlinear methods was found to be at least comparable with linear registration. Mindboggle outperformed every other method, as measured by the agreement between overlapping atlas labels and manually assigned subject labels, with respect to the union or the intersection of voxels. After applying the same procedure that Mindboggle uses to fill a subject's segmented gray matter mask with labels (step 5), the results of the other methods improved. However, after performing a one-way ANOVA (and Tukey's honestly significant difference criterion) in a multiple comparison between the results obtained by the different methods, Mindboggle was still found to be the only nonlinear method whose labeling performance was significantly better than that of linear registration or SPM2. Further advantages to Mindboggle include a high degree of robustness against image artifacts, poor image quality, and incomplete brain data. We tested the latter hypothesis by conducting all of the tests again, this time registering the atlas to an artificially lesioned version of itself, and found that Mindboggle was the only method whose performance did not degrade significantly as the lesion size increased.

Crest lines for surface segmentation and flattening

We present a method for extracting feature curves called crest lines from a triangulated surface. Then, we calculate the geodesic Voronoi diagram of crest lines to segment the surface into several regions. Afterward, barycentric surface flattening using theory from graph embeddings is implemented and, using the Geodesic Voronoi diagram, we develop a faster surface flattening algorithm.

Automated human frontal lobe identification in MR images based on fuzzy-logic encoded expert anatomic knowledge

Identification of human brain structures in MR images comprises an area of increasing interest, which also presents numerous methodological challenges. Here we describe a new knowledge-based automated method designed to identify several major brain sulci and then to define the frontal lobes by using the identified sulci as landmarks. To identify brain sulci, sulcal images were generated by morphologic operations and then separated into different components based on connectivity analysis. Subsequently, the individual anatomic features were evaluated by using fuzzy membership functions. The crisp decisions, i.e., the identification of sulci, were made by taking the maximum of the summation of all the membership functions. The identification was designed in a hierarchical order. The longitudinal fissure was extracted first. The left and right central sulci were then identified based on the left and right hemispheres. Next, the lateral sulci were identified based on the central sulci and hemispheres. Finally, the left and right frontal lobes were defined from the two hemispheres. The method was evaluated by visual inspection, comparison with manual segmentation, and comparison with manually volumetric results in references. The average Jaccard similarities of left and right frontal lobes between the automated and manual segmentation were 0.89 and 0.91, respectively. The average Kappa indices of left and right frontal lobes between the automated and manual segmentation were 0.94 and 0.95, respectively. These results show relatively high accuracy of using this novel method for human frontal lobe identification and segmentation.

Using replicator dynamics for analyzing fMRI data of the human brain

The understanding of brain networks becomes increasingly the focus of current research. In the context of functional magnetic resonance imagery (fMRI) data of the human brain, networks have been mostly detected using standard clustering approaches. In this work, we present a new method of detecting functional networks using fMRI data. The novelty of this method is that these networks have the property that every network member is closely connected with every other member. This definition might to be better suited to model important aspects of brain activity than standard cluster definitions. The algorithm that we present here is based on a concept from theoretical biology called "replicator dynamics."

Automated sulcal segmentation using watersheds on the cortical surface

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 oversegmentation, meaning that a single region is segmented as several pieces. To address this problem, we propose a postprocessing 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.

LIPSIA--a new software system for the evaluation of functional magnetic resonance images of the human brain

This paper describes the non-commercial software system LIPSIA that was developed for the processing of functional magnetic resonance images (fMRI) of the human brain. The analysis of fMRI data comprises various aspects including filtering, spatial transformation, statistical evaluation as well as segmentation and visualization. In LIPSIA, particular emphasis was placed on the development of new visualization and segmentation techniques that support visualizations of individual brain anatomy so that experts can assess the exact location of activation patterns in individual brains. As the amount of data that must be handled is enormous, another important aspect in the development LIPSIA was the efficiency of the software implementation. Well established statistical techniques were used whenever possible.

A probabilistic atlas and reference system for the human brain: International Consortium for Brain Mapping (ICBM)

Motivated by the vast amount of information that is rapidly accumulating about the human brain in digital form, we embarked upon a program in 1992 to develop a four-dimensional probabilistic atlas and reference system for the human brain. Through an International Consortium for Brain Mapping (ICBM) a dataset is being collected that includes 7000 subjects between the ages of eighteen and ninety years and including 342 mono- and dizygotic twins. Data on each subject includes detailed demographic, clinical, behavioural and imaging information. DNA has been collected for genotyping from 5800 subjects. A component of the programme uses post-mortem tissue to determine the probabilistic distribution of microscopic cyto- and chemoarchitectural regions in the human brain. This, combined with macroscopic information about structure and function derived from subjects in vivo, provides the first large scale opportunity to gain meaningful insights into the concordance or discordance in micro- and macroscopic structure and function. The philosophy, strategy, algorithm development, data acquisition techniques and validation methods are described in this report along with database structures. Examples of results are described for the normal adult human brain as well as examples in patients with Alzheimer's disease and multiple sclerosis. The ability to quantify the variance of the human brain as a function of age in a large population of subjects for whom data is also available about their genetic composition and behaviour will allow for the first assessment of cerebral genotype-phenotype-behavioural correlations in humans to take place in a population this large. This approach and its application should provide new insights and opportunities for investigators interested in basic neuroscience, clinical diagnostics and the evaluation of neuropsychiatric disorders in patients.

Building 3D sulcal models using local geometry

This paper presents a series of 3D statistical models of the cortical sulci. They are built from points located automatically over the sulcal fissures, and corresponded automatically using variants on the iterative closest point algorithm. The models are progressively improved by adding in more and more structural and configural information, and the final results are consistent with findings from other anatomical studies. The models can be used to locate and label anatomical features automatically in 3D MR images of the head, for analysis, visualisation, classification, and normalisation.

Automatic labelling of the human cortical surface using sulcal basins

Human brain mapping aims at establishing correspondences between brain function and brain anatomy. One of the most intriguing problems in this field is the high interpersonal variability of human neuroanatomy which makes studies across many subjects very difficult. The cortical folds ('sulci') often serve as landmarks that help to establish correspondences between subjects. In this paper, we will present a method that automatically detects and attributes neuroanatomical names to the cortical folds using image analysis methods applied to magnetic resonance data of human brains. We claim that the cortical folds can be subdivided into a number of substructures which we call sulcal basins. The concept of sulcal basins allows us to establish a complete parcellation of the cortical surface into separate regions. These regions are neuroanatomically meaningful and can be identified from MR data sets across many subjects. Sulcal basins are segmented using a region growing approach. The automatic labelling is achieved by a model matching technique.