Python多波段输出tif图

在Python中，处理多波段的TIFF图像通常需要使用一些地理空间分析库，如 rasterio、gdal 或者 pillow。例如，如果你想要读取一个TIFF文件并将其多个波段合并成一个多通道的图像，可以按照以下步骤操作：

from rasterio import.open

# 打开TIFF文件
with rio.open('your_tif_file.tif') as src:
    # 获取波段数（bands）
    bands = src.count
    
    # 创建一个空的numpy数组来存储所有波段的数据
    multi_band_array = np.empty((src.height, src.width, bands), dtype=src.read(1).dtype)
    
    # 逐个读取并复制每个波段到多通道数组
    for i in range(bands):
        multi_band_array[:, :, i] = src.read(i + 1)  # 加1是因为rasterio默认从1开始计数

# 可以选择将这个数组保存回一个新的TIFF文件，或者用于后续的图像处理
with rio.open('output_multi_band.tif', 'w', driver='GTiff', height=src.height, width=src.width, count=bands, dtype=multi_band_array.dtype, crs=src.crs, transform=src.transform) as dst:
    dst.write(multi_band_array)