[GRASS-user] Data ranges of geotiff not preserved after import to GRASS
Markus Metz
markus.metz.giswork at gmail.com
Sat Aug 13 11:52:59 PDT 2016
On Sat, Aug 13, 2016 at 12:01 AM, William Hudspeth
<bhudspeth at edac.unm.edu> wrote:
> Thank you. That solution worked....
The real problem is that PCI Geomatica writes wrong band statistics to
the metadata. The solution is to request PCI Geomatica to write out
correct values.
>
>
>
> 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