CN110686861A - 一种基于四向高程差的净雨量转栅格格式方法 - Google Patents
一种基于四向高程差的净雨量转栅格格式方法 Download PDFInfo
- Publication number
- CN110686861A CN110686861A CN201911016711.7A CN201911016711A CN110686861A CN 110686861 A CN110686861 A CN 110686861A CN 201911016711 A CN201911016711 A CN 201911016711A CN 110686861 A CN110686861 A CN 110686861A
- Authority
- CN
- China
- Prior art keywords
- grid
- elevation
- net rainfall
- sum
- rainfall
- 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.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M10/00—Hydrodynamic testing; Arrangements in or on ship-testing tanks or water tunnels
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01W—METEOROLOGY
- G01W1/00—Meteorology
- G01W1/14—Rainfall or precipitation gauges
-
- 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
-
- 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
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Geometry (AREA)
- Environmental & Geological Engineering (AREA)
- Software Systems (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Fluid Mechanics (AREA)
- Theoretical Computer Science (AREA)
- Computer Graphics (AREA)
- Hydrology & Water Resources (AREA)
- Remote Sensing (AREA)
- Atmospheric Sciences (AREA)
- Biodiversity & Conservation Biology (AREA)
- Ecology (AREA)
- Environmental Sciences (AREA)
- Sewage (AREA)
Abstract
本发明提出了一种基于四向高程差的净雨量转栅格格式方法,包括:步骤S1,计算每个栅格在四个方向的高程差之和,根据所述高程差之和设置权重;步骤S2,根据计算得到的每个栅格的权重,计算不同栅格净雨量的分配比例;步骤S3,根据计算得到的降雨量的分配比例,将流域的净雨量数据转栅格格式,分配到具体的每个栅格内。本发明可以解决针对水文数据降尺度时面临的地形空间异质性问题。
Description
技术领域
本发明涉及水文数据降尺度技术领域,特别涉及一种基于四向高程差的净雨量转栅格格式方法。
背景技术
水文预报和水动力预报的耦合是水文研究的热点问题之一。耦合过程是先利用水文模型进行流域产汇流预报/模拟,再通过二维水动力学模型预报/模拟流域洪水过程,即利用水文模型模拟产汇流过程,产流计算出流域上一场降雨的净雨量(=降雨量-下渗-蒸散发);再利用水动力模型模拟地表洪水演进过程,水动力模型基于栅格数据可以对地形变化及地表建筑物的特征进行精细描述,因此在应用中需要将水文模型模拟得到的净雨量进行栅格化处理,转化为二维水动力模型可以使用的栅格格式数据。然而,目前尚未有成熟的流域净雨量转栅格格式方法,常用的方式为流域平均分配净雨量,此方法的缺点是未考虑流域地形上的空间异质性,会直接导致二维计算中出现较大误差,并导致计算的稳定性较差。
发明内容
本发明的目的旨在至少解决所述技术缺陷之一。
为此,本发明的目的在于提出一种基于四向高程差的净雨量转栅格格式方法。
为了实现上述目的,本发明的实施例提供一种基于四向高程差的净雨量转栅格格式方法,包括如下步骤:
步骤S1,计算每个栅格在四个方向的高程差之和,根据所述高程差之和设置权重;
步骤S2,根据计算得到的每个栅格的权重,计算不同栅格净雨量的分配比例;
步骤S3,根据计算得到的降雨量的分配比例,将流域的净雨量数据转栅格格式,分配到具体的每个栅格内。
进一步,在所述步骤S1中,采用数字高程模型DEM计算每个栅格在四个方向的高程差之和。
进一步,在所述步骤S1中,所述根据高程差之和计算权重,包括如下步骤:
当计算的高程差之和大于0时,则设置所述权重为高程差之和;
当计算的高程差之和小于0时,则设置所述权重为0。
进一步,在所述步骤S1中,采用下式计算权重,包括:
其中,Wi,j是栅格(i,j)的权重;i,j分别代表栅格x,y方向的序号;z0是当前计算栅格内的高程;zk(k=1,2,3,4)是相邻四个方向栅格的高程。
进一步,在所述步骤S2中,所述计算不同栅格净雨量的分配比例,包括:
其中,Wi,j是栅格(i,j)的权重;i,j分别代表栅格x,y方向的序号;m,n为二维栅格的x,y方向的栅格总数;Pi,j为不同栅格净雨量的分配比例。
进一步,在所述步骤S3中,所述将流域的净雨量数据转栅格格式,分配到具体的每个栅格内,包括:
WGi,j=WE×Pi,j,
其中,WGi,j表示净雨量的栅格值;WE为某一时刻的流域净雨量。
根据本发明实施例的基于四向高程差的净雨量转栅格格式方法,在净雨量分配时充分考虑了四个方向上高程差对汇流过程的影响,充分考虑了流域产汇流机理及水文过程的物理机制,降低了传统的平均分配方法带来的计算误差问题,计算过程简单、好操作、易推广。
本发明附加的方面和优点将在下面的描述中部分给出,部分将从下面的描述中变得明显,或通过本发明的实践了解到。
附图说明
本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解,其中:
图1为根据本发明实施例的基于四向高程差的净雨量转栅格格式方法的流程图;
图2为根据本发明实施例的正方形规则网格示意图。
具体实施方式
下面详细描述本发明的实施例,所述实施例的示例在附图中示出,其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的,旨在用于解释本发明,而不能理解为对本发明的限制。
本发明提出一种基于四向高程差的净雨量转栅格格式方法,可以解决针对水文数据降尺度时面临的地形空间异质性问题。
如图1所示,本发明实施例的基于四向高程差的净雨量转栅格格式方法,包括如下步骤:
步骤S1,计算每个栅格在四个方向的高程差之和,根据所述高程差之和设置权重。
在本步骤中,采用数字高程模型(Digital Elevation Model,简称DEM)计算每个栅格在四个方向的高程差之和。
在本发明的实施例中,基于图2所述正方形规则网格,计算每个栅格在四个方向的高程差之和。
根据高程差之和计算权重,包括如下步骤:
(1)当计算的高程差之和大于0时,则设置所述权重为高程差之和;
(2)当计算的高程差之和小于0时,则设置所述权重为0。
在本发明的实施例中,采用下式计算权重,包括:
其中,Wi,j是栅格(i,j)的权重;i,j分别代表栅格x,y方向的序号;z0是当前计算栅格内的高程;zk(k=1,2,3,4)是相邻四个方向栅格的高程。
步骤S2,根据计算得到的每个栅格的权重,计算不同栅格净雨量的分配比例。
在本步骤中,计算不同栅格净雨量的分配比例,包括:
其中,Wi,j是栅格(i,j)的权重;i,j分别代表栅格x,y方向的序号;m,n为二维栅格的x,y方向的栅格总数;Pi,j为不同栅格净雨量的分配比例。
步骤S3,根据计算得到的降雨量的分配比例,将流域的净雨量数据转栅格格式,分配到具体的每个栅格内。
具体的,将流域的净雨量数据转栅格格式,分配到具体的每个栅格内,包括:
WGi,j=WE×Pi,j,
其中,WGi,j表示净雨量的栅格值;WE为某一时刻的流域净雨量。
根据本发明实施例的基于四向高程差的净雨量转栅格格式方法,在净雨量分配时充分考虑了四个方向上高程差对汇流过程的影响,充分考虑了流域产汇流机理及水文过程的物理机制,降低了传统的平均分配方法带来的计算误差问题,计算过程简单、好操作、易推广。
在本说明书的描述中,参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不一定指的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
尽管上面已经示出和描述了本发明的实施例,可以理解的是,上述实施例是示例性的,不能理解为对本发明的限制,本领域的普通技术人员在不脱离本发明的原理和宗旨的情况下在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。本发明的范围由所附权利要求及其等同限定。
Claims (6)
1.一种基于四向高程差的净雨量转栅格格式方法,其特征在于,包括如下步骤:
步骤S1,计算每个栅格在四个方向的高程差之和,根据所述高程差之和设置权重;
步骤S2,根据计算得到的每个栅格的权重,计算不同栅格净雨量的分配比例;
步骤S3,根据计算得到的降雨量的分配比例,将流域的净雨量数据转栅格格式,分配到具体的每个栅格内。
2.如权利要求1所述的基于四向高程差的净雨量转栅格格式方法,其特征在于,在所述步骤S1中,采用数字高程模型DEM计算每个栅格在四个方向的高程差之和。
3.如权利要求1所述的基于四向高程差的净雨量转栅格格式方法,其特征在于,在所述步骤S1中,所述根据高程差之和计算权重,包括如下步骤:
当计算的高程差之和大于0时,则设置所述权重为高程差之和;
当计算的高程差之和小于0时,则设置所述权重为0。
6.如权利要求1所述的基于四向高程差的净雨量转栅格格式方法,其特征在于,在所述步骤S3中,所述将流域的净雨量数据转栅格格式,分配到具体的每个栅格内,包括:
WGi,j=WE×pi,j,
其中,WGi,j表示净雨量的栅格值;WE为某一时刻的流域净雨量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911016711.7A CN110686861B (zh) | 2019-10-24 | 2019-10-24 | 一种基于四向高程差的净雨量转栅格格式方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911016711.7A CN110686861B (zh) | 2019-10-24 | 2019-10-24 | 一种基于四向高程差的净雨量转栅格格式方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110686861A true CN110686861A (zh) | 2020-01-14 |
CN110686861B CN110686861B (zh) | 2020-09-01 |
Family
ID=69114120
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911016711.7A Active CN110686861B (zh) | 2019-10-24 | 2019-10-24 | 一种基于四向高程差的净雨量转栅格格式方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110686861B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112699336A (zh) * | 2021-01-07 | 2021-04-23 | 中国科学院地理科学与资源研究所 | 基于八向高程差的地表径流流域空间离散化方法及装置 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR100898617B1 (ko) * | 2008-11-20 | 2009-05-27 | 주식회사 범아엔지니어링 | LiDAR 측량값 및 MBES 측량값의 불규칙삼각망(TIN) 보간 성과 검증을 통한 육상부와 수심부가 함께 공존하는 지역의 수치표고모델을 구축하는 방법 |
CN102034001A (zh) * | 2010-12-16 | 2011-04-27 | 南京大学 | 一种以栅格为模拟单元的分布式水文模型设计方法 |
CN106777919A (zh) * | 2016-11-29 | 2017-05-31 | 南京信息工程大学 | 一种基于雨型的暴雨积涝积水动态演进过程预报计算方法 |
CN107194156A (zh) * | 2017-05-03 | 2017-09-22 | 南京信息工程大学 | 一种水量动态分配的有源积水扩散算法 |
CN109033599A (zh) * | 2018-07-18 | 2018-12-18 | 福州大学 | 一种基于随机森林的土壤侵蚀影响因子重要性分析方法 |
-
2019
- 2019-10-24 CN CN201911016711.7A patent/CN110686861B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR100898617B1 (ko) * | 2008-11-20 | 2009-05-27 | 주식회사 범아엔지니어링 | LiDAR 측량값 및 MBES 측량값의 불규칙삼각망(TIN) 보간 성과 검증을 통한 육상부와 수심부가 함께 공존하는 지역의 수치표고모델을 구축하는 방법 |
CN102034001A (zh) * | 2010-12-16 | 2011-04-27 | 南京大学 | 一种以栅格为模拟单元的分布式水文模型设计方法 |
CN106777919A (zh) * | 2016-11-29 | 2017-05-31 | 南京信息工程大学 | 一种基于雨型的暴雨积涝积水动态演进过程预报计算方法 |
CN107194156A (zh) * | 2017-05-03 | 2017-09-22 | 南京信息工程大学 | 一种水量动态分配的有源积水扩散算法 |
CN109033599A (zh) * | 2018-07-18 | 2018-12-18 | 福州大学 | 一种基于随机森林的土壤侵蚀影响因子重要性分析方法 |
Non-Patent Citations (3)
Title |
---|
张镀光等: "基于DEM的地形指数提取方法及应用", 《长江流域资源与环境》 * |
翟家瑞: "网格点计算面平均雨量的方法及其改进", 《人民黄河》 * |
袁定波等: "基于地理空间要素的雅砻江流域面雨量估算", 《水科学进展》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112699336A (zh) * | 2021-01-07 | 2021-04-23 | 中国科学院地理科学与资源研究所 | 基于八向高程差的地表径流流域空间离散化方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN110686861B (zh) | 2020-09-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Berhe et al. | MODSIM-based water allocation modeling of Awash River Basin, Ethiopia | |
Rodríguez‐Iturbe et al. | Energy dissipation, runoff production, and the three‐dimensional structure of river basins | |
Dümenil et al. | A rainfall–runoff scheme for use in the Hamburg climate model | |
Timbadiya et al. | A 1D–2D coupled hydrodynamic model for river flood prediction in a coastal urban floodplain | |
CN106202790B (zh) | 一种分布式河北模型构建方法及其应用 | |
CN108022047A (zh) | 一种海绵城市水文计算方法 | |
Wu et al. | Present and future of urban water balance in the rapidly urbanizing Heihe River Basin, Northwest China | |
CN108182543A (zh) | 一种精细化网格内涝水淹预报方法 | |
Chen et al. | Real-time error correction method combined with combination flood forecasting technique for improving the accuracy of flood forecasting | |
CN105678067A (zh) | 城市河涌防洪排涝动态预报控制方法及*** | |
CN105869070A (zh) | 一种跨流域梯级水电站群效益均衡的合作优化调度方法 | |
CN109101706B (zh) | 集总式水文模型和二维水动力模型的耦合方法 | |
Khan et al. | Back to the future: assessing the damage of 2004 Dhaka flood in the 2050 urban environment | |
CN110686861B (zh) | 一种基于四向高程差的净雨量转栅格格式方法 | |
CN110321611A (zh) | 一种多材料结构拓扑优化方法 | |
CN108763615A (zh) | 基于管网和道路双***swmm对城市涝积水深模拟方法 | |
CN115173445B (zh) | 一种城市配电网与流域水系网灵活性联动运行方法 | |
CN108897940B (zh) | 基于矩形网格的分布式水文模型和二维水动力模型的单向耦合方法 | |
CN107330086A (zh) | 一种提高无资料高海拔流域水文过程模拟精度的方法 | |
Fan et al. | Assessment of land cover resolution impact on flood modeling uncertainty | |
CN117236501A (zh) | 地下水资源管理与海水入侵防控方法及*** | |
CN108896473B (zh) | 耦合不同时空尺度模型的流域侵蚀产沙量预测方法 | |
CN109190213B (zh) | 基于土地利用类型的地表径流流域空间离散化方法 | |
CN110532500B (zh) | 一种区域大气水文模型中水库群参数化方案的构建方法 | |
CN115564238B (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 |