[发明专利]一种基于DEM的汇水区降雨积水深度的计算方法

摘要

本发明公开了一种基于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，其中0h，则计算t=t+hi；令T=T-t，返回计算积水高程为hmax时的积水量R步骤。该方法计算速度快，便于应用GIS（地理信息系统）来实现。

主权项

1.一种基于DEM的汇水区降雨积水深度的计算方法，其特征在于，包括以下步骤：     步骤1、遍历DEM栅格，查找最大高程值h_max，令h=h_max；步骤2、计算所有栅格高程值的和T：T=，h_i为单元栅格高程值；步骤3、计算积水高程为h_max时的积水量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、计算区域单元积水深度H_i：H_i=h-h_i，结束计算；步骤7、令h=h*K，其中0<K<1；；步骤8、设t=0 ；遍历DEM，从1到n计算，n为单元编号，如果 h_i>h，则计算t=t+ h_i；；步骤9、令T=T-t，返回步骤3。