<div dir="auto">How is your region settings? Do they match the DEMs? Maybe importing with r.import will help, as it can override region settings <div dir="auto"><br></div><div dir="auto">Carlos</div><div dir="auto"><br></div><div dir="auto"><br></div><div dir="auto"><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sat, 25 Sep 2021, 11:43 Rich Shepard, <<a href="mailto:rshepard@appl-ecosys.com">rshepard@appl-ecosys.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">I have 77 LiDAR topo files in .tif format. What gdalinfo tells me about a<br>
representative file:<br>
$ gdalinfo columbia_2010_e_dtm_35.tif<br>
Driver: GTiff/GeoTIFF<br>
Files: columbia_2010_e_dtm_35.tif<br>
Size is 10745, 15264<br>
Coordinate System is:<br>
PROJCRS["NAD_1983_HARN_StatePlane_Washington_South_FIPS_4602_Feet",<br>
BASEGEOGCRS["GCS_North_American_1983_HARN",<br>
DATUM["North_American_1983_HARN",<br>
ELLIPSOID["GRS_1980",6378137,298.257222101,<br>
LENGTHUNIT["metre",1,<br>
ID["EPSG",9001]]]],<br>
PRIMEM["Greenwich",0,<br>
ANGLEUNIT["degree",0.0174532925199433,<br>
ID["EPSG",9122]]]],<br>
CONVERSION["Lambert Conic Conformal (2SP)",<br>
METHOD["Lambert Conic Conformal (2SP)",<br>
ID["EPSG",9802]],<br>
PARAMETER["Latitude of false origin",45.3333333333333,<br>
ANGLEUNIT["degree",0.0174532925199433],<br>
ID["EPSG",8821]],<br>
PARAMETER["Longitude of false origin",-120.5,<br>
ANGLEUNIT["degree",0.0174532925199433],<br>
ID["EPSG",8822]],<br>
PARAMETER["Latitude of 1st standard parallel",45.8333333333333,<br>
ANGLEUNIT["degree",0.0174532925199433],<br>
ID["EPSG",8823]],<br>
PARAMETER["Latitude of 2nd standard parallel",47.3333333333333,<br>
ANGLEUNIT["degree",0.0174532925199433],<br>
ID["EPSG",8824]],<br>
PARAMETER["Easting at false origin",500000,<br>
LENGTHUNIT["metre",1],<br>
ID["EPSG",8826]],<br>
PARAMETER["Northing at false origin",0,<br>
LENGTHUNIT["metre",1],<br>
ID["EPSG",8827]]],<br>
CS[Cartesian,2],<br>
AXIS["easting",east,<br>
ORDER[1],<br>
LENGTHUNIT["US survey foot",0.304800609601219,<br>
ID["EPSG",9003]]],<br>
AXIS["northing",north,<br>
ORDER[2],<br>
LENGTHUNIT["US survey foot",0.304800609601219,<br>
ID["EPSG",9003]]]]<br>
Data axis to CRS axis mapping: 1,2<br>
Origin = (1512164.000000000000000,152311.000000000000000)<br>
Pixel Size = (3.000000000000000,-3.000000000000000)<br>
Metadata:<br>
AREA_OR_POINT=Area<br>
DataType=Generic<br>
Image Structure Metadata:<br>
COMPRESSION=LZW<br>
INTERLEAVE=BAND<br>
Corner Coordinates:<br>
Upper Left ( 1512164.000, 152311.000) (121d 0' 8.57"W, 45d44'59.58"N)<br>
Lower Left ( 1512164.000, 106519.000) (121d 0' 4.46"W, 45d37'27.52"N)<br>
Upper Right ( 1544399.000, 152311.000) (120d52'34.01"W, 45d45' 1.35"N)<br>
Lower Right ( 1544399.000, 106519.000) (120d52'30.94"W, 45d37'29.29"N)<br>
Center ( 1528281.500, 129415.000) (120d56'19.49"W, 45d41'14.50"N)<br>
Band 1 Block=512x256 Type=Float32, ColorInterp=Gray<br>
Min=159.316 Max=1037.520<br>
Minimum=159.316, Maximum=1037.520, Mean=262.041, StdDev=157.494<br>
NoData Value=-3.40282306073709653e+38<br>
Overviews: 5373x7632, 2687x3816, 1344x1908, 672x954, 336x477, 168x239<br>
Metadata:<br>
BandName=Band_1<br>
RepresentationType=ATHEMATIC<br>
STATISTICS_COVARIANCES=24804.31297971113<br>
STATISTICS_MAXIMUM=1037.5201416016<br>
STATISTICS_MEAN=262.04051816209<br>
STATISTICS_MINIMUM=159.31565856934<br>
STATISTICS_SKIPFACTORX=1<br>
STATISTICS_SKIPFACTORY=1<br>
STATISTICS_STDDEV=157.49385060919<br>
<br>
I started GRASS with a new location:<br>
grass -c /path/to/project/data/topography/columbia_2010_e_dtm_35.tif /data/grassdata/columbia_river_dtms/topography<br>
<br>
and imported all raw data; e.g.,<br>
<br>
r.in.gdal<br>
in=/path/to/project/data/topography/columbia_2010_e_dtm_35.tif out=dtm_4 --o<br>
<br>
The problem is the imported maps have no elevation data. For example,<br>
r.report dtm_35<br>
100%<br>
+-----------------------------------------------------------------------------+<br>
| RASTER MAP CATEGORY REPORT |<br>
|LOCATION: new_nevins_dock Sat Sep 25 07:28:25 2021|<br>
|-----------------------------------------------------------------------------|<br>
| north: 106729 east: 1544180 |<br>
|REGION south: 60934 west: 1511876 |<br>
| res: 3 res: 3 |<br>
|-----------------------------------------------------------------------------|<br>
|MASK: none |<br>
|-----------------------------------------------------------------------------|<br>
|MAP: (untitled) (dtm_35 in topography) |<br>
|-----------------------------------------------------------------------------|<br>
| Category Information | cell| % |<br>
|#|description | count| cover|<br>
|-----------------------------------------------------------------------------|<br>
|*|no data. . . . . . . . . . . . . . . . . . . . . . . . . .|164373520|100.00|<br>
|-----------------------------------------------------------------------------|<br>
|TOTAL |164373520|100.00|<br>
+-----------------------------------------------------------------------------+<br>
<br>
r.stats -l in=dtm_35 out=-<br>
100%<br>
* no data<br>
<br>
The data file sizes are all in the low 100Ks. I could upload an example to<br>
<a href="http://fileconvoy.com" rel="noreferrer noreferrer" target="_blank">fileconvoy.com</a> if that's helpful.<br>
<br>
What have I done incorrectly that the imported maps have no elevation data<br>
while the raw .tif files do?<br>
<br>
TIA,<br>
<br>
Rich<br>
_______________________________________________<br>
grass-user mailing list<br>
<a href="mailto:grass-user@lists.osgeo.org" target="_blank" rel="noreferrer">grass-user@lists.osgeo.org</a><br>
<a href="https://lists.osgeo.org/mailman/listinfo/grass-user" rel="noreferrer noreferrer" target="_blank">https://lists.osgeo.org/mailman/listinfo/grass-user</a><br>
</blockquote></div>