基于数字高程模型（DEM）进行水文分析，自动提取河网、划分集水区边界，计算流域面积、河长与坡降指标，支持Horton定律验证与Strahler河网分级。
基于D8算法的数字流域提取与河网分析
💡 提示：累积流量阈值越小，提取的支流越多；阈值越大，仅保留主干河道。Strahler分级法下，同级别的两条河流交汇时级别才加1。
出水口（Pour Point）是流域划分的起始点，系统将从该点向上游回溯确定流域边界。
流域提取的理论基础、算法对比与实际应用
基于DEM的流域提取是数字水文分析的基础。完整的分析流程包括：数据预处理（投影转换、拼接裁剪）→ 填洼处理（Fill Sinks）→ 流向分析（Flow Direction）→ 累积流量计算（Flow Accumulation）→ 河网提取（Stream Network）→ 河网分级（Stream Order）→ 流域划分（Watershed/Basin Delineation）。
在ArcGIS中，这些工具集中位于 Spatial Analyst → Hydrology 工具箱下。在QGIS中，可使用 GRASS → r.watershed 或 SAGA → Catchment Area 模块实现相同功能。
| 算法 | 类型 | 适用场景 | 优缺点 |
|---|---|---|---|
| D8 | 单流向 | 陡坡山地、一般流域 | 简单高效；平坦区有偏差 |
| MFD / D∞ | 多流向 | 缓坡、平原区 | 更符合自然；计算量大 |
| FMM | 单流向 | 高精度DEM | 精度高；实现复杂 |
| Lea | 单流向 | 特殊地形 | 处理洼地能力强 |
Strahler分级法：最上游支流为1级；两条同级别的河流交汇时，级别加1；不同级别交汇时取较大者。该方法强调河网的拓扑结构，广泛应用于地貌学和水文学研究。
Shreve分级法：最上游支流为1级；任何两条河流交汇时，级别相加。该方法更强调累积贡献，常用于水量计算和地貌瞬时单位线（GIUH）的参数推导。
实际应用中，Strahler法更为常见，因为它能更好地反映河网的层级结构特征。
洪水模拟与预报：提取的河网和流域边界是分布式水文模型（如SWAT、HEC-HMS、TOPMODEL）的基础输入数据，用于模拟降雨-径流过程和洪水演进。
土壤侵蚀评估：结合坡度、坡长因子（LS因子），可计算流域内的土壤侵蚀量，为水土保持规划提供依据。
水资源管理：划分集水区边界，计算流域面积、河长、坡降等参数，为水库选址、取水口规划和生态流量评估提供支撑。
地质灾害评估：泥石流、滑坡等灾害的敏感性评价需要精确的地形和水系信息。