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

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


On Mon, 1 May 2023, Markus Neteler wrote:

> What does it show in GDAL?

Markus,

I forgot about gdalinfo and the file to which it's applied.

For this data set:
$ gdalinfo hdr.adf
Driver: AIG/Arc/Info Binary Grid
Files: .
        ./prj.adf
        ./w001001.adf
        ./w001001x.adf
        ./sta.adf
        ./metadata.xml
        ./dblbnd.adf
        ./hdr.adf
Size is 11097, 15505
Coordinate System is:
PROJCRS["unnamed",
     BASEGEOGCRS["Unknown datum based upon the GRS 1980 ellipsoid",
         DATUM["Not specified (based on GRS 1980 ellipsoid)",
             ELLIPSOID["GRS 1980",6378137,298.257222101,
                 LENGTHUNIT["metre",1]]],
         PRIMEM["Greenwich",0,
             ANGLEUNIT["degree",0.0174532925199433]],
         ID["EPSG",4019]],
     CONVERSION["unknown",
         METHOD["Lambert Conic Conformal (2SP)",
             ID["EPSG",9802]],
         PARAMETER["Latitude of false origin",41.75,
             ANGLEUNIT["degree",0.0174532925199433],
             ID["EPSG",8821]],
         PARAMETER["Longitude of false origin",-120.5,
             ANGLEUNIT["degree",0.0174532925199433],
             ID["EPSG",8822]],
         PARAMETER["Latitude of 1st standard parallel",43,
             ANGLEUNIT["degree",0.0174532925199433],
             ID["EPSG",8823]],
         PARAMETER["Latitude of 2nd standard parallel",45.5,
             ANGLEUNIT["degree",0.0174532925199433],
             ID["EPSG",8824]],
         PARAMETER["Easting at false origin",1312335.95800525,
             LENGTHUNIT["user-defined",0.3048],
             ID["EPSG",8826]],
         PARAMETER["Northing at false origin",0,
             LENGTHUNIT["user-defined",0.3048],
             ID["EPSG",8827]]],
     CS[Cartesian,2],
         AXIS["(E)",east,
             ORDER[1],
             LENGTHUNIT["user-defined",0.3048]],
         AXIS["(N)",north,
             ORDER[2],
             LENGTHUNIT["user-defined",0.3048]]]
Data axis to CRS axis mapping: 1,2
Origin = (672763.500000000000000,1467820.500000000000000)
Pixel Size = (3.000000000000000,-3.000000000000000)
Corner Coordinates:
Upper Left  (  672763.500, 1467820.500) (123d 0'19.54"W, 45d44'59.57"N)
Lower Left  (  672763.500, 1421305.500) (122d59'59.58"W, 45d37'20.60"N)
Upper Right (  706054.500, 1467820.500) (122d52'30.33"W, 45d45' 9.33"N)
Lower Right (  706054.500, 1421305.500) (122d52'11.40"W, 45d37'30.34"N)
Center      (  689409.000, 1444563.000) (122d56'15.21"W, 45d41'15.03"N)
Band 1 Block=256x16 Type=Float32, ColorInterp=Undefined
   Min=11.370 Max=1163.480
   NoData Value=-3.4028235e+38

Yet ARC/gis apparently has no issues with it.

Regards,

Rich


More information about the grass-user mailing list