Learning a Probabilistic Model for Diffeomorphic Registration

International audienceWe propose to learn a low-dimensional probabilistic deformation model from data which can be used for registration and the analysis of deformations. The latent variable model maps similar deformations close to each other in an encoding space. It enables to compare deformations, generate normal or pathological deformations for any new image or to transport deformations from one image pair to any other image. Our unsupervised method is based on variational inference. In particular, we use a conditional variational autoencoder (CVAE) network and constrain transformations to be symmetric and diffeomorphic by applying a differentiable exponentiation layer with a symmetric loss function. We also present a formulation that includes spatial regularization such as diffusion-based filters. Additionally, our framework provides multi-scale velocity field estimations. We evaluated our method on 3-D intra-subject registration using 334 cardiac cine-MRIs. On this dataset, our method showed state-of-the-art performance with a mean DICE score of 81.2% and a mean Hausdorff distance of 7.3mm using 32 latent dimensions compared to three state-of-the-art methods while also demonstrating more regular deformation fields. The average time per registration was 0.32s. Besides, we visualized the learned latent space and show that the encoded deformations can be used to transport deformations and to cluster diseases with a classification accuracy of 83% after applying a linear projection

Learning a Probabilistic Model for Diffeomorphic Registration

International audienceWe propose to learn a low-dimensional probabilistic deformation model from data which can be used for registration and the analysis of deformations. The latent variable model maps similar deformations close to each other in an encoding space. It enables to compare deformations, generate normal or pathological deformations for any new image or to transport deformations from one image pair to any other image. Our unsupervised method is based on variational inference. In particular, we use a conditional variational autoencoder (CVAE) network and constrain transformations to be symmetric and diffeomorphic by applying a differentiable exponentiation layer with a symmetric loss function. We also present a formulation that includes spatial regularization such as diffusion-based filters. Additionally, our framework provides multi-scale velocity field estimations. We evaluated our method on 3-D intra-subject registration using 334 cardiac cine-MRIs. On this dataset, our method showed state-of-the-art performance with a mean DICE score of 81.2% and a mean Hausdorff distance of 7.3mm using 32 latent dimensions compared to three state-of-the-art methods while also demonstrating more regular deformation fields. The average time per registration was 0.32s. Besides, we visualized the learned latent space and show that the encoded deformations can be used to transport deformations and to cluster diseases with a classification accuracy of 83% after applying a linear projection

Statistical analysis of organs' shapes and deformations: the Riemannian and the affine settings in computational anatomy

International audienceComputational anatomy is an emerging discipline at the interface of geometry, statistics and medicine that aims at analyzing and modeling the biological variability of organs' shapes at the population level. Shapes are equivalence classes of images, surfaces or deformations of a template under rigid body (or more general) transformations. Thus, they belong to non-linear manifolds. In order to deal with multiple samples in non-linear spaces, a consistent statistical framework on Riemannian manifolds has been designed over the last decade. We detail in this chapter the extension of this framework to Lie groups endowed with the affine symmetric connection, a more invariant (and thus more consistent) but non-metric structure on transformation groups. This theory provides strong theoretical bases for the use of one-parameter subgroups and diffeomorphisms parametrized by stationary velocity fields (SVF), for which efficient image registration methods like log-Demons have been developed with a great success from the practical point of view. One can further reduce the complexity with locally affine transformations , leading to parametric diffeomorphisms of low dimension encoding the major shape variability. We illustrate the methodology with the modeling of the evolution of the brain with Alzheimer's disease and the analysis of the cardiac motion from MRI sequences of images

Parallel Transport of Surface Deformations from Pole Ladder to Symmetrical Extension

International audienceCardiac motion contains information underlying disease development , and complements the anatomical information extracted for each subject. However, normalization of temporal trajectories is necessary due to anatomical differences between subjects. In this study, we encode inter-subject shape variations and temporal deformations in a common space of diffeomorphic registration. They are parameterized by stationary velocity fields. Previous normalization algorithms applied in medical imaging were first order approximations of parallel transport. In contrast, pole ladder was recently shown to be a third order scheme in general affine connection spaces and exact in one step in affine symmetric spaces. We further improve this procedure with a more symmetric mapping scheme, which relies on geodesic symmetries around mid-points. We apply the method to analyze cardiac motion among pulmonary hyperten-sion populations. Evaluation is performed on a 4D cardiac database, with meshes of the right-ventricle obtained by commercial speckle-tracking from echo-cardiogram. We assess the stability of the algorithms by computing their numerical inverse error. Our method turns out to be very accurate and efficient in terms of compactness for subspace representation