J'ai téléchargé le jeu de données BigEarthNet et j'ai lu les images TIFF en python avec gdal. Code ci-dessous. Le tableau résultant avait des valeurs bien au-dessus de la plage 0-255 que j'attendais, alors quelle est la plage des valeurs?
[[284 388 554 ... 325 318 325] [211 213 297 ... 319 300 318] [227 206 245 ... 305 318 332] ... [309 612 920 ... 710 643 554] [259 626 862 ... 654 646 536] [260 608 730 ... 501 629 526]]
Voici la sortie
band1 = imagePath + "/" + img + "/" + img + "_B02.tif" band_ds = gdal.Open(band1, gdal.GA_ReadOnly) raster_band = band_ds.GetRasterBand(1) blue = raster_band.ReadAsArray() print(blue)
3 Réponses :
S'il s'agit vraiment d'un tableau Numpy, vous pouvez obtenir la valeur max / min avec:
numpy.amax(raster_band) numpy.amin(raster_band)
... oui mais si ce n'est pas du tout un objet sympathique, ce ne sera pas du tout ... Je voulais juste dire que ça pourrait être quelque chose comme un objet GDAL que numpy ne saurait pas gérer.
Vu que vous avez déjà installé gdal par les looks,
gdalinfo <filename>
vous donnera les statistiques sur votre fichier qui incluront le max et le min des valeurs stockées. Les données d'un tif peuvent représenter beaucoup de choses, par exemple un DEM stocke des valeurs représentant la hauteur dans la bande raster.
La plage de valeurs réelle possible dépend du type de données utilisé pour l'encoder. Par exemple, uint16 a une plage de 0-65535 que je pense avoir utilisé auparavant au format BIL.
Merci, je ne connaissais pas gdalinfo
À partir du guide de l'utilisateur sentinel-2 :
La résolution radiométrique est généralement exprimée en nombre de bits, généralement dans la plage de 8 à 16 bits. La résolution radiométrique de l'instrument MSI est 12 bits, ce qui permet d'acquérir l'image plage de 0 à 4095 valeurs d'intensité lumineuse potentielle.
La plage va de
blue.min ()àblue.max ().