CN102902893A - 一种基于dem的汇水区降雨积水深度的计算方法 - Google Patents
一种基于dem的汇水区降雨积水深度的计算方法 Download PDFInfo
- Publication number
- CN102902893A CN102902893A CN2012104200436A CN201210420043A CN102902893A CN 102902893 A CN102902893 A CN 102902893A CN 2012104200436 A CN2012104200436 A CN 2012104200436A CN 201210420043 A CN201210420043 A CN 201210420043A CN 102902893 A CN102902893 A CN 102902893A
- Authority
- CN
- China
- Prior art keywords
- dem
- calculating
- grid
- ponding
- water
- 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
-
- 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
- Sewage (AREA)
Abstract
本发明公开了一种基于DEM的汇水区降雨积水深度的计算方法,属于城市内涝监测预警技术领域。该方法包括以下步骤:遍历DEM栅格,查找最大高程值hmax,令h=hmax;计算所有栅格高程值的和T,hi为单元栅格高程值;计算积水高程为hmax时的积水量R:R=(n*h-T)*s,n为DEM栅格数量,s为DEM栅格单元面积;求取积水量R区域降水汇流水量Q之间的差值V;判断是否V≤e,其中0≤e≤1,设定的常数,若V≤e,则执行下一步;若V>e,则执行下下一步;计算区域单元积水深度Hi:Hi=h-hi,结束计算;令h=h*K,其中0<K<1;设t=0;遍历DEM,从1到n计算,n为单元编号,如果hi>h,则计算t=t+hi;令T=T-t,返回计算积水高程为hmax时的积水量R步骤。该方法计算速度快,便于应用GIS(地理信息系统)来实现。
Description
技术领域
本发明公开了一种基于DEM的汇水区降雨积水深度的计算方法,属于城市内涝监测预警技术领域。
背景技术
城市内涝监测预警涉及到产流计算、汇流计算、排水计算、汇水区积水深度计算,其中汇水区积水深度是内涝监测预警的重要参数。汇水区是指依据地形、水流方向计算形成的水文计算单元,一个区域往往由若干个汇水区组成。汇水区的总水量由本区降雨径流量和经管网运移来的其他汇水区的排水量组成。汇水区的积水量为汇水区的总水量减去本汇水区的排水量。
DEM(Digital Elevation Model,数字高程模型)是一定范围内规则格网点的平面坐标(X,Y)及其高程(Z)的数据集,它主要是描述区域地貌形态的空间分布。在本发明中,指利用汇水区高程点经过插值形成的汇水区内每个规则网格都包含一个高程数据的栅格数据集,其中网格边长确定后即确定了每个网格面积。
汇水区的积水量在汇水区内要进行淹没分配,淹没分配后的水位高程减去相应点地形高程为该点积水深度。
传统汇水区积水深度计算有两种形式,一种是根据汇水区面积和积水量估算积水深度,这种方法精度低,不能自动化计算;另一种方法是采用二分法进行计算,计算复杂,计算量大,计算速度慢。
发明内容
本发明所要解决的技术问题是针对上述现有技术的不足,而提供一种算法简单计算速度快且精度较高的基于DEM的汇水区降雨积水深度的计算方法。
为解决上述技术问题,本发明采用的技术方案是:
基于DEM的汇水区降雨积水深度的计算方法,其特征在于,包括以下步骤: 步骤1、遍历DEM栅格,查找最大高程值hmax,令h=hmax;
步骤3、计算积水高程为hmax时的积水量R:R=(n*h-T)*s,n为DEM栅格数量,s为DEM栅格单元面积;
步骤4、求取积水量R区域降水汇流水量Q之间的差值V: V=R-Q;
步骤5、判断是否V≤e,其中0≤e≤1,设定的常数,若V≤e,则执行步骤6;若V>e,则执行步骤7;
步骤6、计算区域单元积水深度Hi:Hi=h-hi,结束计算;
步骤7、令h=h*K,其中0<K<1;
V>e说明积水高程应小于hmax,该步骤对最大高程值h(h=hmax)缩小,进而测通过步骤3至步骤5测算h是否就是积水高程;
步骤8、设t=0 ;遍历DEM,从1到n计算,n为单元编号,如果 hi>h,则计算t=t+ hi;
若h是积水高程,则高程大于h的栅格将没有积水,这里t为没有积水的栅格高程值总和;
步骤9、令T=T-t,返回步骤3。
根据步骤2,T为所有栅格高程值之和,因此应将没有积水的栅格高程值之和t去除。
发明效果
1、 本发明计算速度快,从步骤1至步骤9可看出,基本计算为栅格高程值的求和,实现简单,不需进行解算复杂方程;同时可根据不同精度要求调整步骤5的参数e和步骤7的参数k,从而可进一步提高计算速度;
2、 本发明便于利用GIS(地理信息新系统)来实现,栅格计算是GIS的基本功能,步骤1至步骤3、步骤7、步骤8等都可以在GIS属性表中利用GIS统计计算功能直接批量计算得到;
3、 计算结果可基于GIS实现可视化表达,计算得到的各栅格的积水高程可以地图化形式标注于每个地理栅格,也可以对不同积水深度的栅格以不同颜色进行渲染。
附图说明
图1:一种基于DEM的汇水区降雨积水深度的计算方法流程图;
图2:汇水区示意图;
图3:单元积水深度分布图。
实施示例:
设某城市街区一汇水区如下图2所示:
图2示例区域进行了正方形网格划分,网格边长10米,每个单元面积为100m2,整个汇水区面积为7840 m2。网格中所注数字为汇水区域划分为网格后,每个网格单元的高程值,单位为米。该汇水区综合径流系数为0.21,某次降雨强度为每小时56mm,计算得示例区域降水产流水量为950 m3(利用综合径流系数法计算降雨地面产流量,成熟方法,不属于本发明范畴)。
下面结合附图对本发明的计算方法进行详细说明。
如图1所示,包括以下步骤:
步骤1、遍历DEM,查找最大高程值令h=hmax;
此步骤利用遍历算法(成熟算法,不属于本发明范畴),在DEM中查找出汇水区DEM中最大高程值hmax,根据图2所示 DEM,本示例中最大高程值为5.49米;
根据图2所示DEM计算出所有栅格高程值的和,本示例中所有栅格高程值的和为339.78米;
步骤3、计算R=(n*h-T)*s,;
h=5.49米,n=78, s=100m2,计算得R=8844m3;
步骤4、计算V=R-Q,Q为区域降水汇流水量
示例区域综合径流系数为0.21,降雨强度为每小时56mm,计算得示例区域降水产流水量Q为950 m3,计算得到v=7894;
步骤5、判断是否V≤e,若“V≤e”则执行步骤6;若“V>e”则执行步骤7;
这里取e=0.1 m3,判断得到“V>e”,执行步骤7-9;
步骤6、Hi=h-hi,结束计算。
Hi为区域单元积水深度,h为前述积水水位高程,hi为每个栅格单元高程值;
步骤7、计算h=h*K(0<K<1);
K取0.9,;
经过步骤8、步骤9的循环计算得到积水水位高程h为4.22米。
返回到步骤6计算单元积水深度,结束计算。
每个栅格单元高程值与积水水位高程相减得到汇水区每个栅格单元最终积水深度Hi。图3为单元积水深度分布图,表中负值为单元积水深度,正值为所在单元地面高程高于积水水位值,相应单元无积水。
Claims (1)
1.一种基于DEM的汇水区降雨积水深度的计算方法,其特征在于,包括以下步骤: 步骤1、遍历DEM栅格,查找最大高程值hmax,令h=hmax;
步骤3、计算积水高程为hmax时的积水量R:R=(n*h-T)*s,n为DEM栅格数量,s为DEM栅格单元面积;
步骤4、求取积水量R区域降水汇流水量Q之间的差值V: V=R-Q;
步骤5、判断是否V≤e,其中0≤e≤1,设定的常数,若V≤e,则执行步骤6;若V>e,则执行步骤7;
步骤6、计算区域单元积水深度Hi:Hi=h-hi,结束计算;
步骤7、令h=h*K,其中0<K<1;;
步骤8、设t=0 ;遍历DEM,从1到n计算,n为单元编号,如果 hi>h,则计算t=t+ hi;;
步骤9、令T=T-t,返回步骤3。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210420043.6A CN102902893B (zh) | 2012-10-29 | 2012-10-29 | 一种基于dem的汇水区降雨积水深度的计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210420043.6A CN102902893B (zh) | 2012-10-29 | 2012-10-29 | 一种基于dem的汇水区降雨积水深度的计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102902893A true CN102902893A (zh) | 2013-01-30 |
CN102902893B CN102902893B (zh) | 2015-03-11 |
Family
ID=47575121
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210420043.6A Expired - Fee Related CN102902893B (zh) | 2012-10-29 | 2012-10-29 | 一种基于dem的汇水区降雨积水深度的计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102902893B (zh) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104268409A (zh) * | 2014-09-28 | 2015-01-07 | 中国矿业大学 | 一种地面集水时间的计算方法 |
CN105319391A (zh) * | 2015-11-03 | 2016-02-10 | 李爽 | 非线性建筑的雨水流向测量方法及装置 |
CN105760666A (zh) * | 2016-02-04 | 2016-07-13 | 南京信息工程大学 | 一种基于dem的汇水区临界雨量计算方法 |
CN106056851A (zh) * | 2016-05-13 | 2016-10-26 | 国网浙江省电力公司台州供电公司 | 电网设施暴雨预警方法 |
CN106780514A (zh) * | 2016-12-28 | 2017-05-31 | 南京信息工程大学 | 基于监控视频影像的暴雨积涝区积水深度的计算方法 |
CN109712237A (zh) * | 2018-12-21 | 2019-05-03 | 西安恒歌数码科技有限责任公司 | 一种基于高程数据与气象数据的混合设色方法 |
CN111259527A (zh) * | 2020-01-10 | 2020-06-09 | 深圳震有科技股份有限公司 | 基于gis系统的淹没分析方法、智能终端及存储介质 |
CN112329582A (zh) * | 2020-10-29 | 2021-02-05 | 西南大学 | 大数据分析与机理模型协同的道路积水深度监测方法及系统 |
CN112325861A (zh) * | 2020-11-10 | 2021-02-05 | 厦门点云数据科技有限公司 | 一种积水监测和计算方法及终端 |
CN112598910A (zh) * | 2020-11-16 | 2021-04-02 | 河北工程大学 | 一种街道积水预警的引导方法和设备 |
CN112800631A (zh) * | 2021-03-29 | 2021-05-14 | 南京信息工程大学 | 一种城市内涝积水深度计算方法 |
CN112819957A (zh) * | 2021-01-25 | 2021-05-18 | 杭州鲁尔物联科技有限公司 | 一种基于gis技术的城市洼地深度计算方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2001039154A1 (en) * | 1999-11-19 | 2001-05-31 | Deluxe Innovations Pty Ltd | A method for providing up-to-date information on road flooding |
CN101188022A (zh) * | 2007-12-20 | 2008-05-28 | 浙江大学 | 一种面向大城市灾害展示的洪涝淹没分析方法 |
CN101694680A (zh) * | 2009-09-28 | 2010-04-14 | 深圳先进技术研究院 | 城市暴雨水灾害模拟预测方法 |
-
2012
- 2012-10-29 CN CN201210420043.6A patent/CN102902893B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2001039154A1 (en) * | 1999-11-19 | 2001-05-31 | Deluxe Innovations Pty Ltd | A method for providing up-to-date information on road flooding |
CN101188022A (zh) * | 2007-12-20 | 2008-05-28 | 浙江大学 | 一种面向大城市灾害展示的洪涝淹没分析方法 |
CN101694680A (zh) * | 2009-09-28 | 2010-04-14 | 深圳先进技术研究院 | 城市暴雨水灾害模拟预测方法 |
Non-Patent Citations (1)
Title |
---|
XIANFU ZHAO等: "GIS-Based Spatial and Temporal Analysis of Regional Water-logging Confluence", 《2011 FOURTH INTERNATIONAL JOINT CONFERENCE ON COMPUTATIONAL SCIENCES AND OPTIMIZATION》 * |
Cited By (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104268409A (zh) * | 2014-09-28 | 2015-01-07 | 中国矿业大学 | 一种地面集水时间的计算方法 |
CN105319391A (zh) * | 2015-11-03 | 2016-02-10 | 李爽 | 非线性建筑的雨水流向测量方法及装置 |
CN105319391B (zh) * | 2015-11-03 | 2018-05-18 | 李爽 | 非线性建筑的雨水流向测量方法及装置 |
CN105760666A (zh) * | 2016-02-04 | 2016-07-13 | 南京信息工程大学 | 一种基于dem的汇水区临界雨量计算方法 |
CN106056851A (zh) * | 2016-05-13 | 2016-10-26 | 国网浙江省电力公司台州供电公司 | 电网设施暴雨预警方法 |
CN106056851B (zh) * | 2016-05-13 | 2019-03-22 | 国网浙江省电力公司台州供电公司 | 电网设施暴雨预警方法 |
CN106780514A (zh) * | 2016-12-28 | 2017-05-31 | 南京信息工程大学 | 基于监控视频影像的暴雨积涝区积水深度的计算方法 |
CN106780514B (zh) * | 2016-12-28 | 2019-11-19 | 南京信息工程大学 | 基于监控视频影像的暴雨积涝区积水深度的计算方法 |
CN109712237B (zh) * | 2018-12-21 | 2023-03-31 | 西安恒歌数码科技有限责任公司 | 一种基于高程数据与气象数据的混合设色方法 |
CN109712237A (zh) * | 2018-12-21 | 2019-05-03 | 西安恒歌数码科技有限责任公司 | 一种基于高程数据与气象数据的混合设色方法 |
CN111259527A (zh) * | 2020-01-10 | 2020-06-09 | 深圳震有科技股份有限公司 | 基于gis系统的淹没分析方法、智能终端及存储介质 |
CN112329582A (zh) * | 2020-10-29 | 2021-02-05 | 西南大学 | 大数据分析与机理模型协同的道路积水深度监测方法及系统 |
CN112329582B (zh) * | 2020-10-29 | 2022-08-02 | 西南大学 | 大数据分析与机理模型协同的道路积水深度监测方法及系统 |
CN112325861B (zh) * | 2020-11-10 | 2022-05-13 | 厦门点云数据科技有限公司 | 一种积水监测和计算方法及终端 |
CN112325861A (zh) * | 2020-11-10 | 2021-02-05 | 厦门点云数据科技有限公司 | 一种积水监测和计算方法及终端 |
CN112598910A (zh) * | 2020-11-16 | 2021-04-02 | 河北工程大学 | 一种街道积水预警的引导方法和设备 |
CN112819957A (zh) * | 2021-01-25 | 2021-05-18 | 杭州鲁尔物联科技有限公司 | 一种基于gis技术的城市洼地深度计算方法 |
CN112819957B (zh) * | 2021-01-25 | 2024-02-20 | 杭州鲁尔物联科技有限公司 | 一种基于gis技术的城市洼地深度计算方法 |
CN112800631A (zh) * | 2021-03-29 | 2021-05-14 | 南京信息工程大学 | 一种城市内涝积水深度计算方法 |
CN112800631B (zh) * | 2021-03-29 | 2021-07-13 | 南京信息工程大学 | 一种城市内涝积水深度计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102902893B (zh) | 2015-03-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102902893B (zh) | 一种基于dem的汇水区降雨积水深度的计算方法 | |
CN106884405B (zh) | 一种无资料地区溃堤型山洪灾害分析评价方法 | |
CN109815305B (zh) | 一种无资料地区场次洪水径流过程反演的方法 | |
CN104851360B (zh) | 一种洪水风险图的生成方法和系统 | |
CN111339711B (zh) | 一种小流域设计洪水推求方法 | |
CN111507375B (zh) | 一种城市内涝风险快速评估方法及系统 | |
CN104298841A (zh) | 一种基于历史数据的洪水预报方法和系统 | |
CN105760666A (zh) | 一种基于dem的汇水区临界雨量计算方法 | |
CN103886135B (zh) | 基于二维非恒定流数值模型的电力工程选址方法 | |
CN111538798B (zh) | 一种顾及dsm和dlg的城市汇水区精细化提取方法 | |
CN104462774A (zh) | 基于水箱模型的城市道路及低洼地区积水预报方法 | |
Herzog et al. | Least-cost kernel density estimation and interpolation-based density analysis applied to survey data | |
CN103293285A (zh) | 一种在流域或区域尺度上的土壤侵蚀测定方法 | |
Chen et al. | Three-dimensional simulation of regional urban waterlogging based on high-precision DEM model | |
Irvem | Application of GIS to determine storage volume and surface area of reservoirs: the case study of Buyuk Karacay dam | |
Das et al. | Study on urban land suitability assessment using remote sensing and GIS: a case study of Khairagarh in Chhattisgarh | |
CN110796340B (zh) | 一种宏观尺度的小流域综合治理优先级评估方法 | |
CN103559411A (zh) | 一种基于数字高程模型的洪水演进数值模拟计算修正方法 | |
CN104750936A (zh) | 一种基于gis的城市水文径流计算方法及系统 | |
CN105279317A (zh) | 一种基于dem的平地河网水流方向估算方法 | |
CN104123671B (zh) | 气象参证站分区方法 | |
CN106096129A (zh) | 一种基于山地汇水计算的山脚水面规模分析方法 | |
CN106546958A (zh) | 一种优化的雷达数据同化方法 | |
CN110109194A (zh) | 基于栅格陆面模型与向量化流径的河道栅格入流计算方法 | |
Xiao et al. | Impact of digital elevation model resolution on stream network parameters |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150311 Termination date: 20171029 |