【MODIS数据处理#12】例七：基于MOD09Q1数据集合成NDVI

系列文章目录： MODIS数据处理

一、MOD09Q1数据介绍

官网介绍链接：https://lpdaac.usgs.gov/products/mod09q1v006/
user guide:https://lpdaac.usgs.gov/documents/306/MOD09_User_Guide_V6.pdf

MOD09Q1（分辨率8day 250m），相对于MODIS产品系列中已经合成好的MOD13Q1 NDVI（分辨率16day 250m）有更高的时间分辨率16-day。

查看MOD09数据集的user guide可以得知Band1和Band2分别对应红光波段®和近红光波段(NIR)。

---

二、处理步骤

2.1 MRT工具预处理

具体步骤可参考以下两篇专栏，这里不赘述
[MODIS数据处理#1]利用MRT工具预处理MODIS数据——以MOD16、MOD13为例
MRT(MODIS Reprojection Tool)安装、影像批量拼接、重投影和格式转换
``

将栅格数据导入Arcmap中

2.2 去除无效值

根据官方文档的说明，栅格的有效值范围是-100~16000，因此我们需要将在此范围之外的像元设为空值，用到的工具是设为空函数，该工具的详细操作说明可以参考这篇的第四节：
[MODIS数据处理#2]常用的Arcmap内置工具（一）


或者也可以用栅格计算器，地图代数表达式可以参考下式

SetNull(("sur_refl_b01.tif" <  -100) | ("sur_refl_b01.tif" > 16000),"sur_refl_b01.tif")

2.3 计算NDVI

经过上面的步骤，我们已经得到了b01(红色波段)和b02(近红外波段)有效范围-100~16000内的栅格数据

接下来可以用栅格计算器来计算NDVI

地图代数表达式为：

# 注意Float中的F要大写
(Float("A2021121.h28v05.b02.tif") - Float("A2021121.h28v05.b01.tif")) / (Float("A2021121.h28v05.b02.tif") + Float("A2021121.h28v05.b01.tif")) 

NDVI的范围应该是在-1~1之间，这里出现了在此范围之外的值，原因可能是由于云或者海的影响，见下图h28v05区块的位置

可以用栅格计算器对数据的范围进行再次筛选

SetNull(("A2021121.h28v05.NDVI.tif" <  - 1) | ("A2021121.h28v05.NDVI.tif"> 1),"A2021121.h28v05.NDVI.tif")

最终得到的NDVI栅格及局部放大细节见下图：

---