[GRASS-user] Problem importing ESRI grid file (.adf)

[Rich Shepard]

On Sat, 9 Jun 2018, Helmut Kudrnovsky wrote:

> As always recommended, look at the output of gdalinfo, compare it with the
> computional region of your interest and then decide if the data expects
> your need.

Helmut,

   I don't understand what you're suggesting. The coverage of the LiDAR data
is for a specific 1:24000 topographic quad map which contains the project
area of interest. When I import the data (or try to, anyway) I use the '-o'
option to set the computational region to that of the data being imported.

   The LiDAR coordinator for the agency wrote that, "The shapefile is
project[sic] in EPSG 2995 - NAD83(2011) Oregon Lambert Conformal Conic."

   When I run 'gdalinfo --stats' on the Bare_Earth subdirectory I see this:

Driver: AIG/Arc/Info Binary Grid
Files: bh45122c3
        bh45122c3.ovr
        bh45122c3.aux.xml
        bh45122c3/hdr.adf
        bh45122c3/dblbnd.adf
        bh45122c3/sta.adf
        bh45122c3/w001001.adf
        bh45122c3/metadata.xml
        bh45122c3/w001001x.adf
        bh45122c3/prj.adf
Size is 11052, 15426
Coordinate System is:
PROJCS["unnamed",
     GEOGCS["Unknown datum based upon the GRS 1980 ellipsoid",
         DATUM["Not_specified_based_on_GRS_1980_ellipsoid",
             SPHEROID["GRS 1980",6378137,298.257222101,
                 AUTHORITY["EPSG","7019"]],
             AUTHORITY["EPSG","6019"]],
         PRIMEM["Greenwich",0,
             AUTHORITY["EPSG","8901"]],
         UNIT["degree",0.0174532925199433,
             AUTHORITY["EPSG","9122"]],
         AUTHORITY["EPSG","4019"]],
     PROJECTION["Lambert_Conformal_Conic_2SP"],
     PARAMETER["standard_parallel_1",43],
     PARAMETER["standard_parallel_2",45.5],
     PARAMETER["latitude_of_origin",41.75],
     PARAMETER["central_meridian",-120.5],
     PARAMETER["false_easting",1312335.958005249],
     PARAMETER["false_northing",0],
     UNIT["user-defined",0.3048]]
Origin = (829504.500000000000000,1326862.500000000000000)
Pixel Size = (3.000000000000000,-3.000000000000000)
Metadata:
   DataType=Generic
Corner Coordinates:
Upper Left  (  829504.500, 1326862.500) (122d22'44.53"W, 45d22'29.77"N)
Lower Left  (  829504.500, 1280584.500) (122d22'29.72"W, 45d14'52.95"N)
Upper Right (  862660.500, 1326862.500) (122d15' 0.15"W, 45d22'37.00"N)
Lower Right (  862660.500, 1280584.500) (122d14'46.37"W, 45d15' 0.16"N)
Center      (  846082.500, 1303723.500) (122d18'45.19"W, 45d18'45.04"N)
Band 1 Block=256x16 Type=Float32, ColorInterp=Undefined
   Min=206.380 Max=1157.470
   Minimum=206.380, Maximum=1157.470, Mean=526.722, StdDev=308.445
   NoData Value=-3.4028234663852886e+38
   Overviews: 5526x7713, 2763x3857, 1382x1929, 691x965, 346x483, 173x242
   Metadata:
     STATISTICS_MAXIMUM=1157.4699707031
     STATISTICS_MEAN=526.72212274996
     STATISTICS_MINIMUM=206.38000488281
     STATISTICS_STDDEV=308.44451906832

   This looks OK to me and does not help me understand why my r.in.gdal
command line does not produce a rectangular map with the elevations colored
by the default scheme.

Regards,

Rich