摘要: MODIS影象产品(Level-1,2,3,4)均是以一种新的ISIN(Integerized Sinusoidal)投影方法进行投影的,而目前大部分软件包均不支持这种投影格式。另外,MODIS数据文件所采取的HDF-EOS (Hierarchical Dat...
MODIS影象产品(Level-1,2,3,4)均是以一种新的ISIN(Integerized Sinusoidal)投影方法进行投影的,而目前大部分软件包均不支持这种投影格式。另外,MODIS数据文件所采取的HDF-EOS (Hierarchical Data Format – Earth Observing System)存储格式目前也没有得到广泛的支持。MRT就是在这样一种环境背景下产生了,它支持这种新的ISIN投影方式和HDF-EOS存储格式。
MRT是一种针对MODIS数据的处理工具,可以帮助用户把MODIS影象 (Level-2G,Level-3,and Level-4 land data products)重新投影到更为标准的地图投影,而且可以选择影象中的空间子集(spatial subsetting)和波段子集(spectral subsetting)进行投影转换。软件输出格式为raw binary,GeoTIFF(这两种数据格式为大多数软件所支持)和HDF-EOS,而且可以在好几种系统平台上进行运行,包括:Sun Solaris workstations,SGI IRIX workstations,Linux和Microsoft Windows。
MRT可以通过命令行或在MRT图形用户界面(GUI)上进行运行,核心部分便是对影象的重采样和镶嵌。关于MRT的安装及使用说明可以参考MRT用户手册。 在数据处理开始之前,Windows用户需要首先安装Cygwin,并对LINUX语言有一些最基础的语法了解。在确保MRT及Cygwin正确安装好以后,运行Cygwin,利用命令cd work directory (如 cd E:LST)把即将要进行处理的MODIS HDF-EOS文件所在目录设为工作目录。另外,需要注意的是,请确保参数文件(如LST.PRM) 事先利用MRT GUI已经得到并与将要进行大批量处理的HDF-EOS文件保存在同一个工作目录下,以免进行处理的时候无法找到正确的参数文件。在利用cd命令把工作目录设置好以后,可以手工输入或粘贴相应的批处理命令(如果您事先已准备好这些命令)。
关于MRT的使用,这里有几点值得注意的地方:
-
安装MRT的时候,如果用户想利用GUI,则可能需要安装JAVA,否则只能利用命令行的方式进行操作。
-
利用MRT进行数据处理的时候,请确保文件不要保存在太多的子目录下和尽量不要使用中文目录,以免在MRT导入处理数据的时候可能会出现找不到文件的错误。
-
在安装Cygwin的时候,用户可以有选择性的进行安装,只要满足可以利用命令行进行操作的目的就可以了。
-
请确保参数文件事先已经创建并与要进行批处理的MODIS HDF-EOS文件保存在同一工作目录下,只有这样才有可能数据被正确处理。
-
由于LINUX语言对大小写非常敏感,在输入批处理命令的时候,请确保文件名正确输入,包括文件名后缀。代码如下:
day=2001001
while [ $day -le 2001365 ]
do
for trl in MOD11A1.A$day.*.hdf
do
if [ -f $trl ]
then
HDFFILES=$(ls MOD11A1.A$day.*.hdf)
echo $HDFFILES > MOSAICINPUT.TXT
MRTMOSAIC -i MOSAICINPUT.TXT -s "1 0 0 0 0 0 0 0 0 0 0" -o MOSAIC_TMP_$day.hdf
RESAMPLE -p LST.PRM -i MOSAIC_TMP_$day.hdf -o MOD11A1_$day.tif
rm MOSAIC_TMP_$day.hdf
rm MOD11A1.A$day.*.hdf
let "day=day+1"
break
else
let "day=day+1"
break
fi
done
done
关键点:在MODISTOOL成功安装后,用上边的代码在Cygwin中运行可能会不能成功,一方面的原因有可能是参数文件没有配置正确,另一个可能的原因就是Cygwin目录下和MODIStool目录下的cygwin1.dll文件不相同,将其中一个文件拷备到另一个文件目录下使两者一致就没问题了!