1

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

提取波段 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]（作为最小值、最大值、平均值、标准差）。

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

4

1 回答 1

2

很高兴知道即使在文件中指定了比例/偏移量，GDAL 也不会自动应用它们。

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

它也可以从 GDAL 解析的元数据中检索。这可以通过dataset.GetMetadata()which 将返回一个字典来完成。您也可以调用该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}
于 2020-02-06T08:54:24.370 回答