<div dir="ltr"><div dir="ltr">Well, I have gotten a little further in my adventures on the home computer.<div>I found another instance of proj , removed it and cleaned all references out. Just using proj 5.2 now</div><div>compiled libgeotiff from source</div><div>recompiled gdal 2.4.0 - (dynamic link to proj)</div><div>using pdal 1.8</div><div>The following command does a horizontal reprojection but scrubs out the NAVD88 vertical projection and units in the metadata ( as expected), but leaves ( at least the max) vertical values unchanged </div><div><span style="font-family:monospace"><span style="color:rgb(0,0,0)">pdal translate 4699/20130405_las_s23_2510_40.laz 4699/alb_20130405_l</span><br>as_s23_2510_40.laz -f filters.reprojection --filters.reprojection.in_srs="EPSG:4269" --filters.reprojection.out_srs="EPSG:5070"<br></span></div><div><span style="font-family:monospace"><br></span></div><div><font face="monospace">Removing the input projection assignment repeatedly gives:</font></div><div><span style="font-family:monospace"><span style="color:rgb(0,0,0)">(pdal translate filters.reprojection Error) GDAL failure (1) failed to load datum shift file</span><br></span></div><div><span style="font-family:monospace"><span style="color:rgb(0,0,0)"><br></span></span></div><div><font color="#000000" face="monospace">Trying </font><span style="color:rgb(0,0,0);font-family:monospace">pdal translate 4699/20130405_las_s23_2510_40.laz 4699/alb_20130405_l</span></div><span style="font-family:monospace">as_s23_2510_40.laz -f filters.reprojection --filters.reprojection.in_srs="EPSG:4269" --filters.reprojection.out_srs="EPSG:5070+5<br>703"<br></span><div><font face="monospace"><br></font></div><div><font face="monospace">Gives the same error. In an effort to find the right vertical gtx file /usr/local/share/proj now has :</font></div><div><font face="monospace"><br></font></div><div><span style="font-family:monospace"><span style="color:rgb(0,0,0)">alaska epsg GL IGNF nad27 null README.DATUMGRID TN WO
</span><br>BETA2007.gsb esri GL27 ITRF2000 nad83 nzgd2kgrid0005.gsb README.NORTHAMERICA vertconc.gtx world
<br>CH esri.extra gvr2000.gtx ITRF2008 nad.lst other.extra stgeorge vertcone.gtx
<br>conus FL gvr2016.gtx ITRF2014 ntf_r93.gsb proj_def.dat stlrnc vertconw.gtx
<br>egm96_15.gtx g2012ba0.gtx hawaii MD ntv1_can.dat prvi stpaul WI<br>
<br></span></div><div><font face="monospace">not quite stuffing every datum grid into the directory.</font></div><div><font face="monospace"><br></font></div><div><font face="monospace">Still get the same error even with the above projection information in the /usr/local/share/proj directory</font></div><div><font face="monospace">googling around, I found</font></div><div><font face="monospace"><a href="https://www.gadom.ski/2017/05/31/vdatum-with-pdal.html">https://www.gadom.ski/2017/05/31/vdatum-with-pdal.html</a><br></font></div><div><font face="monospace"><br></font></div><div><font face="monospace">So I tried</font></div><div><font face="monospace"><br></font></div><div><span style="font-family:monospace"><span style="color:rgb(0,0,0)">pdal translate 4699/20130405_las_s23_2510_40.laz 4699/alb_20130405_las_s23_2510_40.laz -f filters.reprojection --filters.reprojection.in_srs="EPSG:4269" --filters.reprojection.out_srs="init=EPSG:507</span>0 +geoidgrids=/usr/localshare/proj/g2012ba0.gtx" --writers.las.a_srs=EPSG:5070+5703<br></span></div><div><span style="font-family:monospace"><br></span></div><div><span style="font-family:monospace">and got </span><span style="color:rgb(0,0,0);font-family:monospace">PDAL: Could not import coordinate system 'init=EPSG:5070 +geoidgrids=/usr/localshare/proj/g2012ba0.gtx': (unknown reason).</span></div><div><span style="color:rgb(0,0,0);font-family:monospace"><br></span></div><div><span style="color:rgb(0,0,0);font-family:monospace"><br></span></div><div><span style="color:rgb(0,0,0);font-family:monospace">So now, having thoroughly displayed my ignorance on how to handle gtx files, does anyone have any suggestions?</span></div><div><span style="color:rgb(0,0,0);font-family:monospace"><br></span></div><div><span style="color:rgb(0,0,0);font-family:monospace">Doug </span></div><span style="font-family:monospace">
<br></span><div><span style="font-family:monospace"><br></span></div><div><font face="monospace"><br></font></div><div><font face="monospace"><br></font></div><div><font face="monospace"><br></font></div><div><br></div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Mar 5, 2019 at 12:59 PM Jed Frechette <<a href="mailto:jedfrechette@gmail.com" target="_blank">jedfrechette@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">On Tue, 5 Mar 2019 10:06:08 -0600 Jim Klassen <<a href="mailto:klassen.js@gmail.com" target="_blank">klassen.js@gmail.com</a>> wrote:<br>
> ...clip...<br>
> I'm not sure the "right" answer, but I got around that (first) error by specifying the full path to the geoidgrid file. This required specifying the proj string (instead of the EPSG string) in the reprojection "out_srs" field.<br>
>...clip...<br>
<br>
Not to derail this thread, but since I just ran in to this with the<br>
Windows version installed via conda-forge as well: Is the<br>
filters.reprojection usage shown in Example 2 of the docs expected to<br>
work? I was not able to combine +init +geogrids strings as shown in<br>
that example and needed to specify the full proj string.<br>
<br>
Best,<br>
<br>
-- <br>
Jed Frechette<br>
_______________________________________________<br>
pdal mailing list<br>
<a href="mailto:pdal@lists.osgeo.org" target="_blank">pdal@lists.osgeo.org</a><br>
<a href="https://lists.osgeo.org/mailman/listinfo/pdal" rel="noreferrer" target="_blank">https://lists.osgeo.org/mailman/listinfo/pdal</a></blockquote></div>