一、计算原理
利用NDVI来估算植被覆盖率的原理详见[1] 李苗苗, 吴炳方, 颜长珍, 等.密云水库上游植被覆盖度的遥感估算[J].资源科学,2004 (04): 153-159.
该公式的核心在于确定NDVIsoil和NDVIveg,即裸土地区的NDVI和植被完全覆盖情况下的NDVI。关于这两个值如何选取可以参考下面这篇文章:
[2]甘春英,王兮之,李保生,梁钊雄,李志文,温小浩.连江流域近18年来植被覆盖度变化分析[J].地理科学,2011,31(08):1019-1024.
也有学者根据NDVI累计频率来确定NDVIsoil和NDVIveg的做法,如分别用0.5%、99.5%对应的NDVI值作为NDVIsoil和NDVIveg。
注意:这个0.5%和99.5%完全是经验性的,也有人用5%、95%的。
二、计算植被覆盖率
假设我们经过上面的步骤已经确定NDVIsoil和NDVIveg分别为0.1、0.8,则植被覆盖率FVC的计算公式为:
可以Arcgis中利用栅格计算器工具来实现从NDVI栅格得到植被覆盖率栅格的计算,具体步骤如下:
- 将NDVI栅格添加到Arcmap中
- 执行栅格计算器操作
地图代数表达式:Con((“NDVI.tif”<0.1),0,Con((“NDVI.tif”>=0.8),1,(“NDVI.tif”-0.1)/0.7))
关于地图代数表达式的说明见下图
三、处理前后对比