Computerized Analysis of Magnetic Resonance Images to Study Cerebral Anatomy in Developing Neonates

The study of cerebral anatomy in developing neonates is of great importance for the understanding of brain development during the early period of life. This dissertation therefore focuses on three challenges in the modelling of cerebral anatomy in neonates during brain development. The methods that have been developed all use Magnetic Resonance Images (MRI) as source data. To facilitate study of vascular development in the neonatal period, a set of image analysis algorithms are developed to automatically extract and model cerebral vessel trees. The whole process consists of cerebral vessel tracking from automatically placed seed points, vessel tree generation, and vasculature registration and matching. These algorithms have been tested on clinical Time-of- Flight (TOF) MR angiographic datasets. To facilitate study of the neonatal cortex a complete cerebral cortex segmentation and reconstruction pipeline has been developed. Segmentation of the neonatal cortex is not effectively done by existing algorithms designed for the adult brain because the contrast between grey and white matter is reversed. This causes pixels containing tissue mixtures to be incorrectly labelled by conventional methods. The neonatal cortical segmentation method that has been developed is based on a novel expectation-maximization (EM) method with explicit correction for mislabelled partial volume voxels. Based on the resulting cortical segmentation, an implicit surface evolution technique is adopted for the reconstruction of the cortex in neonates. The performance of the method is investigated by performing a detailed landmark study. To facilitate study of cortical development, a cortical surface registration algorithm for aligning the cortical surface is developed. The method first inflates extracted cortical surfaces and then performs a non-rigid surface registration using free-form deformations (FFDs) to remove residual alignment. Validation experiments using data labelled by an expert observer demonstrate that the method can capture local changes and follow the growth of specific sulcus

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 cuts 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

Cortical Surface Registration and Shape Analysis

A population analysis of human cortical morphometry is critical for insights into brain development or degeneration. Such an analysis allows for investigating sulcal and gyral folding patterns. In general, such a population analysis requires both a well-established cortical correspondence and a well-defined quantification of the cortical morphometry. The highly folded and convoluted structures render a reliable and consistent population analysis challenging. Three key challenges have been identified for such an analysis: 1) consistent sulcal landmark extraction from the cortical surface to guide better cortical correspondence, 2) a correspondence establishment for a reliable and stable population analysis, and 3) quantification of the cortical folding in a more reliable and biologically meaningful fashion. The main focus of this dissertation is to develop a fully automatic pipeline that supports a population analysis of local cortical folding changes. My proposed pipeline consists of three novel components I developed to overcome the challenges in the population analysis: 1) automatic sulcal curve extraction for stable/reliable anatomical landmark selection, 2) group-wise registration for establishing cortical shape correspondence across a population with no template selection bias, and 3) quantification of local cortical folding using a novel cortical-shape-adaptive kernel. To evaluate my methodological contributions, I applied all of them in an application to early postnatal brain development. I studied the human cortical morphological development using the proposed quantification of local cortical folding from neonate age to 1 year and 2 years of age, with quantitative developmental assessments. This study revealed a novel pattern of associations between the cortical gyrification and cognitive development.Doctor of Philosoph

Multi-Manifold Diffeomorphic Metric Mapping for Cortical Registration with its Applications in Brain Structural and Functional Studies

Ph.DDOCTOR OF PHILOSOPH

Connectomics across development:towards mapping brain structure from birth to childhood

The brain is probably the most complex system of the human body, composed of numerous neural units interconnected at dierent scales. This highly structured architecture provides the ability to communicate, synthesize information and perform the analytical tasks of human beings. Its development starts during the transition between the embryonic and fetal periods, from a simple tubular to a highly complex folded structure. It is globally organized as early as birth. This developing process is highly vulnerable to antenatal adverse conditions. Indeed, extreme prematurity and intra uterine growth restriction are major risk factors for long-term morbidities, including developmental ailments such as cerebral palsy, mental retardation and a wide spectrum of learning disabilities and behavior disorders. In this context, the characterization of the brainâs normative wiring pattern is crucial for our understanding of its architecture and workings, as the origin of many neurological and neurobehavioral disorders is found in early structural brain development. Diusion magnetic resonance imaging (dMRI) allows the in vivo assessment of biological tissues at the microstructural level. It has emerged as a powerful tool to study brain connectivity and analyse the underlying substrate of the human brain, comprising its structurally integrated and functionally specialized architecture. dMRI has been widely used in adult studies. Nevertheless, due to technical constraints, this mapping at earlier stages of development has not yet been accomplished. Yet, this time period is of extreme importance to comprehend the structural and functional integrity of the brain. This thesis is motivated by this shortfall, and intends to fill the gap between the clinical and neuroscience demands and the methodological developments needed to fulfill them. In our work, we comprehensibly study the brain structural connectivity of children born extremely prematurely and/or with additional prenatal restriction at school-age. We provide evidence that brain systems that mature early in development are the most vulnerable to antenatal insults. Interestingly, the alterations highlighted in these systems correlate with the neurobehavioral and cognitive impairments seen in these children at school-age. The overall brain organization appear also altered after preterm birth and prenatal restriction. Indeed, these children show dierent brain network modular topology, with a reduction in the overall network capacity. What remains unclear is whether the alterations seen at school age are already present at birth and, if yes, to what extent. In this thesis we set the technical basis to enable the connectome analysis as early as at birth. This task is challenging when dealing with neonatal data. Indeed, most of the assumptions used in adult data processing methods do not hold, due to the inverted image contrast and other MRI artefacts such as motion, partial volume and intensity inhomogeneities. Here, we propose a novel technique for surface reconstruction, and provide a fully-automatic procedure to delineate the newborn cortical surface, opening the way to establish the newborn connectome

Cortical structure: Linking MRI and cytoarchitecture

MRI provides a powerful tool to investigate brain structure in living humans. However a major challenge is interpreting the biological underpinnings of changes at this scale. This dissertation describes investigations into the problem of linking microscale post mortem cortical cytoarchitecture with millimeter-scale measures of cortical anatomy accessible through in vivo MRI. Chapter 1 introduces the problem and previous work done to address it. The following two chapters apply classical atlases of cortical cytoarchitecture to understanding morphological changes both in health (Chapter 2) and in disease (Chapter 3). Chapter 2 demonstrates that sensory processing hierarchies exhibit increasing gradients of cortical thickness, related to changes in cortical cytoarchitecture. In Chapter 3, cytoarchitectonically described differences in gyral and sulcal laminar structure were used to create markers of laminar change from MRI changes in schizophrenia. Classical measurements of histology have limitations; they are observer dependent, two-dimensional with limited coverage of the cortex. To address these issues, Chapters 4-6 document work carried on BigBrain, a 3D 20$\mu$m resolution histological dataset. I created a high-resolution 3D atlas of laminar cytoarchitecture, which was mapped to MRI-compatible cortical surface reconstructions. Chapter 4 records the development of an automated 1D profile-based approach to laminar analysis, revealing basic principles of cortical cytoarchitecture. In Chapter 5 this approach was extended to identify 6 cortical layers throughout the isocortex. These tools can be used to segment 1D cortical intensity profiles derived from any modality. In Chapter 6, the analysis of cortical gradients initially identified using MRI cortical thickness in Chapter 2 was replicated and extended using novel histological data. First histological cortical thicknesses were tested for the same patterns organization measured on in vivo MRI in Chapter 2. These analyses were extended to test which layers contributed most to overall thickness. High-resolution, complete maps of cortical cytoarchitecture mapped to MRI-template cortical surface reconstructions, are a powerful tool and dataset for the neuroimaging community. They offer new possibilities for linking cortical microstructure to in vivo neuroimaging

3D Normal Coordinate Systems for Cortical Areas

A surface-based diffeomorphic algorithm to generate 3D coordinate grids in the cortical ribbon is described. In the grid, normal coordinate lines are generated by the diffeomorphic evolution from the grey/white (inner) surface to the grey/csf (outer) surface. Specifically, the cortical ribbon is described by two triangulated surfaces with open boundaries. Conceptually, the inner surface sits on top of the white matter structure and the outer on top of the gray matter. It is assumed that the cortical ribbon consists of cortical columns which are orthogonal to the white matter surface. This might be viewed as a consequence of the development of the columns in the embryo. It is also assumed that the columns are orthogonal to the outer surface so that the resultant vector field is orthogonal to the evolving surface. Then the distance of the normal lines from the vector field such that the inner surface evolves diffeomorphically towards the outer one can be construed as a measure of thickness. Applications are described for the auditory cortices in human adults and cats with normal hearing or hearing loss. The approach offers great potential for cortical morphometry

Geodesic Active Fields:A Geometric Framework for Image Registration

Image registration is the concept of mapping homologous points in a pair of images. In other words, one is looking for an underlying deformation field that matches one image to a target image. The spectrum of applications of image registration is extremely large: It ranges from bio-medical imaging and computer vision, to remote sensing or geographic information systems, and even involves consumer electronics. Mathematically, image registration is an inverse problem that is ill-posed, which means that the exact solution might not exist or not be unique. In order to render the problem tractable, it is usual to write the problem as an energy minimization, and to introduce additional regularity constraints on the unknown data. In the case of image registration, one often minimizes an image mismatch energy, and adds an additive penalty on the deformation field regularity as smoothness prior. Here, we focus on the registration of the human cerebral cortex. Precise cortical registration is required, for example, in statistical group studies in functional MR imaging, or in the analysis of brain connectivity. In particular, we work with spherical inflations of the extracted hemispherical surface and associated features, such as cortical mean curvature. Spatial mapping between cortical surfaces can then be achieved by registering the respective spherical feature maps. Despite the simplified spherical geometry, inter-subject registration remains a challenging task, mainly due to the complexity and inter-subject variability of the involved brain structures. In this thesis, we therefore present a registration scheme, which takes the peculiarities of the spherical feature maps into particular consideration. First, we realize that we need an appropriate hierarchical representation, so as to coarsely align based on the important structures with greater inter-subject stability, before taking smaller and more variable details into account. Based on arguments from brain morphogenesis, we propose an anisotropic scale-space of mean-curvature maps, built around the Beltrami framework. Second, inspired by concepts from vision-related elements of psycho-physical Gestalt theory, we hypothesize that anisotropic Beltrami regularization better suits the requirements of image registration regularization, compared to traditional Gaussian filtering. Different objects in an image should be allowed to move separately, and regularization should be limited to within the individual Gestalts. We render the regularization feature-preserving by limiting diffusion across edges in the deformation field, which is in clear contrast to the indifferent linear smoothing. We do so by embedding the deformation field as a manifold in higher-dimensional space, and minimize the associated Beltrami energy which represents the hyperarea of this embedded manifold as measure of deformation field regularity. Further, instead of simply adding this regularity penalty to the image mismatch in lieu of the standard penalty, we propose to incorporate the local image mismatch as weighting function into the Beltrami energy. The image registration problem is thus reformulated as a weighted minimal surface problem. This approach has several appealing aspects, including (1) invariance to re-parametrization and ability to work with images defined on non-flat, Riemannian domains (e.g., curved surfaces, scalespaces), and (2) intrinsic modulation of the local regularization strength as a function of the local image mismatch and/or noise level. On a side note, we show that the proposed scheme can easily keep up with recent trends in image registration towards using diffeomorphic and inverse consistent deformation models. The proposed registration scheme, called Geodesic Active Fields (GAF), is non-linear and non-convex. Therefore we propose an efficient optimization scheme, based on splitting. Data-mismatch and deformation field regularity are optimized over two different deformation fields, which are constrained to be equal. The constraint is addressed using an augmented Lagrangian scheme, and the resulting optimization problem is solved efficiently using alternate minimization of simpler sub-problems. In particular, we show that the proposed method can easily compete with state-of-the-art registration methods, such as Demons. Finally, we provide an implementation of the fast GAF method on the sphere, so as to register the triangulated cortical feature maps. We build an automatic parcellation algorithm for the human cerebral cortex, which combines the delineations available on a set of atlas brains in a Bayesian approach, so as to automatically delineate the corresponding regions on a subject brain given its feature map. In a leave-one-out cross-validation study on 39 brain surfaces with 35 manually delineated gyral regions, we show that the pairwise subject-atlas registration with the proposed spherical registration scheme significantly improves the individual alignment of cortical labels between subject and atlas brains, and, consequently, that the estimated automatic parcellations after label fusion are of better quality

Shape analysis of the human brain.

Autism is a complex developmental disability that has dramatically increased in prevalence, having a decisive impact on the health and behavior of children. Methods used to detect and recommend therapies have been much debated in the medical community because of the subjective nature of diagnosing autism. In order to provide an alternative method for understanding autism, the current work has developed a 3-dimensional state-of-the-art shape based analysis of the human brain to aid in creating more accurate diagnostic assessments and guided risk analyses for individuals with neurological conditions, such as autism. Methods: The aim of this work was to assess whether the shape of the human brain can be used as a reliable source of information for determining whether an individual will be diagnosed with autism. The study was conducted using multi-center databases of magnetic resonance images of the human brain. The subjects in the databases were analyzed using a series of algorithms consisting of bias correction, skull stripping, multi-label brain segmentation, 3-dimensional mesh construction, spherical harmonic decomposition, registration, and classification. The software algorithms were developed as an original contribution of this dissertation in collaboration with the BioImaging Laboratory at the University of Louisville Speed School of Engineering. The classification of each subject was used to construct diagnoses and therapeutic risk assessments for each patient. Results: A reliable metric for making neurological diagnoses and constructing therapeutic risk assessment for individuals has been identified. The metric was explored in populations of individuals having autism spectrum disorders, dyslexia, Alzheimers disease, and lung cancer. Conclusion: Currently, the clinical applicability and benefits of the proposed software approach are being discussed by the broader community of doctors, therapists, and parents for use in improving current methods by which autism spectrum disorders are diagnosed and understood