[發(fā)明專利]基于格蘭杰檢驗分析氣候參數之間的因果關系的方法

在審

• 申請?zhí)枺?/td>: 202211128534.3
• 申請日: 2022-09-16
• 公開（公告）號: CN115526469A
• 公開（公告）日: 2022-12-27
• 發(fā)明（設計）人: 劉楊曉月;楊雅萍;陳曉娜
• 申請（專利權）人: 中國科學院地理科學與資源研究所
• 主分類號: G06Q10/06
• 分類號: G06Q10/06;G06Q50/26
• 代理公司: 廣州駿思知識產權代理有限公司 44425
• 代理人: 葉瓊園

1.一種基于格蘭杰檢驗分析氣候參數之間的因果關系的方法，其特征在于，包括以下步驟：

獲取目標區(qū)域的目標時間段內的遙感氣候數據，構建所述目標區(qū)域的遙感氣候時間序列，其中，所述遙感氣候數據包括若干個不同類型的氣候參數數據，所述遙感氣候時間序列包括若干個氣候參數時間序列；

獲取第一類型氣候標識以及第二類型氣候標識，所述第一類型氣候標識以及第二類型氣候標識用于指示兩個待校驗的不同類型的氣候參數數據；從所述目標區(qū)域的遙感氣候時間序列中，提取與所述第一類型氣候標識對應的氣候參數時間序列以及第二類型氣候標識對應的氣候參數時間序列；

獲得所述第一類型氣候標識對應的氣候參數時間序列與所述第二類型氣候標識對應的氣候參數時間序列之間的時間滯后參數，所述時間滯后參數為所述氣候參數時間序列中氣候參數數據的時間單位的索引；根據所述時間滯后參數，對所述第一類型氣候標識對應的氣候參數時間序列以及第二類型氣候標識對應的氣候參數時間序列進行錯位處理，獲得錯位處理后的第一類型氣候標識對應的氣候參數時間序列以及第二類型氣候標識對應的氣候參數時間序列；

對所述錯位處理后的第一類型氣候標識對應的氣候參數時間序列以及第二類型氣候標識對應的氣候參數時間序列進行格蘭杰因果關系校驗，獲得因果關系校驗結果，根據所述因果關系校驗結果，獲得所述第一類型氣候標識與第二類型氣候標識對應的所述兩個待校驗的不同類型的氣候參數數據之間的因果關系。

2.根據權利要求1所述的基于格蘭杰檢驗分析氣候參數之間的因果關系的方法，其特征在于：所述氣候參數數據包括土壤水分數據、溫度數據、降水數據以及植被指數數據。

3.根據權利要求2所述的基于格蘭杰檢驗分析氣候參數之間的因果關系的方法，其特征在于：所述氣候參數時間序列包括土壤水分時間序列、溫度時間序列、降水時間序列以及植被指數時間序列。

4.根據權利要求3所述的基于格蘭杰檢驗分析氣候參數之間的因果關系的方法，其特征在于：所述氣候參數時間序列包括若干個按照時間排序的時間單元對應的氣候參數數據。

5.根據權利要求4所述的基于格蘭杰檢驗分析氣候參數之間的因果關系的方法，其特征在于，所述獲得所述第一類型氣候標識對應的氣候參數時間序列與所述第二類型氣候標識對應的氣候參數時間序列之間的時間滯后參數，包括步驟：

根據所述第一類型氣候標識對應的氣候參數時間序列、第二類型氣候標識對應的氣候參數時間序列以及預設的互相關系數計算算法，獲得所述第一類型氣候標識與第二類型氣候標識，在不同時間滯后單元的互相關系數序列，其中，所述互相關系數序列包括若干個互相關系數，其中，所述互相關系數計算算法為：

Cross correlation＝∑_nX[n+k]×Y^*[n]

式中，Cross correlation為所述互相關系數，X[n]為所述第一類型氣候標識對應的氣候參數時間序列中第n個時間單元對應的氣候參數數據，k為時間滯后單元，Y^*[n]為所述第二類型氣候標識對應的氣候參數時間序列中第n個時間單元對應的氣候參數數據的復共軛；

從所述互相關系數序列中提取最大互相關系數，獲取所述最大互相關系數對應的時間滯后參數，作為所述第一類型氣候標識對應的氣候參數時間序列與所述第二類型氣候標識對應的氣候參數時間序列之間的時間滯后參數。