Communications in Mathematical Sciences

Volume 13 (2015)

Number 5

Topology preservation for image-registration-related deformation fields

Pages: 1135 – 1161

DOI: https://dx.doi.org/10.4310/CMS.2015.v13.n5.a4

Authors

Solène Ozeré (Laboratoire de Mathématiques de l’INSA de Rouen, Normandie Université, Saint-Etienne-du-Rouvray, France)

Carole Le Guyader (Laboratoire de Mathématiques de l’INSA de Rouen, Normandie Université, Saint-Etienne-du-Rouvray, France)

Abstract

In this paper, we address the issue of designing a theoretically well-motivated and computationally efficient method ensuring topology preservation on image-registration-related deformation fields. The model is motivated by a mathematical characterization of topology preservation for a deformation field mapping two subsets of $\mathbb{Z}^2$, namely, positivity of the four approximations to the Jacobian determinant of the deformation on a square patch. The first step of the proposed algorithm thus consists in correcting the gradient vector field of the deformation (that does not comply with the topology preservation criteria) at the discrete level in order to fulfill this positivity condition. Once this step is achieved, it thus remains to reconstruct the deformation field, given its full set of discrete gradient vectors. We propose to decompose the reconstruction problem into independent problems of smaller dimensions, yielding a natural parallelization of the computations and enabling us to reduce drastically the computational time (up to 80 in some applications). For each subdomain, a functional minimization problem under Lagrange interpolation constraints is introduced and its well-posedness is studied: existence/uniqueness of the solution, characterization of the solution, convergence of the method when the number of data increases to infinity, discretization with the Finite Element Method and discussion on the properties of the matrix involved in the linear system. Numerical simulations based on OpenMP parallelization and MKL multi-threading demonstrating the ability of the model to handle large deformations (contrary to classical methods) and the interest of having decomposed the problem into smaller ones are provided.

Keywords

$D^m$-splines, constrained optimization, variational formulation, convergence, finite element method, image registration

2010 Mathematics Subject Classification

41A15, 65D05, 65D07, 65D10, 65K10, 65Y05, 68U10

Published 22 April 2015