[pdal] Trouble reading projection from some 3DEP tiles

[Andrew Bell]

Jim,

pdal info reports the stored spatial reference. I think all the other
information can be found from projinfo. I would ask on the proj
mailing list if you can't figure it out.

On Wed, Sep 11, 2024 at 5:12 PM Andrew Bell <andrew.bell.ia at gmail.com> wrote:
>
> Is this a different issue than the one you reported earlier in the day?
>
> On Wed, Sep 11, 2024 at 3:49 PM Jim Klassen via pdal
> <pdal at lists.osgeo.org> wrote:
> >
> > I've noticed that for some of the 3DEP tiles that PDAL doesn't calculate
> > a 4326 BBOX.  This seems to be related to Proj not liking the how the
> > projection is specified in the LAZ file.  If I force the projection as
> > EPSG:6344+5703 in the reader then "pdal info" works as expected.  My
> > main question is how can I tell if the issue is with the CRS in the LAZ
> > file or the problem is with PDAL/Proj not being flexible enough in
> > parsing the CRS?
> >
> > Steps to reproduce:
> >      wget
> > https://rockyweb.usgs.gov/vdelivery/Datasets/Staged/Elevation/LPC/Projects/MN_UpperMissRiver_B22/MN_UpperMissRiver_2_B22/LAZ/USGS_LPC_MN_UpperMissRiver_B22_334_5271.laz
> >
> >      pdal info --all USGS_LPC_MN_UpperMissRiver_B22_334_5271.laz
> >
> > Observations:
> >      Under "stac" there is a error with message "Failed to create STAC
> > Feature with missing key. 'EPSG:4326'".
> >
> >      Under stats.bbox there is no 4326 bbox calculated.
> >
> > And with --verbose 8 I get:
> >
> > (PDAL Debug) Debugging...
> > (pdal info Debug) Executing pipeline in stream mode.
> > (pdal info filters.stats Debug) GDAL debug:
> > OGRSpatialReference::Validate: Coordinate system of GeographicCRS in the
> > WKT definition is different from the one of the authority. Unsetting the
> > identifier to avoid confusion
> > (pdal info filters.stats Debug) GDAL debug:
> > OGRSpatialReference::Validate: Parsing error : syntax error, unexpected
> > AUTHORITY, expecting number. Error occurred around:
> > UM["North American Vertical Datum
> > 1988",AUTHORITY["EPSG","5103"]],UNIT["Meter",1
> >                                          ^
> > (pdal info filters.hexbin Debug) GDAL debug:
> > OGRSpatialReference::Validate: Coordinate system of GeographicCRS in the
> > WKT definition is different from the one of the authority. Unsetting the
> > identifier to avoid confusion
> > (pdal info filters.hexbin Debug) GDAL debug:
> > OGRSpatialReference::Validate: Parsing error : syntax error, unexpected
> > AUTHORITY, expecting number. Error occurred around:
> > UM["North American Vertical Datum
> > 1988",AUTHORITY["EPSG","5103"]],UNIT["Meter",1
> > _______________________________________________
> > pdal mailing list
> > pdal at lists.osgeo.org
> > https://lists.osgeo.org/mailman/listinfo/pdal
>
>
>
> --
> Andrew Bell
> andrew.bell.ia at gmail.com



-- 
Andrew Bell
andrew.bell.ia at gmail.com