logo资料库

基于MRT工具批量给MODIS拼接和重投影.docx

第1页 / 共4页
第2页 / 共4页
第3页 / 共4页
第4页 / 共4页
资料共4页,全文预览结束
A. 先进行拼接处理 1. 先新建一个 DATA 文件夹,并将 HDF 数据拷贝到该文件夹下 2. 在 DATA 文件夹下再新建一个 Result 文件夹(用于存放拼接后的数 据) 3. 将以下代码内容粘贴到一个 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\data set /a DAY=2011097 rem **batch data start time** set /a DEADLINE=2011273 rem **batch data end time** :start if %DAY% leq %DEADLINE% (goto ORDER) else exit :ORDER rem **save the file name into a notepad** dir *%DAY%.*.hdf/a/b/s > MOSAICINPUT.TXT rem **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%.hdf rem **copy the result to a file and delete the input data** copy MOSAIC_TMP_%DAY%.hdf Result & del MOSAIC_TMP_%DAY%.hdf del *%DAY%.*.hdf set /a DAY= %DAY% + 16 goto start 4. 点击 MODISmosiac.bat,即可进行拼接操作 注意:该操作数据及 BAT 文件须放在一个文件夹(文件夹起名请
用英文,MRT 对中文路径数据会打不开)下。代码未考虑中间天 数间断的情况,比如对 MODIS 时间分辨率为 16 天的数据,在一个 应连续的等差数列中间有数据缺失,存在相差 32 天的情况,可能 会报错,请注意。同时,需要备份好你的输入数据,因为程序运 行中会自动删除拼接好的数据。 进行到此处则得到了初步的拼接数据,下面进行第 二步 B. 再进行重投影,改格式,缩小范围 1. 新建记事本,在记事本中输入下面代码,改 txt 后缀为 bat,定义 modis.bat 文件。 rem Set the MRTDATADIR environmental var to the MRT data directory. set MRTDATADIR=C:\MRT\data for %%i in (*.hdf) do C:\MRT\bin\resample -p my.prm -i %%i -o %%iout.tif 其中,MRTDATADIR 为你的 MRT 安装文件中 data 的路径 2. 利用 MRT 的图形界面工具定义投影文件 选择一期刚刚拼接好的影像,在 UL Corner 和 LR Corner 可以输入 经纬度来缩小你需要的区域,设置你的保存路径以及命名(最好 以日期命名,到时候程序会自动按照你的规则全部按日期命名), 选择 tif 格式,选择你要的投影类型,点击 Save Parameter File 将
投影文件保存。记住:这两个文件必须和刚刚拼接后的数据放在 一个文件夹中 3. 双击 modis.bat 执行批处理即可! 操作如图:
总的来说,就是先进行拼接得到一个 hdf 的临时文 件,然后在进行重投影,转格式,缩范围的操作。整个过 程你只需要安装 MRT 即可进行,不需要再安装 Cygwin 模 拟 lunix 环境,安装 Cygwin 很麻烦。 另外网上有小熊工具箱挺好,但是他编写的只能对 同一个 tile 进行长时序列的批量处理,所以遇到需要多景 影像的时候可以采用此法。 PS:番茄甜土豆就是我,我为他俩代言!
分享到:
收藏