Python遥感图像处理应用篇(二十一)：Python+GDAL 批量计算遥感图像NDVI指数

Posted 2022-12-07 空中旋转篮球

1.NDVI指数计算公式

之前也写过NDVI批量计算的实现方法，采用的是Acrpy python2.7实现的，这里我们采用GDAL实现NDVI指数的批量计算。

使用数据，多波段遥感图像数据。Landsat05 Colection2 Level2数据。

NDVI=(NIR-RED)/(NIR+RED)  就算之前需要搞清楚多光谱数据中的波段对应关系。

影像数据地表反射率修正系数：0.0000275+（-0.2）

2.实现代码

运行环境：运行环境：windows10  pycharm  python3.7.7  GDAL-3.2.3-cp37-cp37m-win_amd64.whl包等。

import pandas as pd
import numpy as np
from osgeo import gdal
from osgeo import gdal_array
import os
import matplotlib.pyplot as plt
from scipy import stats

class getRaterInfo:
    def __init__(self, filepath):
        gdal.UseExceptions()
        self.filepath=filepath

    gdal.UseExceptions()
    def getRsaterDs(self):
        ds = gdal.Open(filepath)
        return ds

    def getRa