功能需求

很多遥感影像的原始数据都是每个波段数据用单独一个文件存储，比如Landsat系列影像数据，而我们在进行遥感解译相关任务时需要多波段数据用一个文件存储，比如rgb波段合并后的真彩色数据，因此对于原始影像需要进行波段合并。

解决思路

思路其实也很简单，只需要将我们想合并的波段依次读取，然后按照自己规定的波段顺序写入同一个文件中即可。主要代码如下所示：

# 写入多波段文件
with rasterio.open(output_path, 'w', **meta) as dst:
    for i, data in enumerate(band_data, start=1):
        dst.write(data, i)  # 按顺序写入波段

未波段合并前如下图所示：

合并后的结果如下图所示：

获取本博客全部代码资料或咨询交流见 博主首页→个人简介
1、本博客相关代码；
2、指导代码运行；
3、同类任务咨询交流；