我们经常会使用研究区的DEM数据进行水文分析，在进行水文分析之前，需要进行洼地的填充，然后进行流向、流量的分析，最后根据最符合研究区现实情况的阈值，提取更加精确的河网，然后可以将河网矢量化，最后根据实际情况，将不符合现实情况的河网删去。

1、加载DEM数据，这里的DEM数据是NASA_DEM_n22e113，具体的现在可以到earthdata免费下载，首先将研究区的DEM按掩膜提取，以中山市为例

2、然后进行一个填洼的处理，将研究区范围内的洼地填充：

填洼之后我们可以看到，研究区的DEM最高点没有变化，但是最低点变高了，说明研究区内存在一定数量的洼地：

3、然后进行流向分析，流向分析可选择三种方法，流向工具支持三种流向建模算法。分别为 D8、多流向 (MFD) 和 D-Infinity (DINF)：

D8 流向法可对每个像元到其最陡下坡邻域的流向进行建模。
以 D8 流向类型运行的流向工具的输出是值范围介于 1 到 255 之间的整型栅格。从中心出发的各个方向值为：

例如，如果最陡下降方向位于当前处理像元的左侧，则将该处理像元的流向编码将为 16，如果最陡下降方向位于当前处理像元的右侧，则将该处理像元的流向编码将为 1，如果最陡下降方向位于当前处理像元的左上，则将该处理像元的流向编码将为 32，以此类推。
这里我们选择D8法：

然后进行一个流量分析：

然后进行栅格计算器，将流量大于300的提取出来（这里的300是一个阈值，阈值的设置不是固定的，需要根据研究区的现实状况，进行合理的选择，也可以进行多次试验进行合理选择，本例选择300）：

得到河网后放大观察，是树状的河流网络：

最后，可以将栅格的河网进行矢量化，得到矢量的河网图层，当然最后的结果实际情况可以叠加到卫星图上进行比对，然后再针对流量的阈值进行调整，最后使用更加精确的阈值得到河网，当然最后的河网还是会和现实情况有出入，这个时候进行矢量化之后手动删除一些不符合现实情况的河网即可。