ArcGIS Pro中基于NHD Flow lines计算野火到湖泊沿水系最短距离

针对线类型水文网络的沿程距离计算，不推荐使用栅格类最小成本路径工具——这类工具基于栅格像元做成本累加，线转栅格时会产生锯齿误差，距离计算偏差通常在10%以上，还容易出现跨支流、跳过拓扑断点的错误路径，完全匹配需求的方案是用ArcGIS Pro的网络分析模块做基于河网拓扑的最短路径计算，具体流程如下：

1. 前置数据预处理 #

• 河网数据裁剪：以每个目标湖泊为中心生成100英里范围的缓冲区，用缓冲区裁剪对应范围的NHD Flow lines数据，剔除范围外的冗余河网，大幅降低后续计算量
• 河网属性补全：给裁剪后的河网线要素新增长度字段，用计算几何工具计算每段河网的英里长度，作为后续路径计算的通行成本
• 源/目标点提取：
  • 野火侧：提取每个野火斑块与河网的所有交点/近邻接入点（如果野火不直接贴河网，可设置500米搜索半径匹配最近河网位置），作为路径计算的源点
  • 湖泊侧：提取每个湖泊岸线与河网的所有交点（即所有入水口、出水口位置），作为路径计算的目标点集，确保覆盖所有连通支流的接入位置

2. 核心计算步骤 #

• 构建河网网络数据集：基于预处理后的NHD河网线创建网络数据集，连通性规则设置为端点连通（匹配NHD原生拓扑规则，避免线段交叉处产生伪连通），将之前计算的河网长度字段设为唯一成本属性，不需要添加单行限制、转弯限制等道路网络规则，河网默认支持双向追溯
• 运行最近设施点分析：在网络分析工具集中新建「最近设施点」分析图层：
  • 将所有湖泊的入/出水口点加载为「设施点」，按湖泊ID分组
  • 将所有野火的河网接入点加载为「事件点」，按野火斑块ID分组
  • 分析参数设置：搜索截止值设为100英里，取消通行方向限制，设置为遍历所有连通河网分支后返回最短路径
  • 运行求解后，工具会自动遍历所有连通支流的可行路径，自动筛除不连通的分支，返回每个野火接入点到对应湖泊所有入/出水口的沿程路径和长度
• 结果聚合：对同一个野火斑块、同一个湖泊对应的所有路径长度取最小值，即为该野火到对应湖泊的最短沿水系距离。

3. 注意事项 #

• NHD部分区域的河网可能存在拓扑断点，如果出现无有效路径的结果，可先对裁剪后的河网运行检查几何、修复几何工具，再重新构建网络数据集
• 可抽取单个野火样例（如CAPCO野火）做目视校验，将导出的最短路径要素叠加原始河网、野火、湖泊图层，确认路径没有偏离实际连通支流
• 80个湖泊的计算量不大，按上述流程裁剪后单批次计算耗时通常在1小时以内，不需要拆分任务

内容的提问来源于stack exchange，提问作者brittbarreto