ArcGIS水文分析实战教程（7）细说流域提取

ArcGIS水文分析实战教程（7）细说流域提取

本章导读：流域是水文领域上是频繁出现的词汇，贯穿整个水文专业。因为流域涉及到的空间范围很是的广，通常不能经过测量的方式去精确划定。并且流域这一词在水利行业上属于管理上的词汇，用做划分水利单位的职能范围。而水文学上的流域就相对来讲更倾向于学术，其指向在地理学上称分水岭，在空间物理上叫集水区。在GIS中对于流域的划分，都是基于DEM数据计算，这样可以比较精准（相对于人为的划定区域）的肯定流域范围，为水文结果提供合理的论证。本章会从一些实际的问题出发去说明流域提取要注意的一些问题。BY 李远祥

流域提取

流域提取，或者叫集水区提取，是根据河流、流向和出水口共同肯定其空间范围。从水文学和地理学上，其区域必须是与河流有对应关系。因此，在流域提取以前，必须指定河流。河流数据能够从DEM数据提取，也能够从已有的矢量河流中转化过来。

若是河流数据是从DEM中提取的，那么对河流的划定就必定要仔细，必需要对流量栅格进行正确的估算，以确保河流的尺度，这样经过符合研究口径的河流才能正确的提取出所对应的流域。流量计算和河流提取，能够分别参考《ArcGIS水文分析实战教程（5）细说流向与流量》和《ArcGIS水文分析实战教程（6）河流提取与河网分级》的内容。

上述的这种状况通常都是比较理想的。所谓的理想，就是河流和流域数据都是出自于同一份DEM，原则上重合得比较好。但若是研究的对象是已经划定的矢量河流，从矢量河流转换成栅格河流去参与计算，两份数据之间会存在差距，每每会致使结果误差很大。

若是是研究小流域，这些小流域通常都是一些不显眼的小河流，甚至不多能在地图上查找出来，那么经过DEM方式直接提取流域，能很好的解决问题。

提取的办法很简单，ArcGIS提供了【分水岭】工具用做划定流域。其帮助说明对于GIS人员和水文人员都很是难理解，甚至不少人都无法搞懂其中的参数。下图是该工具的界面：

【流向栅格数据】比较好理解，但下面的【输入栅格数据或要素倾泻点数据】就有点让人摸不着头脑。

要素倾泻点数据通常就是出水口，在上一章中笔者已经讨论过能够经过河流分级以后提取了和流线，并利用【要素折点转点】工具来提取其终点做为出水口。但若是这种方式提取的出水口必需要使用【捕捉倾泻点】工具来从新校订。

而现实中大多数水文研究者对于这个参数，使用的就是其水文站点数据，由于他们研究的对象通常都是针对水文站的数据肯定其上游流域，以寻找上游流域与水文站数据之间的关联关系。这种模式一样是须要经过【捕捉倾泻点】工具来校订。由于数据不是来自于DEM自己，须要如今DEM中找到最近的流量值最大的点才做为计算用的倾泻点。

捕捉倾泻点

正是其余参与计算的数据，不论是河流仍是倾泻点（出水口、水文站） ，也不论是栅格仍是矢量的，只要他们不是从同一份DEM中分析出来的，都须要跟DEM进行校订，这个校订的过程就是【捕捉倾泻点】。其目的是要找到这个点附近的区域流量栅格中的最大值，确保这个流量值是由上游区域贡献的流量。不然，实际输出的倾泻点数据有可能与之叠加的流量栅格不是该区域的最大值，致使出现很是细小的流域面。

来看【捕捉倾泻点】工具的一些参数设置，以下图

第一个参数【输入栅格数据或者倾泻点数据】，这里能够是倾泻点的栅格数据和矢量数据。若是是矢量数据，工具会自动转化为栅格数据以后进行计算。站位水文分析角度来讲，这个数据通常就是水文站位置、研究的断面位置。

【蓄积栅格数据】指的就是流量统计的数据。
【捕捉距离】是个可选参数，但这个参数必定要注意，默认值为0，意思是找到附近最大的流量的栅格做为该点对应的出水口位置。但若是输入的点恰好在多条河流的附近，就有可能捕捉到其余河流对应的流量栅格。 因此，必须对自身参与计算的数据进行摸底，特别是多个倾泻点距离较近的时候，必须设置合理的搜索距离。

只要数据没有原则上的错误，捕捉到的结果一定在原来倾泻点的附近，以下图，捕捉到的像元点都在倾泻点的边上，是一个栅格数据。

分水岭

【分水岭】工具的使用前面介绍过，有了上一步的捕捉到的倾泻点数据，【分水岭】工具就很是容易使用了。
先来看第一种状况，倾泻点是外部数据，与原始DEM数据没有关系的状况。利用捕捉倾泻点工具先生成新的倾泻点数据，而后才利用分水岭工具进行分析。

结果以下图，流域范围被提取出来。

只要再使用【数据管理】工具箱里面的【栅格转面】就能够将该流域的矢量面提取出来。

这种应用场景通常是水文分析人员为了研究该断面数据与上游流域的关系，也就是说经过检测该断面的水文数据及时发现上游流域的一些变化，反之亦然。而实际上流域提取大多数都是由于这个缘由而作的。

若是这种状况没有作【捕捉倾泻点】，则会出现什么状况？因为提供的出水口位置与DEM数据的偏差，这个点有可能偏在一些很是小的河流的流量最大值上。假如斯特拉勒分级方法的话，那这个点就有可能落在1级河流的出水口位置。以下图，之因此出现不少很是小的面，缘由就在于此。

若是不懂原理，那么永远都不可能知道这里出错的缘由。这也就是笔者为何要花大量的篇幅来讲明D8算法和流量统计等原理。

另外一种状况是比较理想的状况，没有采集回来的倾泻点数据，全部的数据都从DEM数据中提取。
这种状况的研究对象就不是上游流域了，而是整个河流的流域。通常来讲大江大河不会用这种方式去作研究，毕竟国家大江大河的监控是很是专一的，按照目前的状况，大江大河很难出现特大洪涝灾害。洪涝灾害通常多发在中小河流。而对中小河流的流域进行提取，有助于分析其集水区对河流的影响。因此，利用DEM来提取流域（非指定外部倾泻点），通常用做小流域的提取。

其作法仍是要回到【分水岭】工具中。【分水岭】工具的第一个参数，在第一种状况中咱们已经使用了捕捉倾泻点后的栅格数据。使用栅格数据做为输入，除了能够是倾泻点，还能够是执行过河流连接的栅格数据。详细作法能够参考《ArcGIS水文分析实战教程（6）河流提取与河网分级》

只要明白原理，就能够很是清楚了。在水文学上，河流连接以后会生成河流链，其中也包括了河源、节点和出口。虽然ArcGIS工具的河流连接数据是只是包含了节点，但因为在连接以前已经作了河网分级，因此会带有其出口信息（也就是说这个工具自动识别了倾斜点数据）。其设置以下

最终生成的结果与河流的数量是一致，以下图所示

一样，若是须要转为矢量面数据，还须要使用【栅格转面】工具进行转换。以下图

总结

流域提取是一门学问，并非GIS课堂上所说的只是经过DEM提取出对应的流域面就完事。流域做为水文学上基础的研究对象，其意义在于水文的相关性，流域面做为空间范围，能够圈定该空间位置的一切地理要素，时空要素与水系之间的关系。因此，必需要了解其使用的场景和实际的做用，才能够提取出这些基础的数据。

若是读者对水文分析感兴趣的话，能够持续关注CSDN的GIS制图乐园，以及微信公众号【GIS制图乐园】。BY 李远祥