6、MRT(Modis Reprojection Tool)批处理MODIS影象(以海河流域为例)

摘自：http://blog.sina.com.cn/s/blog_4edba2dd01000bji.html#articlecommentname

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

MODIS植被覆盖数据产品的5种分类方法

Class 0 1 2 3 4 5 6 7 8 9 10 11 IGBP (Type 1) 水体 常绿针叶林t 常绿阔叶林 落叶针叶林 落叶阔叶林 混交林 郁闭灌丛 开放灌丛 多树草原 稀树草原 草原 永久湿地 UMD (Type 2) LAI/FPAR (Type 3) NPP (Type 4) 水体 常绿针叶林 常绿阔叶林 落叶针叶林 落叶阔叶林 混交林 郁闭灌丛 开放灌丛 多树草原 稀树草原 草原 水体 草地/谷类作物 灌木丛 阔叶作物 稀树草原 阔叶林 针叶林 无植被 城市 水体 常绿针叶林 常绿阔叶林 落叶针叶林 落叶阔叶林 一年生阔叶植被 一年生草地植被 无植被 城市 PFT(Type 5) 水体 常绿针叶林 常绿阔叶林 落叶针叶林 落叶阔叶林 灌木丛 草地 谷类作物 阔叶作物 城市和建设用地 冰和雪 裸地或低植被覆盖地 12 13 14 作物 城市和建设用地 作物/自然植被的镶嵌体 作物 城市和建设用地 15 永久冰雪 16 裸地或低植被覆盖地 裸地或低植被覆盖地 254 未分类 未分类 未分类 未分类 未分类

图形用户界面（GUI）上进行运行，核心部分便是对影象的重采样和镶嵌。关于MRT的安装及使用说明可以参考MRT用户手册（User Guide）。

MODIS HDF-EOS文件名通常包括.h##v##.。它指定了每一幅影象（刈幅）所在的水平方向和垂直方向的位置。ISIN投影把全球按水平从左到右分别进行编码为h00，h01，?，h35，按从上到下分别编码为v00，v01，?，v17。每一景就按照MODIS HDF-EOS文件名里的.h##v##.来确定具体位置。

由于海河流域面积较大，一次需要4景影象才能完全覆盖整个海河流域，所包含的4景所在的位置分别为：.h26v04.，.h26v05.，.h27v04.，.h27v05.（如上图浅绿色说示）。下面以2001年01月01日的海河流域地表温度数据为例讲述一下如何利用MRT对该日地表数据进行处理。所需要的完全覆盖海河流域的4景MODIS影象为：MOD11A1.A2001001.h26v04.004.2003115130411.hdf MOD11A1.A2001001.h26v05.004.2003115134041.hdf MOD11A1.A2001001.h27v04.004.2003115150352.hdf MOD11A1.A2001001.h27v05.004.2003115125635.hdf

首先运行MRT GUI，导入这4景地表温度数据，在可选择的12项波段数据中，选择第一项LST_Day_1km，然后确定相应的输出文件名称及投影方式并填写相应的投影参数和输出像元的大小（可选），如果需要的话，还可以保存用户所输入的参数，最后运行MRT便可自动完成对该日地表温度数据进行重采样和镶嵌的工作。

当然，如果要处理大量类似的MODIS HDF-EOS数据，如本文需要处理海河流域2001年全年的每日地表温度数据（MOD11A1），16天合成的归一化植被指数产品（MOD13A2），16天合成的地表反照率产品（MOD43B3），地表覆盖产品（MOD12Q1），而且如果依旧选择运行这种MRT GUI方式进行数据操作的话，显然会让人感到无法忍受的。这时，可以选择通过命令行输入的方式来操作MRT对数据进行重采样和镶嵌处理。这就是常说的利用MRT对MODIS HDF-EOS文件进行批处理（Batch Processing）。命令行输入最大的好处就是可以大批量的处理相关的MODIS数据。本文以处理海河流域2001年全年的每日地表温度数据为例来简要介绍如何利用MRT命令行对MODIS HDF-EOS数据进行重采样和镶嵌处理。2001年全年的海河流域每日地表温度数据共有1000多景数据影象。利用MRT对数据进行重采样和镶嵌处理的参数与上述LST.PRM相同。在数据处理开始之前，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命令把工作目录设置好以后，可以手工输入或粘贴相应的批处理命令（如果您事先已准备好这些命令）。

利用这个批处理命令，可以完成海河流域2001年01月01日到12月31日的每日地表温度数据的批处理。关于该命令中每条命令所代表的含义，可参阅LINUX的基本语法及MRT的用户手册，在此不再赘述。其他MODIS HDF-EOS数据（如本文中所用到的归一化植被指数产品（MOD13A2），地表反照率产品（MOD43B3），地表覆盖产品（MOD12Q1））也可以类似的按上述方法进行处理，只要稍作修改便可。 关于MRT的使用，这里有几点值得注意的地方：

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

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

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

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

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