Python处理遥感影像多波段生成图像

刚开始接触Python的学习，计划从最基本的遥感影像读写操作入手。在学习过程中参考了网络上许多先行者的经验分享，期间遇到了一些小问题，特此记录下来以备后续查阅。

这是一幅从USGS平台获取的Landsat8卫星影像数据，拍摄区域位于青海省。

通过ENVI软件对扎陵湖区域进行了图像裁剪，所获取的真彩色合成图像如下所示：

希望借助Python实现将nir、red、green三个波段组合生成假彩色图像，所采用的代码如下：

 import os

    
  
    
 import numpy as np
    
 from osgeo import gdal
    
  
    
 cla