[Insight-users] output region of ResampleImageFilter

[Christoph Niedermayr]

Hi all!

I would like to get an output image containing the *entire* result of an
input image under some transformation T.

Right now i'm doing it this way:
- calculate T' := inverse of T
- calculate physical coordinates of the corners of the input image
- transform corners by T'
- get min/max coordinates of transformed corners
- use those coordinates to tell ResampleImageFilter region and origin of
output image
- run the ResampleImageFilter with transformation T

What i don't like is that i have to explicitly calculate the output
image's extents using the inverse transformation T' because this is not
generic (not all transforms can be inverted).

Anyone know of a better way to do this in itk?

best regards,
Chris