1.一种城市内涝积水深度计算方法,其特征在于,用于计算目标城市内的积水深度,方法包括如下步骤:
步骤1:基于目标城市的数字高程数据,将目标城市划分为I个栅格单元,获取各栅格单元的高程值、边长r;
分别使用预设的积水深度下限初始值、积水深度上限初始值、高程值最小的栅格单元的积水深度初始值初始化目标城市内积水深度的下限值d0、积水深度的上限值d1、目标城市内最小高程值所对应栅格单元的积水深度值d;
然后进入步骤2;
步骤2:获取目标城市的实际积水总量W;然后进入步骤3;
步骤3;按如下步骤3.1至步骤3.2,基于目标城市内最小高程值所对应栅格单元的积水深度值d计算目标城市的估算积水总量W1,根据公式C=|W1‑W|获取积水总量差值C;然后进入步骤4;
步骤3.1:分别针对各栅格单元i,i=1,2...I,依次执行步骤3.1.1至步骤3.1.2以获取各栅格单元i的积水深度d'i,进而获取各栅格单元的积水深度之和Q;然后进入步骤3.2;
步骤3.1.1:计算hmin+d‑hi的值;其中,hi为栅格单元i的高程值,hmin为各栅格单元的高程值中的最小高程值;
步骤3.1.2:判断hmin+d‑hi的值是否大于或等于0,是则将hmin+d‑hi的值作为栅格单元i的积水深度d'i;否则取栅格单元i的积水深度d'i=0;
2
步骤3.2:根据公式W1=r×Q获取目标城市的估算积水总量W1;其中,r为各栅格单元的边长;
步骤4:判断积水总量差值C是否小于积水总量差阈值y,是则进入步骤10,否则进入步骤5;
步骤5:判断d1是否等于步骤1中预设的积水深度上限初始值,是则执行步骤6;否则进入步骤8;
步骤6:判断W1是否小于W,是则将d所对应的值赋值给d0,然后用d所对应的值的2倍对d所对应的值进行更新,返回步骤3;否则进入步骤7;
步骤7:将d所对应的值赋值给d1,返回步骤3;
步骤8:判断W1是否小于W,是则将d所对应的值赋值给d0,然后应用 对d所对应的值进行更新,返回步骤3;否则进入步骤9;
步骤9:将d所对应的值赋值给d1,然后应用 对d所对应的值进行更新,返回步骤3;
步骤10:将目标城市内高程值最小的栅格单元的积水深度值d作为目标城市内的积水深度。
2.根据权利要求1所述的城市内涝积水深度计算方法,其特征在于,在步骤10之后,方法还包括:
步骤11:根据公式hmin+d‑hi获取各栅格单元i的积水深度d'i,其中,hi为栅格单元i的高程值,hmin为各栅格单元的高程值中的最小高程值。
3.根据权利要求1所述的城市内涝积水深度计算方法,其特征在于,在步骤2中,根据如下公式:
2
W=r×I×m
获取目标城市的实际积水总量W;
其中,r为单个栅格单元的边长,m为目标城市在预设时间段内的降水量。
4.根据权利要求1所述的城市内涝积水深度计算方法,其特征在于,在步骤1中,预设的积水深度下限初始值为0mm,预设的积水深度上限初始值为0mm,预设的高程值最小的栅格单元的积水深度为100mm。
5.根据权利要求1所述的城市内涝积水深度计算方法,其特征在于,在步骤4中,所述积水总量差阈值y的取值为目标城市的实际积水总量W的十万分之一。