前言:

本篇文章的出发点是因为之前接触过相关研究,困囧于该系列资料匮乏,想做一个系列。个人道行太浅,不足之处还请见谅。愿与诸君共勉。

数据准备:

MODIS数据产品MOD13Q1—以2010年河南省3、4、5三个月为例:

一、MRT单次数据操作

(1).进入GUI界面操作

1.将所需处理的一个.hdf原始数据加载进来

2.通过左右选项选择所需波段(MOD13Q1已将NDVI提供,只需保留选择就行)

(若用的数据比如MOD12Q1计算NDVI,则Modis算法如公式:NDVI=(Band2-Band1)/(Band2+Band1)。那么就需要将Band2和Band1波段提取出来)

3.Spatial Subset(空间子集):选择Input Lat/Long (输入纬度/经度)

input line/sample (输入行/样本)

output projection (X/Y输出投影X / Y)

4和5在一起

4.选择文件保存路径:

5.输出数据类型:  如我保存的路径(必须同原始数据也就是.hdf在同一文件夹下)就是F:\MODIS\.tif   (注意,直接在MODIS文件夹后加\.tif就行。生成的.tif文件名直接等同于在同一个文件夹下.hdf文件同名(建议相同操作,通常数据多,命名易辨且重要))

6.输出文件类型:GEOTIFF(我们要的就是.tif数据文件)

7.Resampling(重采样类型):选择Nearest Neighbor(最近领域)

8.Output Projection type(输出投影类型):选择原则是下载的MODIS数据(.hdf)是什么投影就选什么投影(百度)

本文选用MOD13Q1,应当采用Sinusoidal投影方式。但大部分研究选用的都是Albers Equal Area。终究道行太浅,本文也不能免俗。

9.Edit Projection Parameters(编辑主动参数):

上图的一些参数STDPR1和2不太清楚,参考同是研究区为河南省的研究输入的

centmer为中央经线(研究区中央经线)

大方框里的可以什么都不输入,默认就行。但WGS84一定要选,点击ok

10.Output Pixel Size(输出像素大小):像素范围越小,精度越高。但输出的.tif文件就越大,会达到1-2个G。选择建议:看自己选择的MODIS产品是什么像素范围。

11.点击Run即可

二、MRT多批次数据操作(重点)

注意:无论你有多少期数据。在进行多批量处理时。用到MRT文件时只输如一个.hdr文件。不用担心,之所以输入一个.hdr文件,是因为在进行多批量处理时,用到MRT工具只是通过一个.hdr文件来建立一个.prm文件作为模板。用CMD来批处理

(1).MRT战场

1-10同上文一致(只输入一个.hdr文件)

11.点击SaveParameter File

在文件夹内将会出现.prm文件。那么接下来咱们转移战场。将MAT给Exit关掉

(2)CMD战场

以下所提及的功能都可以在MRT批量进行,通过基于Dos系统的批处理脚本(后缀.bat)来实现。将待处理的MODIS数据(.hdf)、MRT参数设置文件(.prm)、批处理脚本(.bat)放置在一个文件夹中。运行.bat后,系统就会根据脚本的命令批量执行MRT中的各个功能,从而实现将.hdf文件按设置转换为.tif文件

1.F:\MODIS\.tif 文件路径下就生成了.prm文件

2.用Notepad++打开或记事本打开MOD13Q1NDVI.prm

INPUT_FILENAME = F:\MODIS\MOD13Q1.A2010049.h26v05.006.2015200022349.hdf

SPECTRAL_SUBSET = ( 1 0 0 0 0 0 0 0 0 0 0 0 )  //该MODIS数据产品有12个波段1代表选,0代表不选该波段

SPATIAL_SUBSET_TYPE = INPUT_LAT_LONG

SPATIAL_SUBSET_UL_CORNER = ( 39.999999996 104.432583132 )
SPATIAL_SUBSET_LR_CORNER = ( 29.999999997 103.923048442 ) OUTPUT_FILENAME = F:\MODIS\.tif //输出路径(F:\MODIS\.tif) 以及输出的数据类型(. tif) RESAMPLING_TYPE = NEAREST_NEIGHBOR //重采样类型选择最邻近法 OUTPUT_PROJECTION_TYPE = AEA //输出投影类型选择墨卡托(AEA) OUTPUT_PROJECTION_PARAMETERS = (
0.0 0.0 25.0
47.0 105.0 0.0
0.0 0.0 0.0
0.0 0.0 0.0
0.0 0.0 0.0 ) //输出投影参数建议为默认 DATUM = WGS84 //基准面为WCG84. UTM_ZONE = 49 //平面投影为49N OUTPUT_PIXEL_SIZE = 250 //像元大小为500*500

生成好后的.prm可以通过Notepad++或记事本查看,不动

3.通过Notepad++或记事本另存为一个modis.bat文件(.hdr .tif .bat .prm文件都是在同一个文件夹下)

modis.bat文件代码如下

set MRT_DATA_DIR=C:\MRT\data

for %%i in (*.hdf) do resample -p MOD13Q1NDVI.prm -i %%i  -o  %%iout.tif

pause 

第1行:MODIS批量重投影脚本代码参考上图,使用时需要设置MRT的安装路径(C:/MRT/data)。安装路径请查看自己电脑高级设置里的环境变量。

第3行:文件.prm文件名替换自己的.prm文件名,如果运行出错的话,可以把resample替换为其完整的路径,如D:\software_install\MRT\bin\resample.exe;

4.另存为一个modis.bat文件放在同一个文件夹下。如果出现错误,可参见bat窗口的提示,将安装目录下的datum.txt和spheroid.txt文件复制到DATA文件夹里;

5.双击modis.bat文件。运行完即可。如果运行错误,因为双击modis.bat文件信息一闪而逝,不能看到错误信息。建议打开该目录文件夹下将modis.bat输入

批处理运行过程如下:

运行后关闭即可。

MODIS系列之NDVI(MOD13Q1)四:MRT单次及批次处理数据的更多相关文章

  1. MODIS系列之NDVI(MOD13Q1)七:时间序列S-G滤波之Python

    时间序列S-G滤波之Python 根据上上篇博文(MODIS系列之NDVI(MOD13Q1)五:NDVI处理流程)做出的NDVI.我们求NDVI时间序列图,但该NDVI时序图为地表各土地类型综合的ND ...

  2. MODIS系列之NDVI(MOD13Q1)三:.jdk文件配置+MRT安装

    MRT(MODIS Reprojection Tool)简介: MODIS的全称为中分辨率成像光谱仪(Moderate-Resolution Imaging Spectroradiometer),是搭 ...

  3. MODIS系列之NDVI(MOD13Q1)五:NDVI处理流程

    前言:(个人建议) 1.进行数据处理工作,由于通常数据量比较大.所以在个人电脑中,要将基础数据.不断增长的过程数据.结果数据等分门别类,使得简单易找. 2.将重要数据备份,因为在数据处理过程中,由于出 ...

  4. MODIS系列之NDVI(MOD13Q1)二:modis数据相关信息

    1.MODIS数据的特点 (1)全球免费:NASA对MODIS数据实行全球免费接收的政策(TERRA卫星除MODIS外的其他传感器获取的数据均采取公开有偿接收和有偿使用的政策),这样的数据接收和使用政 ...

  5. MODIS系列之NDVI(MOD13Q1)一:数据下载(一)基于插件

    引言: 写MODIS数据处理这个系列文章的初衷,主要是为了分享本人处理MODIS数据方面的一些经验.鉴于网上对这方面系统性的总结还比较少,我搜集资料时也是走了许多的弯路,因此希望通过此文让初学者能够更 ...

  6. MODIS系列之NDVI(MOD13Q1)一:数据下载(二)基于FTP

    这一篇我们来介绍下MODIS数据的下载方式.当然这边主要是介绍国外网站的下载方式,国内网站的普遍是在地理空间数据云和遥感集市下载.国外网站(NASA官网)下载方式主要介绍两种.本篇主要针对第一种方式, ...

  7. BootStrap 智能表单系列 四 表单布局介绍

    表单的布局分为自动布局和自定义布局两种: 自动布局就是根据配置项中第二级配置项中数组的长度来自动使用不同的bootstrap栅格,通过设置autoLayout为true可以实现自动布局 自动以布局就是 ...

  8. 【HANA系列】【第四篇】SAP HANA XS使用服务器JavaScript Libraries详解

    公众号:SAP Technical 本文作者:matinal 原文出处:http://www.cnblogs.com/SAPmatinal/ 原文链接:[HANA系列][第四篇]SAP HANA XS ...

  9. Spark JDBC系列--取数的四种方式

    Spark JDBC系列--取数的四种方式 一.单分区模式 二.指定Long型column字段的分区模式 三.高自由度的分区模式 四.自定义option参数模式 五.JDBC To Other Dat ...

随机推荐

  1. Python习题集(四)

    每天一习题,提升Python不是问题!!有更简洁的写法请评论告知我! https://www.cnblogs.com/poloyy/category/1676599.html 题目 如果一个 3 位数 ...

  2. annaconda的安装及使用

    一.安装 1.安装包下载: 方式1:官网:https://www.anaconda.com/distribution/ 方式2:清华镜像源:https://mirrors.tuna.tsinghua. ...

  3. 安卓App自动化测试环境

    一.appium安装 1.nodejs依赖 2..net framework 3.Appium桌面程序安装 3.1.安装包获取 3.2.安装过程 二.Java_ jdk安装 1.Java_jdk版本说 ...

  4. python3爬虫爬取金庸小说所有角色

    # -*- coding:utf-8 -*- import requests from bs4 import BeautifulSoup url = 'http://www.jinyongwang.c ...

  5. 毫米波大规模阵列中的AOA估计

    1.AOA估计在毫米波大规模MIMO中的重要性 在毫米波大规模MIMO的CSI估计中,AoA估计具有重要地位,主要原因归纳如下: 毫米波大规模MIMO 的信道具有空域稀疏性,可以简单通过AoA 和路径 ...

  6. 题解 P5663 【加工零件【民间数据】】

    博客园体验更佳 讲讲我的做法 确定做法 首先,看到这道题,我直接想到的是递归,于是复杂度就上天了,考虑最短路. 如何用最短路 首先,看一张图 我们该如何解决问题? 问题:\(3\)做\(5\)阶段的零 ...

  7. nopcommerce4.0 安装步骤

    前言:近期因工作要求接触nopcommerce,最新版本为4.0,以下所有安装都是基于此版本.接下来我可能会写一系列,为了让自己更好的掌握,也希望能帮助到大家 好记性不如烂笔头,新手也可以避免走我的弯 ...

  8. 【Unity游戏开发】跟着马三一起魔改LitJson

    一.引子 在游戏开发中,我们少不了和数据打交道,数据的存储格式可谓是百花齐放,xml.json.csv.bin等等应有尽有.在这其中Json以其小巧轻便.可读性强.兼容性好等优点受到广大程序员的喜爱. ...

  9. OpenCV-Python 直方图-1:查找、绘制和分析 | 二十六

    目标 学会 使用OpenCV和Numpy函数查找直方图 使用OpenCV和Matplotlib函数绘制直方图 你将看到以下函数:cv.calcHist(),np.histogram()等. 理论 那么 ...

  10. 使用PyTorch进行迁移学习

    概述 迁移学习可以改变你建立机器学习和深度学习模型的方式 了解如何使用PyTorch进行迁移学习,以及如何将其与使用预训练的模型联系起来 我们将使用真实世界的数据集,并比较使用卷积神经网络(CNNs) ...