保姆级教程：用ArcGIS Pro的字段计算器，给DEM和地形起伏度分类地貌（附避坑指南）

从DEM到地貌类型：ArcGIS Pro高效分类实战指南

当你手握DEM数据和地形起伏度栅格，却卡在如何将它们转化为直观的地貌类型名称时，这篇教程将成为你的救星。我们将跳过基础操作，直击属性表联合与字段计算器这两个核心环节，解决从数据到图例的"最后一公里"问题。

1. 数据准备与字段创建

在开始地貌分类前，确保已完成DEM和地形起伏度的重分类。假设你已经得到两个栅格图层：DEM_Reclass和Relief_Reclass，并完成了它们的联合操作。

关键检查点：

• 确认联合后的图层属性表包含gridcode和gridcode_1字段
• 检查两个字段的值是否与你的重分类方案一致

注意：ArcGIS Pro会自动将重名字段添加后缀"_1"，gridcode通常对应地形起伏度，gridcode_1对应DEM数据

创建地貌类型字段的操作如下：

1. 右键联合后的图层，选择属性表
2. 点击表选项按钮，选择添加字段
3. 设置字段参数：
   • 名称：Landforms
   • 类型：文本
   • 长度：20（根据实际需要调整）

# 通过ArcPy实现字段创建的等效代码 import arcpy arcpy.AddField_management("联合图层", "Landforms", "TEXT", field_length=20)

2. 构建精准的SQL选择表达式

地貌分类的核心在于正确构建SQL表达式来选择特定组合的区域。以下是详细步骤：

操作流程：

1. 在属性表中点击按属性选择按钮
2. 在表达式构建器中，按以下顺序操作：
   • 双击gridcode字段（地形起伏度）
   • 点击=运算符
   • 点击获取唯一值，选择需要的数值
   • 点击And按钮添加逻辑连接
   • 双击gridcode_1字段（DEM）
   • 重复运算符和值选择过程

常见地貌类型对应值示例：

提示：表达式构建完成后，可以先点击验证检查语法，确认无误再应用选择

3. 字段计算器的精准赋值

选中目标区域后，通过字段计算器为Landforms字段赋值是分类的关键步骤：

1. 右键Landforms字段标题，选择字段计算器
2. 在表达式框中输入地貌类型名称，必须使用英文双引号，例如：
   • "低海拔平原"
   • "中海拔山地"
3. 点击确定完成赋值

# 等效的ArcPy赋值代码 arcpy.CalculateField_management("联合图层", "Landforms", '"低海拔平原"', "PYTHON3")

常见错误及解决方案：

• 错误1：使用中文引号导致语法错误
  解决：确保所有引号为"而非"

• 错误2：字段计算后显示<Null>
  解决：检查是否先做了正确选择，或尝试重新选择

• 错误3：提示字段长度不足
  解决：删除字段后重新创建，设置更大长度

4. 批量处理的效率技巧

手动逐个选择并赋值效率低下，以下是几种提升效率的方法：

方法一：使用预定义函数

1. 在字段计算器中选择Python解析器
2. 使用以下代码模板：

def classify_landform(relief, dem): if relief == 1 and dem == 1: return "低海拔平原" elif relief == 1 and dem == 2: return "中海拔平原" elif relief == 2 and dem == 1: return "低海拔台地" elif relief == 3 and dem == 3: return "高海拔山地" else: return "其他类型"

然后在表达式框中调用：

classify_landform(!gridcode!, !gridcode_1!)

方法二：创建选择集批量处理

1. 将所有分类组合的SQL表达式保存为文本文件
2. 使用以下Python脚本批量处理：

import arcpy layer = "联合图层" field = "Landforms" classifications = [ {"name": "低海拔平原", "sql": "gridcode = 1 AND gridcode_1 = 1"}, {"name": "中海拔平原", "sql": "gridcode = 1 AND gridcode_1 = 2"}, # 添加其他分类... ] for cls in classifications: arcpy.SelectLayerByAttribute_management(layer, "NEW_SELECTION", cls["sql"]) arcpy.CalculateField_management(layer, field, f'"{cls["name"]}"', "PYTHON3")

5. 成果可视化与验证

完成所有赋值后，需要合理展示分类结果：

1. 右键图层选择属性
2. 切换到符号系统选项卡
3. 选择唯一值渲染方式
4. 值字段选择Landforms
5. 点击添加所有值按钮
6. 为不同类型设置合适颜色

质量检查要点：

• 检查属性表中是否有未赋值的记录
• 验证边界区域分类是否合理
• 对比实际地形与分类结果的匹配度

常见问题排查表：

完成所有步骤后，建议导出图层属性表进行备份，同时保存工程文件时勾选保存相对路径，方便后续共享和迁移项目。