MODIS系列之NDVI(MOD13Q1)九:單一土地類型NDVI及時序圖(以耕地類型為例)


單一土地類型耕地NDVI及耕地時序圖

MODIS系列之NDVI(MOD13Q1)五:NDVI處理流程 MODIS系列之NDVI(MOD13Q1)七:時間序列S-G濾波之Python  處理的都是地表各土地類型的NDVI綜合。

 本次博文則以單一土地類型(以耕地類型為例)得出耕地的NDVI及時序圖

一、所需數據

1.研究區域時間(年份)的土地利用類型數據(又叫地表覆蓋數據)下載。下載詳見該系列博文:MODIS系列之NDVI(MOD13Q1)八:土地利用數據下載

本文以2010年河南省地表覆蓋數據為例(.tif數據)

2.研究區矢量數據准備(常見矢量數據下載請見ArcMap 2 (常見矢量數據下載)

本文以河南省省界矢量數據為例(.shp數據)

二、ArcMap戰場

1.將下載好的2010年河南省地表覆蓋數據通過ArcMap打開。(注意:因為網站提供的地表覆蓋數據受投影及條帶的影響,不是按照行政區划分幅。所以下載的地表覆蓋數據要包含研究區域。

 2.打開ArcToolbox

若研究區域較大,需對幾幅地表覆蓋數據“合並”(非術語解釋)。若研究區域較小,一副地表覆蓋數據就能覆蓋研究區域,請省略"合並"步驟。由於研究區為河南,故所需“合並”。

 3.找到“數據管理工具”下的“鑲嵌”和“鑲嵌至新柵格”

 “鑲嵌”和“鑲嵌至新柵格”根據需要,進行二選一。

鑲嵌:將多個輸入柵格鑲嵌到現有柵格數據集。

 鑲嵌至新柵格:將多個柵格數據集鑲嵌到一個新的柵格數據集中。

 4.借助矢量數據,通過“裁剪”工具將“合並”后的地表覆蓋數據進行裁剪

(注:此處用的“裁剪”工具為按掩膜提取。選擇原因請參考ArcMap 3 “裁剪”)

 得到2010年河南省土地分類數據如圖:

 下載的土地分類數據有空隙,始終解決不了。如果知道問題,請通過評論區告知,非常感謝。

5.先打開2010年河南省土地分類數據的屬性表

 右擊河南土地利用分類.tif。打開屬性表。看到Value為10。

6.按屬性提取。因為耕地分類為數字“10”。MODIS系列之NDVI(MOD13Q1)八:土地利用數據(又稱地表覆蓋數據)下載

 7.在點開的按屬性提取頁面,點擊SQL欄

 8.進入“查詢構建器”

 

 9.通過給出的按鍵選項,選出"Value" =10

 

 

 得到如圖:

 10.將之前得到的綜合NDVI.tif數據,通過按掩膜提取工具提取耕地分類數據。得到耕地NDVI.tif數據。

 由於通常所需數據量大,該步驟可以通過python批量提取。也可通過ArcToolbox中按掩膜提取的批處理(右擊按掩膜提取)。

 

 

 11.得到各耕地NDVI.tif數據。后通過參考MODIS系列之NDVI(MOD13Q1)七:時間序列S-G濾波之Python 的代碼。得到耕地NDVI時序圖。(該NDVI時間范圍內,如果種植的是某種作物。那么得到的耕地時序圖,也是作物時序圖。由於該數據選擇的是3、4、5月MODIS數據。在該時間范圍內耕地種植的是小麥。所以本耕地時序圖,又叫小麥時序圖。(不考慮誤差))

 

 至此,NDVI基本系類處理結束。該系列后續會對之前的博文進行不斷更進,以及適當補充。

對本系列(包括所有博文)有建議的,請您留言,與君共勉,非常感謝。

由於最近忙於其它,更新較慢,不足之處,還請見諒。


免責聲明!

本站轉載的文章為個人學習借鑒使用,本站對版權不負任何法律責任。如果侵犯了您的隱私權益,請聯系本站郵箱yoyou2525@163.com刪除。



 
粵ICP備18138465號   © 2018-2025 CODEPRJ.COM