没有合适的资源?快使用搜索试试~ 我知道了~
MRT批处理MODIS数据步骤
4星 · 超过85%的资源 需积分: 47 161 下载量 187 浏览量
2013-05-19
15:16:42
上传
评论 8
收藏 4KB TXT 举报
温馨提示
试读
2页
自己编写的利用MRT批处理MODIS数据的详细步骤
资源推荐
资源详情
资源评论
1.list的生成。
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
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
nwsxer
- 粉丝: 1
- 资源: 2
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
最新资源
- Screenshot_20240504_153703.jpg
- video_20240505_113243.mp4
- 2_GMT20240401-113240_Recording.m4a
- 基于QT+C++开发的实现了网格显示多张带文字的图片在界面中自适应布局+源码
- CS224N_PyTorch_Tutorial.ipynb
- Redis要点和难点实际应用方向和具体案例代码&代码解析
- 566827122690661伏羲X多开安装包.apk
- matlab 图像质量评价经典算法代码
- 基于QT+C++开发的带有标签、下划线的Material Design风格输入框QLineEdit+源码
- PS3111固件,,亲测能搞
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功
- 1
- 2
前往页