01概述

洪水灾害是我国发生频率高、危害范围广、对国民经济影响最为严重的自然灾害，亦是威胁人类生存的十大自然灾害之一。以流域内水库的防洪库容密度作为体现水库的防洪能力的指标，若流域水库防洪库容密度大，则水库防洪能力强，该流域遇到洪涝灾害时，水库越体现其防洪减灾作用。这里，我们介绍一下如何使用ArcGIS计算水库库容量。

02数据获取

本教程所使用的数据是水经注微图下载的DEM数据，除了DEM数据，常见的GIS数据均可以使用水经注微图进行下载。

水经注微图

（关注微信公众号“水经注GIS”,回复“微图4.0”获取微图软件安装包）

03填洼处理

在ArcToolbox内点击“Spatial Analyst工具\水文分析\填洼”，调用填洼工具。

调用填洼工具

在显示的填洼对话框内，输入表面栅格为下载的DEM数据，如下图所示。

填洼设置

点击确定之后可以得到填洼的结果，如下图所示。

填洼结果

04获取流向数据

在ArcToolbox内点击“Spatial Analyst工具\水文分析\流向”，调用流向工具。

调用流向工具

在显示的流向对话框内，输入表面栅格数据为上一步填洼后的栅格。

流向设置

点击确定之后可以得到流向栅格，如下图所示。

流向栅格

05获取流量数据

在ArcToolbox内点击“Spatial Analyst工具\水文分析\流量”，调用流量工具。

调用流量工具

在显示的流量对话框内，输入流向栅格数据为上一步得到的流向栅格。

流量设置

点击确定之后可以得到流量栅格，如下图所示。

流量栅格

06捕捉倾泻点

虽然坝口位置便可以看成是出水口，但也不能直接使用。

因为它与原DEM不是同一份数据，位置上跟DEM对应的出水口还是有差别的，所以，必须要进行捕捉，以确保捕捉到的倾泻点是该点流量最大的栅格。

在ArcToolbox内点击“Spatial Analyst工具\水文分析\捕捉倾泻点”，调用捕捉倾泻点工具。

调用捕捉倾泻点工具

在显示的捕捉倾泻点对话框内，输入栅格数据或要素倾泻点数据为手工选择的倾泻点数据，输入蓄积栅格数据为之前得到的流量数据。

捕捉倾泻点设置

点击确定之后就可以得到捕捉的倾泻点，如下图所示。

捕捉的倾泻点

07计算集水区

在ArcToolbox内点击“Spatial Analyst工具\水文分析\分水岭”，调用分水岭工具。

调用分水岭工具

在显示的分水岭对话框内，输入流向栅格数据为之前获取到的流向栅格，输入栅格数据或要素倾泻点数据为上一步捕捉到的倾泻点数据。

分水岭设置

点击确定之后可以得到集水区，如下图所示。

集水区

08获取集水区范围

在ArcToolbox内点击“转换工具\由栅格转出\栅格转面”，调用栅格转面工具。

调用栅格转面工具

在显示的栅格转面对话框内，输入栅格为上面得到的集水区。

栅格转面设置

点击确定之后可以得到集水区范围，如下图所示。

集水区范围

09裁剪DEM数据

在ArcToolbox内点击“数据管理工具\栅格\栅格处理\裁剪”，调用裁剪工具。

调用裁剪工具

在显示的裁剪对话框内，输入栅格为下载的DEM数据，输出范围为上一步得到的集水面范围，勾选上使用输入要素裁剪几何。

裁剪设置

点击确定之后可以看到裁剪后的DEM数据，如下图所示。

裁剪后的DEM数据

10计算蓄水库容

在ArcToolbox内点击“Spatial Analyst工具\条件分析\设为空函数”，调用设为空函数工具。

调用设为空函数工具

在显示的设为空函数对话框内，输入条件栅格数据为裁剪后的DEM数据，表达式设置为Value >650，表示将高程值大于650的部分设置为NoData，输入条件为false时所取的栅格数据或常量值为裁剪后的DEM。

设为空函数设置

点击确定之后可以得到处理后的DEM数据，如下图所示。

处理后DEM

在ArcToolbox内点击“3D Analyst工具\功能性表面\表面体积”，调用表面体积工具。

调用表面体积工具

在显示的表面体积对话框内，输入表面为处理后的栅格，输出文本文件选择一个合适的位置，参考平面选择BLOW，表示计算下垫面的体积。

表面体积设置

点击确定之后可以得到表面体积信息（即水库库容），单位为立方米，如下图所示。

表面体积信息

11结语

以上就是如何使用ArcGIS计算水库库容量的详细说明，主要包括了数据获取、填挖处理、获取流向数据、获取流量数据、捕捉倾泻点、计算集水区、获取集水区范围、裁剪DEM数据和计算蓄水库容等功能。如果在使用过程中有任何疑问，请联系我们客