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
_______________________________________________
grass-user mailing list
[email protected]
https://lists.osgeo.org/mailman/listinfo/grass-user
