Mehr

Gdal ReadAsArray() andere Werte als gdal-Bandstatistiken

Gdal ReadAsArray() andere Werte als gdal-Bandstatistiken


Ich bin verwirrt darüber, warum ich zum Beispiel bekomme, wenn ich eine Band mit GDAL öffneband.GetMaximun(): 1000, aber wenn ich es macheReadAsArrayIch bekommemyArray.max(): 1200(das gleiche wenn ich es tue)np.max(meinArray)). Weiß jemand warum dieser Unterschied?


Der Unterschied ist höchstwahrscheinlich auf veraltete Statistiken in Ihrem GeoTiff zurückzuführen.band.GetMaximum()liest den Wert aus den vorberechneten Statistiken, die im GeoTiff enthalten sind.

Wenn Sie die Statistiken aktualisieren möchten, können Sie entwedergdalinfo -statsoder tun Sie es mit den Python-Bindungen, vorausgesetzt, Sie haben Schreibzugriff auf die Datei.

aus osgeo importieren gdal ds = gdal.Open("myGeo.tif", gdal.GA_Update) ds.GetRasterBand(1).ComputeStatistics()

Schau das Video: Read and write raster files with GDAL in Python