流域水系自动提取在西苕溪流域的应用.pdf

　第 23 卷 　第 4 期 热 　带 　地 　理 Vol23 , No 4 　 　2003 年 12 月 TROPICAL GEO GRA PH Y Dec , 2003 　 流域水系自动提取在西苕溪流域的应用 李昌峰1 ,2 , 赵 　锐 1 ( 1 中国科学院南京地理与湖泊研究所 , 江苏 南京 　2 10008 ; 2 中国科学院研究生院 , 北京 　100039) 摘 　要 : 论述了如何基于栅格 D EM 自动提取流域自然水系的原理 、方法和流程 。以西苕溪中上游流域为例 , 根 据 D EM 精度 、上游集水区面积阈值和下垫面地形的不同 , 对所提取水系进行了比较 。针对在平均地形坡度小于 3°的平坦区域所提取水系与实际河网偏差较大的问题 , 提出了利用主干河道和平原水系数字化作为约束条件生 成河网的新方法 。 关键词 : 水系 ; D EM ; 自动提取 ; 集水区面积阈值 中图分类号 : P334 + 92 　　　　文献标识码 : A 　　　　文章编号 : 100 1 - 522 1 (2003) 04 - 03 19 - 05 　 　数 字 高 程 模 型 ( Digit al Elevation Model ) 个高程数据点水流方向 , 然后根据高程数据点的水 (D EM) 是描述地面高程值空间分布的一组有序数 流方向数据计算每一个高程数据点的上游集水区 , 组 , 通常有栅格 ( GR ID) 、不规则三角网 ( T IN s) 再根据上游集水区的高程数据 , 用阈值法确定属于 和矢量 (数字划线) (DL Gs) 3 种形式 , 其中栅格 水系的高程数据点 , 最后 , 根据水流方向数据 , 从 D EM 因计算处理方法简单有效 , 且与遥感数据在 水系源头开始将整个水系追索出来[5 ] , 同时进行 结构上容易匹配 , 在流域水文模拟中得到普遍应 子流域的划分 , 并对河网与子流域进行编码 , 构建 用[ 1 ] 。基于栅格 D EM 提取流域水系一般采用 了 河网结构拓扑关系 。由D EM 自动提取流域 自然水 O ’Callaghan 和 Mark 的坡面径流模拟方法[2 ] , 这 系的模块化流程见图 1 。 种算法通常被称为 D8 算法 , 后来 J ense[3 ] 、Gar brecht [4 ] 等在此基础上作了进一步的改进 , 提出了 Rho8 、FRho8 、TA P ES - C 等算法 。不同的算法产 生的结果会有所差别 。由于 D8 是比较传统有效的 算法 , 因而得到了较为广泛的应用 。 目前 , 国内外学者的研究大都侧重于如何利用 各种方法基于栅格 D EM 提取流域 自然水系 , 对于 D EM 精度 、下垫面地形复杂程度以及上游集水区 面积阈值如何影响流域水文模拟的研究则较少 。本 研究中利用 ArcView + Sp atial Analyst 及其相应的 扩展模块和脚本程序 , 基于 D8 算法 , 根据 D EM 精度 、上游集水区面积阈值和下垫面类型的不同 , 对西苕溪中上游流域进行了自然水系的提取试验 , 并对不同条件下生成的流域水系进行了对比分析 。 1 　基于 D EM 提取流域自然水系 12 　D EM 数据的预处理 11 　基本原理 从 D EM 中自动提取流域 自然水系的前提条件 　　从 D EM 中自动提取流域 自然水系是依据水总 是 : 在 D EM 中不能存在小平原和洼地 ,