1. 为什么需要掌握栅格表面分析？

作为一名GIS从业者，我经常遇到这样的场景：客户拿着一堆高程数据问我"这片区域哪里最容易积水？"、"这个山坡的日照时间怎么计算？"、"我们想修的这条路坡度会不会太大？"这些问题看似简单，但如果没有掌握栅格表面分析的核心技能，就只能对着数据干瞪眼。

ArcToolbox中的3D Analyst工具集正是为解决这类问题而生。栅格表面分析作为其核心功能之一，能够将平面的数字高程模型（DEM）转化为具有三维特性的分析对象。不同于简单的二维栅格操作，表面分析考虑的是连续变化的高程场，可以揭示地形中隐藏的空间关系。

提示：很多人误以为3D Analyst只是用来做漂亮的三维可视化，实际上它的分析功能才是真正的价值所在。

2. 栅格表面分析的核心工具解析

2.1 坡度与坡向计算

坡度（Slope）工具是我日常使用频率最高的表面分析工具之一。它的计算原理是基于每个栅格像元与其相邻8个像元的高程差，通过二阶差分算法得出坡度值。在ArcGIS中，坡度结果可以按度数或百分比输出：

• 度数表示法：0°表示完全平坦，90°表示垂直悬崖
• 百分比表示法：100%表示45°斜坡

坡向（Aspect）工具则告诉我们地形朝向哪个方向。它的输出是0-360度的罗盘方位角，其中：

• 0度表示正北
• 90度表示正东
• 180度表示正南
• 270度表示正西
• -1表示平坦区域

在实际项目中，我经常将坡度和坡向结合使用。比如在太阳能板选址时，需要找南向（坡向约180度）且坡度在15-40度之间的区域。

2.2 山体阴影与光照模拟

Hillshade工具可以生成极具视觉冲击力的地形图，但它的价值远不止美观。通过调整太阳方位角（azimuth）和高度角（altitude）参数，我们可以：

1. 模拟不同季节、时间的日照情况
2. 识别地形阴影区域（可能影响植被生长）
3. 增强地形特征的可视化效果

我常用的参数组合是：

• 方位角315度（西北方向）
• 高度角45度
• Z因子1（除非使用非标准高程单位）

注意：山体阴影只是光照模型，不能直接用于定量分析。如需精确计算日照时长，需要使用Solar Radiation工具集。

2.3 等高线生成与曲率计算

Contour工具看似简单，但有几个关键参数常被忽略：

• 等高距（Contour interval）：决定生成等高线的密度
• 起始高程（Base contour）：避免等高线与零值线混淆
• Z因子：当Z值与XY单位不同时需要调整

曲率（Curvature）分析则更为专业，它包括：

• 剖面曲率（影响水流加速度）
• 平面曲率（影响水流汇聚）
• 总曲率（综合指标）

在土壤侵蚀研究中，高曲率区域往往是侵蚀风险区。我曾用曲率分析成功预测了一处工地开挖后的水土流失热点。

3. 高级表面分析技巧

3.1 视线分析与视域分析

Viewshed工具可以回答"从这里能看到什么"的问题。在通信基站选址时，我通常会：

1. 设置观察点高度（如基站天线高度）
2. 考虑地球曲率影响（对长距离分析很重要）
3. 设置最大可视距离（根据实际需求）

Observer Points工具更强大，可以同时计算多个观察点的可视范围。在景区观景台规划中，这种分析能确保每个观景台都有独特视野。

3.2 体积计算与填挖方分析

Cut/Fill工具是工程规划的神器。它通过比较两个表面（通常是设计前后地形）来计算：

• 填方量（需要填充的土方）
• 挖方量（需要挖除的土方）
• 净变化量

我曾用这个工具为一个高尔夫球场项目节省了30%的土方运输成本。关键在于：

1. 精确设置输入基准面
2. 合理定义Z因子
3. 对结果进行分区统计

3.3 表面参数与地形指数

Surface Parameters工具集提供了更专业的指标：

• 粗糙度（Roughness）
• 起伏度（Relief）
• 地形位置指数（TPI）
• 地形湿度指数（TWI）

这些指数在生态研究中非常有用。例如，TPI可以帮助识别山谷和山脊，而TWI可以预测土壤湿度分布。

4. 实战案例：滑坡风险评估

去年我参与了一个山区滑坡风险评估项目，完整的工作流程如下：

1. 数据准备：

   • 获取1米分辨率LiDAR DEM
   • 预处理（填充洼地、去除异常值）

2. 地形参数提取：

   • 坡度（大于30度高风险）
   • 坡向（阳坡更易干燥开裂）
   • 曲率（高曲率区易发生表层滑动）
   • TWI（高湿度区域风险增加）

3. 加权叠加分析：

   # 伪代码示例 risk = 0.4*slope + 0.2*aspect + 0.3*curvature + 0.1*TWI

4. 验证与调整：

   • 与历史滑坡点对比
   • 调整权重系数
   • 设置风险等级阈值

最终成果不仅准确预测了已知风险区，还发现了三处新的潜在风险点。当地政府根据我们的分析调整了防灾预案。

5. 常见问题与性能优化

5.1 大区域处理技巧

处理省级甚至全国范围的DEM数据时，我采用以下策略：

1. 分块处理：

   • 使用Raster Split工具分割数据
   • 并行处理各区块
   • 使用Mosaic合并结果

2. 分辨率选择：

   • 初步分析用中等分辨率（如30米）
   • 重点区域再用高分辨率数据

3. 金字塔构建：

   • 预处理时建立金字塔
   • 选择NEAREST重采样方法保持原始值

5.2 异常值处理心得

DEM数据常有异常值，我的处理步骤是：

1. 使用Raster Calculator识别异常：

   "DEM" < 0 OR "DEM" > 5000

2. 替换方法选择：

   • 邻域均值（适用于小范围异常）
   • 插值（适用于数据缺失）
   • 手动编辑（关键区域）

3. 验证：

   • 检查统计量（均值、标准差）
   • 生成剖面线查看地形连续性

5.3 参数设置陷阱

几个容易出错的参数设置：

1. Z因子：

   • 当XY单位是度（地理坐标系）而Z单位是米时，需要设置适当的Z因子
   • 近似公式：1度≈111km，所以Z因子≈1/111000

2. 输出测量单位：

   • 坡度选择度还是百分比
   • 曲率选择适合的单位

3. 处理范围：

   • 确保所有输入数据范围一致
   • 使用Snap Raster对齐像元

6. 与其他工具的协同应用

6.1 与Spatial Analyst结合

栅格表面分析常需要与Spatial Analyst工具集配合：

1. 重分类（Reclassify）：

   • 将连续坡度分为陡、中、缓三级
   • 为不同等级赋予权重值

2. 栅格计算器（Raster Calculator）：

   • 组合多个地形指数
   • 创建自定义分析模型

3. 区域统计（Zonal Statistics）：

   • 计算各流域的平均坡度
   • 统计不同海拔带的面积

6.2 与3D场景集成

分析结果可以导入ArcScene或ArcGIS Pro的3D场景：

1. 设置基准高度：

   • 使二维栅格浮在三维空间
   • 实现真实地形叠加

2. 垂直夸大：

   • 突出细微地形特征
   • 典型值2-3倍（过大导致失真）

3. 动画制作：

   • 飞行动画展示地形
   • 日照变化模拟

6.3 与Python自动化

对于重复性工作，我编写Python脚本自动处理：

import arcpy from arcpy.sa import * # 设置工作环境 arcpy.env.workspace = "C:/data/terrain" arcpy.env.extent = "study_area.shp" # 批量计算坡度坡向 dem = "elevation.tif" slope = Slope(dem, "DEGREE") aspect = Aspect(dem) # 保存结果 slope.save("slope_deg.tif") aspect.save("aspect.tif")

这个脚本可以扩展为处理整个项目文件夹中的所有DEM数据。