OSGeo / grass

GRASS GIS - free and open-source geospatial processing engine
https://grass.osgeo.org
Other
845 stars 307 forks source link

r.in.pdal fails with LAS files without projection #4157

Open CarlosGrohmann opened 2 months ago

CarlosGrohmann commented 2 months ago

Describe the bug

I'm using GRASS 8.4 from ubuntugis-unstable. When trying to import a lidar dataset I got this error:

terminate called after throwing an instance of
'pdal::pdal_error'
  what():  readers.las: Global encoding WKT flag not set for
point format 6 - 10.

Searching for the error, I found mentions that this can be avoided by using readers.las.nosrs=true but there is no way to add it this option to r.in.pdal.

One field to add extra options to the pdal pipeline could fix this.

System description

external links

https://github.com/PDAL/PDAL/issues/3803 https://gis.stackexchange.com/questions/413191/python-pdal-error-reading-format-1-4-las-file-readers-las-error-global-enco

marisn commented 2 months ago

Could you, please, provide a small file with test data?

CarlosGrohmann commented 2 months ago

Sure, here it is: https://www.dropbox.com/scl/fi/pa6baxi237tqq0o9nu6ou/example.las.zip?rlkey=macixpzd9czt4dgriperrj45b&dl=0