geotiff sanity check

geotiff_info.py

from osgeo import gdal

raster_file = "some_geotiff.tif"

print("gdal.__version__:", gdal.__version__)

print("############ Get STATS using GetStatistics() ############")
ds = gdal.Open(raster_file)
band_count = ds.RasterCount
print("band_count:", band_count)

for band in range(1, band_count + 1):
    stats = ds.GetRasterBand(band).GetStatistics(0, 1)
    print("#############")
    print("band:", band)
    print("min:", stats[0])
    print("max:", stats[1])
    print("mean:", stats[2])
    print("stddev:", stats[3])


print("############ Get STATS using gdal.Info() ############")
result = gdal.Info(raster_file, options=gdal.InfoOptions(stats=True))
print("result:", result)


print("############ Get STATS some other way... ############")
a = ds.ReadAsArray()

for band in range(a.shape[0]):
    band_data = a[band, ...]
    for stat in ['min', 'max', 'mean']:
        print(f"Band {band + 1} {stat}: {getattr(band_data, stat)()}")