23,659 research outputs found
Soliton Dynamics in Computational Anatomy
Computational anatomy (CA) has introduced the idea of anatomical structures
being transformed by geodesic deformations on groups of diffeomorphisms. Among
these geometric structures, landmarks and image outlines in CA are shown to be
singular solutions of a partial differential equation that is called the
geodesic EPDiff equation. A recently discovered momentum map for singular
solutions of EPDiff yields their canonical Hamiltonian formulation, which in
turn provides a complete parameterization of the landmarks by their canonical
positions and momenta. The momentum map provides an isomorphism between
landmarks (and outlines) for images and singular soliton solutions of the
EPDiff equation. This isomorphism suggests a new dynamical paradigm for CA, as
well as new data representation.Comment: published in NeuroImag
Multiple Shape Registration using Constrained Optimal Control
Lagrangian particle formulations of the large deformation diffeomorphic
metric mapping algorithm (LDDMM) only allow for the study of a single shape. In
this paper, we introduce and discuss both a theoretical and practical setting
for the simultaneous study of multiple shapes that are either stitched to one
another or slide along a submanifold. The method is described within the
optimal control formalism, and optimality conditions are given, together with
the equations that are needed to implement augmented Lagrangian methods.
Experimental results are provided for stitched and sliding surfaces
Weighted Diffeomorphic Density Matching with Applications to Thoracic Image Registration
In this article we study the problem of thoracic image registration, in
particular the estimation of complex anatomical deformations associated with
the breathing cycle. Using the intimate link between the Riemannian geometry of
the space of diffeomorphisms and the space of densities, we develop an image
registration framework that incorporates both the fundamental law of
conservation of mass as well as spatially varying tissue compressibility
properties. By exploiting the geometrical structure, the resulting algorithm is
computationally efficient, yet widely general.Comment: Accepted in Proceedings of the 5th MICCAI workshop on Mathematical
Foundations of Computational Anatomy, Munich, Germany, 2015
(http://www-sop.inria.fr/asclepios/events/MFCA15/
Bayesian Estimation of White Matter Atlas from High Angular Resolution Diffusion Imaging
We present a Bayesian probabilistic model to estimate the brain white matter
atlas from high angular resolution diffusion imaging (HARDI) data. This model
incorporates a shape prior of the white matter anatomy and the likelihood of
individual observed HARDI datasets. We first assume that the atlas is generated
from a known hyperatlas through a flow of diffeomorphisms and its shape prior
can be constructed based on the framework of large deformation diffeomorphic
metric mapping (LDDMM). LDDMM characterizes a nonlinear diffeomorphic shape
space in a linear space of initial momentum uniquely determining diffeomorphic
geodesic flows from the hyperatlas. Therefore, the shape prior of the HARDI
atlas can be modeled using a centered Gaussian random field (GRF) model of the
initial momentum. In order to construct the likelihood of observed HARDI
datasets, it is necessary to study the diffeomorphic transformation of
individual observations relative to the atlas and the probabilistic
distribution of orientation distribution functions (ODFs). To this end, we
construct the likelihood related to the transformation using the same
construction as discussed for the shape prior of the atlas. The probabilistic
distribution of ODFs is then constructed based on the ODF Riemannian manifold.
We assume that the observed ODFs are generated by an exponential map of random
tangent vectors at the deformed atlas ODF. Hence, the likelihood of the ODFs
can be modeled using a GRF of their tangent vectors in the ODF Riemannian
manifold. We solve for the maximum a posteriori using the
Expectation-Maximization algorithm and derive the corresponding update
equations. Finally, we illustrate the HARDI atlas constructed based on a
Chinese aging cohort of 94 adults and compare it with that generated by
averaging the coefficients of spherical harmonics of the ODF across subjects
Invariant higher-order variational problems II
Motivated by applications in computational anatomy, we consider a
second-order problem in the calculus of variations on object manifolds that are
acted upon by Lie groups of smooth invertible transformations. This problem
leads to solution curves known as Riemannian cubics on object manifolds that
are endowed with normal metrics. The prime examples of such object manifolds
are the symmetric spaces. We characterize the class of cubics on object
manifolds that can be lifted horizontally to cubics on the group of
transformations. Conversely, we show that certain types of non-horizontal
geodesics on the group of transformations project to cubics. Finally, we apply
second-order Lagrange--Poincar\'e reduction to the problem of Riemannian cubics
on the group of transformations. This leads to a reduced form of the equations
that reveals the obstruction for the projection of a cubic on a transformation
group to again be a cubic on its object manifold.Comment: 40 pages, 1 figure. First version -- comments welcome
- …