python-gdal-处理Landsat8影像

之前处理图像的时候都是用ENVI软件进行的，一个个用鼠标点，这次学习了python-arcgis处理图像，真的是太方便了，相融合哪个波段，改几个字符就可以了。
下面展示一些代码。

#用landsat8 数据的三个波段合成一张图像融合的是5、6、4波段，主要表示陆地和水。
from osgeo import gdal
import os
os.chdir(r'F:\烨叶\work\study\ENVI\LC81160272018129LGN00\LC81160272018129LGN00')
band4_fn='LC08_L1TP_116027_20180509_20180517_01_T1_B4.tif'
band5_fn='LC08_L1TP_116027_20180509_20180517_01_T1_B5.tif'
band6_fn='LC08_L1TP_116027_20180509_20180517_01_T1_B6.tif'

in_ds=gdal.Open(band4_fn)
in_band=in_ds.GetRasterBand(1)  #