对于多个需要拼接的MODIS数据可以采用MRT进行投影,拼接。
所需要的工具有MRT软件,Cygwin软件,以及Java软件。(注意:MRT和Cygwin软件要安装在c盘的根目录下,Java安装在program files文件夹下)
1、首先安装Java软件,安装后再安装MRT软件,安装MRT过程中指定好Java的位置。确认MRT的bin文件夹下的ModisTool能够正常运行,并能够完成具体的投影、拼接。(下图中给定的是XJ流域的MRT设置参数)
![[转载]MODIS <wbr>MRT <wbr>批处理方法 [转载]MODIS <wbr>MRT <wbr>批处理方法](https://s2.51cto.com/images/blog/202410/05205323_670136c38568f79276.gif)
![[转载]MODIS <wbr>MRT <wbr>批处理方法 [转载]MODIS <wbr>MRT <wbr>批处理方法](https://s2.51cto.com/images/blog/202410/05205323_670136c38568f79276.gif)
在这个过程中,将
Xiangjiang.prm文件进行保存。
2、 当MRT软件能够正常运行后,安装Cygwin软件,在安装完Cygwin后。首先检查一下MRT的安装环境是否正确。
![[转载]MODIS <wbr>MRT <wbr>批处理方法 [转载]MODIS <wbr>MRT <wbr>批处理方法](https://s2.51cto.com/images/blog/202410/05205323_670136c38568f79276.gif)
注意: 如果是按照根目录安装的话,这一步可以不再检查。
完成之后,将Cygwin安装文件夹内的bin文件夹的
cygwin1.dill
复制到MRT安装文件夹的bin文件夹中,而把MRT的bin文件夹中的
mrtmosaic.exe
和
resample.exe
两个文件复制到Cygwin安装文件夹内的bin文件夹中。
到这,基本上完成了软件的设置和安装。
3、批处理文件的设置。
用记事本打开保存的
Xiangjiang.prm文件
![[转载]MODIS <wbr>MRT <wbr>批处理方法 [转载]MODIS <wbr>MRT <wbr>批处理方法](https://s2.51cto.com/images/blog/202410/05205323_670136c38568f79276.gif)
注意将红框中的位置及文件信息进行更改。其中可以2000改为2001等进行2001年的数据处理。而(0 1 0 0 0 0 0 0 0 0 0 0 )中1代表的是提取是哪个波段的数据。如果1的位置在第一位置那就提取的是第一个波段。
![[转载]MODIS <wbr>MRT <wbr>批处理方法 [转载]MODIS <wbr>MRT <wbr>批处理方法](https://s2.51cto.com/images/blog/202410/05205323_670136c38568f79276.gif)
注意将红框中的位置及文件信息进行更改。将更改后的文件另存到2000年的需要批处理的文件夹内。如果是需要处理2001年的数据则将2000改为2001年并该文件存入2001年的数据文件夹中,以此类推可以处理其他年份的数据。
4、批处理程序的编写
新建一个txt文本文件,打开文件,将以下代码进行编写。同时注意红线中的文件名和需要修改的地方。
![[转载]MODIS <wbr>MRT <wbr>批处理方法 [转载]MODIS <wbr>MRT <wbr>批处理方法](https://s2.51cto.com/images/blog/202410/05205323_670136c38568f79276.gif)
cd cd f:/Xiangjiang MODIS/2001
注意目录间斜杠号“/”的方向必须如上所述。(该步骤是将采用cygwin软件工作目录指导到需要处理的文件夹位置。)需要特别注意的是如果文件夹命名中带有空格则一定要注意空格中的否则会出现错误。
完成输入回车,再将下面的代码复制到cygwin的命令栏中(注意红色的字体是解释部分,不需要输入)
day=2001001(开始文件的时间)
while [ $day -le 2001353 ](结束文件的时间)
do
for trl in MOD13Q1.A$day.*.hdf
do
if [ -f $trl ]
then
HDFFILES=$(ls MOD13Q1.A$day.*.hdf)
echo $HDFFILES > MOSAICINPUT.TXT
MRTMOSAIC -i MOSAICINPUT.TXT -s "0 1 0 0 0 0 0 0 0 0 0 0" -o MOSAIC_TMP_$day.hdf
RESAMPLE -p Xiangjiang.PRM -i MOSAIC_TMP_$day.hdf -o MOD13Q1_$day.tif (Xiangjiang.PRM 之前保存的prm文件)
rm MOSAIC_TMP_$day.hdf
(该步骤是将运行过程中的中间文件进行删除,如果要删除源文件则用代码rm MOD13Q1.A$day.*.hdf)
let "day=day+16"
(这个16天产品的时间循环间隔为16天,其它modis产品按照实际时间更改)
break
else
let "day=day+16"
break
fi
done
done按回车键即开始运行。
















