MODIS系列之NDVI（MOD13Q1）四：MRT單次及批次處理數據

前言：

本篇文章的出發點是由於以前接觸過相關研究，困囧於該系列資料匱乏，想作一個系列。我的道行太淺，不足之處還請見諒。願與諸君共勉。

數據準備：

MODIS數據產品MOD13Q1—以2010年河南省三、四、5三個月爲例：

 1、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便可

 

2、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輸入

 

 

 

批處理運行過程以下：

 

運行後關閉便可。