植被覆盖度计算

一、实验内容

使用郑州地区的Landsat-8遥感影像数据，完成郑州地区植被覆盖度计算。

二、实验内容

（1）辐射校正：根据个人需求，将DN值转化为辐射亮度、表观反射率或使用ENVI FLAASH模块进行大气校正，转变为真实反射率；

（2）植被指数提取：可选择NDVI、RVI、EVI等植被参数；

（3）根据线性二元模型进行植被覆盖度计算（植被指数可以选择其他指数替代）

（4）输出计算结果

三、实验步骤

由于已经进行了辐射校正，这里直接进行植被指数提取。

（1）植被指数提取

NDVI：

方法一：使用ENVI内置工具

打开ENVI Classic，加载图像zz_lc8_2017118.tif，Transform > NDVI，设置以下参数：

得到的结果如图：

右键NDVI图层 > Quick Stats，如图所示：5.2512%和94.9844%对应的NDVI值分别为0.076429和0.501905

方法二：使用Band Math

打开ENVI Classic，加载图像zz_lc8_2017118.tif，Basic Tools > Band Math，其中b1为近红外波段，b2为红光波段,设置以下参数：

 

输出结果如图：

如图所示：5.2512%和94.9844%对应的NDVI值分别为0.076429和0.501905

RVI:

Basic Tools > Band Math，其中b1为近红外波段，b2为红光波段，进行以下设置：

输出结果如图：

EVI：

Basic Tools > Band Math，其中b1为近红外波段，b2为红光波段，b3为蓝光波段，输入以下公式并进行以下设置：

2.5 * (float(b1) - float(b2)) / (float(b1) + 6 * float(b2) - 7.5 * float(b3) + 1)

输出结果如图：

（2）植被覆盖度计算

Basic Tools > Band Math，这里基于NDVI进行植被覆盖度计算。

输入以下公式：(b1 lt 0.076429) * 0 + (b1 gt 0.501905) * 1 + (b1 ge 0.076429 and b1 le 0.501905) * ((b1 - 0.076429) / 0.425476)

b1对应：NDVI.tif

结果如图：

（3）保存为ArcMap兼容格式

File → Save As → Image File，输出格式为GeoTIFF，如图所示：

（4）ArcMap制图

将保存的.tiff格式数据加载进ArcMap中，给地图添加图名、图例、比例尺等，结果如图：