Python遥感图像处理应用篇(十三)：Arcpy批量计算绿度归一化植被指数GNDVI

Posted 2022-06-30 空中旋转篮球

 1.绿度归一化植被指数GNDVI

公式：GNDVI=TM4-TM2/TM4+TM2（Landst5-7）   

RVI=TM5-TM3/TM5+TM3（Landst8-9）

GNDVI=NIR-G/NIR+G

NIR是近红外波段，R是绿光波段的反射率。

2.Python代码实现

运行环境 arcpy python2.7  ILDE GUI运行。使用数据 Landsat8 OLI 未修正代码

"""
@soderayer
绿度归一化植被指数
GNDVI=TM4-TM2/TM4+TM2   GNDVI=NIR-G/NIR+G
"""
# -*- coding: utf-8 -*-
import os, arcpy, datetime
from arcpy.sa import *

# 存放文件夹的位置
Red_path = r\'G:\\landsat8SI\\B3\'  # 绿光波段存放文件夹
NIR_path = r\'G:\\landsat8SI\\B5\'  # 近红外波段存放文件夹
path = r\'G:\\landsat8SI\\GNDVIOUT\'  # 输出结果存放文件夹

# set the intermediate data folder
intermediateDataPath = path + "\\\\" + "IntermediateData"
# set result data folder
resultDataPath = path + "\\\\" + "Result"

# determine if the folder exists
if os.path.exists(inte