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

Rich Shepard rshepard at appl-ecosys.com
Mon May 1 14:01:25 PDT 2023


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.

I create a location:
grass -c EPSG:6557 /data1/grassdata/Metro_LiDAR

Then import the data set (N.B. a copy of the data source tarball is on the
fileconvoy.com web site; see below):
r.in.gdal in=Bare_Earth/be45122f8/hdr.adf out=topo_2010 mem=2000 --o

GRASS returns:
ERROR: Projection of dataset does not appear to match current location.

        Location PROJ_INFO is:
        name: NAD83(2011) / Oregon GIC Lambert (ft)
        datum: nad83_2011
        ellps: grs80
        proj: lcc
        lat_0: 41.75
        lon_0: -120.5
        lat_1: 43
        lat_2: 45.5
        x_0: 399999.9999984
        y_0: 0
        no_defs: defined

        Dataset PROJ_INFO is:
        name: unnamed
        ellps: grs80
        proj: lcc
        lat_0: 41.75
        lon_0: -120.5
        lat_1: 43
        lat_2: 45.5
        x_0: 400000
        y_0: 0
        no_defs: defined

        Difference in: datum

        In case of no significant differences in the projection definitions,
        use the -o flag to ignore them and use current location definition.
        Consider generating a new location from the input dataset using the
        'location' parameter.

I added the -o option to r.in.gdal and this is the result:
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

If the file loads into ARC/gis why not into grass?

I do not see why this isn't working here. I've uploaded the Bare_Earth
directory to
http://www.fileconvoy.com/dfl.php?id=g4e3d18b830b60bef10004916731ceed9a68e0d28c3.
It's a 474M file.

I need to learn why I have problems importing the data using either EPSG
code, ignoring the data project by override or not.

Regards,

Rich


More information about the grass-user mailing list