基于DEM的ArcGIS流域边界提取

一、操作目标：以已知目标点位作为流域出口，提取上游集水区范围。

二、操作步骤：

        1.准备DEM数据（流向、流量、水系提取基础）

        2.填洼处理：打开填洼工具（Spatial Analyst 工具/水文分析/填洼），输入dem数据，输出Fill.tif。

        3.流向计算：打开流向工具（Spatial Analyst 工具/水文分析/流向）计算栅格流向，输入Fill.tif，输出FlowDir.tif。

        4.流量计算：打开流量工具（Spatial Analyst 工具/水文分析/流量）计算栅格累积流量，输入FlowDir.tif，输出FlowAcc.tif。

        5.冲沟提取：打开大于等于工具（Spatial Analyst 工具/数学分析/逻辑运算/大于等于）针对上一步计算的流量结果提取，提取大于等于常量a的冲沟栅格，输入FlowAcc.tif，输出Greater_Flow.tif。

        6.重分类：打开重分类工具（Spatial Analyst 工具/重分类/重分类）针对上一步提取的冲沟去掉背景值，将背景值设置为Nodata，提大于等于常量a的冲沟栅格，输入Greater_Flow.tif，输出Reclass_Greater.tif。

        7.河流链接：打开河流链接工具（Spatial Analyst 工具/水文分析/河流链接）对重分类好的栅格水系进行附加流向属性。输入Reclass_Greater.tif，输出Stream_Rec.tif。

        8.栅格河网矢量化：打开栅格河网矢量化工具（Spatial Analyst工具/水文分析/栅格河网矢量化）对计算好的具有流向属性的栅格水系转成矢量线图层。输入Stream_Rec.tif，输出Stream.shp。

        9.定义倾泻点（流域出口点），提取流域范围：在任意工作文件夹中新建倾泻点shp文件，命名为“OutPoint.shp”，此点就是指定的流域出口点，根据上一步提取出来的河道水系线，尽量把点标在水系线上。打开集水区工具（Spatial Analyst工具/水文分析/集水区，部分汉化包翻译为分水岭工具），输入流向FlowDir.tif数据、倾泻点OutPoint.shp，输出watershed.shp。

     以上提取的watershed.shp就是所需要的流域面。

注：针对提取流域范围为什么先提取流域水系的问题：限于DEM的数据精度，可能与实际情况出现偏离，因此先提取水系可以有效地确定倾泻点（流域出口）的位置，尽量将倾泻点设置在提取的水系上，因为将倾泻点设置在水系上不容易报错。

以上也是个人学习实践经验，仅供参考！