MODIS系列之NDVI（MOD13Q1）五：NDVI處理流程

本文轉載自查看原文 2020-04-17 22:00 11940 應用

前言：（個人建議）

1.進行數據處理工作，由於通常數據量比較大。所以在個人電腦中，要將基礎數據、不斷增長的過程數據、結果數據等分門別類，使得簡單易找。

2.將重要數據備份，因為在數據處理過程中，由於出錯使重要數據損失，勿讓重要數據成絕版。所以生成或保存壓縮文件是很好的選擇。

3.當數據量大，需要進行批量處理時。模板文檔或代碼要杜絕出現細節錯誤。容易一鍋端，費時間且低效率。

處理工具：ArcMap

處理數據：MOD13Q1的NDVI成品（.tif）文件

研究區矢量文件（以河南省省界為例）

處理流程：

1.導入研究區矢量文件（做空心處理）

2.將.tif文件導入ArcMap中（可以多導，以電腦承受能力為主）。

是否構建金字塔：是

3.批量定義投影（MRT進行批處理時已做過，但雙保險且該步簡單）（ArcMap有很多單次或多批次投影方法，本文只介紹一種）

雙擊圖層

點擊坐標系

點擊投影坐標系

點擊UTM

點擊WGS1984

北半球選北

本文采用WGS 1984 UTM ZONE 49N。

點擊應用，再點擊確定。（如果有常用的坐標，右擊點擊坐標選擇收藏到文件夾）

4.按掩膜提取（裁剪也行，但建議用按掩膜提取）（該步驟只能一個一個處理，Python批處理暫未更新）

單擊ArcToolbox，喚出ArcToolbox，如下圖找到按掩膜提取，雙擊按掩膜提取

點擊確定

為了方便，可以將原始.tif文件移除。結果如下：

該步驟需要代碼的，代碼法如下：

同一掩膜文件，批量裁切：

import arcpy
from arcpy import env
from arcpy.sa import *

#設置環境
arcpy.env.workspace=r"E: FAN InterpolationlKrigO"
rasters=arcpy . ListRasters ("*", "GRID")#獲取文件列表

#輸出路徑
outRasterPath="E: EANInterpolationKrigO_ extractlKrigO_ ."

for ras in rasters:
    print outRasterPathtras

    # Check out any necessary licenses
    arcpy. CheckOutExtension ("Spatial")  
    
    # Local variables
    inRaster ras
    inMaskData="E : FUJIANModelBoundary.shp" #掩膜文件
    outRaster-outRasterPath-ras+".tif"
    
    # Process: Extract by Mask
    arcpy.gp. ExtractByMask_ sa (inRaster, inMaskData, outRaster)

print ("OK")