DEM水文分析中一种有效消除伪水道的简易方法.docVIP

　　文章编号:049420911(20081220053203中图分类号:P208　　　　　　文献标识码:B D EM 水文分析中一种有效消除伪水道的简易方法 董有福 1,2 ,汤国安1,罗明良 3 (1.南京师范大学虚拟地理环境教育部重点实验室,江苏南京210046;2.南京工业大学 土木工程学院,江苏南京210009;3.西华师范大学国土资源学院,四川南充637002 A S i m ple M ethod to Eli m i n a te Pseudo Channels Effecti vely i n Hydrolog i c Ana lysis Ba sed on D E M DONG You 2fu,T ANG Guo 2an,LUO M ing 2liang 摘要:宽广河谷平坦地形区的存在,对基于DE M 提取的水系结构产生了严重影响,目前国内外已提出很多方法来处理这种平坦地 形情况。在对这些方法进行总结回顾的基础上,提出一种简易的伪水道消除方法。该方法是在常用的河网提取算法的基础上,采用缓冲区分析和栅格叠置分析原理实现河网的二次有效提取。实例研究表明,该方法原理简单、易于操作实现、效果明显,使基于 DE M 在宽广河谷区生成的河网结构趋于合理,河网提取结果更加准确。 关键词:伪水道;水文分析;DE M;数字地形分析 　　收稿日期:2008206210 基金项目:国家自然科学基金资助项目40571120 作者简介:董有福(19762,男,河南信阳人,博士生,主要从事数字地形分析和空间数据库方面的研究。 　　 一、引　言数字地形分析技术在数字流域研究中获得了广泛应用,大比例尺系列的数字高程模型数据源的建立同时为水文模拟分析提供了较高精度保障。沟谷网络的提取是相关流域特征参数的基础,而各个栅格单元水流方向的确定是自动生成连续和单宽河网的前提和关键。提取沟谷网络及相关流域特征参数的算法经过20多年的发展,内容不断丰富和相对完善。在众多沟谷网络提取的算法中,目前使用最普遍的确定排水方向的方法是O ’Callaghan 和Mark 提出的最大坡降法(D8,即中心栅格的水流向周围八个栅格单元中高程最低的那个栅格单元。但是,由于受DE M 水平或垂直分辨率的限制、DE M 生成过程中的内插和输出结果的取整以及高程数据误差的影响,或者真实平坦地形的存在,如宽广谷地、河漫滩、河流阶地和人工梯田等,使得无论使用哪一种方法,平坦地区的栅格单元水流方向都难以有效确定,结果是容易出现大量的平行河道和河道变形等问题,当汇流阈值较小时尤为明显。因此,平坦地形的存在严重影响着模拟河网的正确提取。 国内外学者针对DE M 沟谷提取中的平地问题作了大量研究。Jens ont 等提出首先寻找出与具有 流向的网格单元相邻的那些网格单元,然后确定这 些网格单元的流向,重复进行这一步骤直到所有的 平地网格单元都确定了流向为止[1] ;Martz 等在DEDNM 中采用微地貌起伏算法实现了在平坦地形 区赋流向,但它产生的是一种平行模式的流向[2] ;Andrea 提出了一种汇聚模式的赋流向方法,使河流 在平坦地形区内通过的距离最短[3] ;Garbrecht 等在T OP AZ 中改进了微地貌起伏算法,实现了汇聚模式 的赋流向[4] 。以上这些方法没有考虑模拟与实测河网之间的精确拟合,因此,自动生成的河网是不可能与实际河网相吻合的。孔凡哲等通过增加输入地形高程信息避免DE M 中平坦区域生成,从而使由DE M 生成的河网与实际河网能够精确拟合[6] ;谢顺平等提出以出流代价法构建河谷平地栅 格流向和流序等[7] 。 上述这些方法在某种程度上都能够解决平坦地形流向问题,但是由于需要较多的辅助资料,或者算法比较复杂,使得适用性受到了较大限制,因此在很多场合下仍然采用手工编辑伪水道方式,效率相对比较低下。本文则在常用河网提取算法的基础上,采用缓冲区分析和叠置分析等简易方法,通过二次提取河网法实现河谷平坦地形区大量伪水道的消除,从而达到与实际自然河网的有效拟合。 3 52008年　第12期　　　　　　　　　　　　测　绘　通　报 　　 二、原理与方法1.基本原理 正是由于平坦地形区栅格单元的水流方向难 以确定,才使得较宽河谷大量伪水道的出现。因此,该方法的基本思想是,对主干河道一定距离内的区域,在原有高程的基础上有规律地逐步降低其高程值,人为使得平坦地形区栅格单元产生向主干河道的水流方向,从而有效避免平行河道的出现。该区域栅格单元降低的高程值可以用式(1表示。 h =(1-d /d 0?h 0 (1 式中,h 表示每个栅格单元降低的高程值;h 0表示主干河道中心降低的高程值;d 表示栅格单元到主干河道的缓冲距离;d 0表示缓冲半径距离。 从式(1