【ArcGIS教程】基于DEM数据的洪水淹没分析——提取淹没区并计算淹没区洪水量

---

一、案例背景

南方某地位于河口三角洲，属于低洼易涝区。现需要计算一定洪水量可淹没的区域面积。

二、数据说明：

研究区.shp：研究区范围（面状要素）。
DEM.tif：研究区的数字高程模型。

三、分析要求：

1、提取洪水淹没范围。

自行假定洪水倾泻点，提取洪水淹没区范围。

2、计算洪水淹没量。

假定淹没高度，并根据该地高程计算不同高度下的洪水淹没量。

四、分析步骤：

1、精准填洼:

运行ArcToolbox，打开【空间分析工具】，选择【水文分析】里面的【填洼】工具，调出面板进行参数设置。输入栅格选择原始DEM数据，设定输出路径保存为Fill.tif，输出表面栅格。

2、流向计算：

选择【水文分析】里面的【流向】工具，调出面板进行参数设置。输入表面栅格选择填洼好的数据（Fill.tif），保存为FlowDir.tif。

3、流量统计：

选择【水文分析】里面的【流量统计】工具，调出面板进行参数设置。输入流向栅格选择【FlowDir.tif】数据，输出为FlowAcc.tif保存路径即可。

如下，【流量统计】完成：参照这个结果的图例，可以直接从图中读取到一些信息。这个流量统计栅格亮白色代表的是流量越大，其中最大一点共有 44309个像元向其汇入。按照前面章节获取该原始DEM数据的元数据信息可以知道，该数据为像元大小为12.4912.49，那么流量最大点的汇入栅格面积就是12.4912.49 * 44309 平方米。

4、定义最小地表径流

通过对流量统计结果FlowAcc.tif栅格分类按照自然间断点法分类，定义不同值下的最小地表径流，这里假设分为3，可以得出最小地表径流值为4171。

5、流量划分：

通过【条件函数】提取流量大于最小地表径流值为4171的栅格值。输出为Con_FlowAcc.tif。可以看到大于最小地表径流值为4171的栅格值都被提取出来了。

6、河流链接

选择【水文分析】里面的【河流链接】工具，调出面板进行参数设置。输入河流链接栅格选择【Con_FlowAcc.tif】数据；输入流向栅格数据为【FlowDir.tif】，输出为StreamL(HLLJ).tif保存路径即可。


【河流链接】完成：从分析结果可以看到，链接之后离散的像元会被串联起来，形成独立的河流。如上图所示，河网已经被串联成17条，仍然使用栅格的形式表示。

7、河网分级：

选择【水文分析】里面的【河网分级】工具，调出面板进行参数设置。输入河流栅格选择【StreamL(HLLJ).tif】数据；输入流向栅格数据为【FlowDir.tif】，输出为StreamO(HWFJ).tif保存路径即可。

8、分级后河网进行矢量化：

选择【水文分析】里面的【栅格河网矢量化】工具，调出面板进行参数设置。输入河流栅格选择【StreamO(HWFJ).tif】数据；输入流向栅格数据为【FlowDir.tif】，输出为StreamT(HWSLH).shp保存路径即可。

9、要素折点转点工具：

运行ArcToolbox，打开【数据管理】工具，选择【要素】里面的【要素折点转点】工具，执行时点类型选择【end】。输出为StreamTHWSLH_HK.shp保存路径即可。

如下，【要素折点转点】完成：这样提取出来的就是河口的位置。这个点在水文分析中又称之为“出水口”或“倾泻点”但如果这种方式提取的出水口必须要使用【捕捉倾泻点】工具来重新校正。

10、捕捉倾泻点：

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

选择【水文分析】里面的【捕捉倾泻点】工具，调出面板进行参数设置。输入栅格选择【StreamTHWSLH_HK】数据；输入流量栅格数据为【FlowAcc.tif】，输出为SnapPou(BZQXD).tif保存路径即可。

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

11、提取集水区域

选择【水文分析】里面的【分水岭】工具，调出面板进行参数设置。输入流向栅格数据为【FlowDir.tif】；输入栅格选择【SnapPou(BZQXD).tif】数据；输出为Watersh(JSQ).tif保存路径即可。

12、选择目标区流域

选取最大河流流经的几个集水区作为目标区域，执行栅格转面工具。

13、裁剪目标区DEM：

为什么要利用集水区（淹没区）的面去裁剪DEM。这样做的原因是在下一步使用条件函数去过滤DEM一定高度的像元时可以剔除掉不相关的数据。这部分是非常重要的。
运行ArcToolbox，打开【空间分析】工具，选择【提取】里面的【按掩膜提取】工具，调出面板进行参数设置，输入栅格选择原始DEM；输入掩膜数据选择提取出的矢量。

14、设定淹没高度

假设掩膜高度为10米，通过件函数】提取DEM中小于10的栅格值。输出为YMQ10tif保存路径即可。

假设掩膜高度为5，通过件函数】提取DEM中小于5栅格值。输出为YMQ5M.tif保存路径即可。

15、淹没洪水量计算：

运行ArcToolbox，打开【3D空间分析】工具，选择【表面功能】里面的【表面体积】工具，调出面板进行参数设置，输入表面数据为【YMQ10M.tif或YMQ5.tif】；参考Plane选择【BELOW】保存路径即可。


可以看到淹没该区域10m时，洪水量为106969914.76324立方米，淹没面积为6701788.01075平方米。

可以看到淹没该区域5时，洪水量为73751844.942729立方米，淹没面积为6664295.1152349平方米。

16、专题图制作：

五、总结

以上就是如何使用ArcGIS开展洪水淹没区分析、提取淹没区并计算淹没区洪水量的详细说明，如果想获取教程中实验数据，请关注微信公众号【GIS小能手】然后在后台留言。