A bi-variant variational model for diffeomorphic image registration with relaxed Jacobian determinant constraints.

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 ⁡ (∇ φ (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 ⁡ (∇ φ (x)) = f (x) > 0. This allows local deformations to shrink and grow within a flexible range 0 < κ m < det ⁡ (∇ φ (x)) < κ M. The Jacobian determinant of transformation is explicitly controlled by optimizing the relaxation function f (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 (x) , and a regularizer is used to ensure the smoothness of f (x). Furthermore, the positivity constraint ensures that f (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. • A novel bi-variant diffeomorphic image registration model with relaxed Jacobian determinant constraints is proposed. • The existence of the minimizer for the variational model is proved. • A penalty splitting algorithm with a multilevel strategy is designed to solve the new model. • Numerical experiments show that the proposed algorithm is convergent and the new model has good performance. [ABSTRACT FROM AUTHOR]