CN113688755B - 基于六边形格网的多流向流域特征提取方法 - Google Patents
基于六边形格网的多流向流域特征提取方法 Download PDFInfo
- Publication number
- CN113688755B CN113688755B CN202111004964.XA CN202111004964A CN113688755B CN 113688755 B CN113688755 B CN 113688755B CN 202111004964 A CN202111004964 A CN 202111004964A CN 113688755 B CN113688755 B CN 113688755B
- Authority
- CN
- China
- Prior art keywords
- grid
- hexagonal
- flow
- hexagonal grid
- water flow
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000605 extraction Methods 0.000 title claims abstract description 19
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 87
- 238000000034 method Methods 0.000 claims abstract description 7
- 238000004364 calculation method Methods 0.000 claims description 17
- 238000012876 topography Methods 0.000 claims description 4
- 238000010276 construction Methods 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 5
- 238000012545 processing Methods 0.000 description 3
- 101000667821 Homo sapiens Rho-related GTP-binding protein RhoE Proteins 0.000 description 2
- 102100039640 Rho-related GTP-binding protein RhoE Human genes 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 1
- 238000004162 soil erosion Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/40—Scaling of whole images or parts thereof, e.g. expanding or contracting
- G06T3/4007—Scaling of whole images or parts thereof, e.g. expanding or contracting based on interpolation, e.g. bilinear interpolation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/05—Geographic models
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/10—Constructive solid geometry [CSG] using solid primitives, e.g. cylinders, cubes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
- G06T17/205—Re-meshing
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Graphics (AREA)
- Remote Sensing (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种基于六边形格网的多流向流域特征提取方法,该方法首先按双线性插法将规则格网DEM转换成球面六边形格网表示,分别计算得到六边形格网的顶点高程和中心点高程;然后依次连接顶点和中心点,将每个六边形格网划分成三角形,计算每个三角形中心点高程,根据三角形与邻近六边形格网高程差确定水流流向;根据与邻近六边形格网高程差比例,计算格网流量;最后遍历格网,跟踪流域路径。本发明利用六边形格网作为基础格网,解决了基于规则格网DEM流域特征提取中,因边邻近、角邻近格网距离不一致引起的地形表达失衡性问题;利用三角形与邻近六边形格网高程差比例,确定流向和流量,解决了原六边形单流向算法无法真实描述水流特征的问题。
Description
技术领域
本发明涉及流域水文特征提取技术领域,特别涉及一种基于六边形格网的多流向流域特征提取方法。
背景技术
流域水文特征提取与分析是水文水资源的重要组成部分,为水文建模、地形地貌建模、水土侵蚀评估、水资源管理等提供重要的基础信息。平向汇流的转移路径通常以水流方向、径流节点、径流网络和流域边界等水文要素表达。其中,水流方向是流域水文特征提取的基础。
随着DEM数据获取能力的提升,基于DEM的流域水文特征提取得到广泛研究。其中基于规则格网的DEM因其数据结构简单、计算处理方便等特点成为流域水文特征提取的主要研究对象。基于规则格网DEM的水流方向算法主要的方法包括单流向算法、多流向算法和混合流向算法。
单流向算法假设每个格网有且仅有一个水流方向。典型的算法有D8算法:根据与8个邻近格网间的高程差确定水流方向,水流向高程差最大的格网。Rho8算法:引入一个均匀分布的随机变量计算坡度权重,在D8算法的基础上,增加了水流方向选择的随机性。
多流向算法假设每个格网有多个水流方向,即水流可以流向比它低的所有相邻格网。典型的算法有FD8算法:在D8算法基础上,根据高程差和比例分配流向和水流量。FRho8算法:在Rho8算法基础上添加多方向权重因子。多流向算法避免了水流偏向某格网,能较好的模拟处于坡面的水域特征。
混合流向算法将单流向和多流向算法结合,根据流域地形特征,分别采用不同的流向算法。在水流辐合区域,使用单流向算法,在坡面-水流发散区域,使用多流向算法。
基于六边形格网的算法:因为六边形格网只有边邻近格网,与邻近格网中心点距离相等,在处理最近邻域、移动路径等邻域处理问题中具有一定的优势。将基于规则格网的DEM数据转换至六边形格网,利用与邻近格网间的高程差,计算水流方向。
现有方法存在如下不足:
(1)正方形格网与其边邻近和角邻近格网间距离之比为流域水文特征提取中水流量的计算依赖于格网间距离,不一致的格网间距导致边邻近方向和角邻近方向上地形表达的失衡性,影响流域水文特征提取精度。
(2)基于六边形格网的地形特征提取研究中,水流方向基于单流向计算,而真实流域水流特征较复杂,在水流辐合区,单个格网水流全部流向高程更低格网,而在坡面,水流呈发散状。因此,单流向算法不能反映坡面的发散水流状况,无法真实描述水流特征。
发明内容
本发明针对现有技术的缺陷,提供了一种基于六边形格网的多流向流域特征提取方法。
为了实现以上发明目的,本发明采取的技术方案如下:
一种基于六边形格网的多流向流域特征提取方法,以球面六边形格网作为基础格网表达地形,分解每个六边形格网为三角形,根据与邻近六边形格网的高程差确定水流流向,根据高程差比例计算流量。
进一步地,基于六边形格网的多流向流域特征提取方法,包括以下步骤:
S1:基于六边形格网的DEM构建,包括以下子步骤:
S11:根据规则格网DEM分辨率,确定球面六边形格网剖分层次,按照Snyder等积投影,生成正二十面体六边形格网;
S12:根据规则格网DEM范围,确定参与计算的六边形格网;
S13:利用双线性插值法,计算六边形格网顶点高程;
S14:根据六边形格网顶点高程均值获得格网中心点高程;
S15:基于六边形格网的DEM的区域预处理,即对区域洼地进行填洼处理。
S2:水流多流向计算,包括以下子步骤:
S21:连接每个格网的中心点和顶点,形成三角形格网;
S22:根据顶点插值,计算每个三角形格网的高程值;
S23:计算每个三角形与邻近六边形格网的高程差,记为Δi=E△-EHex;
S24:当Δi<0时,水流方向为当前三角形单元指向邻近格网。
S3:流量分配;
假设当六边形格网与其邻近格网间存在高程差时,六边形格网中的水流量全部流向其高程值较小的邻格网中。设每个格网单元的水流量为Fo,流入格网单元的水流量为Fin,每个格网单元水流出阈值为ε,流出格网单元的水流总量为Fout=Fo+Fin-ε。根据每个格网单元的流向和高程差分配流入邻近格网的水流量。
S31:根据多流向计算方法,统计Δi<0的个数,即为六边形格网的出水口总数,记为N;
S32:当N=1时,表示仅有一个出水口,则六边形格网水流总量Fout全部流入对应邻近格网;
S33:当N≠1时,表示六边形格网存在多个出水口,根据每个出水口高程差Δi确定水流量。
首先,归一化Δi,计算公式如式(1)所示:
其次,按比例分配水流量,计算公式如式(2)所示:
Fout_i=Fout×λi (2)
S34:更新邻近六边形格网水量总量,Fin=Fin+Fout_i。
S4流域路径跟踪;
使用广度优先遍历法则,逐格网遍历,判定水流方向,计算水流量,跟踪流域路径。
与现有技术相比,本发明的优点在于:
1.利用六边形格网作为基础格网,解决了基于规则格网DEM流域特征提取中,因边邻近、角邻近格网距离不一致引起的地形表达失衡性问题。
2.利用三角形与邻近六边形格网高程差比例,确定流向和流量,解决了原六边形单流向算法无法真实描述水流特征的问题。
附图说明
图1为本发明实施例基于六边形格网的多流向流域特征提取方法流程图;
图2为本发明实施例六边形格网DEM结果。其中图2a为球面正二十面体等积六边形格网;图2b为规则格网DEM与规则六边形DEM叠加及插值示意图;
图3为本发明实施例基于六边形格网DEM的多流向示意图。图3a为六边形各顶点形成的不等高六棱台,顶点与中心点连接构成三角形;图3b单个六边形格网确定的多个流向示意图;
图4为本发明实施例流向及流域特征曲线图,其中图4a为六边形流向图,图4b为流域特征曲线图。
具体实施方式
为使本发明的目的、技术方案及优点更加清楚明白,以下根据附图并列举实施例,对本发明做进一步详细说明。
如图1所示,一种基于六边形格网的多流向流域特征提取方法,包括以下步骤:
S1:基于六边形格网的DEM构建,包括以下子步骤:
S11:根据规则格网DEM分辨率,确定球面六边形格网剖分层次,按照Snyder等积投影,生成正二十面体六边形格网,如图2a所示;
S12:根据规则格网DEM范围,确定参与计算的六边形格网;
S13:利用双线性插值法,计算六边形格网顶点高程,如图2b所示;
S14:根据六边形格网顶点高程均值获得格网中心点高程;
S15:基于六边形格网的DEM的区域预处理,即对区域洼地进行填洼处理。
S2:水流多流向计算,包括以下子步骤:
S21:连接每个格网的中心点和顶点,形成三角形格网,如图3a所示;
S22:根据顶点插值,计算每个三角形格网的高程值;
S23:计算每个三角形与邻近六边形格网的高程差,记为Δi=E△-EHex;
S24:当Δi<0时,水流方向为当前三角形单元指向邻近格网,如图3b所示。
S3:流量分配;
如图3所示。假设当六边形格网与其邻近格网间存在高程差时,六边形格网中的水流量全部流向其高程值较小的邻格网中。设每个格网单元的水流量为Fo,流入格网单元的水流量为Fin,每个格网单元水流出阈值为ε,流出格网单元的水流总量为Fout=Fo+Fin-ε。水流向各邻近单元的概率是相等的,高程差越大,水流速度越快,流量越大,根据每个格网单元的流向和高程差分配流入邻近格网的水流量。
S31:根据多流向计算方法,统计Δi<0的个数,即为六边形格网的出水口总数,记为N;
S32:当N=1时,表示仅有一个出水口,则六边形格网水流总量Fout全部流入对应邻近格网;
S33:当N≠1时,表示六边形格网存在多个出水口,根据每个出水口高程差Δi确定水流量。
首先,归一化Δi,计算公式如式(1)所示:
其次,按比例分配水流量,计算公式如式(2)所示:
Fout_i=Fout×λi (2)
S34:更新邻近六边形格网水量总量,Fin=Fin+Fout_i。
S4流域路径跟踪;
使用广度优先遍历法则,逐格网遍历,判定水流方向,计算水流量,跟踪流域路径。
与基于栅格提取流域特征曲线结果对比,如图4所示本发明提取的流域特征曲线能更好地表达流域的细节特征,对流域方向描述更符合客观现实。
本领域的普通技术人员将会意识到,这里所述的实施例是为了帮助读者理解本发明的实施方法,应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各种具体变形和组合,这些变形和组合仍然在本发明的保护范围内。
Claims (1)
1.一种基于六边形格网的多流向流域特征提取方法,其特征在于:以球面六边形格网作为基础格网表达地形,分解每个六边形格网为三角形,根据与邻近六边形格网的高程差确定水流流向,根据高程差比例计算流量;
所述多流向流域特征提取方法,包括以下步骤:
S1:基于六边形格网的DEM构建,包括以下子步骤:
S11:根据规则格网DEM分辨率,确定球面六边形格网剖分层次,按照Snyder等积投影,生成正二十面体六边形格网;
S12:根据规则格网DEM范围,确定参与计算的六边形格网;
S13:利用双线性插值法,计算六边形格网顶点高程;
S14:根据六边形格网顶点高程均值获得格网中心点高程;
S15:基于六边形格网的DEM的区域预处理,即对区域洼地进行填洼处理;
S2:水流多流向计算,包括以下子步骤:
S21:连接每个格网的中心点和顶点,形成三角形格网;
S22:根据顶点插值,计算每个三角形格网的高程值;
S23:计算每个三角形与邻近六边形格网的高程差,记为Δi=E△-EHex;
S24:当Δi<0时,水流方向为当前三角形单元指向邻近格网;
S3:流量分配;
假设当六边形格网与其邻近格网间存在高程差时,六边形格网中的水流量全部流向其高程值较小的邻格网中;设每个格网单元的水流量为Fo,流入格网单元的水流量为Fin,每个格网单元水流出阈值为ε,流出格网单元的水流总量为Fout=Fo+Fin-ε;根据每个格网单元的流向和高程差分配流入邻近格网的水流量;
S31:根据多流向计算方法,统计Δi<0的个数,即为六边形格网的出水口总数,记为N;
S32:当N=1时,表示仅有一个出水口,则六边形格网水流总量Fout全部流入对应邻近格网;
S33:当N≠1时,表示六边形格网存在多个出水口,根据每个出水口高程差Δi确定水流量;
首先,归一化Δi,计算公式如式(1)所示:
其次,按比例分配水流量,计算公式如式(2)所示:
Fout_i=Fout×λi (2)
S34:更新邻近六边形格网水量总量,Fin=Fin+Fout_i;
S4流域路径跟踪;
使用广度优先遍历法则,逐格网遍历,判定水流方向,计算水流量,跟踪流域路径。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111004964.XA CN113688755B (zh) | 2021-08-30 | 2021-08-30 | 基于六边形格网的多流向流域特征提取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111004964.XA CN113688755B (zh) | 2021-08-30 | 2021-08-30 | 基于六边形格网的多流向流域特征提取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113688755A CN113688755A (zh) | 2021-11-23 |
CN113688755B true CN113688755B (zh) | 2023-08-08 |
Family
ID=78584051
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111004964.XA Active CN113688755B (zh) | 2021-08-30 | 2021-08-30 | 基于六边形格网的多流向流域特征提取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113688755B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115546443B (zh) * | 2022-11-25 | 2023-03-24 | 西南林业大学 | 一种球面六边形格网的局部等距性优化方法和系统 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2010126176A1 (ko) * | 2009-04-28 | 2010-11-04 | Choi Joonsoo | 등고선이 표현하는 지역을 다각형 세부 영역들로 분해하고 수치표고모델 데이터와 지형표면구성 데이터를 산출하는 방법과, 그 방법을 실행하기 위한 프로그램이 기록된 기록 매체 |
CN102915227A (zh) * | 2012-09-03 | 2013-02-06 | 南京师范大学 | 面向大区域流域提取的并行方法 |
CN104331870A (zh) * | 2014-11-25 | 2015-02-04 | 武汉大学 | 一种适于规则网格dem的结构化多尺度重构方法 |
CN105138722A (zh) * | 2015-07-14 | 2015-12-09 | 南京师范大学 | 基于数字河湖网络的平原河网区流域集水单元划分方法 |
CN105303612A (zh) * | 2014-12-03 | 2016-02-03 | 河南理工大学 | 一种基于不规则三角网模型的数字河网提取方法 |
CN106845074A (zh) * | 2016-12-19 | 2017-06-13 | 中国人民解放军信息工程大学 | 建立六边形格网模型的方法、洪水推演模拟方法及其系统 |
CN111062128A (zh) * | 2019-12-16 | 2020-04-24 | 南京大学 | 一种流域汇流模拟估测方法 |
CN112084643A (zh) * | 2020-08-31 | 2020-12-15 | 西湖大学 | 基于数字高程和土壤参数的流域提取方法 |
-
2021
- 2021-08-30 CN CN202111004964.XA patent/CN113688755B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2010126176A1 (ko) * | 2009-04-28 | 2010-11-04 | Choi Joonsoo | 등고선이 표현하는 지역을 다각형 세부 영역들로 분해하고 수치표고모델 데이터와 지형표면구성 데이터를 산출하는 방법과, 그 방법을 실행하기 위한 프로그램이 기록된 기록 매체 |
CN102915227A (zh) * | 2012-09-03 | 2013-02-06 | 南京师范大学 | 面向大区域流域提取的并行方法 |
CN104331870A (zh) * | 2014-11-25 | 2015-02-04 | 武汉大学 | 一种适于规则网格dem的结构化多尺度重构方法 |
CN105303612A (zh) * | 2014-12-03 | 2016-02-03 | 河南理工大学 | 一种基于不规则三角网模型的数字河网提取方法 |
CN105138722A (zh) * | 2015-07-14 | 2015-12-09 | 南京师范大学 | 基于数字河湖网络的平原河网区流域集水单元划分方法 |
CN106845074A (zh) * | 2016-12-19 | 2017-06-13 | 中国人民解放军信息工程大学 | 建立六边形格网模型的方法、洪水推演模拟方法及其系统 |
CN111062128A (zh) * | 2019-12-16 | 2020-04-24 | 南京大学 | 一种流域汇流模拟估测方法 |
CN112084643A (zh) * | 2020-08-31 | 2020-12-15 | 西湖大学 | 基于数字高程和土壤参数的流域提取方法 |
Non-Patent Citations (1)
Title |
---|
王政 ; 赵学胜 ; 罗富丽 ; 李亚路 ; 孙中昶.球面三角格网单元面积的分形计算.《武汉大学学报(信息科学版)》.2019,1541-1546+1554. * |
Also Published As
Publication number | Publication date |
---|---|
CN113688755A (zh) | 2021-11-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113742910A (zh) | 基于中小流域洪水预报的水库来水量预警预报方法及系统 | |
CN108345735A (zh) | 一种流域水文模型参数校准方法 | |
CN110543984B (zh) | 一种基于水安全的城市水面率规划方法 | |
CN108022047A (zh) | 一种海绵城市水文计算方法 | |
CN108755565B (zh) | 一种多空间尺度流域产流产沙预测方法及装置 | |
CN108021780A (zh) | 一种基于无规则非结构网格模型的山洪动态仿真方法 | |
CN110929443B (zh) | 一种基于高精度地形概化的二维洪水模拟方法 | |
CN107038308B (zh) | 一种基于线性内插的规则格网地形建模方法 | |
CN115130396A (zh) | 一种河道型水库库区的分布式水文模型建模方法 | |
CN113688755B (zh) | 基于六边形格网的多流向流域特征提取方法 | |
CN112632871B (zh) | 基于遥感的无资料自由溢流水库出流过程动态估算方法 | |
CN114580283B (zh) | 山溪性强潮河口分汊段中长期动力地貌演变数值模拟方法 | |
CN110990926B (zh) | 一种基于面积修正率的城市地表建筑水动力学仿真方法 | |
CN111985129A (zh) | 一种城市暴雨内涝精细化模拟方法 | |
CN110287595B (zh) | 一种城市不同下垫面减灾效应分析方法 | |
CN115905963A (zh) | 一种基于支持向量机模型的洪水预报方法及系统 | |
CN113379828A (zh) | 一种融合地表形态特征的坡长提取方法 | |
CN110457772B (zh) | 一种结合平面曲率和最陡下坡方向的dem流向估计方法 | |
CN115564238B (zh) | 一种基于水动力机制的小水库洪水预报调洪分析方法 | |
CN112182814A (zh) | 一种基于稀疏断面点据的河道水下地形建模方法 | |
CN109635317B (zh) | 一种高海拔地区cfd仿真地形拓扑方法 | |
CN113869804B (zh) | 一种洪涝灾害下的电网设备风险预警方法及系统 | |
CN114896909A (zh) | 一种基于水位高度的明渠流量计算方法 | |
CN105740812B (zh) | 一种基于数字表面模型的城市汇水区提取方法 | |
CN111861014B (zh) | 一种分布式单位线洪水预报方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |