[發明專利]一種城市內澇積水深度計算方法

有效

• 申請號: 202110330233.8
• 申請日: 2021-03-29
• 公開（公告）號: CN112800631B
• 公開（公告）日: 2021-07-13
• 發明（設計）人: 薛豐昌;韓宇
• 申請（專利權）人: 南京信息工程大學;南京宸象空間信息技術有限公司
• 主分類號: G06F30/20
• 分類號: G06F30/20;G06T17/05
• 代理公司: 南京經緯專利商標代理有限公司 32200
• 代理人: 田凌濤

1.一種城市內澇積水深度計算方法，其特征在于，用于計算目標城市內的積水深度，方法包括如下步驟：

步驟1：基于目標城市的數字高程數據，將目標城市劃分為I個柵格單元，獲取各柵格單元的高程值、邊長r；

分別使用預設的積水深度下限初始值、積水深度上限初始值、高程值最小的柵格單元的積水深度初始值初始化目標城市內積水深度的下限值d₀、積水深度的上限值d₁、目標城市內最小高程值所對應柵格單元的積水深度值d；

然后進入步驟2；

步驟2：獲取目標城市的實際積水總量W；然后進入步驟3；

步驟3；按如下步驟3.1至步驟3.2，基于目標城市內最小高程值所對應柵格單元的積水深度值d計算目標城市的估算積水總量W₁，根據公式C＝|W₁-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：計算h_min+d-h_i的值；其中，h_i為柵格單元i的高程值，h_min為各柵格單元的高程值中的最小高程值；

步驟3.1.2：判斷h_min+d-h_i的值是否大于或等于0，是則將h_min+d-h_i的值作為柵格單元i的積水深度d'_i；否則取柵格單元i的積水深度d'_i＝0；

步驟3.2：根據公式W₁＝r²×Q獲取目標城市的估算積水總量W₁；其中，r為各柵格單元的邊長；

步驟4：判斷積水總量差值C是否小于積水總量差閾值y，是則進入步驟10，否則進入步驟5；

步驟5：判斷d₁是否等于步驟1中預設的積水深度上限初始值，是則執行步驟6；否則進入步驟8；

步驟6：判斷W₁是否小于W，是則將d所對應的值賦值給d₀，然后用d所對應的值的2倍對d所對應的值進行更新，返回步驟3；否則進入步驟7；

步驟7：將d所對應的值賦值給d₁，返回步驟3；

步驟8：判斷W₁是否小于W，是則將d所對應的值賦值給d₀，然后應用對d所對應的值進行更新，返回步驟3；否則進入步驟9；

步驟9：將d所對應的值賦值給d₁，然后應用對d所對應的值進行更新，返回步驟3；

步驟10：將目標城市內高程值最小的柵格單元的積水深度值d作為目標城市內的積水深度。

2.根據權利要求1所述的城市內澇積水深度計算方法，其特征在于，在步驟10之后，方法還包括：

步驟11：根據公式h_min+d-h_i獲取各柵格單元i的積水深度d'_i，其中，h_i為柵格單元i的高程值，h_min為各柵格單元的高程值中的最小高程值。

3.根據權利要求1所述的城市內澇積水深度計算方法，其特征在于，在步驟2中，根據如下公式：

W＝r²×I×m

獲取目標城市的實際積水總量W；

其中，r為單個柵格單元的邊長，m為目標城市在預設時間段內的降水量。

4.根據權利要求1所述的城市內澇積水深度計算方法，其特征在于，在步驟1中，預設的積水深度下限初始值為0mm，預設的積水深度上限初始值為0mm，預設的高程值最小的柵格單元的積水深度為100mm。

5.根據權利要求1所述的城市內澇積水深度計算方法，其特征在于，在步驟4中，所述積水總量差閾值y的取值為目標城市的實際積水總量W的十萬分之一。