MODIS数据下载与处理_new
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
MODIS数据(1B、产品、常用处理软件)免费下载
1、MODIS 1B数据下载
L1B数据下载地址:/data/search.html
打开网页,选中需要的数据源类型
注意:其中,MOD03数据是用于对1KM,QKM,HKM数据进行几何纠正所用。
别忘记下载。
其中,日期类型为:月/日/年时:分:秒
其中,网页中显示的时间为UTC时间,换算为北京时间为:UTC时间=北京时间-8小时。
因此,要获得1月16日的数据则范围为:01/15/2003/ 16:00:00~01/16/2003/ 16:00:00
在‘spatial selection’选项中选择“latitude/longtitude”,按经纬度形式选择影像范围。
点击’search’查到需要的数据:
勾选需要的数据,点击‘order files now’,输入你接收信息的邮箱,点’order’开始订购该数据。
(如果要搜索多天数据,可以选‘add files to shopping cart’继续搜索其他日期的数据。
所订购数据的存放位置信息:
点击‘Data->Track Orders ’可以查看所有已订购的数据的状态。
如果’state’显示‘avalable’即可开始下载。
ftp:
username: anonymous
password: *****************(自己申请数据的邮箱)一般使用FTP下载工具下载。
首先新建站点:
点击‘连接’。
则在右边的框中会显示所有数据,找到自己数据所在的文件夹,并拖到左边的框中,开始下载数据。
OK!
2、MODIS产品下载地址
https:///api//u515127/u515368/u515406/u515421/u515529/u516066/u516091/u
516110/u516150/u516307/u517021/u517029(要装IE7支持)
或已知所下载文件名的modis09~17数据可以从FTP上下载:
ftp:/
tarra的数据放在MOLT里。
Aqua 数据放在MOLA。
上述的搜索网站(https://)(如果遇到能查到数据但是下载不了的情况,也可以在FTP中通过查找与该数据相同日期同轨道号来下载。
)在网站查到该数据的好处是能查看数据的覆
盖区域。
可以查看覆盖范围。
假如:MOD11A1.A2008288.h28v06.005.2008290030125.hdf 从:
ftp:///MOLT/MOD11A1.005/2008.10.14/查找到2008-10-14 轨道号为H26V05,就从里面的文件夹里选出来下载就可以了。
一般使用FTP下载工具下载。
注:数据的轨道号:
从左到右为H,从上到下为V,可根据具体省份范围自己来数。
例如江西省的为(H27V05、H27v06, H28V05、H28v06四块组成)
3、常用的MODIS软件:
常用的MODIS查看软件还有:
a、MODIS explorer(推荐使用)其下载地址:
/software/HDFExplorer/或:
/,该软件可以方便查看HDF格式的MODIS元数据或信息。
b、modis swath tool(推荐使用)
可用NASA网站提供的modis swath tool 对HDF格式的1B数据进行几何精纠正(下载地址:/pub/imswelcome/),该软件使用MOD03数据对影像进行纠正,处理速度快且使用简单方便。
比直接用ENVI的Georeference 更好用。
(好像也可进行影像拼接)
C、MRT(推荐使用)
可以用来对下载的hdf格式的MODIS产品进行重投影等,下载地址:/pub/imswelcome/)。
更多软件下载:/software.php#HDFExplorer
包括:
•Binary_Dumper -- dumps HDF-EOS objects in binary format
•CV -- reads the HDF and HDF-EOS files in native ISIS cube format
•EOSVIEW --displays and verifies HDF and HDF-EOS files
•HDFExplorer -- visualization program that reads HDF, HDF-EOS and HDF5 files
•HDFView -- a visual tool for browsing and editing NCSA HDF4 and HDF5 files
•HDF-EOS plug-in for HDFView -- extends HDFView functionality for browsing any HDFEOS file
•HE5View -- a file viewing tool for examining and verifying HDF-EOS 5.x files
•HEG -- converts HDF-EOS to GeoTIFF
•HEMU -- a tool to update metadata inside an HDF-EOS file
•HMR -- reads metadata from HDF-EOS file and writes to a text file
•HMU -- updates the metadata embedded in an HDF-EOS file
•JEB -- a Java EOS Browse tool
•MS2GT -- reads HDF-EOS files containing MODIS swath data and produces flat binary files containing gridded data
•MTD -- allows EOSDIS extended data providers to format their products in ECS standard formats
•MetCheck -- validates the inventory meta-data in a ".met" file
•ODL2XML -- converts any ECS ODL file to an ECS XML file
•PGS_Toolkit -- a set of tools to manage the metadata that are generated with each EOS product
•PHDIS -- reads any HDF-EOS file containing data gridded in the Lambert Azimuthal Equal Area projection
•TRMM_Data_Reader -- reads in a TRMM HDF-EOS data file, and writes out SDS and Vdata as binary files
•UFM -- accepts an ODL file as input, and generates a simple HTML representation of the ODL
•hdf2bin -- converts HDF or HDFEOS file into plain binary file
•hdfeos-netcdf -- converts HDF-EOS4/5 files to netCDF
•hdfeos52netcdf -- converts HDF-EOS5 files to netCDF
•he2he5_lib -- a library of HDF-EOS 4/5 Compatibility Layer
•he5cc -- a shell script that does for HDF-EOS5 what the h5cc script does for HDF5
•he5ex -- extracts HDF-EOS5 data to binary or ASCII data format in HTML or XML index
•he5mu -- updates metadata inside an HDF-EOS5 file
•he5v -- a Profile Validator for HDF-EOS 5
•heconvert -- converts HDF-EOS Version 2.x files to HDF-EOS Version 5.x
•heex -- extracts HDF-EOS data to binary or ASCII data format in HTML or XML index •view_hdf -- a visualizationtool for accessing data stored in HDF and HDF-EOS files •xml2he -- takes an XML document and creates an HDF-EOS4/5 file matching the document
主要网站:浏览,格式转换,地图投影等工具。
美国:
/hdfeos/softwarelist.cfm
/hdfeos/Info/index.cfm (HDF/ HDF-EOS 资料大全:从入门到专业级别都有)。
/tools.html
/hdfeos/forum/categories.cfm?catid=15 ( 各种HDF 工具使用的问题和讨论)
/hdfeoss.html
/pub/imswelcome/ (MODIS Swath Reprojection Tool (MRT Swath))
/pub/imswelcome/ (工具清单。
如:HEG, )
俄罗斯:
The International MODIS/AIRS Processing Package (IMAPP)
http://www.scanex.ru/eosdb/eosdb.htm(俄罗斯:ENVI modistools_envi.zip 的来源)
http://www.scanex.ru/eosdb/imapp.htm (俄罗斯:把IMAPP 转化为标准的HDF-EOS:
imapp2daac.zip)
http://www.scanex.ru/downloads/imapp/imappw_um.htm#faq (有关IMAPP 的问答)Format Converters
HDF-EOS Binary Dumper
bindmp creates an output file which contains binary data for a specified
HDF-EOS object.
Platform(s): DEC, SGI
Data Application(s): All Data Applications, gridded data, point data, swath
data
License: SESDA
10/27/2001
This tool will convert many EOS products written in HDF-EOS to GeoTIFF, native
binary or HDF-EOS Grid. It also has re-projection, resampling, subsetting,
stitching (mosaicing), and metadata creation capabilities. New features
implemented since v0.6
Platform(s): Linux, SGI, Sun, Windows
Data Application(s): ASTER, MISR, MODIS
License: ECS
02/12/2003
hdf2bin
converts HDF or HDFEOS file into binary format. Information about the
extracted science data sets is printed to standard output (screen) or can be
saved to an ASCII file.
Platform(s): SGI, Sun, Win98, Win2000, WinNT
Data Application(s): All Data Applications
License: DAAC
05/03/2002
hdfeos2bin
This tool extracts data arrays from an HDF-EOS file and writes the arrays to
separate flat binary data files.
Platform(s): Win98, Win2000, WinNT, Unix
Data Application(s): All Data Applications
License: DAAC
5/6/2003
Not rated
HDF Explorer v1.3
HDF Explorer is a data visualization program that reads Hierarchical Data Format files (HDF, HDF-EOS and HDF5) and netCDF files.
Platform(s): Win95, Win98, Win2000, WinNT
Data Application(s): All Data Applications
License: Commercial
09/21/2004
Not rated
MODIS Swath-to-Grid Toolbox (MS2GT) v0.5
The MODIS Swath-to-Grid Toolbox (MS2GT) is a set of software tools that can be used to read HDF-EOS files containing MODIS swath data and produce flat
binary files containing gridded data in a variety of map projections.
Platform(s): Linux, Unix
Data Application(s): MODIS, swath data
License: DAAC
4、MODIS 1B数据需处理的过程(几何纠正,DN转反射率,拼接)
1、几何纠正
1、安装从网站上下载MODIS Swath tool 软件(用于1B数据的几何纠正,下载地址:
/pub/imswelcome/需要注册(免费)才能下载) 以及java 软件。
并安装,安装方法与MODIS tools相同。
2、然后从安装路径中ModisSwathTool.bat,打开MODIS swath tool,’Open input file …’,打
开下载的1B数据。
3、从菜单‘file’->’open geolocation file’打开与该数据想配套的03级数据。
4、设置投影参数:
点击‘edit Projection Parameters …’修改参数:然后点击”run”,几何纠正开始。
二、DN值转反射率(scale_offset改正)
1、几何纠正完成后即可用ENVI打开:
2、从basic tools->preprocessing->Data-specific Utilities->view HDF attributes,从未经几何纠正前的数据查看该数据scale_offset值。
3、纠正gain_offset值: 从basic tools->preprocessing->General Purpose Utilities->Apply Gain and Offset.
输入相应的scale 和offset值,并保存图像。
5、拼接
在弹出来的对话框中,import->”import files and proporties…”选择要进行拼接的影像,
在“data value to ignore”中输入要忽略的背景的值。
这样做能使两影像之间的拼接时影像外的区域为透明显示。
显示结果,如下,如果选择了import->”import files”可以通过右击影像列表中的影像,选”edit entry…”来修改要忽略的背景色的值,效果是一样的。
然后选file->apply输出拼接的结果。
5、MODIS产品处理之MRT重投影参数设置
如果下载的是多幅图像该软件可以进行自动拼接:方法是,输入多幅影像如下图:
6、MRT(Modis Reprojection Tool)批处理MODIS影象(以海河流域为例)
摘自:/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存储格式。
免费用户注册及下载地址:/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种分类方法
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的情况及下载请参阅:
/。
在确保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(),否则只能利用命令行的方式进行操作。
(2)利用MRT进行数据处理的时候,请确保文件不要保存在太多的子目录下和尽量不要使用中文目录,以免在MRT导入处理数据的时候可能会出现找不到文件的错误。
(3)在安装Cygwin的时候,用户可以有选择性的进行安装,只要满足可以利用命令行进行操作的目的就可以了。
(4)请确保参数文件事先已经创建并与要进行批处理的MODIS HDF-EOS文件保存在同一工作目录下,只有这样才有可能数据被正确处理。
(5)由于LINUX语言对大小写非常敏感,在输入批处理命令的时候,请确保文件名正确输入,包括文件名后缀。
MRT(MODIS Reprojection Tool)批处理MOD11A1的代码(2007-10-14 10:46:50)
标签:mrt批处理modis
以处理MOD11A1数据为例,处理其他MODIS产品类似:
代码:
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
let "day=day+1"
break
else
let "day=day+1"
break
fi
done
done
LST.PRM参数文件:
#The "INPUT_FILENAMES" field would be commented. If you want to load multiple input files please uncomment the "INPUT_FILENAMES" field and comment the"INPUT_FILENAME" field.
#Also the "ORIG_SPECTRAL_SUBSET" field needs to be uncommented and changed to "SPECTRAL_SUBSET". The initial "SPECTRAL_SUBSET" field should be deleted.
#INPUT_FILENAMES = ( E:\LST\MOD11A1.A2001131.h26v04.004.2003115130411.hdf
E:\LST\MOD11A1.A2001131.h26v05.004.2003115134041.hdf
E:\LST\MOD11A1.A2001131.h27v04.004.2003115150352.hdf
E:\LST\MOD11A1.A2001131.h27v05.004.2003115125635.hdf )
INPUT_FILENAME = E:\LST\TmpMosaic.hdf
SPECTRAL_SUBSET = ( 1 )
#ORIG_SPECTRAL_SUBSET = ( 1 0 0 0 0 0 0 0 0 0 0 0 )
SPATIAL_SUBSET_TYPE = INPUT_LAT_LONG
SPATIAL_SUBSET_UL_CORNER = ( 50.0 124.457906149 )
SPATIAL_SUBSET_LR_CORNER = ( 30.0 115.470053838 )
OUTPUT_FILENAME = E:\LST\LST.TIF
RESAMPLING_TYPE = NEAREST_NEIGHBOR
OUTPUT_PROJECTION_TYPE = AEA
OUTPUT_PROJECTION_PARAMETERS = (
0.0 0.0 25.0
47.0 105.0 0.0
0.0 0.0 0.0
0.0 0.0 0.0
0.0 0.0 0.0 )
DATUM = WGS84
OUTPUT_PIXEL_SIZE = 1000
最近很多朋友反映我BLOG上的利用MRT处理MODIS的代码通不过的问题,问题大概就是如下图所显示的\mrt\Modis\bin下的cygwin1.dll版本(1005.24.0.0)与\cygwin\bin目录下的cygwin1.dll版本(1005.25.0.0)不相容的错误:
解决方法:
各位可以手动把上面提到的两个目录下的cygwin1.dll版本调成一致就可以正常运行我BLOG上运行的代码了。
(拷贝)
MRT SWATH TOOL批处理MODIS影象代码(以处理MOD021KM为例)
本代码为处理MOD021KM的程序,处理其他MODIS数据产品类似:
代码:
day=2003121
time=100
while [ $day -le 2003151 ]
do
while [ $time -le 500 ]
for trl in /cygdrive/I/tmp/mod02/MOD021KM.A$day.0$time.005.*.hdf
do
for geo in /cygdrive/i/tmp/MOD03.A$day.0$time.005.*.hdf
do
if [ -f $trl ]
then
swath2grid -if=$trl -gf=$geo -of=MOD021KM.A$day.0$time -off=GEOTIFF_FMT
-sds="EV_1KM_Emissive, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1; EV_500_Aggr1km_RefSB, 1, 1, 1, 1, 1" -kk=NN -oproj=ALBERS -oprm=0.0,0.0,25.0,47.0,105.0 -osp=8
-opsz=1000.0
let "time=time+1"
break
else
let "time=time+1"
break
fi
done
done
done
let "day=day+1"
let "time=100"
done
本程序稍微涉及到LINUX语言,并且要安装CYGWIN,代码含义请参看
MRTSwath_Users_Manual。
代码的含义:
-if 输入文件名
-of 输出文件名
-off 输出格式
-sds 输出哪个波段的数据,0表示不输出,1代表输出。
在MODIS021KM中
EV_1KM_Emissive有16个波段的数据。
如果不需要输出EV_1KM_Emissive,可以不写,也可以全部写0
-kk 重采样方法 NN为最近邻采样
-oproj 投影参数
-osp=8 WGS1984
具体参数可以参考mrtswath使用手册
在安装完cygwin,MRTswath以后,运行cygwin,
在$下设置数据所在的目录如:cd I:/tmp
在打开的TXT文本中复制上述代码,然后粘贴:
就开始运行了
7、 cygwin安装下载setup.exe
安装Cygwin,您需要下载setup.exe,此程序根据您选择的软件包,再从互联网下载所有组件。
点击此处下载setup.exe
安装
1.双击setup.exe
2.按照安装提示“下一步”
3.选择从“互联网”安装:
4.选择安装的路径,一般是C:\Cygwin,选择本地临时目录,用于存放下载的软件;选择您连接互联网的方式,一般选“直接连接”就可以了。
5.选择下载源,您可以在下载列表里选择:或者,直接在URL里输入/pub/
6.根据安装提示,选择需要安装的软件包,完成安装。
如果您在国内,是您所能使用的最快的镜像。
使用
•双击桌面的Cygwin图标,即可启动Cygwin里的Bash。
•或者,你可将C:/Cygwin/bin加到%PATH%当中,您就可以直接在cmd.exe里面使用Linux命令了,比如less, cat, wc , wget
•如果您有一个Bash脚本,您可以直接用C:/Cygwin/bin/bash.exe
X:/scripts/Monitor.sh的形式调用。