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