[GRASS-user] Data ranges of geotiff not preserved after import to GRASS
William Hudspeth
bhudspeth at edac.unm.edu
Fri Aug 12 15:01:56 PDT 2016
Thank you. That solution worked....
On 08/11/2016 01:34 AM, Markus Metz wrote:
> On Wed, Aug 10, 2016 at 10:29 PM, Bill Hudspeth <bhudspeth at edac.unm.edu> wrote:
>> Hello
>>
>> I am trying to import a multiband geotiff into GRASS 7.0.3 (running on Ubuntu 16.01). The gdalinfo output shows the 4 bands, and their data ranges. Both my geotiff and my GRASS location are in the projection UTM Zone 13N, WGS84 (epsg:32613). My GRASS region is 0.5 meters to import a panchromatic band, but the 4 bands in the multiband file are 2.0 m resolution. After I import, the data ranges of the 4 new grass rasters are very different from those listed in the gdalinfo output below.
> The data ranges reported by GRASS exceed those reported by GDAL. Most
> probably the data ranges reported by GDAL are approximations and thus
> a subset of the real ranges. Try gdalinfo -mm
> o10APR11180709-M2AS-055552078010_01_P001.tif and check if min/max
> values change.
>
> Markus M
>
>> gdalinfo o10APR11180709-M2AS-055552078010_01_P001.tif
>> Driver: GTiff/GeoTIFF
>> Files: o10APR11180709-M2AS-055552078010_01_P001.tif
>> o10APR11180709-M2AS-055552078010_01_P001.tif.aux.xml
>> Size is 3534, 2684
>> Coordinate System is:
>> PROJCS["WGS 84 / UTM zone 13N",
>> GEOGCS["WGS 84",
>> DATUM["WGS_1984",
>> SPHEROID["WGS 84",6378137,298.257223563,
>> AUTHORITY["EPSG","7030"]],
>> AUTHORITY["EPSG","6326"]],
>> PRIMEM["Greenwich",0],
>> UNIT["degree",0.0174532925199433],
>> AUTHORITY["EPSG","4326"]],
>> PROJECTION["Transverse_Mercator"],
>> PARAMETER["latitude_of_origin",0],
>> PARAMETER["central_meridian",-105],
>> PARAMETER["scale_factor",0.9996],
>> PARAMETER["false_easting",500000],
>> PARAMETER["false_northing",0],
>> UNIT["metre",1,
>> AUTHORITY["EPSG","9001"]],
>> AUTHORITY["EPSG","32613"]]
>> Origin = (318626.000000000000000,3848648.000000000000000)
>> Pixel Size = (2.000000000000000,-2.000000000000000)
>> Metadata:
>> AREA_OR_POINT=Area
>> TIFFTAG_ARTIST=PCI
>> TIFFTAG_DATETIME=2010:04:11 18:07:09
>> TIFFTAG_SOFTWARE=PCI Geomatica
>> Image Structure Metadata:
>> INTERLEAVE=PIXEL
>> Corner Coordinates:
>> Upper Left ( 318626.000, 3848648.000) (106d58'54.51"W, 34d45'49.98"N)
>> Lower Left ( 318626.000, 3843280.000) (106d58'50.35"W, 34d42'55.81"N)
>> Upper Right ( 325694.000, 3848648.000) (106d54'16.62"W, 34d45'54.41"N)
>> Lower Right ( 325694.000, 3843280.000) (106d54'12.63"W, 34d43' 0.24"N)
>> Center ( 322160.000, 3845964.000) (106d56'33.53"W, 34d44'25.13"N)
>> Band 1 Block=3534x2 Type=UInt16, ColorInterp=Red
>> Min=145.000 Max=473.000
>> Minimum=145.000, Maximum=473.000, Mean=286.443, StdDev=29.355
>> NoData Value=0
>> Metadata:
>> STATISTICS_MAXIMUM=473
>> STATISTICS_MEAN=286.44285428343
>> STATISTICS_MINIMUM=145
>> STATISTICS_STDDEV=29.354503358529
>> Band 2 Block=3534x2 Type=UInt16, ColorInterp=Green
>> Min=150.000 Max=771.000
>> Minimum=150.000, Maximum=771.000, Mean=447.144, StdDev=61.481
>> NoData Value=0
>> Metadata:
>> STATISTICS_MAXIMUM=771
>> STATISTICS_MEAN=447.14438693812
>> STATISTICS_MINIMUM=150
>> STATISTICS_STDDEV=61.481469436348
>> Band 3 Block=3534x2 Type=UInt16, ColorInterp=Blue
>> Min=116.000 Max=940.000
>> Minimum=116.000, Maximum=940.000, Mean=566.192, StdDev=92.234
>> NoData Value=0
>> Metadata:
>> STATISTICS_MAXIMUM=940
>> STATISTICS_MEAN=566.19172957414
>> STATISTICS_MINIMUM=116
>> STATISTICS_STDDEV=92.234459276031
>> Band 4 Block=3534x2 Type=UInt16, ColorInterp=Undefined
>> NoData Value=0
>>
>>
>> My outputs from r.info for the 4 imported bands are below...Why are the data ranges not correct??
>>
>> +----------------------------------------------------------------------------+
>> | Map: o10APR11180709_M2AS_055552078 Date: Wed Aug 10 14:21:20 2016 |
>> | Mapset: worldview2 Login of Creator: wilbur |
>> | Location: PotteryMound |
>> | DataBase: /home/wilbur/grass_data |
>> | Title: ( o10APR11180709_M2AS_055552078010_01_P001.1 ) |
>> | Timestamp: none |
>> |----------------------------------------------------------------------------|
>> | |
>> | Type of Map: raster Number of Categories: 0 |
>> | Data Type: CELL |
>> | Rows: 2684 |
>> | Columns: 3534 |
>> | Total Cells: 9485256 |
>> | Projection: UTM (zone 13) |
>> | N: 3848648 S: 3843280 Res: 2 |
>> | E: 325694 W: 318626 Res: 2 |
>> | Range of data: min = 133 max = 780 |
>> | |
>> | Data Description: |
>> | generated by r.in.gdal |
>> | |
>> | Comments: |
>> | r.in.gdal -k input="/home/wilbur/Rio_Abajo_Active/isleta_pottery_mou\ |
>> | nd/WorldView/WorldView-2_Image_Data/GeoTIFF_16-Bit/o10APR11180709-M2\ |
>> | AS-055552078010_01_P001.tif" output="o10APR11180709_M2AS_05555207801\ |
>> | 0_01_P001" memory=300 offset=0 |
>> | |
>>
>>
>>
>>
>>
>> +----------------------------------------------------------------------------+
>> | Map: o10APR11180709_M2AS_055552078 Date: Wed Aug 10 14:21:20 2016 |
>> | Mapset: worldview2 Login of Creator: wilbur |
>> | Location: PotteryMound |
>> | DataBase: /home/wilbur/grass_data |
>> | Title: ( o10APR11180709_M2AS_055552078010_01_P001.2 ) |
>> | Timestamp: none |
>> |----------------------------------------------------------------------------|
>> | |
>> | Type of Map: raster Number of Categories: 0 |
>> | Data Type: CELL |
>> | Rows: 2684 |
>> | Columns: 3534 |
>> | Total Cells: 9485256 |
>> | Projection: UTM (zone 13) |
>> | N: 3848648 S: 3843280 Res: 2 |
>> | E: 325694 W: 318626 Res: 2 |
>> | Range of data: min = 119 max = 1193 |
>> | |
>> | Data Description: |
>> | generated by r.in.gdal |
>> | |
>> | Comments: |
>> | r.in.gdal -k input="/home/wilbur/Rio_Abajo_Active/isleta_pottery_mou\ |
>> | nd/WorldView/WorldView-2_Image_Data/GeoTIFF_16-Bit/o10APR11180709-M2\ |
>> | AS-055552078010_01_P001.tif" output="o10APR11180709_M2AS_05555207801\ |
>> | 0_01_P001" memory=300 offset=0 |
>> | |
>>
>> +----------------------------------------------------------------------------+
>> | Map: o10APR11180709_M2AS_055552078 Date: Wed Aug 10 14:21:20 2016 |
>> | Mapset: worldview2 Login of Creator: wilbur |
>> | Location: PotteryMound |
>> | DataBase: /home/wilbur/grass_data |
>> | Title: ( o10APR11180709_M2AS_055552078010_01_P001.3 ) |
>> | Timestamp: none |
>> |----------------------------------------------------------------------------|
>> | |
>> | Type of Map: raster Number of Categories: 0 |
>> | Data Type: CELL |
>> | Rows: 2684 |
>> | Columns: 3534 |
>> | Total Cells: 9485256 |
>> | Projection: UTM (zone 13) |
>> | N: 3848648 S: 3843280 Res: 2 |
>> | E: 325694 W: 318626 Res: 2 |
>> | Range of data: min = 59 max = 1176 |
>> | |
>> | Data Description: |
>> | generated by r.in.gdal |
>> | |
>> | Comments: |
>> | r.in.gdal -k input="/home/wilbur/Rio_Abajo_Active/isleta_pottery_mou\ |
>> | nd/WorldView/WorldView-2_Image_Data/GeoTIFF_16-Bit/o10APR11180709-M2\ |
>> | AS-055552078010_01_P001.tif" output="o10APR11180709_M2AS_05555207801\ |
>> | 0_01_P001" memory=300 offset=0 |
>> | |
>> +----------------------------------------------------------------------------+
>>
>>
>> +----------------------------------------------------------------------------+
>> | Map: o10APR11180709_M2AS_055552078 Date: Wed Aug 10 14:21:20 2016 |
>> | Mapset: worldview2 Login of Creator: wilbur |
>> | Location: PotteryMound |
>> | DataBase: /home/wilbur/grass_data |
>> | Title: ( o10APR11180709_M2AS_055552078010_01_P001.4 ) |
>> | Timestamp: none |
>> |----------------------------------------------------------------------------|
>> | |
>> | Type of Map: raster Number of Categories: 0 |
>> | Data Type: CELL |
>> | Rows: 2684 |
>> | Columns: 3534 |
>> | Total Cells: 9485256 |
>> | Projection: UTM (zone 13) |
>> | N: 3848648 S: 3843280 Res: 2 |
>> | E: 325694 W: 318626 Res: 2 |
>> | Range of data: min = 50 max = 1341 |
>> | |
>> | Data Description: |
>> | generated by r.in.gdal |
>> | |
>> | Comments: |
>> | r.in.gdal -k input="/home/wilbur/Rio_Abajo_Active/isleta_pottery_mou\ |
>> | nd/WorldView/WorldView-2_Image_Data/GeoTIFF_16-Bit/o10APR11180709-M2\ |
>> | AS-055552078010_01_P001.tif" output="o10APR11180709_M2AS_05555207801\ |
>> | 0_01_P001" memory=300 offset=0 |
>> | |
>>
>>
>> Thanks!!
>> _______________________________________________
>> grass-user mailing list
>> grass-user at lists.osgeo.org
>> http://lists.osgeo.org/mailman/listinfo/grass-user
More information about the grass-user
mailing list