![a670bff3388e89d3ef8d2b8d8e33b8cc.png](https://i-blog.csdnimg.cn/blog_migrate/9a49a3d96b04678359f0939833cc7b23.jpeg)
引言
在遥感应用领域,植被指数已广泛用于定性和定量评价植被覆盖及其长势。由于植被光谱表现为植被、土壤亮度、环境影响、阴影、土壤颜色和湿度复杂混合反应,而且受大气时空间变化的影响,因此植被指数没有一个普遍的值,其研究经常表明不同的结果。而归一化植被指数(NDVI)是遥感领域中最常用的植被指数之一。
本文提供了一种批量化计算NDVI的Python计算程序,将图像置于某一文件夹后,设定计算结果目录,运行该代码即可。注意:本程序无法处理含有NoData值的遥感图像,解决方法是先利用ENVI中的波段运算方法,去除NoData值。
数据介绍
4波段遥感图像,分别为蓝、绿、红、近红外波段,数据类型为uint16。下图显示了本实验所用的真彩色图像:
![640a22e8495d6e5e565a881d5422093c.png](https://i-blog.csdnimg.cn/blog_migrate/c61d142bd38d79d7d8d3ecb72172f8c6.jpeg)
Python代码实现
import numpy as np
from osgeo import gdal
import os
#定义图像的打开方式
def image_open(img):