[pdal] Support for proj4 pipeline string in filter.reprojection

[Aurelien Vila]

Hello dear PDAL community,

I would like to reproject a point cloud using a proj4 pipeline string such
as the following :
"+proj=pipeline
+step +proj=pipeline +step +proj=merc +step +proj=axisswap +order=2,1
+step +proj=unitconvert +xy_in=m +xy_out=us-ft" (see
https://proj.org/operations/pipeline.html for more details)
Since version 3.0, GDAL support this feature in this CLI tools gdalwarp and
ogr2ogr. It would be nice to have the same feature in PDAL.
The OSR class CoordinateTransformation that is used by PDAL to accomplish
coordinates transformation, support this through the method
SetCoordinateOperation(). Thus, it should be easy to implement by modifying
the SrsTransform class and providing minor API modifications in the
re-projection filters (add of coordinates_operation option and out_srs will
be no longer mandatory when specifying a coordinates operation).
I would like to know if someone in the community is currently working on
that. If not, I could handle it myself and submit a PR soon.

Thanks in advance

Aurélien Vila
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/pdal/attachments/20190716/29e21db2/attachment.html>