以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文件。直接点击,就会运行批处理代码