使用MRT批量处理MODIS数据教程
"MRT批处理MODIS数据步骤" 这篇内容是关于如何使用MRT(MODIS Reprojection Tool)进行批处理MODIS数据的详细教程。MRT是一款由NASA开发的工具,用于重投影、裁剪、组合以及转换MODIS(Moderate Resolution Imaging Spectroradiometer)卫星数据。MODIS数据广泛应用于气候变化、环境监测等领域。 1. **生成文件列表**: 首先,你需要在指定目录下找到特定日期(如2000年的第49天)的MODIS数据文件,并创建一个文本文件(如lt2000049.txt),列出所有相关的HDF文件。通过命令行工具,可以使用`dir /bH:\MODISData\nasa\2000\*A2000049*.hdf>>H:\MODISData\nasa\2000\lt2000049.txt`来完成这个操作。 2. **批处理设置**: 创建一个批处理脚本(.bat文件),该脚本将包含MRT的操作指令。在这个过程中,你需要指定输入文件名,这可以通过取消注释“INPUT_FILENAMES”字段并删除或更改“INPUT_FILENAME”字段来实现。同时,你需要调整“ORIG_SPECTRAL_SUBSET”字段,改为“SPECTRAL_SUBSET”,初始的“SPECTRAL_SUBSET”字段应被删除。 3. **参数设置**: 在MRT中,你需要定义输出文件的参数。这包括选择要使用的可用波段,以及选择特定的子集,如果需要。此外,还需要指定重采样类型(例如,最近邻插值),输出投影类型(如Albers Equal Area),标准纬度(如STDPR125和STDPR247),中央经度(如105度),地理坐标系(WGS84),输出像素大小(默认为250*250米,可能需要调整为其他尺寸,如500*500米),以及分辨率(可能需要保持与原始MODIS数据相同,即231*231,或者自定义)。 4. **保存参数文件**: 完成参数设置后,保存为 prm 文件,例如 MOD16DAY.prm。这个文件将记录所有的处理设置,以便于后续的批处理任务。 5. **执行批处理**: 使用保存的 prm 文件运行MRT,它会根据文件列表中的每个文件进行处理,将MODIS数据转换为你所设定的新格式和投影。输出文件通常为.tif格式,可以直接用于GIS分析或其他应用。 批处理MODIS数据是一个复杂的过程,涉及到文件管理、参数配置和地理空间数据处理等多个环节。正确地执行这些步骤可以有效地处理大量MODIS数据,提高工作效率,同时也确保了数据的一致性和准确性。对于进行气候建模、遥感分析等工作的科研人员来说,掌握这种方法至关重要。
dir /b H:\MODISData\nasa\2000\*A2000049*.hdf >> H:\MODISData\nasa\2000\lt2000049.txt
本语句用于把你的所有的.HDF文件生成.txt的list。
*A2000049*.hdf 意思是含有“A2000049”字段的所有文件,用于合成。这里是200年第49天的数据为例。
其他隔天按照此格式重复。换掉“A2000049”就行。如第65天用“*A2000065*.hdf”
>>后面部分是存贮路径和名称,为了便于理解和记忆,建议按照天数存贮。
写完之后,保存为.bat文件。
以上书写在txt中进行。
2..prm的生成
点开MRT,opern input file,打开几景要合并的影像;available band,选择你要输出的波段,右边的是要输出的波段,不需要的送到左边;下面的坐标设置,打开文件后显示的是文件文件的坐标,你可以不修改,就是后面的输出合成的是所有影响数据;你也可以修改,例如你的研究区域比较小,远小于几张影像的合成区,你可修改四角坐标,这样数据处理里速度快一些,保存的的数据库量也小一些;specify output file设置输出路径和名称,记住一定要在名字后面加.tif什么的,这个表示输出格式;RESAMPL Type选择一个重采样的方法,一般选择nearest neighber;output projection type,选择一种投影方式,如果选择Albers Equal Area,那么在Edit projection parameters时STDPR1填25,STDPR2填47,CentMer填105,这个全中国一样;下面的Datum选择一个坐标系,如WGS84;output pixel size 输出栅格大小,可以默认,也可以修改,如MODIS16day数据网站标明250*250,实际是231*231,如果默认就是231*231,你也可以修改为500*500,直接填入500就可以了。然后run.之后save parameters file,这一步最重要了,哥们,做这么多,就是为了这一步。保存为一个你习惯的名字。当然你也可以自己写。但是对于变成能力差和对于投影变换头痛的我们大多数人而言,这种方法获取.prm文件显然是最容易的。
生成之后的文件可以直接用。
下面是我生成的,prm文件(我保存的文件名字为MOD16DAY.prm,后面MRT命令中可以看到)
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 = ( H:\MODISData\nasa\2000\MOD13Q1.A2000049.h26v05.005.2006269215315.hdf H:\MODISData\nasa\2000\MOD13Q1.A2000049.h27v05.005.2006269214353.hdf )
INPUT_FILENAME = H:\MODISData\nasa\2000\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 = ( 39.6 107.1 )
SPATIAL_SUBSET_LR_CORNER = ( 34.8 111.3 )
OUTPUT_FILENAME = H:\MODISData\nasa\2000\MOD13Q1.A2000049.h26v05.005.2006269215315.tif
下载后可阅读完整内容,剩余1页未读,立即下载
- 粉丝: 1
- 资源: 2
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
最新资源
- WebLogic集群配置与管理实战指南
- AIX5.3上安装Weblogic 9.2详细步骤
- 面向对象编程模拟试题详解与解析
- Flex+FMS2.0中文教程:开发流媒体应用的实践指南
- PID调节深入解析:从入门到精通
- 数字水印技术:保护版权的新防线
- 8位数码管显示24小时制数字电子钟程序设计
- Mhdd免费版详细使用教程:硬盘检测与坏道屏蔽
- 操作系统期末复习指南:进程、线程与系统调用详解
- Cognos8性能优化指南:软件参数与报表设计调优
- Cognos8开发入门:从Transformer到ReportStudio
- Cisco 6509交换机配置全面指南
- C#入门:XML基础教程与实例解析
- Matlab振动分析详解:从单自由度到6自由度模型
- Eclipse JDT中的ASTParser详解与核心类介绍
- Java程序员必备资源网站大全