[發明專利]無資料地區水文模型參數率定方法、裝置和終端設備

有效

• 申請號: 202110695306.3
• 申請日: 2021-06-23
• 公開（公告）號: CN113254861B
• 公開（公告）日: 2021-09-24
• 發明（設計）人: 陳曉宏;于澤興;吳杰峰
• 申請（專利權）人: 中山大學
• 主分類號: G06F17/11
• 分類號: G06F17/11;G01C13/00;G01N21/359
• 代理公司: 北京集佳知識產權代理有限公司 11227
• 代理人: 陳嘉雯

1.一種無資料地區水文模型參數率定方法，其特征在于，所述方法包括：

獲取待研究區域率定期內的至少一個近紅外波段影像和根據水文模型隨機生成的至少一組初始參數；

計算各個所述近紅外波段影像中參考區域和淹沒區域的多像素比例指數，以得到一組多像素比例指數數據；

將每一組所述初始參數輸入所述水文模型生成每一組徑流數據；

計算所述多像素比例指數數據與每一組所述徑流數據的相關系數；

選擇所述相關系數排名在預設閾值內的所述徑流數據對應的初始參數，記為有效參數；

對所述有效參數采用水量平衡濾波器進行過濾，以得到過濾參數；

根據所述過濾參數得到所述水文模型的率定參數。

2.根據權利要求1所述的方法，其特征在于，所述計算各個所述近紅外波段影像中參考區域和淹沒區域的多像素比例指數，包括：

在每一個所述近紅外波段影像中選擇參考區域和淹沒區域；

分別計算所述參考區域和所述淹沒區域的反射率；

根據所述參考區域的反射率和所述淹沒區域的反射率計算所述多像素比例指數。

3.根據權利要求2所述的方法，其特征在于，所述分別計算所述參考區域和所述淹沒區域的反射率，包括：

分別獲取所述參考區域和所述淹沒區域中每一個像素的反射率；

根據參考區域中每一個像素的反射率來計算所述參考區域的反射率；

根據淹沒區域中每一個像素的反射率、所述淹沒區域的長度和所述淹沒區域的寬度來計算所述淹沒區域的反射率。

4.根據權利要求1所述的方法，其特征在于，所述計算所述多像素比例指數數據與每一組所述徑流數據的相關系數，包括：

計算所述多像素比例指數數據與每一組所述徑流數據的斯皮爾曼等級相關系數。

5.根據權利要求1所述的方法，其特征在于，所述對所述有效參數采用水量平衡濾波器進行過濾，以得到過濾參數，包括：

基于Budyko框架構建Budyko水量平衡濾波器；

對所述有效參數采用所述Budyko水量平衡濾波器進行過濾，以得到過濾參數。

6.根據權利要求1-5任一項所述的方法，其特征在于，所述水文模型包括CemaNeige-GR4J水文模型；所述獲取根據水文模型隨機生成的至少一組初始參數，包括：

采用超拉丁采樣方法隨機生成至少一組CemaNeige-GR4J水文模型的初始參數。

7.根據權利要求2所述的方法，其特征在于，所述根據所述參考區域的反射率和所述淹沒區域的反射率計算所述多像素比例指數，包括：

計算所述參考區域的反射率與所述淹沒區域的反射率的比值，記為所述多像素比例指數。

8.根據權利要求2所述的方法，其特征在于，所述在每一個所述近紅外波段影像中選擇參考區域和淹沒區域，包括：

在所述近紅外波段影像中選擇山坡區域或城市區域作為所述參考區域，其中所述參考區域包括至少一個像素；

在所述近紅外波段影像中選擇河流區域作為所述淹沒區域，其中所述淹沒區域的位置與水文站的位置之間的距離大于或等于預設距離，所述淹沒區域的長度在預設范圍內，所述淹沒區域的寬度小于或等于河道寬度最大值。