Skip to content

VectorDataReprojection app does not output VectorData in physical space when using a reference image in sensor geometry

Description

When using the VectorDataReprojection application on a vector data to project it on a reference sensor image coordinate system, the ouptut vector data is in ContinuousIndex reference (rows and columns index), not in physical space attached to image (i.e. using spacing and origin). This usually not the desired behaviour, since origin and spacing account for extract and possible subsampling of image, and GIS systems will make use of it to superimpose data.

This is due to those two lines :

https://gitlab.orfeo-toolbox.org/orfeotoolbox/otb/blob/develop/Modules/Applications/AppProjection/app/otbVectorDataReprojection.cxx#L154 https://gitlab.orfeo-toolbox.org/orfeotoolbox/otb/blob/develop/Modules/Applications/AppProjection/app/otbVectorDataReprojection.cxx#L155

If those two lines are set, the output origin and spacing will be passed down to the GenericRSTransform, which will apply it to get back to the ContinuousIndex reference.

We can add a mode to enable / disable those two lines or (preferabily for me) remove them completely.

Even the support of input / output origin and spacing is questionable I think.

Steps to reproduce

Reproject a GIS vector data using the application on an sensor geometry image for each spacing != (1,1) or origin !=(0,0).

Configuration information

OTB 6.6, all OS.