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

前言：（个人建议）

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")