李在田1. 认识ArcToolbox中的3D Analyst工具箱
作为一名长期使用ArcGIS进行空间分析的老用户，我深刻体会到3D Analyst工具箱在地形分析和栅格处理中的独特价值。这个工具箱最初是作为ArcGIS的扩展模块出现的，经过多年发展已经成为三维空间分析的核心工具集。特别是在处理数字高程模型（DEM）、地形可视化和表面分析等场景时，3D Analyst提供的工具链几乎无可替代。
3D Analyst工具箱主要包含以下几类功能：
- 表面创建与插值（如Topo to Raster、IDW等）
- 表面分析（坡度、坡向、山体阴影等）
- 体积计算与视线分析
- 栅格重分类与计算
- TIN（不规则三角网）的创建与编辑
其中，栅格计算功能是我们日常工作中使用频率最高的工具之一。它允许我们通过数学表达式对多个栅格图层进行代数运算，实现复杂的地形指标计算。比如计算地形湿度指数（TWI）、进行多因子叠加分析等场景都离不开栅格计算。
提示：使用3D Analyst功能前，请确保已在ArcMap中激活该扩展模块（自定义→扩展模块→勾选3D Analyst）
需要模型API调用？ 免费领10W Token，多模型网关一键接入 Claude、DeepSeek 等主流模型。
2. 栅格计算器的核心功能解析
2.1 栅格计算器的基本界面
栅格计算器（Raster Calculator）位于3D Analyst工具箱的"栅格计算"分类下。打开后会看到一个简洁但功能强大的界面，主要包含以下区域：
- 表达式输入框：用于输入栅格计算的数学表达式
- 栅格图层列表：显示当前地图文档中可用的栅格数据集
- 运算符面板：提供算术、逻辑、关系等各类运算符
- 函数列表：包含Con、IsNull等常用的栅格函数
- 输出设置：指定结果栅格的保存位置和名称
2.2 栅格计算的基本语法
栅格计算使用的表达式语法类似于常规的数学表达式，但有一些特殊规则需要注意：
- 栅格图层引用：使用方括号包裹图层名称，如"[dem]"
- 运算符优先级：与常规数学运算一致（括号优先，然后是乘除，最后加减）
- 函数调用：函数名后跟括号内的参数，如"Con([slope] > 30, 1, 0)"
一个典型的栅格计算表达式示例：
code复制([dem] - FocalStatistics([dem], NbrRectangle(3,3), "MINIMUM")) /
(FocalStatistics([dem], NbrRectangle(3,3), "MAXIMUM") -
FocalStatistics([dem], NbrRectangle(3,3), "MINIMUM"))
这个表达式计算了DEM的地形位置指数（TPI），展示了如何使用栅格计算器进行复杂的邻域分析。
3. 栅格计算的实战应用案例
3.1 案例一：地形湿度指数（TWI）计算
地形湿度指数是描述地表水文特征的重要指标，计算公式为：
code复制TWI = Ln(As / tan(Slope))
其中As为上游集水面积，Slope为坡度。
在ArcGIS中实现这一计算的完整步骤如下：
-
首先计算流向（Flow Direction）：
- 使用"流向"工具（Flow Direction），输入DEM数据
- 输出结果表示每个像元水流向的8个可能方向
-
计算汇流累积量（Flow Accumulation）：
- 使用"流量"工具（Flow Accumulation），输入上一步的流向栅格
- 结果表示流入每个像元的上游像元总数，即As
-
计算坡度（Slope）：
- 使用"坡度"工具，输入原始DEM
- 注意选择输出单位为"度"而非"百分比"
-
栅格计算器表达式：
code复制Ln( ( [FlowAcc] + 1 ) / Tan( [Slope] * 3.1415926 / 180 ) )
注意：这里+1是为了避免对0取对数；角度转换为弧度是因为Tan函数使用弧度制
3.2 案例二：多因子叠加分析
假设我们需要评估某区域的地质灾害风险，考虑坡度、坡向和高程三个因子：
- 对各因子进行归一化处理（0-1范围）：
code复制// 坡度归一化
([Slope] - 0) / (90 - 0)
// 高程归一化
([DEM] - 500) / (2000 - 500) // 假设研究区高程范围500-2000米
- 坡向重分类（假设北坡最危险）：
code复制Con( [Aspect] >= 315 || [Aspect] < 45, 1,
Con( [Aspect] >= 45 && [Aspect] < 135, 0.5,
Con( [Aspect] >= 135 && [Aspect] < 225, 0.2, 0.8 ) ) )
- 加权叠加计算风险指数：
code复制0.5 * [Slope_Normalized] + 0.3 * [Aspect_Reclass] + 0.2 * [Elevation_Normalized]
4. 栅格计算的高级技巧与问题排查
4.1 处理NoData值的常见策略
栅格计算中NoData值的处理经常导致意外结果，以下是几种应对方法：
- 使用IsNull函数检测并替换NoData：
code复制Con(IsNull([input]), 0, [input])
-
设置计算环境中的NoData处理选项：
- 在"环境设置"→"栅格分析"中
- 可以设置"忽略NoData"或指定替代值
-
使用SetNull函数有条件地设置NoData：
code复制SetNull([slope] > 30, [slope]) // 坡度大于30度的区域设为NoData
4.2 性能优化技巧
处理大范围高分辨率栅格时，计算可能非常耗时。以下方法可以提升效率：
-
设置合适的处理范围（Processing Extent）和像元大小：
- 在环境设置中限定处理范围
- 统一所有输入栅格的像元大小
-
使用临时内存栅格进行中间计算：
- 将中间结果保存为"内存栅格"（如"in_memory/temp"）
- 减少磁盘I/O开销
-
分块处理超大栅格：
- 使用"栅格分块"工具（Raster Tile）分割输入数据
- 分别计算后使用"栅格拼接"（Mosaic）合并结果
4.3 常见错误与解决方案
-
错误：表达式无效
- 检查括号是否匹配
- 确认所有栅格名称拼写正确（包括大小写）
- 确保运算符两侧有空格（如"a + b"而非"a+b"）
-
错误：像元大小或范围不匹配
- 在环境设置中统一像元大小和范围
- 使用"重采样"工具调整栅格分辨率
-
错误：输出结果全为NoData
- 检查输入栅格是否有有效值
- 确认表达式逻辑正确（如条件语句的阈值设置）
5. 栅格计算与其他工具的协同应用
5.1 与ModelBuilder集成
栅格计算器可以无缝集成到ModelBuilder中，实现自动化处理：
- 在ModelBuilder中添加"栅格计算器"工具
- 将表达式中的栅格参数设置为模型变量
- 使用"计算值"工具动态生成复杂表达式
技巧：在模型中使用"%变量名%"引用变量值，实现动态表达式构建
5.2 与Python脚本结合
通过arcpy模块可以在Python脚本中调用栅格计算：
python复制import arcpy
from arcpy.sa import *
# 设置工作环境
arcpy.env.workspace = "C:/data"
arcpy.CheckOutExtension("3D")
# 执行栅格计算
out_raster = Raster("dem") * 3.28 # 米转英尺
out_raster.save("dem_feet.tif")
# 更复杂的表达式
out_twi = Ln((FlowAccumulation(FlowDirection("dem")) + 1) /
Tan(Slope("dem") * 3.1415926 / 180))
out_twi.save("twi.tif")
5.3 与空间统计工具链结合
栅格计算结果常作为空间统计的输入，例如：
- 先用栅格计算生成地形指标
- 使用"栅格转点"工具转换为样本点
- 应用地统计工具（如克里金插值）进行分析
- 结果再通过栅格计算与其他因子叠加
这种工作流在生态环境评估、地质灾害预警等应用中非常常见。