MRT批处理MODIS数据(详细教程)
MRT 批处理命令详解(附加代码)
利用MRT进行MODIS产品的批处理步骤：1.首先利用MRT生成参数文件.prm（parameters file）生成参数文件主要分为三个部分，一是选择要合成的文件，选择需要的波段，二是确定空间范围；三是确定投影信息，空间分辨。
确定好这三个部分，把参数文件保存起来。
2.MRT批处理命令第一条命令：dir /b/s D:\modistools\bin\NDVI\*A2013001*.hdf >> D:\modistools\bin\NDVI\list.txt该命令的意思是找到文件存储路径下所有的*A2013001*.hdf文件，生成这些文件的list.txt。
比如以东北地区，MOD123 产品为例，2013年001天的文件共有12个。
第二条命令：D:MRTMOSAIC -i D:\modistools\bin\NDVI\list.txt -o D:\NDVI\MOD123.hdf -s "1 0 0 0 0 0 0 0 0 0 0 "该命令首先是利用list.txt 文件了的文件进行合成，-o 是输出的意思，-o D:\NDVI\MOD123.hdf 是输出的临时文件，名字随便起，这个临时文件里存储的是想要的数据集。
-s "1 0 0 0 0 0 0 0 0 0 0 " 代表输出的是第一数据集，用1和0 进行选择，需要的数据集设为1，不需要的设为0，以MOD123为例，一个MOD123.hdf 里共有11个数据集（data fields），如果只需要NDVI，则是"1 0 0 0 0 0 0 0 0 0 0 "根据自己选取的数据集，MOD123.hdf 临时文件就存储成如下结构：第三条命令：RESAMPLE -p D:\modistools\bin\NDVI.prm -i D:\NDVI\MOD123.hdf-o D:\modistools\bin\NDVI\NDVI.tif该命令是重采样，根据第一步保存的参数文件和临时文件，进行重采样，然后输出，需要注意的是，输出的文件名是又两部分组成自己起的名字+临时文件数据集的名字，如在这里自己起的名字是NDVI.tif，临时文件数据集的名字是 1 km monthly NDVI, 那么存储到路径下的文件的名字是NDVI. 1_km monthly_NDVI.tif，临时文件数据集名字里的空格在输出时保存为下划线。
MODIS数据下载与处理_new
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: anonymouspassword: *****************（自己申请数据的邮箱）一般使用FTP下载工具下载。
首先新建站点：点击‘连接’。
则在右边的框中会显示所有数据，找到自己数据所在的文件夹，并拖到左边的框中，开始下载数据。
OK！2、MODIS产品下载地址https:///api//u515127/u515368/u515406/u515421/u515529/u516066/u516091/u516110/u516150/u516307/u517021/u517029(要装IE7支持)或已知所下载文件名的modis09~17数据可以从FTP上下载：ftp:/tarra的数据放在MOLT里。
modis数据处理流程
modis数据处理流程下载温馨提示:该文档是我店铺精心编制而成，希望大家下载以后，能够帮助大家解决实际的问题。
文档下载后可定制随意修改，请根据实际需要进行相应的调整和使用，谢谢!并且，本店铺为大家提供各种各样类型的实用资料，如教育随笔、日记赏析、句子摘抄、古诗大全、经典美文、话题作文、工作总结、词语解析、文案摘录、其他资料等等，如想了解不同资料格式和写法，敬请关注!Download tips: This document is carefully compiled by theeditor. I hope that after you download them,they can help yousolve practical problems. The document can be customized andmodified after downloading,please adjust and use it according toactual needs, thank you!In addition, our shop provides you with various types ofpractical materials,such as educational essays, diaryappreciation,sentence excerpts,ancient poems,classic articles,topic composition,work summary,word parsing,copy excerpts,other materials and so on,want to know different data formats andwriting methods,please pay attention!MODIS 数据处理流程一般包括以下几个步骤：1. 数据下载：确定需要处理的 MODIS 数据产品和时间范围。
MODIS数据处理过程
1.MCD12Q1数据下载下载地址：/data/search.html打开网页，出现如下图：1.1选中需要的数据源类型：1.2 选择所要下载数据的时间，通道，X围等：1.3点击search查到需要的数据：1.4出现如下页面：1.5勾选需要的数据，点击‘order files now’，输入你接收信息的，点’order’开始订购该数据。
(如果要搜索多天数据，可以选‘add files to shopping cart’继续搜索其他日期的数据。
出现下面的页面，填写你的订购，qq不能用，其他的都可以，最好是你能够进去的，点击下面的Order。
1.6然后出现如下页面：点击红色标记部分能看到如下你订购的信息：（注意：红色部分是你要下载数据的文件名字）1.7 一般使用FTP下载工具下载。
（FTP在网上可以下载得到）首先新建站点：ftp:username: anonymouspassword: .（自己申请数据的）1.7点击‘连接’。
在右边的框中会显示所有数据，找到自己数据所在的文件夹（我这里的文件夹名字跟上面的不是同一个，这里仅供参考，以上面的那个文件夹为准），并拖到左边的框中（选择你要保存数据的目录），开始下载数据。
OK！！2.处理数据第一步：数据读取在ENVI中file中第一个选项打开hdf格式的文件：第二步：数据地理定位MODIS三级数据产品土地覆盖数据是已经经过投影了的，其投影类型为正弦曲线投影，但是很多遥感软件在打开HDF数据的时候，会有投影丢失现象，ENVI 也是，可以通过hdf属性查看工具，查看ENVI打开hdf时未识别的投影信息，以及投影参数，利用该参数在ENVI中自定义正弦曲线投影，为MCD12Q1数据添加投影。
（1）用hdf属性查看工具查看该hdf数据的地理定位信息，选择你刚才打开的那副影像，可得图像左上角、右下角坐标值，投影类型，投影参数，地球半径，起始位置。
2）在波段列表中，单击右键打开edit Header，打开Edit Attributes->Map Info。
MODIS影像处理流程——ok
1、同时期影像拼接、重投影、转格式（MRT）黄土高原地区Projected coordinate system name: Beijing_1954_3_Degree_GK_CM_108EGeographic coordinate system name: GCS_Beijing_1954Map Projection Name: Transverse MercatorScale Factor at Central Meridian: 1.000000Longitude of Central Meridian: 108.000000Latitude of Projection Origin: 0.000000False Easting: 500000.000000False Northing: 0.000000北京54参考的椭圆：椭球坐标参数：长半轴a=6378245m；短半轴=6356863.0188m2、在Erdas中裁剪，得到研究区的原值NDVI -2000~9663方法：aoi格式文件裁剪：在ERDAS图标面板工具条中单击Data Prep图标，Subset，打开Subset对话框。
在Subset对话框中需要设置下列参数：①输入文件名（input file）②输出文件名（output file）③单击aoi选中裁剪文件名：在choose aoi选项中选中文件名④其余的我选择默认，然后ok等结果3、得到范围在-1~1的NDVI方法：ENVI — Basic Tools — Band Math（波段运算，把负值去掉）在弹出的对话框Band Math中，Enter an expression： (b1 lt 0)*0+(b1 ge 0)*(b1*0.0001)。
这个公式意思就是：要是值小于0 就乘以0，使其变为0；同时，值大于等于0的话就乘以0.0001这个系数。
这样就OK了。
关于B1 代表所选的那个影像。
MODIS数据处理方法(整理)
一般处理步骤：
1几何校正
Map----georeference MODIS--------选择波段--------选择投影方式（经纬度）------OK(选择辐射率数据，波段1-7)
2.拼接
3.大气校正(flaash)
用envi自带的flaash模块做大气校正，注意，输入文件的要求：1）数据是经过定标后的辐射率数据，2）数据要有中心波长值，可通过edit header输入，3）支持四种数据类型：浮点型等等，4）数据储存类型：标准栅格文件，且必须是BIP或BIL，5)适用的数据光谱范围是0.4-2500μm。
1.1打开经过几何校正的数据，输入中心波段，直接导入输好的txt文件；
1.2修改数据类型，点开basic tools---convert data把数据类型改成BIL，选择
in place直接覆盖原文件；
1.3点开Spectral-----Flaash做大气校正，输入辐射率文件，
2创建ROI区域
3导出所需数据。
Envi调用MODIS Reprojection Tool(MRT)对MODIS产品进行批处理拼接
Envi调用MODISReprojectionTool（MRT）对MODIS产品进行批处理拼接、重投影、裁切1熟悉MRTMODIS产品的类型不同，一景HDF格式的影像包含的波段也各不相同。
MRT处理时需要选择处理波段，0表示不作处理，1表示处理，首先要确定影像的波段数。
1、拼接Mrtmosaic.exe程序用来拼接影像。
调用方式为：mrtmosaic-i"g:\n%1.txt"-s"0000000000000000000000000001110000"-o"g:\te mp.hdf"命令行帮助如下Usage:mrtmosaic-iinput_filenames_file-t-h-ooutput_filename-sspectral_subset"b1b2...bN"-gfilenameforthelogfile需要指定三个参数：-I指定输入文件，可以采取两种方式1、1.hdf2.hdf3.hdf2、把影像的完整路径保存到txt文件中，作为input参数。
如-IMOD092008001.txt。
-s指定需要处理的波段，同样可以采用两种方式1、直接给出，如–s“000000000000”，注意英文引号；2、指定一个txt路径，让程序读取；-o指定输出路径，一般直接给出1、直接给出，如–og:\tmp.hdf。
注意直接存为HDF格式，便于后续处理。
2、重投影、裁切Resample.exe用来重投影、裁切是MRT程序的核心。
调用方式为：命令行帮助如下：Usage:resample-pparameter_file[options] Optionsthatoverrideparameterfilespecifications:-iinput_file_name-ooutput_file_name-rresampling_type[NNBICCNONE]-tprojection_type[AEAERGEOHAMIGHISINLALCCMERCATMOLPSSINTMUTM] -jprojection_parameter_list"p1p2 (15)-sspectral_subset"b1b2...bN"Ifusingthe-sswitch,theSDSsshouldberepresentedasanarrayof0sand1s.A'1'specifiestoprocessthatSDS;'0'specifiestoskipthatSDS.UnspecifiedSDSswillnotbeprocessed.Ifthe-sswitchisnotspecified,thenallSDSswillbeprocessed.-aspatial_subset_type[INPUT_LAT_LONGINPUT_LINE_SAMPLEOUTPUT_PROJ_C OORDS]-lspatial_subset"ULlatULlongLRlatLRlong"-or-"ULlineULsampleLRlineLRsample(0-based)"-or-"ULprojxULprojyLRprojxLRprojy"NOTE:line/samplemustbespecifiedforthehighestresolution ofallSDSsspecifiedtobeprocessedintheproduct.-uUTM_zone-xpixel_size-gfilenameforthelogfile可以只指定1个参数：用–p读入prm参数文件，进行处理。
批量给MODIS拼接和重投影
资料范本本资料为word版本，可以直接编辑和打印，感谢您的下载批量给MODIS拼接和重投影地点：__________________时间：__________________说明：本资料适用于约定双方经过谈判，协商而共同承认，共同遵守的责任与义务，仅供参考，文档可直接下载或修改，不需要的部分可直接删除，使用时请详细阅读内容先进行拼接处理先新建一个DATA文件夹，并将HDF数据拷贝到该文件夹下在DATA文件夹下再新建一个Result文件夹（用于存放拼接后的数据）将以下代码内容粘贴到一个txt文件中，更改txt文件扩展名为bat即可，假设你命名该bat文件名为MODISmosaic.bat；程序导读：rem 表示开始的为注释；MOSAICINPUT.TXT这是程序运行时自动生成的，不用管；MRTDATADIR为你的MRT安装文件中data的路径，“c:/MRT/bin/mrtmosaic.exe”改成你的mrtmosaic.exe的安装路径。
"set /a DAY= %DAY% + 16 "则是因为输入数据是16天间隔的，根据你的数据进行相应修改。
rem Set the MRTDATADIR environmental var to the MRT data directory.set MRTDATADIR=C:\MRT\dataset /a DAY=2011097 rem **batch data start time**set /a DEADLINE=2011273 rem **batch data end time**:startif %DAY% leq %DEADLINE% (goto ORDER) else exit:ORDERrem **save the file name into a notepad**dir *%DAY%.*.hdf/a/b/s > MOSAICINPUT.TXTrem **execute mosaic **rem Set the mrtmosaic.exe directory.c:/MRT/bin/mrtmosaic.exe -i MOSAICINPUT.TXT -s "1 0 0 0 0 0 0 0 0 0 0 0" -o MOSAIC_TMP_%DAY%.hdfrem **copy the result to a file and delete the input data**copy MOSAIC_TMP_%DAY%.hdf Result & del MOSAIC_TMP_%DAY%.hdfdel *%DAY%.*.hdfset /a DAY= %DAY% + 16goto start点击MODISmosiac.bat，即可进行拼接操作注意：该操作数据及BAT文件须放在一个文件夹（文件夹起名请用英文，MRT对中文路径数据会打不开）下。