0 配置及环境
- 系统: win7 ×64
- python: 3.6.4 ×64
- 包: numpy、GDAL（安装方法可查看GDAL的安装）
1 载入数据并劈开
利用GDAL读取栅格数据(遥感数据.tif):
image = gdal.Open(tifpath) # 读取整个遥感影像，多行多列多层
band = image.GetRasterBand(n) # 得到第n层的数据
band_array = band.ReadAsArray(0,0,lie,hang) # 转化为数组形式
2 叠加数据图层
nir、red 、blue为植被标准假彩色合成，能够目视植被的长势状态，为此我们需要读取并叠加这三个图层：
layer3 = np.array((nir, red, blue)) # nir是第一层， red是第二层，blue是第三层
run.tif_write(out_path, layer3, p, g) # 输出即可
run.tif_write是我自己写的输出遥感图像的包，具体请查阅GDAL的使用