601 research outputs found
Diffeomorphic density matching by optimal information transport
We address the following problem: given two smooth densities on a manifold,
find an optimal diffeomorphism that transforms one density into the other. Our
framework builds on connections between the Fisher-Rao information metric on
the space of probability densities and right-invariant metrics on the
infinite-dimensional manifold of diffeomorphisms. This optimal information
transport, and modifications thereof, allows us to construct numerical
algorithms for density matching. The algorithms are inherently more efficient
than those based on optimal mass transport or diffeomorphic registration. Our
methods have applications in medical image registration, texture mapping, image
morphing, non-uniform random sampling, and mesh adaptivity. Some of these
applications are illustrated in examples.Comment: 35 page
Diffeomorphic random sampling using optimal information transport
In this article we explore an algorithm for diffeomorphic random sampling of
nonuniform probability distributions on Riemannian manifolds. The algorithm is
based on optimal information transport (OIT)---an analogue of optimal mass
transport (OMT). Our framework uses the deep geometric connections between the
Fisher-Rao metric on the space of probability densities and the right-invariant
information metric on the group of diffeomorphisms. The resulting sampling
algorithm is a promising alternative to OMT, in particular as our formulation
is semi-explicit, free of the nonlinear Monge--Ampere equation. Compared to
Markov Chain Monte Carlo methods, we expect our algorithm to stand up well when
a large number of samples from a low dimensional nonuniform distribution is
needed.Comment: 8 pages, 3 figure
Diffeomorphic density registration
In this book chapter we study the Riemannian Geometry of the density
registration problem: Given two densities (not necessarily probability
densities) defined on a smooth finite dimensional manifold find a
diffeomorphism which transforms one to the other. This problem is motivated by
the medical imaging application of tracking organ motion due to respiration in
Thoracic CT imaging where the fundamental physical property of conservation of
mass naturally leads to modeling CT attenuation as a density. We will study the
intimate link between the Riemannian metrics on the space of diffeomorphisms
and those on the space of densities. We finally develop novel computationally
efficient algorithms and demonstrate there applicability for registering RCCT
thoracic imaging.Comment: 23 pages, 6 Figures, Chapter for a Book on Medical Image Analysi
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/
Indirect Image Registration with Large Diffeomorphic Deformations
The paper adapts the large deformation diffeomorphic metric mapping framework
for image registration to the indirect setting where a template is registered
against a target that is given through indirect noisy observations. The
registration uses diffeomorphisms that transform the template through a (group)
action. These diffeomorphisms are generated by solving a flow equation that is
defined by a velocity field with certain regularity. The theoretical analysis
includes a proof that indirect image registration has solutions (existence)
that are stable and that converge as the data error tends so zero, so it
becomes a well-defined regularization method. The paper concludes with examples
of indirect image registration in 2D tomography with very sparse and/or highly
noisy data.Comment: 43 pages, 4 figures, 1 table; revise
- …