Large Deformation Diffeomorphic Metric Mapping Framework Research Articles

String Methods for Stochastic Image and Shape Matching

Matching of images and analysis of shape differences is traditionally pursued by energy minimization of paths of deformations acting to match the shape objects. In the Large Deformation Diffeomorphic Metric Mapping (LDDMM) framework, iterative gradient descents on the matching functional lead to matching algorithms informally known as Beg algorithms. When stochasticity is introduced to model stochastic variability of shapes and to provide more realistic models of observed shape data, the corresponding matching problem can be solved with a stochastic Beg algorithm, similar to the finite temperature string method used in rare event sampling. In this paper, we apply a stochastic model compatible with the geometry of the LDDMM framework to obtain a stochastic model of images and we derive the stochastic version of the Beg algorithm which we compare with the string method and an expectation-maximization optimization of posterior likelihoods. The algorithm and its use for statistical inference is tested on stochastic LDDMM landmarks and images.

Multiscale Frame-Based Kernels for Large Deformation Diffeomorphic Metric Mapping

We present a set of multiscale frame-based kernels that can be used to construct diffeomorphic transformation in the large deformation diffeomorphic metric mapping (LDDMM) framework. We construct multiscale kernels via compact wavelet frames that are equipped with the hierarchical multiresolution analysis. We show that these kernels under certain conditions can form reproducing kernel Hilbert spaces of smooth velocity fields and hence can be used to generate multiscale diffeomorphic transformation for LDDMM. As a proof of concept, we incorporate these kernels in the LDDMM framework. We show the improvement of whole brain mapping accuracy using the LDDMM with frame-based kernels in comparison to that obtained using the LDDMM with Gaussian kernels. Moreover, we evaluate whole brain mapping accuracy of the LDDMM with frame-based kernels against that obtained from the 14 brain mapping methods given by Klein et al.. Our results suggest that the LDDMM with frame-based kernels has the potential to outperform the 14 brain mapping methods for whole brain mapping.

Diffeomorphic Surface Registration with Atrophy Constraints.

Diffeomorphic registration using optimal control on the diffeomorphism group and on shape spaces has become widely used since the development of the large deformation diffeomorphic metric mapping (LDDMM) algorithm. More recently, a series of algorithms involving sub-Riemannian constraints have been introduced in which the velocity fields that control the shapes in the LDDMM framework are constrained in accordance with a specific deformation model. Here, we extend this setting by considering, for the first time, inequality constraints in order to estimate surface deformations that only allow for atrophy, introducing for this purpose an algorithm that uses the augmented Lagrangian method. We prove the existence of solutions of the associated optimal control problem and the consistency of our approximation scheme. These developments are illustrated by numerical experiments on simulated and real data.

Computing Diffeomorphic Paths for Large Motion Interpolation.

In this paper, we introduce a novel framework for computing a path of diffeomorphisms between a pair of input diffeomorphisms. Direct computation of a geodesic path on the space of diffeomorphisms Diff(Ω) is difficult, and it can be attributed mainly to the infinite dimensionality of Diff(Ω). Our proposed framework, to some degree, bypasses this difficulty using the quotient map of Diff(Ω) to the quotient space Diff(M)/Diff(M) μ obtained by quotienting out the subgroup of volume-preserving diffeomorphisms Diff(M) μ . This quotient space was recently identified as the unit sphere in a Hilbert space in mathematics literature, a space with well-known geometric properties. Our framework leverages this recent result by computing the diffeomorphic path in two stages. First, we project the given diffeomorphism pair onto this sphere and then compute the geodesic path between these projected points. Second, we lift the geodesic on the sphere back to the space of diffeomerphisms, by solving a quadratic programming problem with bilinear constraints using the augmented Lagrangian technique with penalty terms. In this way, we can estimate the path of diffeomorphisms, first, staying in the space of diffeomorphisms, and second, preserving shapes/volumes in the deformed images along the path as much as possible. We have applied our framework to interpolate intermediate frames of frame-sub-sampled video sequences. In the reported experiments, our approach compares favorably with the popular Large Deformation Diffeomorphic Metric Mapping framework (LDDMM).

A large deformation diffeomorphic metric mapping solution for diffusion spectrum imaging datasets

Spatial transformation for diffusion spectrum imaging (DSI) is an important step for group analyses of DSI datasets. In this study, we developed a transformation method for DSI datasets under the framework of large deformation diffeomorphic metric mapping (LDDMM), which is termed LDDMM-DSI. The proposed method made use of the fact that a DSI dataset is 6D, and generalized the original 2D/3D LDDMM algorithm to the 6D case with some modifications made for the DSI datasets. In this manner, the conventional reorientation problem that arises from transforming diffusion-weighted datasets was avoided by making the DSI datasets capable of being freely deformed in the q-space. The algorithm treated the data-matching task as a variational problem under the LDDMM framework and sought optimal velocity fields from which the generated transformations were diffeomorphic and the transformation curve was a geodesic. The mathematical materials and numerical implementation are detailed in the paper, and experiments were performed to analyze the proposed method on real brain DSI datasets. The results showed that the method was capable of registering different DSI datasets in both global structural shapes and local diffusion profiles. In conclusion, the proposed method can facilitate group analyses of DSI datasets and the generation of a DSI template.

Simultaneous Multi-scale Registration Using Large Deformation Diffeomorphic Metric Mapping

In the framework of large deformation diffeomorphic metric mapping (LDDMM), we present a practical methodology to integrate prior knowledge about the registered shapes in the regularizing metric. Our goal is to perform rich anatomical shape comparisons from volumetric images with the mathematical properties offered by the LDDMM framework. We first present the notion of characteristic scale at which image features are deformed. We then propose a methodology to compare anatomical shape variations in a multi-scale fashion, i.e., at several characteristic scales simultaneously. In this context, we propose a strategy to quantitatively measure the feature differences observed at each characteristic scale separately. After describing our methodology, we illustrate the performance of the method on phantom data. We then compare the ability of our method to segregate a group of subjects having Alzheimer's disease and a group of controls with a classical coarse to fine approach, on standard 3D MR longitudinal brain images. We finally apply the approach to quantify the anatomical development of the human brain from 3D MR longitudinal images of pre-term babies. Results show that our method registers accurately volumetric images containing feature differences at several scales simultaneously with smooth deformations.

Registration of Anatomical Images Using Paths of Diffeomorphisms Parameterized with Stationary Vector Field Flows

Computational Anatomy aims for the study of variability in anatomical structures from images. Variability is encoded by the spatial transformations existing between anatomical images and a template selected as reference. In the absence of a more justified model for inter-subject variability, transformations are considered to belong to a convenient family of diffeomorphisms which provides a suitable mathematical setting for the analysis of anatomical variability. One of the proposed paradigms for diffeomorphic registration is the Large Deformation Diffeomorphic Metric Mapping (LDDMM). In this framework, transformations are characterized as end points of paths parameterized by time-varying flows of vector fields defined on the tangent space of a Riemannian manifold of diffeomorphisms and computed from the solution of the non-stationary transport equation associated to these flows. With this characterization, optimization in LDDMM is performed on the space of non-stationary vector field flows resulting into a time and memory consuming algorithm. Recently, an alternative characterization of paths of diffeomorphisms based on constant-time flows of vector fields has been proposed in the literature. With this parameterization, diffeomorphisms constitute solutions of stationary ODEs. In this article, the stationary parameterization is included for diffeomorphic registration in the LDDMM framework. We formulate the variational problem related to this registration scenario and derive the associated Euler-Lagrange equations. Moreover, the performance of the non-stationary vs the stationary parameterizations in real and simulated 3D-MRI brain datasets is evaluated. Compared to the non-stationary parameterization, our proposal provides similar results in terms of image matching and local differences between the diffeomorphic transformations while drastically reducing memory and time requirements.