MODIS数据下载与处理

本文格式为Word版，下载可任意编辑——MODIS数据下载与处理MODIS数据(1B、产品、常用处理软件)免费下载

1、MODIS1B数据下载

L1B数据下载地址：/data/search.html开启网页，选中需要的数据源类型

注意：其中，MOD03数据是用于对1KM,QKM,HKM数据进行几何改正所用。别忘掉下载。其中，日期类型为：月/日/年时:分:秒

其中，网页中显示的时间为UTC时间，换算为北京时间为：UTC时间=北京时间-8小时。因此，要获得1月16日的数据则范围为：01/15/2023/16:00:00~01/16/2023/16:00:00

在‘spatialselection’选项中选择“latitude/longtitude〞,按经纬度形式选择影像范围。

点击’search’查到需要的数据：

勾选需要的数据，点击‘orderfilesnow’，输入你接收信息的邮箱，点’order’开始订购该数据。(假使要探寻多天数据，可以选‘addfilestoshoppingcart’继续探寻其他日期的数据。

所订购数据的存放位置信息：

点击‘Data->TrackOrders’可以查看所有已订购的数据的状态。

假使’state’显示‘avalable’即可开始下载。

ftp:

username:anonymous

password:rainaleaf@163.com（自己申请数据的邮箱）

一般使用FTP下载工具下载。首先新建站点：

点击‘连接’。则在右边的框中会显示所有数据，找到自己数据所在的文件夹，并拖到左边的框中，开始下载数据。

OK！

2、MODIS产品下载地址

/api//u515127/u515368/u515406/u515421/u515529/u516066/u516091/u516110/u516150/u516307/u517021/u517029(要装IE7支持)

或已知所下载文件名的modis09~17数据可以从FTP上下载：ftp:/

tarra的数据放在MOLT里。Aqua数据放在MOLA。上述的探寻网站（）（假使遇到能查到数据但是下载不了的状况，也可以在FTP中通过查找与该数据一致日期同轨道号来下载。）在网站查到该数据的好处是能查看数据的覆

盖区域。可以查看覆盖范围。

假使：MOD11A1.A2023288.h28v06.005.2023290030125.hdf从：

/MOLT/MOD11A1.005/2023.10.14/查找到2023-10-14轨道号为H26V05,就从里面的文件夹里选出来下载就可以了。一般使用FTP下载工具下载。

注：数据的轨道号：

’opengeolocationfile’开启与该数据想配套的03级数据。

4、设置投影参数：

点击‘editProjectionParameters?’修改参数：

然后点击〞run〞,几何改正开始。

二、DN值转反射率（scale_offset改正）

1、几何改正完成后即可用ENVI开启：

2、从basictools->preprocessing->Data-specificUtilities->viewHDFattributes,从未经几何改正前的数据查看该数据scale_offset值。

3、改正gain_offset值:从basictools->preprocessing->GeneralPurposeUtilities->ApplyGainandOffset.

输入相应的scale和offset值，并保存图像。5、拼接

在弹出来的对话框中，import->〞importfilesandproporties…〞选择要进行拼接的影像,

在“datavaluetoignore〞中输入要忽略的背景的值。这样做能使两影像之间的拼接时影像外的区域为透明显示。显示结果，如下，假使选择了import->〞importfiles〞可以通过右击影像列表中的影像，选〞editentry…〞来修改要忽略的背风景的值，效果是一样的。

然后选file->apply输出拼接的结果。

5、MODIS产品处理之MRT重投影参数设置

假使下载的是多幅图像该软件可以进行自动拼接：方法是，输入多幅影像如下图：

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

摘自：/s/blog_4edba2dd01000bji.html#articlecommentname

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

免费用户注册及下载地址：/landdaac/tools/modis/index.asp。MRT是一种针对MODIS数据的处理工具，可以帮助用户把MODIS影象(Level-2G，Level-3，andLevel-4landdataproducts)重新投影到更为标准的地图投影，而且可以选择影象中的空间子集（spatialsubsetting）和波段子集（spectralsubsetting）进行投影转换。软件输出格式为rawbinary，GeoTIFF（这两种数据格式为大多数软件所支持）和HDF-EOS，而且可以在好几种系统平台上进行运行，包括：SunSolarisworkstations，SGIIRIXworkstations，Linux和MicrosoftWindows。MRT可以通过命令行或在MRT

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

Class01234567891011IGBP(Type1)水体常绿针叶林t常绿阔叶林落叶针叶林落叶阔叶林混交林郁闭灌丛开放灌丛多树草原稀树草原草原永久湿地UMD(Type2)LAI/FPAR(Type3)NPP(Type4)水体常绿针叶林常绿阔叶林落叶针叶林落叶阔叶林混交林郁闭灌丛开放灌丛多树草原稀树草原草原水体草地/谷类作物灌木丛阔叶作物稀树草原阔叶林针叶林无植被城市水体常绿针叶林常绿阔叶林落叶针叶林落叶阔叶林一年生阔叶植被一年生草地植被无植被城市PFT(Type5)水体常绿针叶林常绿阔叶林落叶针叶林落叶阔叶林灌木丛草地谷类作物阔叶作物城市和建设用地冰和雪裸地或低植被覆盖地121314作物城市和建设用地作物/自然植被的镶嵌体作物城市和建设用地

15永久冰雪16裸地或低植被覆盖地裸地或低植被覆盖地254未分类未分类未分类未分类未分类

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

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

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

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

当然，假使要处理大量类似的MODISHDF-EOS数据，如本文需要处理海河流域2023年全年的每日地表温度数据（MOD11A1），16天合成的归一化植被指数产品（MOD13A2），16天合成的地表反照率产品（MOD43B3），地表覆盖产品（MOD12Q1），而且假使依旧选择运行这种MRTGUI方式进行数据操作的话，显然会让人感到无法忍受的。这时，可以选择通过命令行输入的方式来操作MRT对数据进行重采样和镶嵌处理。这就是常说的利用MRT对MODISHDF-EOS文件进行批处理（BatchProcessing）。命令行输入最大的好处就是可以大批量的处理相关的MODIS数据。本文以处理海河流域2023年全年的每日地表温度数据为例来简要介绍如何利用MRT命令行对MODISHDF-EOS数据进行重采样和镶嵌处理。2023年全年的海河流域每日地表温度数据共有1000多景数据影象。利用MRT对数据进行重采样和镶嵌处理的参数与上述LST.PRM一致。在数据处理开始之前，Windows用户需要首先安装Cygwin，并对LINUX语言有一些最基础的语法了解。关于Cygwin的状况及下载请参阅：

http://./。在确保MRT及Cygwin正确安装好以后，运行Cygwin，利用命令cdworkdirectory（如cdE:\\LST）把即将要进行处理的MODISHDF-EOS文件所在目录设为工作目录。另外，需要注意的是，请确保参数文件（如LST.PRM）事先利用MRTGUI已经得到并与将要进行大批量处理的HDF-EOS文件保存在同一个工作目录下，以免进行处理的时候无法找到正确的参数文件。在利用cd命令把工作目录设置好以后，可以手工输入或粘贴相应的批处理命令（假使您事先已准备好这些命令）。

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

（1）安装MRT的时候，假使用户想利用GUI，则可能需要安装JAVA（），否则只能利用命令行的方式进行操作。

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

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

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

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

MRT(MODISReprojectionTool)批处理MOD11A1的代码

(2023-10-1410:46:50)

标签：mrt批处理modis

以处理MOD11A1数据为例,处理其他MODIS产品类似:代码:

day=2023001

while[$day-le2023365]do

fortrlinMOD11A1.A$day.*.hdfdo

if[-f$trl]then

HDFFILES=$(lsMOD11A1.A$day.*.hdf)echo$HDFFILES>MOSAICINPUT.TXT

MRTMOSAIC-iMOSAICINPUT.TXT-s\RESAMPLE-pLST.PRM-iMOSAIC_TMP_$day.hdf-oMOD11A1_$day.tifrmMOSAIC_TMP_$day.hdflet\breakelse

let\breakfidonedone

LST.PRM参数文件:

#The\fieldwouldbecommented.Ifyouwanttoloadmultipleinputfilespleaseuncommentthe\fieldandcommentthe\field.

#Alsothe\\

#INPUT_FILENAMES=(E:\\LST\\MOD11A1.A2023131.h26v04.004.2023115130411.hdfE:\\LST\\MOD11A1.A2023131.h26v05.004.2023115134041.hdfE:\\LST\\MOD11A1.A2023131.h27v04.004.2023115150352.hdfE:\\LST\\MOD11A1.A2023131.h27v05.004.2023115125635.hdf)INPUT_FILENAME=E:\\LST\\TmpMosaic.hdf

SPECTRAL_SUBSET=(1)

#ORIG_SPECTRAL_SUBSET=(100000000000)SPATIAL_SUBSET_TYPE=INPUT_LAT_LONG

SPATIAL_SUBSET_UL_CORNER=(50.0124.457906149)SPATIAL_SUBSET_LR_CORNER=(30.0115.470053838)OUTPUT_FILENAME=E:\\LST\\LST.TIFRESAMPLING_TYPE=NEAREST_NEIGHBOROUTPUT_PROJECTION_TYPE=AEAOUTPUT_PROJECTION_PARAMETERS=(0.00.025.047.0105.00.00.00.00.00.00.00.00.00.00.0)DATUM=WGS84

OUTPUT_PIXEL_SIZE=1000

最近好多朋友反映我BLOG上的利用MRT处理MODIS的代码通不过的问题，问题大约就是如下图所显示的\\mrt\\Modis\\bin下的cygwin1.dll版本（100）与\\cygwin\\bin目录下的cygwin1.dll版本（100）不相容的错误：

解决方法：

各位可以手动把上面提到的两个目录下的cygwin1.dll版本调成一致就可以正常运行我BLOG上运行的代码了。（拷贝）

MRTSWATHTOOL批处理MODIS影象代码(以处理MOD021KM为例)

本代码为处理MOD021KM的程序，处理其他MODIS数据产品类似:代码:

day=2023121time=100

while[$day-le2023151]do

while[$time-le500]

do

fortrlin/cygdrive/I/tmp/mod02/MOD021KM.A$day.0$time.005.*.hdfdo

forgeoin/cygdrive/i/tmp/MOD03.A$day.0$time.005.*.hdfdo

if[-f$trl]then

swath2grid-if=$trl-gf=$geo-of=MOD021KM.A$day.0$time-off=GEOTIFF_FMT-sds=\0,0,0,0,0,0,1,1,1,1,1,1,1,1;EV_500_Aggr1km_RefSB,1,1,1,1,1\-opsz=1000.0

let\breakelse

let\breakfidonedonedone

let\let\done

本程序稍微涉及到LINUX语言,并且要安装CYGWIN,代码含义请参看MRTSwath_Users_Manual。

代码的含义：

-if输入文件名-of输出文件名-off输出格式

-sds输出哪个波段的数据，0表示不输出，1代表输出。在MODIS021