使用ArcGIS的Python模块ArcPy对栅格图按属性批量导出的代码

使用ArcGIS的Python模块ArcPy对栅格图按属性批量导出的代码

在Arcgis中，想对栅格的某一部分按属性提取，可以用“按属性导出”工具；需要同时提取多个属性时，可以用“按属性导出”工具的批处理。
但是，当你有一个包含多个字段的栅格图，需要对某一字段的全部属性值批量导出，数据量很大时（以千计），批处理也很麻烦，特别是如果需要导出的字段中每个属性值都包含了多个栅格值时，如何解决？

准备：

（1）一个包含要导出字段所有值的CSV文件，您可以使用这个文件来确保您的脚本处理所有唯一的字段值。
（2）需要操作的栅格数据
（3）工作路径
（4）ArcGIS的Python模块ArcPy

步骤：

（1）读取CSV文件：首先，您需要使用Python来读取CSV文件，并提取出所有唯一的值。
（2）使用读取的value值：然后，您可以使用这些从CSV文件中提取的value值来替换原本直接从栅格数据中获取的value值。

代码：

import arcpy
import csv

# 设置工作空间
arcpy.env.workspace = "D:/output"
inRaster = "10km.tif"
csv_file = "value.csv"
# 从CSV文件中读取value值
uniquevalues = []
with open(csv_file, mode='r') as infile:
    reader = csv.reader(infile)
    next(reader, None)  # 跳过表头
    for row in reader:
        if row:  # 确保行不是空的
            id_value = row[0]  # 假设ID100在第一列
            if id_value not in uniqueIDs:
                uniqueIDs.append(id_value)

# 检查读取的value值
print(f"Unique values from CSV: {uniquevalues}")

# 使用CSV中的value值进行处理
for id in uniquevalues:
    try:
        print(f"Processing value={id}")
        queryString = f"value = {id}"
        outRaster = ExtractByAttributes(inRaster, queryString)
        output_path = f"D:/output/value_{id}.tif"
        outRaster.save(output_path)
        print(f"Successfully processed value={id}")
    except Exception as e:
        print(f"Error processing value={id}: {e}")