从 gdal ReadAsArray 返回的 numpy 数组中的值范围是多少？答案

【问题标题】：What is the range of values in the numpy array returned from gdal ReadAsArray?从 gdal ReadAsArray 返回的 numpy 数组中的值范围是多少？
【发布时间】：2019-09-28 08:19:30
【问题描述】：

我下载了 BigEarthNet 数据集，并使用 gdal 在 python 中读取了 TIFF 图像。代码如下。结果数组的值远高于我预期的 0-255 范围，那么值的范围是多少？

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)

这是输出

[[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]]

【问题讨论】：

范围从blue.min()到blue.max()。

标签： python numpy gdal

【解决方案1】：

从外观上看你已经安装了gdal，

gdalinfo <filename>

将为您提供文件上的统计信息，其中包括存储值的最大值和最小值。 tif 中的数据可以表示很多东西，例如 DEM 存储表示栅格波段中高度的值。

实际可能的值范围取决于用于对其进行编码的数据类型。例如 uint16 的范围为 0-65535，我认为我以前在 BIL 格式中使用过。

【讨论】：

谢谢，我不知道 gdalinfo

【解决方案2】：

如果它确实是一个 Numpy 数组，您可以通过以下方式获取最大值/最小值：

numpy.amax(raster_band)
numpy.amin(raster_band)

【讨论】：

...是的，但如果它根本不是一个 numpy 友好的对象，它就不会......我只是说它可能是一个 GDAL 对象，numpy 不知道如何处理与。

【解决方案3】：

From the sentinel-2 user guide:

辐射分辨率通常表示为位数，通常在 8 到 16 位的范围内。辐射分辨率为 MSI 仪器是 12 位的，可以通过范围为 0 到 4095 个潜在光强度值。

【讨论】：