GLASS数据预处理/MRT工具批量处理MODIS数据

以GLASS17E01为例
NPP
时间分辨率：8天
空间分辨率：500m
001-361

1 数据准备

2 Prm文件制作

因为批处理，所以不要点运行

3 批处理代码

使用MODIS数据的处理过程有问题，所以不得逐年处理

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

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

::设定数据存放文件夹(修改4)
set workdir =H:\NPP\2001

::设定结果文件存放文件夹(修改5)
set savedir= H:\NPP_result\

: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
C:\MRT\bin\resample.exe -p %workdir%GLASS2001.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文件。直接点击，就会运行批处理代码

参考博客:
[1] GLASS文件处理
[2]MODIS批处理