ArcGIS水文分析提取河网：步骤和填洼失败排查

ArcGIS水文分析提取河网：步骤和填洼失败排查

用 DEM 做流域边界、汇流路径和沟谷线识别时，ArcGIS水文分析最常见的任务就是提取河网。新手通常能找到 Hydrology 工具箱，却容易卡在两个地方：一是提取出来的河网太密、太碎或断断续续；二是前置的 Fill 填洼工具运行失败，后面的流向和汇流累积都无法继续。

本文按真实项目流程讲清楚ArcGIS水文分析提取河网的完整步骤，并重点排查ArcGIS水文分析填洼失败的原因。你可以把它当作一份从 DEM 检查、填洼、流向、汇流累积到矢量河网输出的操作清单。

问题背景：为什么河网提取不是直接把 DEM 转成线

DEM 只是连续高程栅格，每个像元记录地表高度。河网提取要做的是从高程变化中推断水从哪里流向哪里、哪些位置汇集了足够大的上游面积，再把这些高汇流位置识别为沟道或河网。因此，ArcGIS 的水文处理不是一个单独按钮，而是一条依赖关系很强的栅格处理链。

在这条链里，Fill 填洼通常排在前面。DEM 中可能存在真实洼地，也可能存在由噪声、数据拼接、采样误差造成的微小坑洼。如果不处理，水流会停在这些像元里，导致 Flow Direction 结果异常，Flow Accumulation 不能形成连续汇流路径，最后的河网就会断裂或偏移。

但填洼也不是越彻底越好。山地、水库、矿坑、湖泊、岩溶洼地等场景中，部分洼地可能是真实地形。正式项目里要结合 DEM 来源、研究区水文背景和成果用途判断，而不是把 Fill 当成无条件修复工具。

核心原理：流向、汇流累积和阈值决定河网形态

提取河网的底层逻辑可以拆成三步。第一步是让每个像元有明确的下游方向；第二步是计算每个像元上游有多少像元汇入；第三步是设定一个阈值，只保留汇流累积量足够大的像元作为河网。

Flow Direction 常见输出是方向编码栅格，用来表示每个像元向邻近哪个方向流动。Flow Accumulation 则在流向栅格基础上计算上游贡献面积，数值越大，说明上游汇入的像元越多，越可能形成沟道或河流。

阈值是河网提取里最关键也最容易被忽略的参数。阈值过小，山坡上的临时汇流路径都会被提取出来，河网会像毛细血管一样过密；阈值过大，小支流会被删除，主河道也可能断裂。一般建议把阈值换算成最小汇水面积，而不是凭感觉输入一个整数。

河网提取的结果不是软件“自动识别真实河流”，而是 DEM、填洼方式、流向算法、汇流阈值和后处理规则共同生成的水文推断结果。

ArcGIS水文分析步骤：从 DEM 到矢量河网

下面这套ArcGIS水文分析步骤适合 ArcGIS Pro，也适合按相同逻辑在 ArcMap 中操作。工具名称可能因界面语言略有不同，但处理顺序不要随意打乱。

1. 准备 DEM。优先使用覆盖完整研究区的单波段高程栅格，确认没有明显条带、空洞、拼接错位和异常高程值。DEM 最好放在 File Geodatabase 或路径简短的英文目录中。
2. 统一坐标系和单位。水文分析可以在栅格上运行，但面积阈值、像元大小和距离解释更依赖投影坐标系。正式分析建议使用适合研究区的米制投影，而不是直接在经纬度坐标下解释汇水面积。
3. 设置环境参数。在地理处理环境中明确设置 Processing Extent、Cell Size、Snap Raster、Mask 和 Output Coordinate System。若有研究区边界，建议先用一致的边界裁剪 DEM，再继续后续处理。
4. 运行 Fill 填洼。在 Spatial Analyst 的 Hydrology 工具中运行 Fill，输入 DEM，输出填洼后的 DEM。初学阶段可以先不设置 Z limit，观察结果是否正常；正式项目要根据地形和业务判断是否限制填洼深度。
5. 运行 Flow Direction。用填洼后的 DEM 生成流向栅格。后续 Flow Accumulation、Stream Link、Watershed 等工具都依赖这个结果。
6. 运行 Flow Accumulation。输入流向栅格，生成汇流累积栅格。这个栅格通常显示为大范围低值和少量高值河谷，需要用合适的拉伸或分类方式查看。
7. 确定河网阈值。根据像元大小和希望识别的最小汇水面积计算阈值。例如 30 米 DEM，一个像元面积是 900 平方米；若希望最小汇水面积约 0.9 平方公里，阈值可从 900000 / 900 = 1000 附近试起。
8. 用条件表达式提取河网栅格。可用 Raster Calculator 或 Con 工具把高于阈值的汇流累积像元转成河网栅格，例如 Con("flow_acc" >= 1000, 1)。阈值需要多轮对比，不建议一次定死。
9. 生成河段和等级。可选运行 Stream Link 划分河段，再用 Stream Order 生成 Strahler 或 Shreve 河流等级。需要做流域结构分析时，这一步很有用。
10. 转成矢量线。用 Stream To Feature 把河网栅格转为线要素，输入河网栅格和流向栅格。输出后检查线方向、断点、悬挂线和与已知河流的吻合度。

这套流程的重点不是记住工具名字，而是理解每一步输出都服务于下一步。只要 Fill、Flow Direction 或 Flow Accumulation 中任何一步异常，后面的河网线都会放大这个错误。

河网提取阈值应该怎么选

河网结果质量很大程度取决于汇流累积阈值。很多教程会直接给一个 500、1000 或 5000 的数字，但这个数字只有和 DEM 分辨率、研究区地形、目标比例尺结合起来才有意义。

更稳妥的做法是先确定“最小汇水面积”，再换算为像元数量。公式很简单：阈值等于最小汇水面积除以单个像元面积。若 DEM 像元大小是 10 米，单个像元面积为 100 平方米；若希望提取 0.1 平方公里以上的沟道，阈值约为 100000 / 100 = 1000。

阈值调试时建议同时叠加影像、已有水系、等高线或实测河道。山地丘陵区可能需要较小阈值才能保留支沟；平原区 DEM 起伏弱，阈值过低会产生很多不符合地貌逻辑的线。最终阈值应服务于你的制图比例尺和分析目的，而不是追求线条越多越好。

ArcGIS水文分析填洼失败：先查输入、环境和许可

遇到ArcGIS水文分析填洼失败，不要马上改输出文件名反复运行。Fill 对输入栅格、许可、路径、环境设置和 NoData 边界都比较敏感，排查应从基础条件开始。

• Spatial Analyst 扩展未启用。Fill、Flow Direction、Flow Accumulation 等水文工具通常依赖 Spatial Analyst。若工具不可用或运行报许可错误，先检查扩展许可是否可用并已启用。
• 输入 DEM 不是有效单波段高程栅格。多波段影像、分类栅格、空栅格或错误单位的高程数据都不适合作为 Fill 输入。先查看栅格属性、统计值、像元大小和 NoData 值。
• NoData 边界或内部空洞太多。DEM 内部有大块 NoData，会让流向链路断开。先检查是否需要重新裁剪、镶嵌、填补空洞或换数据源。
• 环境 Mask 或 Extent 继承错误。如果之前做过裁剪、重采样或其他栅格分析，ArcGIS 可能继承了旧的 Mask、Extent 或 Snap Raster，导致 Fill 输出被意外限制或为空。
• 路径、权限或文件名问题。超长路径、网络盘不稳定、目录无写入权限、输出名含特殊字符，都可能造成工具失败。建议先输出到本地 File Geodatabase 或简短英文目录。
• 栅格过大导致资源不足。高分辨率 DEM 覆盖大范围时，Fill 需要较多内存和临时空间。可以先裁剪研究区、分块处理，或确认临时目录所在磁盘空间充足。
• Z limit 设置不合理。Z limit 过小可能不能处理关键洼地，过大则可能改变真实地形。若不确定，先不设或用小范围样区比较不同参数。

如果工具报错信息不清楚，建议新建一个空白工程，只加载一小块 DEM，清空环境参数后试运行 Fill。小样区能成功，完整数据失败，通常说明问题在数据范围、资源、路径或环境设置；小样区也失败，则优先怀疑输入 DEM、许可或栅格格式。

填洼失败后的具体修复流程

填洼失败常常不是一个原因造成的。下面按排查成本从低到高给出修复流程，建议逐项记录，避免反复试错后忘记改过哪些参数。

1. 复制 DEM 到本地工作库。把输入 DEM 放入本地 File Geodatabase 或短路径文件夹，输出也放在同一工作目录，避免网络路径和权限问题。
2. 查看 DEM 属性。检查像元大小、行列数、坐标系、高程范围、NoData 值和像素类型。若最小值最大值明显异常，先处理高程异常。
3. 清空地理处理环境。在 Environments 中检查 Extent、Mask、Cell Size、Snap Raster、Output Coordinate System。不能确认时，新建工程重新跑一次。
4. 裁剪一个小样区测试。选择包含山脊、沟谷和少量边界的区域运行 Fill。若样区正常，再扩大到完整研究区。
5. 处理 NoData 空洞。若 DEM 内部有孔洞，应先用合理方法补齐或重新获取 DEM。不要把 NoData 简单改成 0，因为 0 可能表示真实海拔，会制造人工洼地。
6. 检查镶嵌边界。多幅 DEM 拼接后可能出现接边台阶或条带，Fill 会把这些异常当作地形处理。必要时先重采样、镶嵌到统一网格并做边界质检。
7. 控制处理范围。大范围高分辨率数据可以按流域或分区先处理，再检查边界衔接。分块时要保留足够缓冲区，否则边界附近流向会受影响。
8. 比较不同 Z limit。对存在真实洼地的地区，不要盲目全量填平。用样区比较不设 Z limit、小 Z limit 和较大 Z limit 的差异，再决定正式参数。

修复完成后，不要只看 Fill 是否成功。还要把原 DEM 和填洼 DEM 做差值检查，查看哪些地方被改变、改变幅度是否合理。若大片湖泊、矿坑或城市低洼区被过度填平，后续河网可能会失真。

常见坑点：河网能出图不等于水文正确

很多水文处理问题表面上是工具失败，实际是结果看起来成功但不可靠。下面这些坑点在课程作业和项目制图中都很常见。

• 用经纬度 DEM 直接解释面积阈值。经纬度坐标的单位是度，不适合直接把像元大小当作米来计算汇水面积。正式阈值换算应基于投影后的像元面积。
• 只调阈值，不检查 DEM 质量。DEM 有条带、空洞、异常高程或道路桥涵阻断时，调阈值只能改变线的多少，不能真正修复流向错误。
• 把所有洼地都当成错误。自然湖泊、封闭盆地和工程洼地可能是真实地貌。若研究目标包含这些水文单元，填洼策略要谨慎。
• 在边界附近直接判断河网断裂。研究区裁剪边界会截断上游汇流面积，边界附近的河网通常不适合直接评价完整性。必要时应使用比研究区更大的 DEM 缓冲范围。
• Stream To Feature 后不做线要素检查。矢量化后仍要检查线方向、重复线、短碎线、悬挂线和属性字段。需要制图交付时，还可能要做平滑、裁剪和拓扑整理。
• 把提取河网当作真实水系替代品。DEM 提取河网是地形推断结果，和实测河道、影像解译水系、地名水系数据之间可能存在差异。正式成果应说明数据来源和提取方法。

工具和方法对比：图形界面、ModelBuilder 和 ArcPy 怎么选

同样是从 DEM 提取河网，可以用 ArcGIS Pro 图形界面、ModelBuilder 或 ArcPy。选择哪一种，取决于你是学习单次流程、做批量生产，还是要长期复现实验参数。

初学者建议先用图形界面跑通完整流程，理解每一步输出的含义；当阈值、环境和质检规则稳定后，再用 ModelBuilder 或 ArcPy 固化流程。这样既能降低出错成本，也方便以后复现同一套河网提取参数。

实用检查清单：提交河网成果前看这些

完成河网提取后，建议按下面清单做一次质检。它能帮助你判断结果是否适合进入流域划分、坡面分析、水土保持评价或制图表达。

FAQ：水文分析提取河网和填洼失败

ArcGIS水文分析提取河网的基本步骤是什么？

基本流程是准备 DEM、设置环境参数、运行 Fill 填洼、生成 Flow Direction、计算 Flow Accumulation、用阈值提取河网栅格，最后用 Stream To Feature 转成矢量线。规范流程还应包含 DEM 质量检查、阈值换算和结果质检。

ArcGIS水文分析填洼失败最常见原因是什么？

常见原因包括 Spatial Analyst 扩展未启用、输入 DEM 无效、DEM 内部 NoData 过多、环境 Mask 或 Extent 设置错误、输出路径没有权限、栅格范围太大导致资源不足，以及 Z limit 参数不合理。建议先用小样区清空环境测试，再回到完整数据处理。

提取河网时阈值越小越好吗？

不是。阈值越小，提取出的线越多，但很多可能只是坡面汇流痕迹或 DEM 噪声，并不适合作为河网。设置阈值时应把它换算成最小汇水面积，再结合影像、已有水系和目标比例尺调整。

Fill 填洼成功后，为什么河网还是断的？

填洼只解决 DEM 中部分局部洼地问题。河网断裂还可能来自 DEM 内部 NoData、研究区裁剪范围太小、流向异常、汇流阈值过大、平原区高程起伏不足，或道路桥涵等人工地物造成的地形阻断。需要继续检查 Flow Direction 和 Flow Accumulation，而不是只重复 Fill。

ArcGIS水文分析步骤中必须先裁剪 DEM 吗？

不一定，但建议使用比研究区略大的 DEM 范围来处理。直接按行政边界裁剪可能切断上游汇流，导致边界附近河网不完整。比较稳妥的做法是先用缓冲后的范围裁剪 DEM，完成水文分析后，再按正式研究区裁剪成果。

提取出的河网和已有水系不一致怎么办？

先确认 DEM 分辨率、投影、填洼方式和阈值是否合理，再判断已有水系的来源和比例尺。DEM 河网反映的是地形推断，已有水系可能来自测绘、影像解译或人工修编。两者不一致时，要结合影像、等高线和实地资料判断，不要简单认为某一方一定正确。

填洼失败可以跳过 Fill 吗？

不建议直接跳过。未填洼 DEM 可能导致流向停在局部洼地，汇流累积无法形成连续下游路径。如果确认洼地是真实水文单元，应选择合适的预处理策略或限制填洼深度，而不是完全忽略填洼失败。

总结

从 DEM 提取河网的可靠流程，是从数据质量检查开始，依次完成填洼、流向、汇流累积、阈值提取、河段分级和矢量化。真正决定结果质量的不是某一个工具按钮，而是 DEM 来源、环境参数、阈值换算和后续质检是否清楚。

遇到填洼失败时，先查许可、输入 DEM、NoData、路径权限、环境 Mask、处理范围和资源占用；遇到河网过密、断裂或偏移时，再回到阈值、流向和 DEM 质量中排查。按这套流程处理，提取出的河网才更容易解释、复现和用于后续空间分析。