保姆级教程：用Arcmap 10.0水文分析工具，从DEM到流域边界一步不落

从DEM到流域边界：Arcmap 10.0水文分析全流程实战指南

当你第一次拿到数字高程模型（DEM）数据时，可能会被那些密密麻麻的等高线和复杂的水文分析工具搞得一头雾水。别担心，这篇文章就是为你准备的——无论你是地理信息系统（GIS）专业的学生，还是刚接触水文分析的工程师，甚至是自学GIS的爱好者，都能通过这篇教程一步步完成从DEM数据到闭合流域边界的完整流程。我们将以Arcmap 10.0为操作平台，用最通俗的语言解释每个步骤背后的原理，同时提供那些只有老手才知道的实用技巧和避坑指南。

1. 准备工作与环境配置

在开始正式的水文分析之前，我们需要确保Arcmap 10.0的环境配置正确，并准备好所需的DEM数据。打开Arcmap后，建议先创建一个新的空白地图文档，并设置好工作空间路径。工作空间路径最好选择英文目录，避免中文路径可能导致的潜在问题。

推荐的基础配置步骤：

• 在"地理处理"菜单中勾选"空间分析扩展模块"
• 在"自定义"→"扩展模块"中激活"Spatial Analyst"和"3D Analyst"
• 设置地图文档的坐标系与DEM数据保持一致（通常为WGS84或UTM）

提示：DEM数据的质量直接影响最终分析结果。理想情况下，DEM应该没有数据缺失区域，且分辨率适合你的研究尺度（水文分析常用30米分辨率）。

2. DEM预处理：填洼与流向分析

2.1 为什么必须先填洼？

原始DEM中常存在"洼地"——这些是周围高程高于中心点的区域，会导致水流方向计算错误。填洼（Fill）操作能够消除这些人为误差，确保水流能够连续流出流域。

在Arcmap中执行填洼：

1. 打开"Spatial Analyst Tools"→"水文分析"→"填洼"
2. 输入DEM数据
3. 设置Z限制（通常从默认值1.0开始尝试）
4. 指定输出位置和名称

# 填洼操作的ArcPy代码示例 import arcpy from arcpy.sa import * filled_dem = Fill("input_dem.tif") filled_dem.save("filled_dem.tif")

2.2 流向分析与流量累积

填洼后的DEM可以计算流向（Flow Direction）和流量累积（Flow Accumulation）：

流向分析常见问题解决：

• 如果结果出现明显的直线条带，可能是原始DEM存在质量问题
• 流量累积值异常大时，检查DEM是否进行了正确的填洼处理

3. 流域边界提取的关键步骤

3.1 确定倾泻点位置

倾泻点（Pour Point）是流域的出水口位置，可以通过以下几种方式确定：

1. 已知的测量站点坐标
2. 流量累积值最大的河流交汇处
3. 通过"栅格计算器"筛选高累积值区域

# 使用Con函数筛选高流量累积区域 high_flow = Con(FlowAccumulation("flow_dir") > 1000, 1) high_flow.save("high_flow_areas.tif")

3.2 捕捉倾泻点与流域划分

手动选择的倾泻点可能不完全准确，需要使用"捕捉倾泻点"工具进行微调：

1. 设置合理的捕捉距离（通常为DEM分辨率的2-5倍）
2. 运行"分水岭"工具生成初步流域边界
3. 使用"栅格转多边形"工具将结果转为矢量格式

注意：捕捉距离过大会导致流域范围失真，过小则可能无法正确捕捉到河流网络。

4. 结果验证与常见问题排查

4.1 流域边界闭合性检查

一个合格的流域边界应该是完全闭合的多边形。在Arcmap中可以通过以下步骤验证：

1. 打开属性表，检查是否有多个要素
2. 使用"要素转线"→"要素转面"确保闭合性
3. 检查拓扑错误（重叠、缝隙等）

4.2 常见错误与解决方案

在实际项目中，我遇到过DEM分辨率与捕捉距离不匹配导致流域边界断裂的情况。经过多次测试发现，将捕捉距离设为DEM分辨率的3倍，同时确保所有中间数据使用相同的坐标系，能够有效避免这类问题。