基于ArcGIS水文流域分析

（1）数据基础：无洼地的DEM

（2）关键步骤：流向分析

（3）计算流水累积量

（4）提取河流网络

（5）提取河流网络矢量数据。

（6）平滑处理河流网络。

（7）进行流域分析。

（1）数据基础：无洼地的DEM

在ArcMap中加载 DEM数据，右击DEM图层，点击缩放至图层，显示全部。

在【ArcToolbox】中，（要打开扩展模块）执行命令[SpatialAnalyst工具]——>[水文分析]——> [填洼]，按下图所示指定各参数，其中Z限制——填充阈值，当设置一个值后，在洼地填充过程中，那些洼地深度大于阈值的地方将作为真实地形保留，不予填充；系统默认情况是不设阈值，也就是所有的洼地区域都将被填平。之后点击确定即可。

确定后执行结果得到无洼地的DEM数据[Fill_dem1]

（2）关键步骤：流向分析

在上一步的基础上进行，在【ArcToolbox】中，执行命令[SpatialAnalyst工具]——>[水文分析]——>[流向]，按下图所示指定各参数：

确定后执行完成后得到流向栅格[Flowdir_fill1]

（3）计算流水累积量

在上一步的基础上进行，在【ArcToolbox】中，执行命令[SpatialAnalyst工具]——>[水文分析]——>[流量]，按下图所示指定各参数：

确定后执行完成得到流水累积量栅格[flowacc_flow1]，如图：

（4）提取河流网络

首先，提取河流网络栅格。

在上一步的基础上进行，打开【Arctoolbox】，运行工具[Spatial Analyst 工具]——>[地图代数]——>[栅格计算器]，在[地图代数表达式]中输入公式：Con(Flow Accumulation1>800,1)，（这里的Flow Accumulation1要以上一步得到的文件名为准,注意是Con,不是con,大写第一个字母，不然出错）如图：

[输出栅格]指定为:StreamNet保存路径和文件名任意）

执行后，结果为（关闭除[Streamnet]之外的其它图层）：

结果说明：通过此操作将流水累积量栅格[Flow Accumulation1]中栅格单元值（流水累积量）大于800的栅格赋值为1，从而得到河流网络栅格[StreamNet]

（5）提取河流网络矢量数据。

在上一步的基础上进行，打开【Arctoolbox】工具箱，运行工具[SpatialAnalyst 工具]——>[水文分析]——>[栅格河网矢量化]，按下图设置参数：

后得到河流网络矢量数据，如图

（6）平滑处理河流网络。

打开[编辑器]工具栏，执行工具栏中的命令[编辑器]——>[开始编辑]，确保目标图层为河流网络图层[StreamT_StreamN1],通过打开[StreamT_StreamN1]属性表，并选择属性表的所有行选择图层[StreamT_StreamN1]中的所有要素，也可以通过[要素选择按钮]选择图层中所有要素。执行[编辑器]工具栏中的命令[编辑器]——>[更多的编辑工具]——>[高级编辑]打开工具条：[高级编辑]，点击其上的[平滑]按钮：

在[平滑]处理对话框中输入参数[允许最大偏移]为4，得到平滑后的河流网络矢量图层，执行命令: [编辑器]——>[停止编辑]，保存所做修改。（然后，比较平滑处理后的数据与没有进行处理过的数据）

（7）进行流域分析。

在上一步的基础上进行，打开【Arctoolbox】，运行工具[Spatial Analyst工具]——>[水文分析]——>[盆域分析]，如图：

然后在打开【Arctoolbox】中，打开[空间分析]工具栏，执行命令:[转换工具]——>[由栅格转出]——>[栅格转面]，将流域栅格转换成为矢量图层。如图：