刚刚开始想着手学习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 = da