准备工作

1.数据准备：进行网站批量下载MRT数据，存储方式为数据–年份–天（例如GLASSFVC数据为8天分辨率，001、009–361）
下载步骤见博文批量下载GLASS数据
2.数据处理用到的是MRT工具，下载安装步骤见博文：MRT工具安装

处理步骤

1.Prm文件制作。将MRT工具打开，导入一景要拼接的全部图幅，调整好参数，如下图所示：
注意：要选择FVC波段；将prm文件输入保存路径后，点run进行运行

运行结果如图所示：
等待finish完成后，关闭。

2.批处理代码撰写。
①新建一个TXT文件，将以下代码复制粘贴进去。

::设定MRT的data文件存放路径，即MRT安装路径(修改1)
set  MRTDATADIR=D:\software\MRT\data

::设定所处理的文件的起始日期(修改2)
set /a DAY=2019001
::设定所处理的文件的结束日期(修改3)
::若同时处理多年可以将多年影像放文件夹改变起始与结束日期，如2001001-2019361
set /a DEADLINE=2019361

::设定数据存放文件夹(修改4)
set workdir =D:\glassdata\NM2019

::设定结果文件存放文件夹(修改5)
set savedir= D:\glassdata\NM2019result\

:start
::判断日期是否超限，否则进行循环计算
if %DAY% leq %DEADLINE% (goto ORDER) else exit

:ORDER


::扫描数据文件夹，将modis数据存放进临时文档中
dir %workdir%*%DAY%.*.hdf/a/b/s > %workdir%MOSAICINPUT.TXT

::启动mrtmosaic，将modis数据进行拼接处理
D:\software\MRT\bin\mrtmosaic.exe -i %workdir%MOSAICINPUT.TXT -s "1 0" -o %workdir%MOSAIC_TMP_%DAY%.hdf 

::启动resample，将拼接后modis数据进行采样及转投影操作
::prm文件要拷贝进所需处理的文件夹之中，注意名称要与下面代码中一致(修改6),*.prm文件名为GLASS2019.prm
D:\software\MRT\bin\resample.exe -p %workdir%GLASS2019.prm -i %workdir%MOSAIC_TMP_%DAY%.hdf -o %savedir%GLASS_%DAY%.tif

::将生成的tif结果剪切到后面的路径中


::删除生成的临时hdf文件
del %workdir%MOSAIC_TMP_%DAY%.hdf

rem del *%DAY%.*.hdf
rem del MOSAIC_TMP_%DAY%.hdf
::下一个数据是加几天？时间分辨率(修改7)
set /a DAY= %DAY% + 8

goto start
cmd.exe

②将TXT文件后缀改成.bat文件。直接点击，就会运行批处理代码。

拼接每一景有进度提醒，如下图所示：

等待处理结束就可以了，觉得有帮助的伙伴可以点赞支持下啊！！！

更多精彩请关注GeoLab 219公众号