Abstract:We propose numerical algorithms for solving large deformation diffeomorphic image registration problems. We formulate the nonrigid image registration problem as a problem of optimal control. This leads to an infinite-dimensional partial differential equation (PDE) constrained optimization problem. The PDE constraint consists, in its simplest form, of a hyperbolic transport equation for the evolution of the image intensity. The control variable is the velocity field. Tikhonov regularization on the control ensures well-posedness. We consider standard smoothness regularization based on H1- or H2-seminorms. We augment this regularization scheme with a constraint on the divergence of the velocity field (control variable) rendering the deformation incompressible (Stokes regularization scheme) and thus ensuring that the determinant of the deformation gradient is equal to one, up to the numerical error. We use a Fourier pseudospectral discretization in space and a Chebyshev pseudospectral discretization in time. The latter allows us to reduce the number of unknowns and enables the time-adaptive inversion for nonstationary velocity fields. We use a preconditioned, globalized, matrix-free, inexact Newton-Krylov method for numerical optimization. A parameter continuation is designed to estimate an optimal regularization parameter. Regularity is ensured by controlling the geometric properties of the deformation field. Overall, we arrive at a black-box solver that exploits computational tools that are precisely tailored for solving the optimality system. We study spectral properties of the Hessian, grid convergence, numerical accuracy, computational efficiency, and deformation regularity of our scheme. We compare the designed Newton-Krylov methods with a globalized Picard method (preconditioned gradient descent). We study the influence of a varying number of unknowns in time. The reported results demonstrate excellent numerical accuracy, guaranteed local deformation regularity, and computational efficiency with an optional control on local mass conservation. The Newton-Krylov methods clearly outperform the Picard method if high accuracy of the inversion is required. Our method provides equally good results for stationary and nonstationary velocity fields for two-image registration problems.

On Finite Difference Jacobian Computation in Deformable Image Registration

New method of averaging diffeomorphisms based on Jacobian determinant and curl vector

Diffeomorphic Image Registration with an Optimal Control Relaxation and Its Implementation

A Bi-variant Variational Model for Diffeomorphic Image Registration with Relaxed Jacobian Determinant Constraints

A Novel Diffeomorphic Model for Image Registration and Its Algorithm

Towards Positive Jacobian: Learn to Postprocess Diffeomorphic Image Registration with Matrix Exponential

An Inexact Newton-Krylov Algorithm for Constrained Diffeomorphic Image Registration

Deformable Image Registration Using Functions of Bounded Deformation

An Iterative SIFT Based on Intensity and Spatial Information for Remote Sensing Image Registration

Frequency Diffeomorphisms for Efficient Image Registration

On Reducing Negative Jacobian Determinant of the Deformation Predicted by Deep Registration Networks

Distributed-memory large deformation diffeomorphic 3D image registration

Adaptive Diffeomorphic Multiresolution Demons and Their Application to Same Modality Medical Image Registration with Large Deformation

Explicit–implicit symmetric diffeomorphic deformable image registration with convolutional neural network

A Hybrid Scheme for Efficient Deformable Image Registration.

Landmark and Intensity Based Registration with Large Deformations via Quasi-conformal Maps

A Novel Unsupervised Learning Model for Diffeomorphic Image Registration

Applications of Automatic Differentiation in Image Registration

A Deep Learning Framework for Diffeomorphic Mapping Problems via Quasi-conformal Geometry Applied to Imaging

Diff-Reg v1: Diffusion Matching Model for Registration Problem