文章目录

• 一、数据下载
• • 1.数据下载基本操作
  • 2.如何选择区域li>
  • 3.下载得到的数据是什么样子li>
• 二、数据处理
• • 1.处理分类
  • 2.预处理
  • 3.计算处理

一、数据下载

1.数据下载基本操作

数据下载 址：https://ladsweb.modaps.eosdis.nasa.gov/
也可以通过代码实现批量下载，这里没有实验，仅在官 上按照年份进行下载因此下载速度相对较慢。
打开 址，点击“Find data”，进入界面，按照顺序依次选择：products、time、location、files最后即可获取数据。

2.如何选择区域h2>

这里由于研究需要，我们选择的中国区域，MODIS在中国区域的行列 如下所示。
引自：https://blog.csdn.net/weixin_48807049/article/details/124711496

二、数据处理

1.处理分类

对于MODIS数据的处理大致分为两类：预处理及后期计算处理。我在处理时，预处理主要采用MRT工具，后期计算时采用ArcGIS软件。

2.预处理

（1）MRT软件安装
参考博客：https://blog.csdn.net/gisboygogogo/article/details/75784080，即可完成安装操作。
（2）批量对HDF实现重投影、采样、拼接处理（生成tiff格式）
批处理主要有两种方式：代码与直接在cmd中输入，这里采用第二种。在批处理之前，要先采用MRT软件生成prm文件，prm文件相当于包含了你所设置的投影、采样及拼接等操作。在GUI界面上设置好参数后，点击Sava Parameter frofile设置输出prm文件位置及文件名即可，不需要run。

（具体操作可参考：https://blog.csdn.net/weixin_44588773/article/details/122410487）

完成后关闭MRT，直接在ModisTool.bat所在的bin文件夹上方的地址栏中输入cmd，用到的命令如下：

（命令中hdf_dictionary为需要处理的modis数据所在目录，prm_dictionary为刚才保存的prm参数文件，output_dictionary为处理后的影像的输出目录）

上一行命令输入完成后回车，输入MRTBatch.bat命令，回车即开始批处理工作。

代码如下（示例）：

（2）设空值
由于MOD16数据集对于城市、冻土冰雪、裸地、水体这些区域仅用特殊值标记（32762—32766），不计算ET值，在此对这些区域进行排除。32762—32766代表：

代码如下（示例）：

（3）时间分辨率转化（将8天转为1个月）
由于MODIS时间分辨率为8天，有时研究需要月尺度的时间分辨率数据，因此需要转化，具体的批量设置空值的操作步骤可以参考：https://blog.csdn.net/qq_37948866/article/details/100098259
除此之外，也要进行单位换算，Y=A*Dn+B，其中Dn为原始影像像元值、A为比例系数、B为偏移量、Y为真实值。以MOD16A2数据为例：，Dn单位为kg/m2/8d，乘于比例因子A（0.1）后变为mm/8d。因此，对于8-day分辨率的ET栅格数据，ET(单位mm/d) = 像元值 * 比例因子 / 8。也可以换算为月尺度的数据，但这里要注意的是最后一景（平年时第361天）的数据代表的不是8天的ET数据，而是5天的ET数据。

以上就是处理MOD蒸散数据的大致流程。

文章知识点与官方知识档案匹配，可进一步学习相关知识Python入门技能树首页概览215681 人正在系统学习中