[pdal] filters reprojection problem

[Doug Newcomb]

Well, I have gotten a little further in my adventures on the home computer.
I found another instance of proj , removed it and cleaned all references
out. Just using proj 5.2 now
compiled libgeotiff from source
recompiled gdal 2.4.0 - (dynamic link to proj)
using pdal 1.8
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
pdal translate 4699/20130405_las_s23_2510_40.laz 4699/alb_20130405_l
as_s23_2510_40.laz -f filters.reprojection
--filters.reprojection.in_srs="EPSG:4269"
--filters.reprojection.out_srs="EPSG:5070"

Removing the input projection assignment repeatedly gives:
(pdal translate filters.reprojection Error) GDAL failure (1) failed to load
datum shift file

Trying pdal translate 4699/20130405_las_s23_2510_40.laz 4699/alb_20130405_l
as_s23_2510_40.laz -f filters.reprojection
--filters.reprojection.in_srs="EPSG:4269"
--filters.reprojection.out_srs="EPSG:5070+5
703"

Gives the same error. In an effort to find the right vertical gtx file
/usr/local/share/proj now has :

alaska        epsg          GL           IGNF      nad27         null
               README.DATUMGRID     TN            WO
BETA2007.gsb  esri          GL27         ITRF2000  nad83
        nzgd2kgrid0005.gsb  README.NORTHAMERICA  vertconc.gtx  world
CH            esri.extra    gvr2000.gtx  ITRF2008  nad.lst
      other.extra         stgeorge             vertcone.gtx
conus         FL            gvr2016.gtx  ITRF2014  ntf_r93.gsb
  proj_def.dat        stlrnc               vertconw.gtx
egm96_15.gtx  g2012ba0.gtx  hawaii       MD        ntv1_can.dat  prvi
               stpaul               WI

not quite stuffing every datum grid into the directory.

Still get the same error even with the above projection information in the
/usr/local/share/proj directory
googling around, I found
https://www.gadom.ski/2017/05/31/vdatum-with-pdal.html

So I tried

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:5070
+geoidgrids=/usr/localshare/proj/g2012ba0.gtx"
--writers.las.a_srs=EPSG:5070+5703

and got PDAL: Could not import coordinate system 'init=EPSG:5070
+geoidgrids=/usr/localshare/proj/g2012ba0.gtx': (unknown reason).


So now, having thoroughly displayed my ignorance on how to handle gtx
files, does anyone have any suggestions?

Doug







On Tue, Mar 5, 2019 at 12:59 PM Jed Frechette <jedfrechette at gmail.com>
wrote:

> On Tue, 5 Mar 2019 10:06:08 -0600 Jim Klassen <klassen.js at gmail.com>
> wrote:
> > ...clip...
> > 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.
> >...clip...
>
> Not to derail this thread, but since I just ran in to this with the
> Windows version installed via conda-forge as well: Is the
> filters.reprojection usage shown in Example 2 of the docs expected to
> work? I was not able to combine +init +geogrids strings as shown in
> that example and needed to specify the full proj string.
>
> Best,
>
> --
> Jed Frechette
> _______________________________________________
> pdal mailing list
> pdal at lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/pdal
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/pdal/attachments/20190322/1e7a1b9e/attachment.html>