【MODIS数据处理#10】例五:用NDVI估算植被覆盖度
文章目录一、功能介绍二、脚本代码三、工具参数四、工具界面系列文章目录:一、功能介绍Con((“A2000M06.NDVI.tif”<=0.05),0,Con((“A2000M06.NDVI.tif”>=0.7),1,(“A2000M06.NDVI.tif”-0.05)/0.65))二、脚本代码#!/usr/bin/python# -*- coding: UTF-8 -*-三、工具参数四
一、计算原理
利用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))
关于地图代数表达式的说明见下图
三、处理前后对比


更多推荐



所有评论(0)