<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:0px 0px 0px 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_withGeoid.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."  When running the above pipeline in newer versions of PDAL, I get the above warning, and the output grid does not have the geoid applied.<div><br></div><div>I tried adjusting 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_withGeoid.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 again, the output grid does not 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><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><br></div></div></div></div>