[GRASS-user] Need help importong LiDAR file with r.in.gdal

[Markus Neteler]

Rich,

On Mon, May 1, 2023 at 11:01 PM Rich Shepard <rshepard at appl-ecosys.com> wrote:
>
> I'm trying to import a LiDAR data set. I've used EPSG:6559 while the data
> source manager tells me it imports into ARC/gis using EPSG:6557. I've tried
> both and grass doesn't work for me regardless of which code I use.

[...]
> r.in.gdal -o in=$HOME/projects/oregon/northside-rock/data/LDQ-45122F8/2014-OLC-Metro/Bare_Earth/be45122f8/hdr.adf out=topo_2010 mem=2000 --o
> Over-riding projection check
> Importing raster map <topo_2010>...
> ERROR 1: Block is corrupt, block size was 21532, but expected to be 16384.
> ERROR 1: /home/rshepard/projects/oregon/northside-rock/data/LDQ-45122F8/2014-OLC-Metro/Bare_Earth/be45122f8, band 1: IReadBlock failed at X offset 15, Y offset 308: Block is corrupt, block size was 21532, but expected to be 16384.
> ERROR: Unable to read row 4928

What does it show in GDAL?

gdalinfo  $HOME/projects/oregon/northside-rock/data/LDQ-45122F8/2014-OLC-Metro/Bare_Earth/be45122f8/hdr.adf

The point is that GDAL needs to be able to read it (since r.in.gdal
depends on it).

Markus