浮石分享 http://blog.sciencenet.cn/u/floatrock

博文

[转载]MODIS产品MRTtool批处理

已有 5689 次阅读 2010-11-10 16:40 |个人分类:未分类|系统分类:科研笔记|文章来源:转载

MODIS影象产品(Level-1,2,3,4)均是以一种新的ISIN(Integerized Sinusoidal)投影方法进行投影的,而目前大部分软件包均不支持这种投影格式。另外,MODIS数据文件所采取的HDF-EOS(Hierarchical Data Format – Earth Observing System)存储格式目前也没有得到广泛的支持。MRT就是在这样一种环境背景下产生了,它支持这种新的ISIN投影方式和HDF-EOS存储格式。

免费用户注册及下载地址:http://edcdaac.usgs.gov/landdaac/tools/modis/index.asp

   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用户手册(User Guide)。

   在数据处理开始之前,Windows用户需要首先安装Cygwin,并对LINUX语言有一些最基础的语法了解。关于Cygwin的情况及下载请参阅:http://www.cygwin.com/。在确保MRT及Cygwin正确安装好以后,运行Cygwin,利用命令cd work directory(如 cd E:LST)把即将要进行处理的MODIS HDF-EOS文件所在目录设为工作目录。另外,需要注意的是,请确保参数文件(如LST.PRM)事先利用MRT GUI已经得到并与将要进行大批量处理的HDF-EOS文件保存在同一个工作目录下,以免进行处理的时候无法找到正确的参数文件。在利用cd命令把工作目录设置好以后,可以手工输入或粘贴相应的批处理命令(如果您事先已准备好这些命令)。

关于MRT的使用,这里有几点值得注意的地方:

(1)安装MRT的时候,如果用户想利用GUI,则可能需要安装JAVA(http://java.sun.com),否则只能利用命令行的方式进行操作。

(2)利用MRT进行数据处理的时候,请确保文件不要保存在太多的子目录下和尽量不要使用中文目录,以免在MRT导入处理数据的时候可能会出现找不到文件的错误。

(3)在安装Cygwin的时候,用户可以有选择性的进行安装,只要满足可以利用命令行进行操作的目的就可以了。

(4) 请确保参数文件事先已经创建并与要进行批处理的MODIS HDF-EOS文件保存在同一工作目录下,只有这样才有可能数据被正确处理。

(5)由于LINUX语言对大小写非常敏感,在输入批处理命令的时候,请确保文件名正确输入,包括文件名后缀。

http://www.sciencenet.cn/m/user_content.aspx?id=277277 该地址讲述了如何成功的安装MRTtool。

代码如下:

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文件不相同,将其中一个文件拷备到另一个文件目录下使两者一致就没问题了!至于第二个问题,鄙人研究了半天,差点在最后一刻放弃,最终惊人点播才得以成功实现。



http://blog.sciencenet.cn/blog-499915-382391.html

上一篇:[转载]大气辐射传输模型
下一篇:IDL 地形数据阴影影像自动生成

0

发表评论 评论 (0 个评论)

数据加载中...
扫一扫,分享此博文

Archiver|手机版|科学网 ( 京ICP备14006957 )

GMT+8, 2019-11-15 01:00

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部