利用python计算NDVI

一、代码如下

import osgeo
import gdal
from matplotlib import pyplot as plt
import numpy as np
import cv2
from PIL import Image
import cv2


#--1)读取遥感影像：尺寸（宽度和长度）、波段数、投影信息、数组

def NDVI_tif(filename):
    dataset = gdal.Open(filename)  # 打开文件
    im_width = dataset.RasterXSize  # 宽度：栅格影像的列数
    im_height = dataset.RasterYSize  # 长度：栅格影像的行数
    im_bands = dataset.RasterCount  # 波段数
    im_geotrans = dataset.GetGeoTransform()  # 仿射矩阵，左上角