《批量给MODIS拼接和重投影.pdf》由会员分享，可在线阅读，更多相关《批量给MODIS拼接和重投影.pdf（5页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、 批量给 MODIS 拼接和重投影 The Standardization Office was revised on the afternoon of December 13,2020 A.先进行拼接处理 1.先新建一个 DATA 文件夹，并将 HDF 数据拷贝到该文件夹下 2.在 DATA 文件夹下再新建一个 Result 文件夹（用于存放拼接后的数据）3.将以下代码内容粘贴到一个 txt 文件中，更改 txt 文件扩展名为bat 即可，假设你命名该 bat 文件名为；程序导读：rem 表示开始的 为注释 ；这是 程序运行 时自动 生成的，不用管；MRTDATADIR为 你 的MRT安 装
2、 文 件 中data的 路 径，“c:/MRT/bin/”改成你的的安装路径。set/a DAY=%DAY%+16 则是因为输入数据是 16 天间隔的，根据你的数据进行相应修改。rem Set the MRTDATADIR environmental var to the MRT data directory.set MRTDATADIR=C:MRTdata set/a DAY=2011097 rem*batch data start time*set/a DEADLINE=2011273 rem*batch data end time*:start if%DAY%leq%DEADLINE%(
3、goto ORDER)else exit:ORDER rem*save the file name into a notepad*dir*%DAY%.*.hdf/a/b/s rem*execute mosaic*rem Set the directory.c:/MRT/bin/-i -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
4、%.hdf del*%DAY%.*.hdf set/a DAY=%DAY%+16 goto start 4.点击，即可进行拼接操作 注意：该操作数据及 BAT 文件须放在一个文件夹（文件夹起名请用英文，MRT 对中文路径数据会打不开）下。代码未考虑中间天数间断的情况，比如对 MODIS 时间分辨率为 16 天的数据，在一个应连续的等差数列中间有数据缺失，存在相差 32 天的情况，可能会报错，请注意。同时，需要备份好你的输入数据，因为程序运行中会自动删除拼接好的数据。进行到此处则得到了初步的拼接数据，下面进行第二步 B.再进行重投影，改格式，缩小范围 1.新建记事本，在记事本中输入下面代码，改
5、txt 后缀为 bat，定义文件。rem Set the MRTDATADIR environmental var to the MRT data directory.set MRTDATADIR=C:MRTdata for%i in(*.hdf)do C:MRTbinresample-p -i%i-o%其中，MRTDATADIR 为你的 MRT 安装文件中 data 的路径 2.利用 MRT 的图形界面工具定义投影文件 选择一期刚刚拼接好的影像，在 UL Corner 和 LR Corner 可以输入经纬度来缩小你需要的区域，设置你的保存路径以及命名（最好以日期命名，到时候程序会自动按照你的规则全部按日期命名），选择tif 格式，选择你要的投影类型，点击 Save Parameter File 将投影文件保存。记住：这两个文件必须和刚刚拼接后的数据放在一个文件夹中 3.双击执行批处理即可！操作如图：总的来说，就是先进行拼接得到一个 hdf 的临时文件，然后在进行重投影，转格式，缩范围的操作。整个过程你只需要安装 MRT 即可进行，不需要再安装Cygwin 模拟 lunix 环境，安装 Cygwin 很麻烦。另外网上有小熊工具箱挺好，但是他编写的只能对同一个 tile 进行长时序列的批量处理，所以遇到需要多景影像的时候可以采用此法。PS:番茄甜土豆就是我，我为他俩代言！