<div dir="ltr">Hello,<div><br></div><div>In older versions of PDAL I used to be able to run a pipeline like:</div><div><br></div><blockquote style="margin:0 0 0 40px;border:none;padding:0px"><div><i>{"pipeline": [{"type" : "readers.las",</i></div><div><i>                  "filename": "input_Ellipsoidal.laz"},</i></div><div><i>                 {"type":"filters.reprojection",</i></div><div><i>                  "in_srs":"EPSG:32611",</i></div><div><i>                  "out_srs":"+init=EPSG:32611 +geoidgrids=g2003conus.gtx"},</i></div><div><i>                 {"type" : "writers.las",</i></div><div><i>                  "filename": "output_Geoid.laz",</i></div><div><i>                  "compression": "laszip",</i></div><div><i>                  "a_srs": "EPSG:32611+5703"</i></div><div><i>                 }]}</i></div></blockquote><div><br></div><div>This would run fine, and the output would have the geoid properly applied.</div><div><br></div><div>However, newer versions of PDAL (e.g. 2.5.3) uses newer version of PROJ, which no longer likes its string in the form, "+init=epsg:XXXX".  It will give the warning:</div>"Warning 1: +init=epsg:XXXX syntax is deprecated. It might return a CRS with a non-EPSG compliant axis order."





<div><br></div><div>I adjusted the PROJ string to remove this convention:</div><div><br></div><div><div><i>{"pipeline": [{"type" : "readers.las",</i></div><div><i>                  "filename": "input_Ellipsoidal.laz"},</i></div><div><i>                 {"type":"filters.reprojection",</i></div><div><i>                  "in_srs":"EPSG:32611",</i></div><div><i>                  "out_srs":"</i>+proj=utm +zone=11 +type=crs<i> +geoidgrids=g2003conus.gtx"},</i></div><div><i>                 {"type" : "writers.las",</i></div><div><i>                  "filename": "output_Geoid.laz",</i></div><div><i>                  "compression": "laszip",</i></div><div><i>                  "a_srs": "EPSG:32611+5703"</i></div><div><i>                 }]}</i></div></div><div><br></div><div>This runs without warning or errors, but the output grid <b>does not</b> have the geoid applied.  Most of the documentation on this type of operation I've seen still refers to the older PROJ syntax of, "+init=epsg:XXXX".</div><div><br></div><div>Is there a way to apply a geoid with the new PROJ syntax, and hopefully I am just missing something with the proper syntax to do this operation on the newer versions of PDAL?</div><div><br></div><div>Thanks!</div></div>