That doesn't quite produce the desired result. The output extents of course form a rectangle, but
in the input coordinate system, those extents from a parallelogram. In only being able to specify the min,max X and Y dimensions, I must set a rectangle that encloses the output extents' parallelogram but which will include some data
that will not really be in the output scene. So this approach might reduce the amount of extra data, but it's not really
going to produce an exact match to an existing file.
In one of my searches on this issue I saw a suggestion of perhaps using theenvi_layer_stacking_doit with the /exclusive flag.
It will cause the output to include the 'sample' band in the output data which could subsequently be removed.
Seems to me there must be a better, albeit undocumented, way to do this since the interactive GUI can do it and the
layer stacking process does something very similar.
|