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

[Peter Lim]

Hi,

 

In addition, when I tried to show the info 1 of the file using:

pdal info las\file.las

 

It gave this error:

PDAL: Geometry::transform() failed.  NULL source SRS.

 

 

Pls. advise.

 

Thank you.

 

Best Regards,

 

Peter Lim

GPS Lands (S) Pte Ltd

 

From: Peter Lim [mailto:peter.lim at gpslands.com] 
Sent: Wednesday, 28 April 2021 9:34 am
To: 'Andrew Bell'; 'Kirk Waters - NOAA Federal'
Cc: 'pdal'
Subject: RE: [pdal] (pdal pipeline Error) GDAL failure (6) Cannot find coordinate operations from

 

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



-- 
This email has been checked for viruses by Avast antivirus software.
https://www.avast.com/antivirus
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/pdal/attachments/20210428/cb257156/attachment-0001.html>