MODIS植被覆盖度

最近在做植被覆盖度，网上查了很多的相关资料，经过一段时间的摸索，终于明白了其中的一些逻辑思路，以下是根据网上的资料和个人的理解做的一些整理和总结，不知道是否正确，如有不对的地方还请多多指出。

一、数据下载

MODIS数据下载详情：https://blog.csdn.net/wz_xiaowangzi/article/details/89143038

二、MRT投影转换（以MOD13Q1为例）

1、下载和安装工具MRT：可去百度搜索下载

下载的MRT文件包括：

安转MRT工具的前提是首先安转java，关于MRT的具体安装过程网上有详细教程。

2、使用MRT对MODIS进行投影转换

原因：由于MODIS使用的正弦曲线投影，所以需要投影到平时比较常用的坐标系统。

步骤：①打开MRT界面

②输入要进行投影转换的MODIS文件

③选择提取的波段 

④输出投影设置 

Albers Equal Area投影参数的设置 

（关于投影参数的设置可在网上查找资料，有相关的介绍） 

⑤设置输出位置、名字和分辨率，其他默认，最后点击run即可 （输出的结果为你所提取波段的文件）

三、植被覆盖度计算

利用NDVI计算植被覆盖度

①裁剪出研究区的NDVI（subset  data from rois）

②将NDVI恢复到正常区间（MRT在投影转换的过程中将NDVI扩大了10000倍）：使用band math输入公式b1*0.0001即可运算，（或者是将负值去掉，(b1 lt 0)*0+(b1 ge 0)*(b1*0.0001)，此时NDVI的区间为（-1,1）之间）

③计算置信度（compute statistics）

选择计算的NDVI图层，点击ok

选择掩膜（研究区范围掩膜，目的是只计算研究区范围内的值） 

点击OK，出现NDVI统计视图：

④根据置信度植被覆盖度的计算公式：

(b1 lt NDVI最小值)*0+(b1 gt NDVI最大值)*1+(b1 ge NDVI最小值 and b1 le NDVI最大值)*((b1-NDVI最小值)/(NDVI最大值-NDVI最小值))

即可计算植被覆盖度。 

参考：

MODIS NDVI数据处理相关问题：https://blog.csdn.net/lucky51222/article/details/24939913

MODIS数据计算各年植被覆盖度：http://blog.sina.com.cn/s/blog_935d88e40102vgo7.html