arcgis批量计算TIF数据1和2波段NDVI的python代码

以下是可以用于在ArcGIS中批量计算TIF数据中1和2波段的NDVI的Python代码：

import arcpy

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

# 获取所有TIF文件
rasters = arcpy.ListRasters("*", "TIF")

# 循环处理每个TIF文件
for raster in rasters:
    # 定义输出文件名
    out_ndvi = raster[:-4] + "_NDVI.tif"

    # 定义NDVI计算公式
    ndvi_formula = "(Float(\"{0}\\{1}\\{2}\") - Float(\"{0}\\{1}\\{3}\")) / (Float(\"{0}\\{1}\\{2}\") + Float(\"{0}\\{1}\\{3}\"))".format(arcpy.env.workspace, raster, "Band_2", "Band_1")

    # 执行计算
    arcpy.gp.SingleOutputMapAlgebra_sa(ndvi_formula, out_ndvi)

print("NDVI计算完成")

该代码假设TIF文件中的第1波段为红色和第2波段为近红外色，并使用NDVI公式计算每个像元的NDVI值。计算结果将保存在与原始文件相同的位置，文件名为“原始文件名_NDVI.tif”。

相关问题

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

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

import arcpy
import os

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

# 遍历文件夹中的所有tif文件
for file in arcpy.ListFiles("*.tif"):
    # 计算NDVI
    out_raster = os.path.splitext(file)[0] + "_NDVI.tif"
    arcpy.gp.SingleOutputMapAlgebra_sa("(float(\"" + file + "\" - \"" + os.path.splitext(file)[0] + "_1.tif\") / (\"" + file + "\" + \"" + os.path.splitext(file)[0] + "_1.tif\"))", out_raster)

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

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`文件夹中。你需要根据自己的数据路径和文件名做出相应的修改。