Add `jacobian` method to `_GeometricTransform`

The _GeometricTransform could have the method

@abstractmethod
def jacobian(self, coords)

which returns the Jacobian matrix Jacobian matrix and determinant - Wikipedia for the given points. This can for example be used for optimisation methods as least_squares — SciPy v1.14.0 Manual.

I think it is rather straight forward to implement for the different transformations.

Do you think this is a valuable contribution or is it too specific? If wanted I can take care of the implementation.

1 Like

Interesting idea! It should help with convergence in model fitting, so would be worth trying.