Implement Jacobian attenuation in the target space

[effigies]

What would you like to see added in this software?

SDCFlows should provide a method to generate a Jacobian weight image to multiply unwarped images by.

Like the fieldmap estimate, the Jacobian should be sampled into the target space, so that it is a simple post-resampling multiplication.

To reduce copies and compounding floating-point error (for example, if we resample as float64 and save as float32), ideally our resampler would also accept a Jacobian and multiply before writing to disk.

Do you have any interest in helping implement the feature?

Yes

Additional information / screenshots

No response