Modis 植被数据的批量处理
Windows 下 MRT 批 处 理 MODIS 数 据 — — 拼 接 和 重 投 影 等
注:本文承接《 MRT(MODIS Reprojection Tool) 下载及安装教程》一文。
1.MRT 安装
LP DAAC 最新发布的 MRT 支持 windows 平台下处理 MODIS 影像,能够完全不需要通过 cygwin 来模
拟 LINUX 系统。MRT 下载地址如下:
https://lpdaac.usgs.gov/tools/modis_reprojection_tool
在成功注册之后,就可以点击下载 Windows NT+32-bit
解压之后得到 MRT_Win.zip、mrt_install.bat、reg_set.exe、unzip.exe 四个文件,运行 mrt_install.bat
文件,然后按照提示输入 MRT 和 JAVA 安装路径就可以正确安装并运行 MRT。
ps:安装前,需要确定电脑是否安装了 java,若没有需要去 http://java.sun.com 下载最新版的 java,然后记
住安装路径,以便之后的输入。
2.数据准备
1)本人的 MRT 安装在 C:\MRTmodis\MRT 目录下,
运行 C:\MRTmodis\MRT\bin\Modistool.bat 进入 MRT GUI 界面,
选择一副影像制作批量处理需要的*.prm 文件。
cmd 命令行方式启动 MRT:我的安装目录为 C:\MRTmodis\MRT\bin
2)在界面输入数据,设置所需参数,并保存成 .prm 文件
选择想要提取的波段
声明输出文件的路径,指定输出目录和命名(如果要在 ENVI 中打开,建议保存成.hdr 文件);
设置地理坐标或者投影坐标等参数,本人以 UTM(通用横轴墨卡托)投影为例。(通用
横轴墨卡托是横轴墨卡托的改进,其与横轴墨卡托或高斯-克吕格的主要区别在与 UTM 的
比例因子为 0.9996。在 MRT 中 UTM 的参数设置比较简单:
Lon/Z:区域内任意一点的经度坐标,也可以设置为零,如果设置为 0,那么 UTM Zone
中必须填入一个分区代码。
Lat/Z:区域内任意一点的纬度坐标,也可以设置为 0,如果设置为零,那么 UTM Zone 中
必须填入一个分区代码);
设置输出数据的经纬度范围, UL CORNER = ( 37.5
89 )
LR CORNER = (31 103 )——初步裁剪。
保存 prm 文件,并执行 Run 操作。
3.cmd 命令实现 MRT 批处理
运行 cmd 命令,将工作目录设置到 C:\MRTmodis\MRT\bin 中,即 MRT 安装目录中的 bin 文
件夹中。
输入 java -jar MRTBatch.jar
-d G:\modis2000-2018-05-10\2018-05-10
-p G:\modis2000-2018-05-10\2018-05-10\2000.prm -o G:\modis2000-2018-05-10\2018-05-10
其中,-d 表示的是影像数据存储的目录,-p 表示经过 MRT GUI 处理的 prm 文件路径,
-o 表示输出路径。这串命令表示的是对所有的影像数据批处理得到每个影像的拼接和重投
影的 prm 文件。
运行成功并得到所有影像的 prm 文件后,继续在输入 MRTBatch.bat 命令,执行这个 bat
文件,即进行影像的批处理。
案例具体代码如下:
C:\Windows\System32> CD C:\MRTmodis\MRT\bin
C:\MRTmodis\MRT\bin> java -jar MRTBatch.jar
-d G:\modis2000-2018-05-10\2018-05-10
-p G:\modis2000-2018-05-10\2018-05-10\2000.prm -o G:\modis2000-2018-05-10\2018-05-10
C:\MRTmodis\MRT\bin> mrtbatch.bat
得到的其中结果如下图所示:
ArcGIS Python 实 现 Modis NDVI 批 量 化 月 最 大 合 成
1. 批量处理过程及代码
最大合成法(MVC)可以在 Envi 中的 Band Math 中进行,式子是 B1>B2,但是无法批量化;本文
实现在 ArcGIS 中利用 Python 代码批量进行,如下:用到的 Modis NDVI 数据是在 MRT 中进行拼
接与转投影后的月数据,一个月有两期,数据格式是.tif,文件名的格式如:20000501_NDVI.tif,