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

Abstract

Diffeomorphic registration is a widely used technique for finding a smooth and invertible transformation between two coordinate systems, which are measured using template and reference images. The point-wise volume-preserving constraint $\det (\nabla \mathbf{\phi }(\mathbf{x}))=1$ is effective in some cases, but may be too restrictive in others, especially when local deformations are relatively large. This can result in poor matching when enforcing large local deformations. In this paper, we propose a new bi-variant diffeomorphic image registration model that introduces a soft constraint on the Jacobian equation $\det (\nabla \mathbf{\phi }(\mathbf{x}))=f(\mathbf{x})>0$. This allows local deformations to shrink and grow within a flexible range $0<{\kappa }_m<\det (\nabla \mathbf{\phi }(\mathbf{x}))<{\kappa }_M$. The Jacobian determinant of transformation is explicitly controlled by optimizing the relaxation function $f(\mathbf{x})$. To prevent deformation folding and improve the smoothness of the transformation, a positive constraint is imposed on the optimization of the relaxation function $f(\mathbf{x})$, and a regularizer is used to ensure the smoothness of $f(\mathbf{x})$. Furthermore, the positivity constraint ensures that $f(\mathbf{x})$ is as close to one as possible, which helps to achieve a volume-preserving transformation on average. We also analyze the existence of the minimizer for the variational model and propose a penalty-splitting algorithm with a multilevel strategy to solve this model. Numerical experiments demonstrate the convergence of the proposed algorithm and show that the positivity constraint can effectively control the range of relative volume without compromising the accuracy of the registration. Moreover, the proposed model generates diffeomorphic maps for large local deformations and outperforms several existing registration models in terms of performance.

arXiv comment42 pages, Accepted to Applied Mathematical Modelling,130(2024)

References in corpus (8)

• The varifold representation of non-oriented shapes for diffeomorphic registration
• Estimation of Large Motion in Lung CT by Integrating Regularized Keypoint Correspondences into Dense Deformable Registration
• A Lagrangian Gauss-Newton-Krylov Solver for Mass- and Intensity-Preserving Diffeomorphic Image Registration
• Indirect Image Registration with Large Diffeomorphic Deformations
• Diffeomorphic density matching by optimal information transport
• A New Variational Model for Joint Image Reconstruction and Motion Estimation in Spatiotemporal Imaging
• A unifying framework for $n$-dimensional quasi-conformal mappings
• Diffeomorphic Image Registration with An Optimal Control Relaxation and Its Implementation