arcgis批量计算Modis数据NDVI的python代码

以下是一个基于ArcPy的Python代码，用于批量计算Modis数据的NDVI值：

import arcpy
import os

# 设置工作空间
arcpy.env.workspace = r"C:\ModisData"

# 遍历文件夹中的所有hdf文件
for file in arcpy.ListFiles("*.hdf"):
    # 解包hdf文件
    arcpy.ExtractSubDataset_management(file, "NDVI", r"C:\NDVI")
    
    # 设置NDVI文件夹为工作空间
    arcpy.env.workspace = r"C:\NDVI"
    
    # 遍历文件夹中的所有tif文件
    for tif_file in arcpy.ListFiles("*.tif"):
        # 计算NDVI
        out_raster = os.path.splitext(tif_file)[0] + "_NDVI.tif"
        arcpy.gp.SingleOutputMapAlgebra_sa("(float(\"" + tif_file + "\" - \"" + os.path.splitext(tif_file)[0] + "_1.tif\") / (\"" + tif_file + "\" + \"" + os.path.splitext(tif_file)[0] + "_1.tif\"))", out_raster)

这段代码假定Modis数据存储在`C:\ModisData`文件夹中，NDVI值将被计算并存储在`C:\NDVI`文件夹中。你需要根据自己的数据路径和文件名做出相应的修改。