Python-GDAL 栅格数据处理

将独立的栅格波段合并为一张TIF影像

#objective：将独立的栅格数据合并为一张影像，并保存输出。
#导入一个可以操作系统的库
import os
#导入GDAl库
from osgeo import gdal
#从系统中读取文件的目录，分别读取要和合并的波段
ds = os.chdir(r'G:\TestData22')
b1 = 'LC81270392015294LGN00_B1.TIF'
b2 = 'LC81270392015294LGN00_B2.tif'
b3 = 'LC81270392015294LGN00_B3.tif'
#使用gdal库中的Open（）方法打开一个波段，随后将打开的波段定义为输入的第一个波段
in_ds = gdal.Open(b1)
in_band = in_ds.GetRasterBand(1)
#GetDriverByName官网解释 Fetch a driver based on the short name.
#其实就是根据输入的GTiff确定类型，返回所支持类型的初始化结果（实例）。