Python遥感图像处理应用篇(二十六)：Python+GDAL Sentinel-2数据波段组合

Posted 2022-12-13 空中旋转篮球

1.使用场景描述

批量读取原始压缩包解压后的数据，然后进行波段组合运算。文件列表如下：

 之前写过的是基于Landsat8数据，将所有tif文件放置在一个文件夹种实现的波段组合。

Python遥感图像处理应用篇(二十五)：Python+GDAL 波段组合

文件路径很多，我们需要先读取到需要组合波段的完整路径。

E:\\S2AB_2020_hubei\\S2A_MSIL2A_20200817T025551_N0214_R032_T49RGP_20200817T051650.SAFE\\GRANULE\\L2A_T49RGP_A026913_20200817T030350\\IMG_DATA\\R20m

 选取以上9个波段进行波段组合。

首先读取影像的名称，即合成后的波段名称，我们就以文件夹名称命名，保存在image_list列表中:

image_list=[]  #获取带输出的影像名称
for file in os.listdir(inputpa