水文模型往往必须应用在一个闭合流域中,否则流域边界处的水会流向流域外部,造成多个流域出口和水量不平衡的问题。Arcmap中的水文分析工具是很好用的提取流域的工具,可以为各种水文模型提供流域范围,本文利用Arcmap10.0以上版本,简单介绍如何通过DEM数据提取出DHSVM模型所需要的流域边界。
一、流域边界提取
1.载入DEM
载入DEM,裁剪至合适大小,保证整个流域都在范围之内。
2.填洼
打开ArcToolBox——Spatial Analysis工具——水文分析,打开填洼工具,输入DEM,生成填洼后的DEM。
3.计算流向
打开流向工具,输入填洼后的DEM,生成流向数据direct。
4.计算流量
打开流量工具,输入流向数据,计算汇流累积量amount。得到的流量数据中,值越大的地方(白色的地方)代表越处在下游流域出口处。
5.捕捉倾泻点
这一步是寻找流域出口的位置,需要用到已知的水文站的位置。
1.在当前数据库中新建点数据station
2.编辑器——开始编辑,选择新建的点要素。
3.在水文站的大概位置上,添加一个点,点的位置尽量靠近流量大的单元格(也就是流量数据中白色的区域)
4.停止编辑保存编辑内容。
5.打开捕捉倾泻点工具,输入的要素倾泻点数据是刚才的点数据,蓄积栅格数据是流量数据,输出流域出水口outpoint,捕捉距离可以自定义,一般一到两个单元格尺寸(比如dem分辨率为100,捕捉距离可以设置为200)。
生成之后给outpoint改个显眼的颜色,放大检查,看到station点旁边生成了一个小单元格,且这个单元格处在流量数据白色的区域(值比较高,流量大的地方),这个就是outpoint的位置。如果生成的outpoint不在流量数据白色的区域,说明捕捉距离太小了,需要放大捕捉距离重做。
6.提取流域边界
打开分水岭工具,输入流向数据,输入刚才的流域出水口outpoint,得到流域范围watershed。
至此可以得到只包含一个出水口的流域边界数据,DHSVM里面并不需要转成面的流域边界数据,所以这样就可以了,后续处理方法可以参考前面DHSVM栅格数据制作方法。
二、数据裁剪
有的人可能会发现得到的流域范围数据比DEM小很多,行列号也不一致,为了不必要的麻烦,可以对dem继续裁剪。
1.根据流域范围裁剪DEM
打开ArcToolBox——数据管理工具——栅格——栅格处理——裁剪,输入DEM,输出范围选择刚才的流域边界数据。下方保持裁剪范围勾选上,这样裁剪后的DEM行列号会和流域边界数据一致
其他DHSVM需要的土壤类型、植被类型数据,也可以通过这种方法裁剪至与流域边界数据一致。