Suite

Gdal ReadAsArray() des valeurs différentes des statistiques de bande gdal

Gdal ReadAsArray() des valeurs différentes des statistiques de bande gdal


Je ne comprends pas pourquoi si j'ouvre un groupe avec GDAL, j'obtiens, par exemple,bande.GetMaximun() : 1000, mais si je le faisReadAsArrayJe reçoismonTableau.max() : 1200(le même si je faisnp.max(monTableau)). Est-ce que quelqu'un sait pourquoi cette différence ?


La différence est probablement due à des statistiques obsolètes dans votre GeoTiff.bande.GetMaximum()lira la valeur des statistiques précalculées, qui sont incluses dans le GeoTiff.

Si vous souhaitez mettre à jour les statistiques, vous pouvez soit utilisergdalinfo -statsou faites-le avec les liaisons Python, à condition que vous ayez un accès en écriture au fichier.

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

Voir la vidéo: install gdal, rasterio and fiona in python