基于dem提取流域河网的算法研究.docxVIP

基于dem提取流域河网的算法研究 0 分布式参数模型 水文模型分为总参数模型（lumpedpare模型s）和分布式参数模型（distribut参数模型s）。集总式参数模型把整个流域看成一个整体,不考虑影响水文过程的气候和下垫面条件空间差异,流域参数取其平均值。分布式参数模型通常将流域划分为一些小的地域单元,这些地域单元可以是格网也可以是子流域,通常假定这些地域单元内部是均一的,各有其对应参数,地域单元之间有一定拓扑关系,通过这种拓扑关系能够说明物质的传输方向。分布式参数模型考虑了流域内部气候和下垫面条件空间差异性,比集总式参数模型具有更高的空间分辨率。目前水文模型的研究重点已从集总式参数模型转移到GIS与分布式参数模型相结合的研究方面,分布式参数模型的建立需要能够确定流域内部排水结构等地表水文特征资料(如流向、河网、子流域等),而现有绝大多数分布式参数模型都是从DEM中自动提取这些信息的。 1 基于dem提取流域河网的算法的总结 1.1 dem格网阵列的扫描 Peuker和Douglas(1975)提出了一种简单的移动窗口算法,算法思路是使用一个2×2的窗口对DEM格网阵列进行扫描,将在窗口中具有最低高程值的点做标志,而剩余未做标志的点将表示山脊线上的点。类似地,对在窗口中具有最高高程值的点做标志,而剩余未做标志的点将表示山谷线上的点。通过这样计算,特征点便被提取出来了,这些点表示了潜在的沟谷线(或分水线)所在的位置。 1.2 坡面径流模拟算法 Mark(1984)及O’Callaghan和Mark(1984)注意到Peuker和Douglas方法不能提取出连续河网,因此他们提出了一种能提取连续河网的新算法,即坡面径流模拟算法。该算法以水流在地表沿着坡度最大方向流动为基础。坡面径流模拟算法的实现分四步:(1)确定DEM中每一栅格单元的水流方向;(2)计算每一栅格单元的上游汇水面积;(3)给定一个最小水道集水面积阀值;(4)提取流域边界及子流域。 1.3 线搜索算法 在O’Callaghan和Mark(1984)提出坡面径流模拟算法的同时,Yoeli(1984)提出了谷线搜索算法。该算法的计算步骤是:(1)首先用一个样条曲线函数找到DEM中的所有最低点;(2)然后从这些高程最低点生成谷线,每条谷线起始于这些高程最低点中没有参与生成谷线的高程最大值,将其延伸到与其最近的高程低值点,重复这样的操作直到该谷线碰到别的谷线,如湖泊或海洋;(3)重复步骤(2)直到遍历所有的高程最低点。 2 坡面径流模拟算法中存在的问题 从DEM中提取连续河网是大多数研究者一生致力追求的结果。移动窗口算法最大缺点就是提取的河网不连续,现在已经少有使用了。谷线搜索算法虽说也能提取出连续河网,但其只能应用于很小的DEM,更重要的是其不是基于水文学概念来提取河网的,对其他水文特征,因而很难与分布式参数模型相结合,这些都影响了该算法的广泛应用。坡面径流模拟算法由于其水文学基础,能产生连续河网,能与分布式参数模型紧密结合,且算法成熟稳定,易于程序实现,是目前应用最广泛、最普遍的算法。此处仅对坡面径流模拟算法中存在的问题进行探讨。该算法中主要存在下列问题:(1)水道起始点位置的确定;(2)凹陷与平坦区域水流方向的确定;(3)模拟河网与自然河网之间的偏差。 2.1 初始阀值的确定 对于水道起始点位置,坡面径流模拟算法采用最小水道集水面积阀值的概念,通常认为如果某处的上游集水面积等于某一给定阀值,则该点为河道的起始点,而大于该阀值的栅格点被定义为水道。但该阀值的给定具有很大的随意性和主观性,且值的大小将直接影响生成河网的疏密程度。事实上现在还无法找到一个确切的阀值来确定水道起始点位置,因为即使在很小的流域内,每一条沟谷的阀值都可能各不相同。 Vogt(2003)首先根据决定河网密度大小的环境特性,如气候、土地覆盖、地形、岩性、土壤等信息对流域内的景观地貌进行分类,然后分别提取不同景观地貌类型的坡度-面积关系并分别确定其水道集水面积阀值,最后对不同景观地貌取不同的水道集水面积阀值来提取河网,在意大利进行试验表明所提取的河网与实际河网在位置和河网密度方面都很吻合。 2.2 网向的排水系统 Tribe(1992)认为坡面径流模拟算法能产生连续河网的先决条件是DEM中每个栅格单元都必须被赋予一个确定的水流方向。D8算法在三种情况下确定栅格流向会出问题:(1)栅格有两个或多个可能的流向,(2)在平坦区域处,(3)在闭合洼地处。 对于栅格有两个或多个可能的流向,O’Callaghan和Mark(1984),Skidmore(1990)建议主观上选择顺时针从正北方向开始最先碰到的那个方向;Jenson和Domingue(1988)采用了类似的主观方法,然而对于有两种以上可能流向的栅格