如何获取geotiff堆栈中的乐队名称？

Question

从像 'NDVI_TS.tif' 这样的文件中的 geotiff（NDVI 时间序列）堆栈中，我想获取单独的波段名称。例如：“频段 086：20190803T004719”。例如，当我将堆栈加载到 QGIS 中时，我可以看到这一点。我需要追溯名称中的日期，在上面的示例中为 2019-08-03。但是，我找不到从 Python 访问它的方法。我可以通过索引访问这些乐队，但这并不能帮助我找到它们来自哪个日期。

from osgeo import gdal

NDVI = gdal.Open('NDVI_TS.tif', gdal.GA_ReadOnly) 
#for example I can get data from band 86 as and array with:
band86 = NDVI.GetRasterBand(86).ReadAsArray()

我觉得应该有一些简单的解决方案，但未能找到。

Answer 1

我不知道这是否有效，但这是我解决这个问题的方法（仅用小图像进行测试）：

bands = {NDVI.GetRasterBand(i).GetDescription(): i for i in range(1, NDVI.RasterCount + 1)}

您会得到一本带有乐队名称（您在 QGIS 中看到的）及其各自索引的字典。

{ ...,
 'Band 086: 20190803T004719': 86,
 ...
}

您可以创建一个小函数来解析乐队名称并根据需要获取日期（没有正确测试）：

import datetime
import re

def string2date(string):
   date_string = re.match('.* (.*)T', string).group(1)
   return datetime.datetime.strptime(date_string, '%Y%m%d').strftime('%Y-%m-%d')

然后您可以应用到之前的字典：

bands = {string2date(NDVI.GetRasterBand(i).GetDescription()): i for i in range(1, NDVI.RasterCount + 1)}

希望它有帮助