Python进行.tif格式数据波段组合；多波段tif数据波段分离

阔别已久，我又来了；啊啊啊……
读了研究生，忙到天昏地暗（哈哈哈，这都是借口，主要是花了一定时间适应研究生生活）

这方面的内容我也主要是应用在遥感领域，我们知道在遥感应用中，大多数的图像数据格式为.tif格式的。进行波段组合使用ArcGIS里面的Composite Bands工具是可以实现的，波段分离也可以使用一些工具或巧妙的方法分开。但当你使用Python程序的话，就可以简化这个软件点击和等待的过程，只要修改2个路径就可以了，岂不乐载？

.tif格式数据波段组合

import os
 
import numpy as np
from osgeo import gdal
from skimage._shared.utils import safe_as_int
 
class IMAGE:
 
    # 读图像文件
    def read_img(self,filename):
 
        dataset = gdal.Open(filename)  # 打开文件
 
        im_width = dataset.RasterXSize  # 栅格矩阵的列数
        im_height = dataset.RasterYSize  # 栅格矩阵的行数
        # im_bands = dataset.RasterCount  # 波段数
        im_geotrans = dataset.GetGeoTransform()  # 仿射矩阵，左上角像素的大地坐标和像素分辨率
        im_proj = dataset.GetProjection()  # 地图投影信息，字符串表示
        im_data = dataset.ReadAsArray(0, 0, im_width, im_height)
 
        del dataset   #关闭对象dataset，释放内存
        # return im_width, im_height, im_proj, im_geotrans, im_data,im_bands
        return  im_proj, im_geotrans, im_data, im_width,im_height
 
    # 遥感影像的存储
    # 写GeoTiff文件
    def write_img(self,filename, im_proj, im_geotrans, im_data):
        # 判断栅格数据的数据类型
        if 'int8' in im_data.dtype.name:
            datatype = gdal.GDT_Byte
        elif 'int16' in im_data.dtype