在ArcMap中加载 DEM数据，右击DEM图层，点击缩放至图层，显示全部。
在ArcMap中，需要打开【扩展模块】功能，如下图：
弹出的对话框中勾选所有选项，如下图：
点击【关闭】。
然后，在【ArcToolbox】中执行命令[SpatialAnalyst工具]——>[水文分析]——> [填洼]，按下图所示指定各参数，其中Z限制——填充阈值，当设置一个值后，在洼地填充过程中，那些洼地深度大于阈值的地方将作为真实地形保留，不予填充；系统默认情况是不设阈值，也就是所有的洼地区域都将被填平。之后点击确定即可。
(注意：上图中的“Output surface raster”点击的时候，不要用默认目录，自己新建一个目录，然后命名保存的文件,以后每一步的保存文件都保存在这里面)。
确定后执行结果得到无洼地的DEM数据[Fill_dem1]
(2)关键步骤：流向分析
在上一步的基础上进行，在【ArcToolbox】中，执行命令[SpatialAnalyst工具]——>[水文分析]——>[流向]，按下图所示指定各参数：
确定后执行完成后得到流向栅格[Flowdir_fill1]，如下图：
此时的流向就满足D8算法，分别是1、2、4、8、16、32、64、128
(3)计算流水累积量
在上一步的基础上进行，在【ArcToolbox】中，执行命令[SpatialAnalyst工具]——>[水文分析]——>[流量]，按下图所示指定各参数：
确定后执行完成得到流水累积量栅格[flowacc_flow1]，如图：
(4)提取河流网络
首先，提取河流网络栅格。
在上一步的基础上进行，打开【Arctoolbox】，运行工具[Spatial Analyst 工具]——>[地图代数]——>[栅格计算器]，在[地图代数表达式]中输入公式：Con(Flow Accumulation1>800,1)，(这里的Flow Accumulation1要以上一步得到的文件名为准,注意是Con,不是con,大写第一个字母，>符号两边都有空格，尽量不要键盘输入，用鼠标点面板上的公式，不然容易出错,)如图：
注意：设定阈值(800)。不同级别的沟物对应不同的阈值，不同区域相同级别的沟谷对应的阈值也是不同。实际工作中设定阈值应通过不断实验、结合土壤植被等资料辅助检验方法确定。 可设置任意数字。
[输出栅格]指定为:StreamNet保存路径和文件名任意)
执行后，结果为(关闭除[Streamnet]之外的其它图层)：
结果说明：通过此操作将流水累积量栅格[Flow Accumulation1]中栅格单元值(流水累积量)大于800的栅格赋值为1，从而得到河流网络栅格[StreamNet]
然后，提取河流网络矢量数据。
在上一步的基础上进行，打开【Arctoolbox】工具箱，运行工具[SpatialAnalyst 工具]——>[水文分析]——>[栅格河网矢量化]，按下图设置参数：
后得到河流网络矢量数据，如图：
接着,平滑处理河流网络。
打开[编辑器]工具栏，执行工具栏中的命令[编辑器]——>[开始编辑]，确保目标图层为河流网络图层[StreamT_StreamN1],通过打开[StreamT_StreamN1]属性表，并选择属性表的所有行选择图层[StreamT_StreamN1]中的所有要素，也可以通过[要素选择按钮]选择图层中所有要素
执行[编辑器]工具栏中的命令[编辑器]——>[更多的编辑工具]——>[高级编辑]打开工具条：[高级编辑]，点击其上的[平滑]按钮：
在[平滑]处理对话框中输入参数[允许最大偏移]为4，得到平滑后的河流网络矢量图层，执行命令: [编辑器]——>[停止编辑]，保存所做修改。(然后，比较平滑处理后的数据与没有进行处理过的数据)
最后，进行流域分析。
在上一步的基础上进行，打开【Arctoolbox】，运行工具[Spatial Analyst工具]——>[水文分析]——>[盆域分析]，如图：
然后在打开【Arctoolbox】中，打开[空间分析]工具栏，执行命令:[转换工具]——>[由栅格转出]——>[栅格转面]，将流域栅格转换成为矢量图层。如图：
私信【dem】可免费获取arcgis软件、随机地区的12米分辨率dem数据。