Doctor of Philosophy in Computing

dissertationAn important area of medical imaging research is studying anatomical diffeomorphic shape changes and detecting their relationship to disease processes. For example, neurodegenerative disorders change the shape of the brain, thus identifying differences between the healthy control subjects and patients affected by these diseases can help with understanding the disease processes. Previous research proposed a variety of mathematical approaches for statistical analysis of geometrical brain structure in three-dimensional (3D) medical imaging, including atlas building, brain variability quantification, regression, etc. The critical component in these statistical models is that the geometrical structure is represented by transformations rather than the actual image data. Despite the fact that such statistical models effectively provide a way for analyzing shape variation, none of them have a truly probabilistic interpretation. This dissertation contributes a novel Bayesian framework of statistical shape analysis for generic manifold data and its application to shape variability and brain magnetic resonance imaging (MRI). After we carefully define the distributions on manifolds, we then build Bayesian models for analyzing the intrinsic variability of manifold data, involving the mean point, principal modes, and parameter estimation. Because there is no closed-form solution for Bayesian inference of these models on manifolds, we develop a Markov Chain Monte Carlo method to sample the hidden variables from the distribution. The main advantages of these Bayesian approaches are that they provide parameter estimation and automatic dimensionality reduction for analyzing generic manifold-valued data, such as diffeomorphisms. Modeling the mean point of a group of images in a Bayesian manner allows for learning the regularity parameter from data directly rather than having to set it manually, which eliminates the effort of cross validation for parameter selection. In population studies, our Bayesian model of principal modes analysis (1) automatically extracts a low-dimensional, second-order statistics of manifold data variability and (2) gives a better geometric data fit than nonprobabilistic models. To make this Bayesian framework computationally more efficient for high-dimensional diffeomorphisms, this dissertation presents an algorithm, FLASH (finite-dimensional Lie algebras for shooting), that hugely speeds up the diffeomorphic image registration. Instead of formulating diffeomorphisms in a continuous variational problem, Flash defines a completely new discrete reparameterization of diffeomorphisms in a low-dimensional bandlimited velocity space, which results in the Bayesian inference via sampling on the space of diffeomorphisms being more feasible in time. Our entire Bayesian framework in this dissertation is used for statistical analysis of shape data and brain MRIs. It has the potential to improve hypothesis testing, classification, and mixture models

Elastic shape matching of parameterized surfaces using square root normal fields.

In this paper we define a new methodology for shape analysis of parameterized surfaces, where the main issues are: (1) choice of metric for shape comparisons and (2) invariance to reparameterization. We begin by defining a general elastic metric on the space of parameterized surfaces. The main advantages of this metric are twofold. First, it provides a natural interpretation of elastic shape deformations that are being quantified. Second, this metric is invariant under the action of the reparameterization group. We also introduce a novel representation of surfaces termed square root normal fields or SRNFs. This representation is convenient for shape analysis because, under this representation, a reduced version of the general elastic metric becomes the simple \ensuremathL2\ensuremathL2 metric. Thus, this transformation greatly simplifies the implementation of our framework. We validate our approach using multiple shape analysis examples for quadrilateral and spherical surfaces. We also compare the current results with those of Kurtek et al. [1]. We show that the proposed method results in more natural shape matchings, and furthermore, has some theoretical advantages over previous methods

Doctor of Philosophy

dissertationAn important aspect of medical research is the understanding of anatomy and its relation to function in the human body. For instance, identifying changes in the brain associated with cognitive decline helps in understanding the process of aging and age-related neurological disorders. The field of computational anatomy provides a rich mathematical setting for statistical analysis of complex geometrical structures seen in 3D medical images. At its core, computational anatomy is based on the representation of anatomical shape and its variability as elements of nonflat manifold of diffeomorphisms with an associated Riemannian structure. Although such manifolds effectively represent natural biological variability, intrinsic methods of statistical analysis within these spaces remain deficient at large. This dissertation contributes two critical missing pieces for statistics in diffeomorphisms: (1) multivariate regression models for cross-sectional study of shapes, and (2) generalization of classical Euclidean, mixed-effects models to manifolds for longitudinal studies. These models are based on the principle that statistics on manifold-valued information must respect the intrinsic geometry of that space. The multivariate regression methods provide statistical descriptors of the relationships of anatomy with clinical indicators. The novel theory of hierarchical geodesic models (HGMs) is developed as a natural generalization of hierarchical linear models (HLMs) to describe longitudinal data on curved manifolds. Using a hierarchy of geodesics, the HGMs address the challenge of modeling the shape-data with unbalanced designs typically arising as a result of follow-up medical studies. More generally, this research establishes a mathematical foundation to study dynamics of changes in anatomy and the associated clinical progression with time. This dissertation also provides efficient algorithms that utilize state-of-the-art high performance computing architectures to solve models on large-scale, longitudinal imaging data. These manifold-based methods are applied to predictive modeling of neurological disorders such as Alzheimer's disease. Overall, this dissertation enables clinicians and researchers to better utilize the structural information available in medical images

Singular geodesic coordinates for representing diffeomorphic maps in computational anatomy, with application to the morphometry of early Alzheimer's disease in the medial temporal lobe

In this work we develop novel algorithms for building one to one correspondences between anatomical forms by providing a sparse representation of dense registration information. These sparse parameterizations of complex high dimensional data allow robustness in the face of noise and anomalies, and a platform for inference that is effective in the face of multiple comparisons. We review background in the theory of generating smooth, invertible transformations (the diffeomorphism group), and build our parameterization as a function supported on surfaces bounding anatomical structures of interest. We show how dimensionality can be reduced even further and still provide a rich family of mappings using principal component analysis or Laplace Beltrami eigenfunctions supported on the surface. We develop algorithms for surface matching and image matching within this model, and demonstrate the desired robustness by working with published large neuroimaging datasets that include many low quality examples. Finally we turn to addressing challenges associated with some specific data types: images with multiple labels, and longitudinal data. We use the mapping tools developed to draw conclusions about the progression of early Alzheimer's disease in the medial temporal lobe

Visual analytics methods for shape analysis of biomedical images exemplified on rodent skull morphology

In morphometrics and its application fields like medicine and biology experts are interested in causal relations of variation in organismic shape to phylogenetic, ecological, geographical, epidemiological or disease factors - or put more succinctly by Fred L. Bookstein, morphometrics is "the study of covariances of biological form". In order to reveal causes for shape variability, targeted statistical analysis correlating shape features against external and internal factors is necessary but due to the complexity of the problem often not feasible in an automated way. Therefore, a visual analytics approach is proposed in this thesis that couples interactive visualizations with automated statistical analyses in order to stimulate generation and qualitative assessment of hypotheses on relevant shape features and their potentially affecting factors. To this end long established morphometric techniques are combined with recent shape modeling approaches from geometry processing and medical imaging, leading to novel visual analytics methods for shape analysis. When used in concert these methods facilitate targeted analysis of characteristic shape differences between groups, co-variation between different structures on the same anatomy and correlation of shape to extrinsic attributes. Here a special focus is put on accurate modeling and interactive rendering of image deformations at high spatial resolution, because that allows for faithful representation and communication of diminutive shape features, large shape differences and volumetric structures. The utility of the presented methods is demonstrated in case studies conducted together with a collaborating morphometrics expert. As exemplary model structure serves the rodent skull and its mandible that are assessed via computed tomography scans

Geodesic shape regression with multiple geometries and sparse parameters

International audienceMany problems in medicine are inherently dynamic processes which include the aspect of change over time, such as childhood development, aging, and disease progression. From medical images, numerous geometric structures can be extracted with various representations, such as landmarks, point clouds, curves, and surfaces. Different sources of geometry may characterize different aspects of the anatomy, such as fiber tracts from DTI and subcortical shapes from structural MRI, and therefore require a modeling scheme which can include various shape representations in any combination. In this paper, we present a geodesic regression model in the large deformation (LDDMM) framework applicable to multi-object complexes in a variety of shape representations. Our model decouples the deformation parameters from the specific shape representations, allowing the complexity of the model to reflect the nature of the shape changes, rather than the sampling of the data. As a consequence, the sparse representation of diffeomorphic flow allows for the straightforward embedding of a variety of geometry in different combinations, which all contribute towards the estimation of a single deformation of the ambient space. Additionally, the sparse representation along with the geodesic constraint results in a compact statistical model of shape change by a small number of parameters defined by the user. Experimental validation on multi-object complexes demonstrate robust model estimation across a variety of parameter settings. We further demonstrate the utility of our method to support the analysis of derived shape features, such as volume, and explore shape model extrapolation. Our method is freely available in the software package deformetrica which can be downloaded at www.deformetrica.org