[pdal] (pdal pipeline Error) GDAL failure (6) Cannot find coordinate operations from
Peter Lim
peter.lim at gpslands.com
Tue Apr 27 19:35:05 PDT 2021
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>
More information about the pdal
mailing list