[pdal] pdal pipeline filters.reprojection

[Howard Butler]

You need the datum files added to your docker image. I shut that off because they were quite heavy and we were trying to put our docker image on a diet.

Make a new Dockerfile with the following (commented out at the moment in https://github.com/PDAL/PDAL/blob/master/scripts/docker/dependencies/Dockerfile#L248)

> FROM pdal/pdal:1.4
> MAINTAINER Howard Butler <howard at hobu.co>
> 
> 
> RUN apt-get update && apt-get install -y --fix-missing --no-install-recommends \
>         wget \
>     && rm -rf /var/lib/apt/lists/*
> 
> RUN mkdir /vdatum \
>     && cd /vdatum \
>     && wget http://download.osgeo.org/proj/vdatum/usa_geoid2012.zip && unzip -j -u usa_geoid2012.zip -d /usr/share/proj \
>     && wget http://download.osgeo.org/proj/vdatum/usa_geoid2009.zip && unzip -j -u usa_geoid2009.zip -d /usr/share/proj \
>     && wget http://download.osgeo.org/proj/vdatum/usa_geoid2003.zip && unzip -j -u usa_geoid2003.zip -d /usr/share/proj \
>     && wget http://download.osgeo.org/proj/vdatum/usa_geoid1999.zip && unzip -j -u usa_geoid1999.zip -d /usr/share/proj \
>     && wget http://download.osgeo.org/proj/vdatum/vertcon/vertconc.gtx && mv vertconc.gtx /usr/share/proj \
>     && wget http://download.osgeo.org/proj/vdatum/vertcon/vertcone.gtx && mv vertcone.gtx /usr/share/proj \
>     && wget http://download.osgeo.org/proj/vdatum/vertcon/vertconw.gtx && mv vertconw.gtx /usr/share/proj \
>     && wget http://download.osgeo.org/proj/vdatum/egm96_15/egm96_15.gtx && mv egm96_15.gtx /usr/share/proj \
>     && wget http://download.osgeo.org/proj/vdatum/egm08_25/egm08_25.gtx && mv egm08_25.gtx /usr/share/proj \
>     && rm -rf /vdatum
> 


Then go build a new container:

> docker build -t pdal/datum  - < Dockerfile


The GeoTIFF info in that file does look ok, but if you wanted to hammer it, you might do so like this:

> {
>   "pipeline":[
>     {
>       "type":"readers.las",
>       "spatialreference":"EPSG:4269+5703",
>       "filename":"/data/20110802_2012_NCMP_OH_Erie_24_GeoClassified_LAS_NAVD88.laz"
>     },
>     {
>       "type":"filters.reprojection",
>       "out_srs":"EPSG:3734"
>     },
>     {
>       "type":"filters.sample",
>       "radius":"5"
>     },
>     {
>       "type":"writers.las",
>       "filename":"/data/poisson_8-16.las"
>     }
>   ]
> }

Finally, run your docker environment with the noise turned up to see GDAL/Proj.4 doing your datum transform:

> docker run --env CPL_DEBUG=ON --env PROJ_DEBUG=ON -v `pwd`:/data pdal/datum pdal pipeline /data/erie.json

Hope this helps,

Howard

> On Feb 17, 2017, at 10:25 PM, Stephen Mather <stephen at smathermather.com> wrote:
> 
> Trying to process data from 2011-2012 USACE NCMP Topobathy Lidar: Lake Erie (specifically this tile).
> 
> JSON:
> {
>   "pipeline":[
>     "/data/20001800PAS.las",
>     {
>       "type":"readers.las",
>       "spatialreference":"EPSG:4269",
>       "filename":"/data/20001800PAS.las"
>     },
>     {
>       "type":"filters.reprojection",
>       "out_srs":"EPSG:3734"
>     },
>     {
>       "type":"filters.sample",
>       "radius":"5"
>     },
>     {
>       "type":"writers.las",
>       "filename":"/data/poisson_8-16.las"
>     }
>   ]
> }
> 
> docker run -v $(pwd)/data:/data pdal/pdal:1.4 pdal pipeline /data/pipeline.json --readers.las.filename=/data/lakefront.laz --writers.las.filename=/data/filtered.laz
> 
> I get the following error:
> 
> (pdal pipeline filters.reprojection Error: 0): GDAL failure (1) failed to load datum shift file
> 
> I was hoping by specifying the input spatialreference I could dodge this error, but alas... .
> 
> AFAIK, these are the same datums (data?):
> http://spatialreference.org/ref/epsg/4269/
> http://spatialreference.org/ref/epsg/3734/
> 
> What am I missing?
> Thanks in advance,
> Best,
> Steve
> _______________________________________________
> pdal mailing list
> pdal at lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/pdal