MRT(MODIS Reprojection Tool)安装、影像批量拼接、重投影和格式转换

一、安装MRT(MODIS Reprojection Tool)

  安装准备:检查是否安装java.exe。Java版本至少为Java 2 Runtime Environment version 1.5或者是Java 2 SDK version 1.5或者更高的版本。在Windows平台中可以通过开始/搜索并输入java.exe来检查自己的电脑是否安装有java,也可以不检查直接重新安装一份。

  解压安装包MRT_download_<platform>.zip。包含四个文件:mrt_install.bat,MRT_<platform>.zip,reg_set.exe, 以及unzip.exe。

  双击mrt_install.bat,并按提示完成安装。

二、使用MRT(MODIS Reprojection Tool)

(一)GUI界面操作   

  安装完成后,打开MRT的安装路径文件夹,打开bin文件夹,双击ModisTool.bat或者ModisTool.jar即可进 入MRT的GUI界面。

  输入MODIS数据,选择输出NDVI和EV,指定输出路径输出文件格式为GEOTIFF,采样方式为最近邻,投影类型为Albers Equal Erea,输出像元大小为250m,然后直接点击run就可以了。其中,投影参数设置见下图(STDPR1第一条标准纬线27,STDPR2第二条标准纬线为45,CentMer中央经线为105,datum基准面WGS84):

(二)批处理拼接之一年以内的数据

  备注:该操作数据及BAT文件须放在一个文件夹下(文件夹起名请用英文,MRT对中文路径数据会打不开)。

  操作步骤:

  1.新建一个DATA文件夹,将HDF数据拷到DATA文件夹下;

  2.在DATA文件夹下再新建个Result文件夹(存放拼接后数据);

  3.将以下代码内容粘到一个TXT文件中,更改TXT文件扩展名为BAT即可,假设你命名该BAT文件名为MOSAIC.BAT;

  4:点击MOSAIC.BAT,即可进行拼接操作。

  注意:第一行代码里,把MRTDATADIR替换为你自己的MRTDATADIR名称,MRTDATADIR名称可以在安装目录里,找到ModisTool.bat,然后右键选择编辑中查看;

此外以下代码未考虑中间天数间断的情况,比如对MODIS时间分辨率为16天的数据,在一个应连续的等差数列中间有数据缺失,存在相差32天的情况,可能会报错,请注意。

 1 set MRTDATADIR=f:\MOD_HDE_SOURCE 2
 3 set /a DAY=2000049        rem ***开始日期***
 4
 5 set /a DEADLINE=2000353  rem ***结束日期***
 6
 7 :start
 8
 9 if %DAY% leq %DEADLINE% (goto ORDER) else exit
10
11 :ORDER
12
13 rem **将当天的图幅数据文件名放在一个TXT文件中**
14
15 dir *%DAY%.*.hdf/a/b/s > MOSAICINPUT.TXT
16
17 rem **拼接HDF12种数据中的第一种**
18
19 d:/software_install/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
20
21 rem **将拼接后数据复制到已建Result 文件夹中并删除当天的数据**
22
23 copy MOSAIC_TMP_%DAY%.hdf  Result & del MOSAIC_TMP_%DAY%.hdf
24
25 del *%DAY%.*.hdf
26
27 set /a DAY= %DAY% + 16   rem ***拼接下一个时相的数据***
28
29 goto start

(三)批处理拼接之多年的数据

  以上代码只能处理一年之内或者每隔一年的数据,如果要一次性处理多年的数据,则参见以下的代码

 1 set MRT_DATA_DIR=f:\MOD_HDF_SOURCE
 2
 3 set /a DAY=2015337
 4
 5 set /a DAY1=337
 6
 7 set /a DEADLINE=2016353
 8
 9 :start
10
11 if %DAY% leq %DEADLINE% (goto ORDER1) else exit
12
13 :ORDER1
14
15 if %DAY1% leq 353 (goto ORDER)
16
17 :ORDER
18
19 dir *%DAY%.*.hdf/a/b/s > MOSAICINPUT.TXT
20
21 d:/software_install/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
22
23 copy MOSAIC_TMP_%DAY%.hdf  Result & del MOSAIC_TMP_%DAY%.hdf
24
25 del *%DAY%.*.hdf
26
27 set /a DAY1= %DAY1% + 16
28
29 set /a DAY= %DAY% + 16
30
31 if %DAY1% GTR 357 (goto ORDER2)
32
35 goto start
36
38 :ORDER2
39
41 set /a DAY=%DAY% + 1000 - 368
42
43 set /a DAY1=1
44
45 pause
46
47 goto start 

(四)批处理重投影及格式转换

  新建一个文件夹,假设起名为DATA,把所有要处理的MODIS HDF文件放到该文件夹下;

  在GUI界面中,设定一个HDF文件投影参数后,点击“Save Parameter File…”保存自己设定的参数文件到DATA文件夹下,假设为test.prm;

  在DATA文件夹里,新建一个记事本文件,将下面的内容粘贴进去,需更改几处。更改完毕之后,将记事本文件另存为“.bat”文件,假设为“modis.bat”;

1 set MRT_DATA_DIR=f:\MOD13Q1_REPROJECTION
2
3 for %%i in (*.hdf) do resample.exe -p test.prm -i %%i -o %%iout.tif
4
5 pause

  其中,

  第一行代码里,把MRTDATADIR替换为你自己的MRTDATADIR名称,MRTDATADIR名称可以在安装目录里,找到ModisTool.bat,然后右键选择编辑中找到;

  第二行代码里,把test.prm替换为自己的投影参数文件,如果运行出错的话,可以把resample替换为其完整的路径,如D:\software_install\MRT\bin\resample.exe;

  第三行代码为运行完成后保留bat运行窗口,如果出现错误,可参见bat窗口的提示,将安装目录下的datum.txt和spheroid.txt文件复制到DATA文件夹里;

  

  双击modis.bat,便可看到命令窗口的批出理过程,如发生异常,检查prm的参数是否设置正确(可在GUI中运行看是否报错)或见上一段的注意事项。

时间: 2024-10-22 03:19:24

MRT(MODIS Reprojection Tool)安装、影像批量拼接、重投影和格式转换的相关文章

ffmpeg下载安装和简单应用(C#音频格式转换)

ffmpeg下载安装和简单应用 先介绍一下ffmpeg:FFmpeg是一个自由软件,可以运行音频和视频多种格式的录影.转换.流功能,包含了libavcodec -这是一个用于多个项目中音频和视频的解码器库,以及libavformat--一个音频与视频格式转换库.名字由来:"FFmpeg"这个单词中的"FF"指的是"Fast Forward"[2].有些新手写信给"FFmpeg"的项目负责人,询问FF是不是代表"Fas

MODIS-MRT(MODIS Reprojection Tool)

1.MRT_Downloads:https://files.cnblogs.com/files/haijian/MRT_download_Win.zip(NASA已经不再提供软件支持,官网已下架) 2.MRT_install: ①Install java-jre:见解 ②Install MRT:解压 -> 双击mrt_install.bat -> 输入安装目录 -> "y"确认创建目录 -> “1”win -> 输入jre\bin路径 -> &quo

SuperMapPy 批量拼接 GeoTiff影像

影像拼接工具使用说明 一.影像像素位深检查 1.采用开源库GDAL的gdalinfo.exe读取GeoTiff文件的信息,如类型.投影,范围等. 2.采用DOS脚本遍历所有GeoTiff文件,输出各个GeoTiff文件的信息到日志. 3.统计GeoTiff文件个数,Type=Byte的文件的个数.数目一致表明全部GeoTiff像素类型一致. 运行示例: for /f %i in ('dir /b /s d:\tifdata\*.tif ') do@( e:\gdal1.4.2\bin\gdali

ArcGIS Server 10.2 实战(五)spatial etl tool 格式转换服务

上不同的地图服务平台对地图文件格式的要求多种多样,arcgis使用的文件很难应用于其他平台上,因此需要有格式转换的服务来克服这种使用不同平台带来的麻烦,下面以TIFF格式转GEOTIFF格式为例. 首先需要准备几件事: 1.确保安装了arcgis data interoperability for desktop 2.在自定义菜单>扩展模块中勾选Data Interoperability 3.在ArcToolbox中新建一个自己的工具箱,右键工具箱>新建>选择spatial etl to

批量创建用户与批量文件重命名

title: 批量创建用户与批量文件重命名 tags: 批量,用户,创建,重命名 批量创建用户与批量文件重命名 1. 批量创建用户 如果需要批量创建用户,且用户名为stu01-stu10,用户密码为随机的10位密码. 1.1 分析解答 如果用户名为stu01的话,我们可以采用for循环,直接变量范围为01-10,密码随机的话我们可以采用系统默认的随机变量RANDOM,因为默认的RANDOM只有5位数字,因此我们可以给它md5加密,然后取其中的8位字符当做密码.当然我们也可以不用循环,直接使用命令

网页版批量文件重命名

功能:这是一个网页版的文件批量重命名 操作步骤: 输入文件所在的绝对路径, 输入你想要替换的关键字, 输入你想要用什么字符来替换他(可为空), 点击确定,文件将会批量被重命名, 点击撤回就会还原出起初的文件名字, 点击清空主要是一键清除文本框里面的内容. 总的有1个jsp,四个Java程序共同完成这个功能: start.jsp(显示界面)   CodeFilter.java(过滤器,主要负责编码)   ReceiveStartServlet.java(接受启动,进行文件操作) RevokeSer

使用Spire.Office自制一个进行批量格式转换(Pdf转换为Word)的工具

使用Spire.Office自制一个进行批量格式转换(Pdf转换为Word)的工具 本文文末提供一个可将PDF转换为Word的小工具供下载. +BIT祝威+悄悄在此留下版了个权的信说: Spire.Office 我一直想找一个工具能够把我收集的一些PDF资料转换为Word版.最近偶然发现一个超级强大的工具Spire.Office for .NET可以做到这一点. 而且Spire.Office提供C#.VB的DLL供调用,还提供了丰富简洁的示例项目可供学习. 如下图所示,Spire.Office

asp.net 文件批量移动重命名

最近闲时写了个批量移动重命名文件的工具 点击下载工具

用Python写的批量文件重命名

  有些时候下载图片或其他文件,文件名都怪怪的,可选的办法是下载一个文件批量重命名的软件.当然,如果想自己'DIY'一把的话编个Python脚本最好不过了. 下面的代码实现的对指定类型的文件进行批量重命名.拷贝下面的代码到待批量命名的文件夹下,保存为xx.py直接运行,程序会提示需要批量命名的扩展名,以及重命名时的文件前缀. # -*- coding: cp936 -*- """ Created on Wed Jun 25 16:24:23 2014 @author: Adm