How to calculate Jacobian matrix of a transformation with simpleITK?


How to calculate Jacobian matrix at each pixel with simpleITK after I have estimated a B-spline transformation between fixed_image and moving_image?

Do I need to convert the transformation into displacement field to compute the Jacobian matrix? Many thanks.

Hello @dada,

The only functionality provided by SimpleITK with respect to Jacobians is to compute the determinant of a deformation field. To do that, first get a deformation field representation of your transformation via the TransformToDisplacementFieldFilter and then run the DisplacementFieldJacobianDeterminantFilter filter.

1 Like