如何在 GDAL Python 中获取带数组的正确比例?

Question

我是Python中使用GDAL的新手，我正在尝试使用它从Sentinel 2 SAFE-products中检索波段数据。我设法提取了波段数组，但我无法正确缩放它。

这将提取波段 4 的未缩放数组：

import gdal
product_path = "S2B_MSIL2A_20200124T101219_N0213_R022_T33UUU_20200124T121752.SAFE"
dataset = gdal.Open(product_path + "MTD_MSIL2A.xml")
bands10m_path = dataset.GetSubDatasets()[0][0]
bands10m_dataset = gdal.Open(bands10m_path)
b4_band = bands10m_dataset.GetRasterBand(1)
b4_array = b4_band.ReadArray()

到目前为止一切顺利，但是数组的数据类型是uint16，值范围从0到16896.

b4_band.GetMinimum()和b4_band.GetMaximum()都返回None。

b4_band.GetStatistics(True,True)返回[0.0, 2829.0, 347.05880000000104, 334.8397839901348](作为最小值，最大值，平均值，标准dev(。

这是否有助于我以某种方式提取正确的比例？我一无所知...

Answer 1

请注意，即使在文件中指定了缩放/偏移，GDAL 也不会自动应用它们。

对于 Sentinel 2，它们不在文件中指定，而是在元数据 (xml( 中指定。您可以使用文本编辑器查看您在示例中使用的 xml 文件。并按照@Val的建议搜索"QUANTIFICATION_VALUE"。

也可以从 GDAL 分析的元数据中检索它。这可以通过返回字典的dataset.GetMetadata()来完成。您也可以调用gdal.Info实用程序，这两种方法如下所示。

import gdal
archive = 'S2A_MSIL2A_20200126T082221_N0213_R121_T34HDH_20200126T110103.zip'
# Use a dataset
ds = gdal.Open(archive)
meta = ds.GetMetadata()
ds = None
# Alternatively use gdal.Info
r = gdal.Info(archive, format='json')
meta = r['metadata']['']

您可以过滤掉相关值，并将它们从字符串转换为浮点数

，如下所示：

{k: float(v) for k,v in meta.items() if k.endswith('_QUANTIFICATION_VALUE')}
# Result in:
{'AOT_QUANTIFICATION_VALUE': 1000.0,
'BOA_QUANTIFICATION_VALUE': 10000.0,
'WVP_QUANTIFICATION_VALUE': 1000.0}