利用Python把遥感影像的某几个波段合成

Python合成遥感影像波段实践与问题解析

最新推荐文章于 2023-12-02 21:16:17 发布

原创

最新推荐文章于 2023-12-02 21:16:17 发布 · 5.4k 阅读

56 ·

CC 4.0 BY-SA版权

文章标签：

#python

本文介绍了作者在使用Python处理Landsat8遥感影像时，尝试将nir, red, green三个波段合成假彩色所遇到的问题及解决过程。在尝试过程中，遇到了ValueError和AttributeError，通过调试和理解代码，成功解决了这两个问题，强调了独立理解代码的重要性。" 133947985,7783893,EdgeX v3.0.0 Minnesota 部署实战,"['EdgeX Foundry', 'Docker', 'golang', '后端', '设备管理']

刚刚开始想着手学习Python，决定从最基础的读写遥感影像开始。学习借鉴了网上很多前辈们的经验，自己出现了一些小问题写在这里，算是记录一下。

这是USGS上下载的一景Landsat8影像，地点在青海。

用ENVI截取了扎陵湖的一部分，真彩色合成图如下：

想利用Python把nir,red,green三个波段合成假彩色，代码如下：

import os

import numpy as np
from osgeo import gdal

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.GetProj