This is a response to a query. What to do about a file with no statistics. Using gdalinfo on this SPOT4 scene we get the following:
$ gdalinfo SPOT4_HRVIR2_046-209_2_080911_110453.tif
Driver: GTiff/GeoTIFF
Files: SPOT4_HRVIR2_046-209_2_080911_110453.tif
SPOT4_HRVIR2_046-209_2_080911_110453.tfw
Size is 4440, 4180
Coordinate System is:
PROJCS["SWEREF99 TM",
GEOGCS["SWEREF99",
DATUM["SWEREF99",
SPHEROID["GRS 1980",6378137,298.2572221010002,
AUTHORITY["EPSG","7019"]],
AUTHORITY["EPSG","6619"]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433],
AUTHORITY["EPSG","4619"]],
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",0],
PARAMETER["central_meridian",15],
PARAMETER["scale_factor",0.9996],
PARAMETER["false_easting",500000],
PARAMETER["false_northing",0],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]],
AUTHORITY["EPSG","3006"]]
Origin = (610600.000000000000000,7560400.000000000000000)
Pixel Size = (20.000000000000000,-20.000000000000000)
Metadata:
AREA_OR_POINT=Point
TIFFTAG_COPYRIGHT=Lantm?teriet
TIFFTAG_RESOLUTIONUNIT=1 (unitless)
TIFFTAG_XRESOLUTION=1
TIFFTAG_YRESOLUTION=1
Image Structure Metadata:
INTERLEAVE=BAND
Corner Coordinates:
Upper Left ( 610600.000, 7560400.000) ( 17d39'42.97"E, 68d 8' 9.21"N)
Lower Left ( 610600.000, 7476800.000) ( 17d34'41.45"E, 67d23'12.54"N)
Upper Right ( 699400.000, 7560400.000) ( 19d47'31.12"E, 68d 5'16.24"N)
Lower Right ( 699400.000, 7476800.000) ( 19d38'30.00"E, 67d20'25.87"N)
Center ( 655000.000, 7518600.000) ( 18d40' 6.86"E, 67d44'28.15"N)
Band 1 Block=4440x128 Type=Byte, ColorInterp=Red
Band 2 Block=4440x128 Type=Byte, ColorInterp=Green
Band 3 Block=4440x128 Type=Byte, ColorInterp=Blue
Band 4 Block=4440x128 Type=Byte, ColorInterp=Undefined
Pretty good but no stats. We can force minimum and maximum calculation the the -mm flag
$ gdalinfo -mm SPOT4_HRVIR2_046-209_2_080911_110453.tif
Driver: GTiff/GeoTIFF
Files: SPOT4_HRVIR2_046-209_2_080911_110453.tif
SPOT4_HRVIR2_046-209_2_080911_110453.tfw
Size is 4440, 4180
Coordinate System is:
PROJCS["SWEREF99 TM",
GEOGCS["SWEREF99",
DATUM["SWEREF99",
SPHEROID["GRS 1980",6378137,298.2572221010002,
AUTHORITY["EPSG","7019"]],
AUTHORITY["EPSG","6619"]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433],
AUTHORITY["EPSG","4619"]],
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",0],
PARAMETER["central_meridian",15],
PARAMETER["scale_factor",0.9996],
PARAMETER["false_easting",500000],
PARAMETER["false_northing",0],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]],
AUTHORITY["EPSG","3006"]]
Origin = (610600.000000000000000,7560400.000000000000000)
Pixel Size = (20.000000000000000,-20.000000000000000)
Metadata:
AREA_OR_POINT=Point
TIFFTAG_COPYRIGHT=Lantm?teriet
TIFFTAG_RESOLUTIONUNIT=1 (unitless)
TIFFTAG_XRESOLUTION=1
TIFFTAG_YRESOLUTION=1
Image Structure Metadata:
INTERLEAVE=BAND
Corner Coordinates:
Upper Left ( 610600.000, 7560400.000) ( 17d39'42.97"E, 68d 8' 9.21"N)
Lower Left ( 610600.000, 7476800.000) ( 17d34'41.45"E, 67d23'12.54"N)
Upper Right ( 699400.000, 7560400.000) ( 19d47'31.12"E, 68d 5'16.24"N)
Lower Right ( 699400.000, 7476800.000) ( 19d38'30.00"E, 67d20'25.87"N)
Center ( 655000.000, 7518600.000) ( 18d40' 6.86"E, 67d44'28.15"N)
Band 1 Block=4440x128 Type=Byte, ColorInterp=Red
Computed Min/Max=0.000,255.000
Band 2 Block=4440x128 Type=Byte, ColorInterp=Green
Computed Min/Max=0.000,255.000
Band 3 Block=4440x128 Type=Byte, ColorInterp=Blue
Computed Min/Max=0.000,255.000
Band 4 Block=4440x128 Type=Byte, ColorInterp=Undefined
Computed Min/Max=0.000,255.000
But that is done on the fly every time and isn’t available to any other programme. To get permanently attached statistics we must create a new file with gdal_translate and the -stats flag.
$ gdal_translate -stats SPOT4_HRVIR2_046-209_2_080911_110453.tif output.tif
Input file size is 4440, 4180
0...10...20...30...40...50...60...70...80...90...100 - done.
Then we get this from gdalinfo.
$ gdalinfo output.tif
Driver: GTiff/GeoTIFF
Files: output.tif
Size is 4440, 4180
Coordinate System is:
PROJCS["SWEREF99 TM",
GEOGCS["SWEREF99",
DATUM["SWEREF99",
SPHEROID["GRS 1980",6378137,298.2572221010002,
AUTHORITY["EPSG","7019"]],
AUTHORITY["EPSG","6619"]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433],
AUTHORITY["EPSG","4619"]],
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",0],
PARAMETER["central_meridian",15],
PARAMETER["scale_factor",0.9996],
PARAMETER["false_easting",500000],
PARAMETER["false_northing",0],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]],
AUTHORITY["EPSG","3006"]]
Origin = (610600.000000000000000,7560400.000000000000000)
Pixel Size = (20.000000000000000,-20.000000000000000)
Metadata:
AREA_OR_POINT=Point
TIFFTAG_COPYRIGHT=Lantm?teriet
TIFFTAG_RESOLUTIONUNIT=1 (unitless)
TIFFTAG_XRESOLUTION=1
TIFFTAG_YRESOLUTION=1
Image Structure Metadata:
INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left ( 610600.000, 7560400.000) ( 17d39'42.97"E, 68d 8' 9.21"N)
Lower Left ( 610600.000, 7476800.000) ( 17d34'41.45"E, 67d23'12.54"N)
Upper Right ( 699400.000, 7560400.000) ( 19d47'31.12"E, 68d 5'16.24"N)
Lower Right ( 699400.000, 7476800.000) ( 19d38'30.00"E, 67d20'25.87"N)
Center ( 655000.000, 7518600.000) ( 18d40' 6.86"E, 67d44'28.15"N)
Band 1 Block=4440x1 Type=Byte, ColorInterp=Red
Min=0.000 Max=255.000
Minimum=0.000, Maximum=255.000, Mean=34.493, StdDev=40.536
Mask Flags: PER_DATASET ALPHA
Metadata:
STATISTICS_MAXIMUM=255
STATISTICS_MEAN=34.493035637312
STATISTICS_MINIMUM=0
STATISTICS_STDDEV=40.53643151195
Band 2 Block=4440x1 Type=Byte, ColorInterp=Green
Min=0.000 Max=255.000
Minimum=0.000, Maximum=255.000, Mean=34.832, StdDev=42.388
Mask Flags: PER_DATASET ALPHA
Metadata:
STATISTICS_MAXIMUM=255
STATISTICS_MEAN=34.832346814518
STATISTICS_MINIMUM=0
STATISTICS_STDDEV=42.388441124655
Band 3 Block=4440x1 Type=Byte, ColorInterp=Blue
Min=0.000 Max=255.000
Minimum=0.000, Maximum=255.000, Mean=44.688, StdDev=49.705
Mask Flags: PER_DATASET ALPHA
Metadata:
STATISTICS_MAXIMUM=255
STATISTICS_MEAN=44.688375576533
STATISTICS_MINIMUM=0
STATISTICS_STDDEV=49.70528306057
Band 4 Block=4440x1 Type=Byte, ColorInterp=Alpha
Min=0.000 Max=255.000
Minimum=0.000, Maximum=255.000, Mean=54.075, StdDev=60.230
Metadata:
STATISTICS_MAXIMUM=255
STATISTICS_MEAN=54.075123227294
STATISTICS_MINIMUM=0
STATISTICS_STDDEV=60.229994798131
No comments:
Post a Comment