[pdal] (pdal pipeline Error) GDAL failure (6) Cannot find coordinate operations from

[Andrew Bell]

Peter,

If you have a single file that will cause this behavior I can investigate.
I doubt anything is wrong, but it would be good to understand in detail.

Thanks,

On Tue, Apr 27, 2021 at 9:33 PM Peter Lim <peter.lim at gpslands.com> wrote:

> Hi ,
>
>
>
> When I removed all the srs from the pipeline, I got this error:
>
>
>
> (pdal pipeline writers.las Error) GDAL failure (1) PROJ:
> proj_get_ellipsoid: CRS has no geodetic CRS
>
> (pdal pipeline writers.las Error) GDAL failure (1) PROJ:
> proj_get_ellipsoid: Object is not a CRS or GeodeticReferenceFrame
>
> (pdal pipeline writers.las Error) GDAL failure (1) PROJ:
> proj_get_ellipsoid: CRS has no geodetic CRS
>
> (pdal pipeline writers.las Error) GDAL failure (1) PROJ:
> proj_get_ellipsoid: Object is not a CRS or GeodeticReferenceFrame
>
> (pdal pipeline Error) GDAL failure (6) Cannot find coordinate operations
> from
> `ENGCRS["unnamed",EDATUM[""],CS[Cartesian,2],AXIS["(E)",east,ORDER[1],LENGTHUNIT["metre",1,ID["EPSG",9001]]],AXIS["(N)",north,ORDER[2],LENGTHUNIT["metre",1,ID["EPSG",9001]]]]'
> to
> `ENGCRS["unnamed",EDATUM[""],CS[Cartesian,2],AXIS["(E)",east,ORDER[1],LENGTHUNIT["metre",1,ID["EPSG",9001]]],AXIS["(N)",north,ORDER[2],LENGTHUNIT["metre",1,ID["EPSG",9001]]]]'
>
> PDAL: filters.crop: Unable to reproject bounds.
>
>
>
>
>
> Modified pipeline:
>
> [
>
>     {
>
>         "type":"readers.las",
>
>         "filename":"*.las"
>
>     },
>
>     {
>
>         "type": "filters.merge"
>
>     },
>
>     {
>
>         "type":"filters.crop",
>
>         "bounds":"([364035,364085],[142446,142486],[-1000, 1000])"
>
>     },
>
>     {
>
>         "type":"writers.las",
>
>         "filename":"output.las"
>
>     }
>
> ]
>
>
>
> Please advise.
>
>
>
> Thank you.
>
>
>
> Best Regards,
>
>
>
> Peter Lim
>
> *GPS Lands (S) Pte Ltd*
>
>
>
> *From:* Andrew Bell [mailto:andrew.bell.ia at gmail.com]
> *Sent:* Tuesday, 27 April 2021 9:36 pm
> *To:* Kirk Waters - NOAA Federal
> *Cc:* Peter Lim; pdal
> *Subject:* Re: [pdal] (pdal pipeline Error) GDAL failure (6) Cannot find
> coordinate operations from
>
>
>
>
>
> On Tue, Apr 27, 2021 at 8:59 AM Kirk Waters - NOAA Federal <
> kirk.waters at noaa.gov> wrote:
>
> Perhaps I'm missing it, but it doesn't look like there is a transform to
> be done. The input and output codes are the same. Couldn't he just skip
> specifying the SRS?
>
>
>
> No, the input SRS is some ENGCRS. The option "default_srs" that's being
> used in readers.las only applies if the data doesn't already have an SRS.
> You can override the data SRS by using "override_srs", but I'm not sure if
> that's correct or what's intended in this case.
>
>
>
> Your data seems to be in a custom engineering CRS. PDAL can't convert your
> crop boundary from EPSG:32648 to your custom CRS because no transformation
> exists in PROJ, the library used for this purpose. You might ask on the
> PROJ list or perhaps someone else on this list will know how to deal
> with this.
>
>
>
>  I’m using the following pipeline:
>
>
>
> [
>
>     {
>
>         "type":"readers.las",
>
>         "filename":"*.las",
>
>
> "default_srs":"EPSG:32648"
>
>
>     },
>
> ... snipped ...
>
>
>
> --
>
> Andrew Bell
> andrew.bell.ia at gmail.com
>
>
> <https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=emailclient> Virus-free.
> www.avast.com
> <https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=emailclient>
> <#m_-658076429685796966_DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2>
>


-- 
Andrew Bell
andrew.bell.ia at gmail.com
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/pdal/attachments/20210428/16f22211/attachment.html>