一、数据获取
DEM数据来源于地理空间数据云，该网站是由中科院计算机网络信息中心于2008年创立的地学大数据平台。
或直接进入小程序下载（含全国12.5米-90米DEM下载）：各种数据、插件下载
二、填洼
将DEM数据中凹陷的区域填充至与倾斜点同样高度，这里的【Z限制】说的是设定一个特定的值，凹陷区域的最低点高程与倾斜点高程之间的距离大于这个特定的值视为凹陷区与为有效凹陷，该工具会填充所有差值小于该特定值的凹陷区域。
工具位置：【Spatial Analyst】→【水文分析】→【填洼】
三、流向
使用【流向】工具计算每个像元的水流方向，这里选用【D8】单向流向分配法来计算每个像元最陡下降方向。
工具位置：【Spatial Analyst】→【水文分析】→【流向】
四、流量
使用【流量】工具创建每个像元累积流量的栅格，输出数据类型根据需要选择是【FLOAT】浮点型还是【INTEGER】整型或是【DOUBLE】双精度型。
工具位置：【Spatial Analyst】→【水文分析】→【流量】
五、捕捉倾泻点
这里倾泻点的位置是人工选择的，一般是大坝的出水口，但不能直接进行使用，因为大坝的出水口不代表就是DEM数据的出水口，所以需要捕捉，以确保捕捉到的倾泻点是该点流量最大的栅格。
工具位置：【Spatial Analyst】→【水文分析】→【捕捉倾泻点】
六、集水区
使用【集水区】工具获取集水区。
工具位置：【Spatial Analyst】→【水文分析】→【集水区】
七、栅格转面
使用【要素转面】工具将获取的范围转为SHP格式的文件。
工具位置：【转换工具】→【由栅格转出】→【栅格转面】
八、提取DEM
使用【按掩膜提取】工具获取集水区范围内DEM数据。
工具位置：【Spatial Analyst】→【提取分析】→【按掩膜提取】
九、设置空函数
假定该区域内最高的蓄水高程为300米，利用【栅格计算器】与【SetNull】函数将高于300米的像元设置为NotData。
工具位置：【Spatial Analyst】→【地图代数】→【栅格计算器】
十、计算容量
使用【表面体积】工具来计算水库容量，这里【参考平面】选用的是【BELOW】。
工具位置：【3D Analyst】→【功能性表面】→【表面体积】
最终结果为112849109.828立方米。
- END -
