[pdal] Change vertical datum

[Howard Butler]

> On Jun 7, 2022, at 9:07 AM, Ulises Ibarra <ulisesmartinibarra at gmail.com> wrote:
> 
> Hello everyone:
> 
> I'm not sure that the topics that I expose correspond to the topics that are covered in this list, if not, I wish you could tell me where to go.  Maybe on the PROJ list?
> 
> When I do a vertical datum shift of a point cloud, using the following pipeline:
> 
> [
> "input.las",
> {
> "type":"filters.reprojection",
> "in_srs":"EPSG:32614",
> "out_srs":"+init=epsg:32614 +geoidgrids=C:/nubes/geoides/GGM10_GEO.gtx"
> },
> {
> "type":"writers.las",
> "filename":"output_with_Vertical_Datum_GGM2010.las"
> }
> ]
> 
> Questions:
> 
> 1.- After running the command I get an alert:
> 
> " Warning 1: +init=epsg:XXXX syntax is deprecated. It might return a CRS with a non-EPSG compliant axis order. "
> 
> Could you tell me what I'm doing wrong at this point?
> 
> 2.- What argument would I have to use to see in the metadata something like "Vertical Datum": "GGM10" or "Mexican Gravimetric Geoid 2010"? Of course, in case the datum change has been made.

Ulises,

PDAL supports whatever version of PROJ you have linked.

See https://pdal.io/tutorial/grid-shift.html for some inspiration on how facts might apply to your specific geodetic configuration. 

Additionally, you might find https://pdal.io/stages/filters.projpipeline.html useful for specific application of PROJ transformation steps to data.

Howard