MODIS批量下载及基于MRT的批量格式转换、重采样和拼接处理

denghaijun

MODIS数据为我们提供大量的数据源，特别对于缺资料区域的研究非常重要。最近正好在处理MODIS数据，想把MODIS数据的下载、格式转换、重采样及拼接的一点心得体会和大家分享下。
1 关于下载
网上有多种方法，我给大家分享一下，由中科院遥感所刘金平博士提供的脚本，很好用。进入下面链接
http://jpliu.top/dMODIS.html

图1

a. 数据下载网站如下
https://e4ftl01.cr.usgs.gov/MOLT/
b. 时间尺度为数字，(e.g. 1 or 8)
c. 点击生成脚本，把内容复制到txt中,,如名为**txt“”，修改后缀为“**.bat”
d.把这个bat文件，放到wget目录下，双击，即可。例子：给一个参数填充的例图，用户名和密码要记得填充好。


注意：i)行列号：h24v05，其中5前面的0不能少；ii)最后一个空：MOD13A2.A，这个".A"一定要加。

2.关于格式转换、重采样及拼接
数据下载完后，我们拿到的是hdf格式的数据，涉及到格式转为tif格式（当然hdf格式也是可以直接读取的）；不同的数据分辨率不一样，有的数据500米的，有的数据时1000米的，如果涉及不同分辨率数据的建模，比如基于MODIS的LST数据和EVI/NDVI数据反演土壤水分，就要通过重采样把分辨率统一，这是就涉及到重采样，如何批量处理；又如研究区域较大，会涉及到影像的拼接。一般的常用工具是MRT，但是MODIS官网已经停止更新MRT，推出了新的工具HEG，但是HEG现在在批量处理方面还不是很稳定，正在研究这个HEG的批量拼接，研究明白了，再和大家分享心得体会，这里我们谈谈MRT如果来做这些工作。
a.计算环境，需要安装java，下载链接https://www.java.com/zh_CN/，安装好。
b.安装MRT，官网停止更新，我这里有一个，链接：https://pan.baidu.com/s/1r7kpAAZ7c00UbXmK6RV3-Q ；提取码：86cp
安装步骤，网站的帖子较多，这里就不详细介绍，推荐这个链接https://blog.csdn.net/gisboygogogo/article/details/75784080
c.安装后，即可进行批量处理，

3.关于具体处理步骤
a.先到MRT的安装目录下找到bin文件夹，进入这个文件，找到“Modis Tool”这个快捷方式（安装时会提到这个快捷方式的生成）；
b.打开后，进入到界面，如下图

图2

c.先拼接两幅影响看看效果，选择拼接的变量，输出格式为tif，重采样的方式为“Bilinear”，投影方式根据自己的需要选择，这里介绍下“Albera Equal Area”等积投影参数设置，如下图。

图9

设置后，点击OK

d.output pixel size，即为重采样到多大分辨率，比如500米或1000米

e.把上面设置好的，保存为一个prm格式文件，如下图，方便接下来的批处理。

图3

f.点击run，运行完后，可以查看下处理后的数据，如果满意，即可进行批处理，如果不满意，再重复上面操作修改参数。

注意：批量处理之前，需要修改prm，不管是多个变量还是单个变量，建议都要修改。如下图
修改前：

图4

修改后：

图5

修改“SPECTRAL_SUBSET ”，这个是表示我们需要拼接的变量在HDF文件中位置，即第几个变量，因为一个hdf文件往往有多个变量，若拼接第一个变量，不修改没问题，但是我们可能拼接的不是第一个变量，也可能是多个变量，因此需要修改。

4.运行批处理代码
进入DOS界面，
把工作目录设置为MRT的安装目录，
输入下面代码：

图6

input_directory: 为下载好的MODIS的HDF文件的路径；parameter_directory:为上一步中保存好的prm文件路径（该prm文件一般保存在MRT安装目录的bin文件夹下）。
注意：这个-o parameter_file_directory部分，不要输入。
详见MRT说明文档第40页。

这一行代码输入完成后，回车，输入“matbatch.bat”,回车后，进行批量处理。

图7

看到如下界面即大功告成。

图8

把MRT的说明文档附上。
mrt41_usermanual_032811.pdf (911.44 KB, 下载次数: 66)

2018-11-17 12:33 上传

点击文件名下载附件

叶七

请教老师：做tvdi模型时地表温度数据是只选取LST-day数据吗？

denghaijun

叶七 发表于 2020-3-11 09:40
请教老师：做tvdi模型时地表温度数据是只选取LST-day数据吗？

您好，对这个模型不熟悉哦，不能解答您的问题，抱歉哦。

ZSQ977

每次尝试会出现找不到TepMosaic.bat文件，这个文件不是处理数据的中间文件吗？

小松子

第一个链接打不开了

RS-yiwen

可以请老师分享一下第一个链接，网络上找的都失效了