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

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

1.使用场景描述

之前使用arcpy python2.7写了一篇进行遥感数据波段组合计算的博文，使用场景是将所有单波段数据放在一个文件夹中（如下图），文件名称前缀一样，后缀波段编号用来区分不同的波段，数据按照名称默认排序，读取的时候也是按照排列顺序读取。

Python遥感图像处理基础篇(一)：arcpy遥感图像波段组合批量处理_空中旋转篮球的博客-CSDN博客_python处理遥感

 2.实现思路

（1）获取待输出波段组合的影像列表，如上图，有14个文件，实际上是2景影像，我们需要获取这2景影像的前缀名称，并且存入列表中。实现代码如下：

file_list = os.listdir(inputpath)
#先获取数据列表中的文件名称，每一景数据（7个波段）为一个数据
image_list=[]
for readPath in file_list:
    # 获取文件名称相同的前缀作为输出文件的名称
    filename=readPath[:-7]
    print(filename)
    if filename not in image_list:
        image_list.append(filename)
print(image_list)

输出结果：

[\'LC08_L