ArcMap计可见光算植被指数（VDVI）

NDVI  （Normalized Difference Vegetation Index,归一化植被指数，标准差异植被指数)

表达式：NDVI=(p(nir)-p(red))/(p(nir)+p(red)) 其中p（nir）表示近红外波段，p（red）表示可见光红波段，它和植物的蒸腾作用、太阳光的截取、光合作用以及地表净初生产力（NPP）等密切相关。

1、NDVI的应用：检测植被生长状态、植被覆盖度和消除部分辐射误差等；

2、-1<=NDVI<=1，负值表示地面覆盖为云、水、雪等，对可见光高反射；0表示有岩石或裸土等，NIR和R近似相等；正值，表示有植被覆盖，且随覆盖度增大而增大；

3、NDVI的局限性表现在，用非线性拉伸的方式增强了NIR和R的反射率的对比度。对于同一幅图象，分别求RVI和NDVI时会发现，RVI值增加的速度高于NDVI增加速度，即NDVI对高植被区具有较低的灵敏度；

4、NDVI能反映出植物冠层的背景影响，如土壤、潮湿地面、雪、枯叶、粗超度等，且与植被覆盖有关。

以TM为例：

NDVI=（B4-B3）/(B4+B3)

1.arcgis打开TM，B4，B3两波段数据

2.空间分析里面有个栅格计算器的工具（spatial analyst-raster calculator）

输入（B4-B3）/(B4+B3)

3.点击evaluate（执行）

 

4.查看结果

这个时候你会发现NDVI的值范围是正确的（-1,1），但是它是不连续变化的。根据研究区以及经验判断这个结果是错的。

分析原因，寻找解决的办法。

——数据类型的为整形

需要修改栅格数据的类型，改为float pointing 浮点型！

修改办法：

栅格计算器里面有float工具如下图：

转换4,3两个波段的数据type为float point

再次运算。

查看结果：

可见光VDVI

(2 * "DSC01878.JPG - Band_2" / ("DSC01878.JPG - Band_1" + "DSC01878.JPG - Band_2" + "DSC01878.JPG - Band_3") -  "DSC01878.JPG - Band_1" / ("DSC01878.JPG - Band_1" + "DSC01878.JPG - Band_2" + "DSC01878.JPG - Band_3") -
 "DSC01878.JPG - Band_3" / ("DSC01878.JPG - Band_1" + "DSC01878.JPG - Band_2" + "DSC01878.JPG - Band_3") )-(1.4*"DSC01878.JPG - Band_1" / ("DSC01878.JPG - Band_1" + "DSC01878.JPG - Band_2" + "DSC01878.JPG - Band_3")
 -"DSC01878.JPG - Band_2" / ("DSC01878.JPG - Band_1" + "DSC01878.JPG - Band_2" + "DSC01878.JPG - Band_3") )