Invertibility in image registration is usually enforced using the flow of diffeomorphism approach: the transformation is the solution at t=1 of the nonlinear transport equation. This approach prevent the use of a parametric model to model the transformation. Consequently, when using B-Spline transform, soft or hard constraints on the transformation Jacobian are usually imposed during the registration to penalize non invertible transform as in