一种大区域洪水淹没范围快速提取的分块种子蔓延算法.docxVIP

第49卷第4期 第49卷第4期 华中师范大学学报(自然科学版) V01．49 NO．4 2015年8月 JOURNAL 0F CENTRAL CHINA NORMAL UNIVERSITY(Nat．Sci．) Aug．2015 文章编号：1000—1190(2015)04—0603—05 一种大区域洪水淹没范围快速提取的 分块种子蔓延算法 杨启贵+，王汉东 (长江水利委员会长江勘测规划设计研究有限责任公司，武汉430010) 摘 要：快速、准确的获取洪水淹没范围是洪灾损失评估和防洪决策的核心环节．常规种子蔓延算 法由于存在大量的递归运算，效率不高且无法处理大区域海量地形数据．针对上述问题，提出了一 种大区域洪水淹没范围快速提取的分块种子蔓延算法，将研究区域分成预先设定大小的块，种子 蔓延以块为单位进行组织处理，块的数量远小于原始DEM格网单元的数量．该算法有效解决了常 规种子蔓延算法递归运算层次过深的问题，可实现大区域海量数据在给定淹没洪水水位条件下的 淹没区快速提取． 关键词：数字高程模型；种子蔓延算法；数据分块；淹没分析 中图分类号：TP391 文献标识码：A 洪水淹没分析是水文预测预报、防洪调度以及 种子蔓延算法是一种基于区域空间特征的扩散探测 洪涝灾害评估的一项重要内容，快速、准确、科学地 算法，其核心思想是将给定的种子点赋予特定的属 获取淹没区范围是洪灾损失评估及防洪决策的核 性，然后在区域上4方向或8方向扩散，求取满足给 心环节口]．将GIS技术与水文水动力学模型相结 定条件且具有连通特性的点．利用种子蔓延算法计 合进行洪水淹没灾害评估，已成为防洪减灾领域的 算淹没区，就是按给定水位条件，求取满足精度、连 一个重要研究课题[2。1 ．基于DEM数据的洪水淹 通性要求的点的集合，该集合给出的连续平面就是 没分析包括两类方法：一类是结合水文水动力学模 淹没区范围．种子蔓延算法主要采用递归的方式进 型构建洪水演进模型[6。2|，另一类是根据设定的静 行扩展探测，算法的效率受制于递归的层次，递归 态水位面计算淹没范围[2。4’14‘1 ．洪水演进模型可 层次太深会导致算法效率降低，甚至导致程序崩 以较为精确地模拟洪水淹没过程，但是需要较多的 溃，因此，无法提取大区域的洪水淹没区．降低递归 水文水力学参数条件，建模过程复杂．基于给定静 层次是种子蔓延算法应用于大区域海量数据洪水 水位的洪水淹没分析是洪水淹没达到最终平衡状 淹没需要解决的关键问题．刘仁义等人采用“跳步 态的近似模拟，因其计算简单，能快速确定淹没范 探测”的方式，通过增加种子点探测步距来提高计 围和水深分布，因此得到广泛应用．静水位下洪水 算效率[33；丁志雄等人预先划定大致的淹没范围， 淹没分析主要分为两种情形[2。4]：第1种情形称为 通过缩小研究区域来提高计算效率H1；张煜、沈定 无源淹没，凡是高程值低于给定水位的点，皆计入 涛等人针对海量地形数据，提出了基于数据压缩的 淹没区，相当于整个区域均匀降水；第2种情形称 淹没范围提取算法[14。1 5|．增加探测步距其实是降 为有源淹没，需要考虑水流连通性的问题，洪水只 低空间分辨率，即增加格网单元的尺寸，会降低分 能淹没到它能流到的地方，相当于高发洪水向邻域 析结果的精度；预先判断洪水的大致淹没范围常常 泛滥．对于环形山洼地，第1种淹没计算可能导致 是比较困难的，特别是对于地势相对平坦的研究区 环形山内外都形成淹没区，而在第2种情形中，外 域．本文从降低传统种子蔓延算法递归层次的问题 来的洪水若未及山顶，则只能在环形山外围形成环 出发，提出基于分块的种子蔓延算法，将研究区域 形淹没区，而环形山内部洼地不会被淹没． 分成小的块，在各个块上分别运用种子蔓延算法提 有源淹没常用的算法是种子蔓延算法[2。4’141 ， 取淹没区域，可有效降低算法执行过程中的递归层 收稿日期：2015-04—03． 基金项目：水利部科技推广计划项目(TGl418)；长江勘测规划设计研究有限责任公司自主科研项目(CX2013209) *E-mail：yangqigui@cjwsjy．corn．cn． 万方数据 华中师范大学学报(自然科学版) 华中师范大学学报(自然科学版) 第49卷 次．算法不受研究区域大小的限制，可用于提取大 { 区域洪水淹没范围． int BlockID=0；／／块编号 int SRow一0；／／块起始行号 1分块种子蔓延算法 int SColumn=0；／／块起始列号 1．1算法基本思想 int ERow=0；／／块末行号 种子蔓延算法的递归层次与地形数据量有很 int EColumn=0；／／块末列号 大关系，随着D