本申請涉及災害誘發因子探測，特別是涉及一種時空尺度滑坡災害誘發因子探測方法及相關裝置。

背景技術：

1、近年發展起來的合成孔徑雷達干涉(interferometric?syntheticapertureradar，insar)技術能夠全天時、全天候、高精度對大范圍研究區域實施連續監測，在許多領域中表現出巨大的優勢和價值，但是在后期對監測結果的解譯、分析過程中存在災害誘發因子(也可以稱為致災因子)探測困難的問題，尤其是在單一los向(line-of-sight，雷達視線方向)形變結果下。隨著insar技術應用場景越來越廣泛，所監測研究區域的地理、地形情況也越來越復雜，單一los向形變結果在部分場景下無法滿足災害誘發因子探測的要求，需要反演地表三維形變以獲得更加準確的災害信息，因此，需要一種顧及三維形變的定量分析來完成災害誘發因子探測的方法。

2、國內外研究學者對災害誘發因子探測做了大量研究，開展了包括統計指數(statistical?index，si)、邏輯回歸(logistic?regression，lr)模型、現場地質觀測、隨機森林模型等不同的滑坡災害誘發因子探測方法的研究，但結果不是特別理想，且操作復雜度較高。地理探測器是探測空間異質性的一種新的空間統計學方法，對于空間尺度下因子分異性可以表現出良好的量化效果，并能夠準確判斷各因子兩兩之間線性與非線性相互作用，但是地理探測器難以量化時間尺度下潛在災害誘發因子影響程度，需要補充時間尺度下災害誘發因子探測方法。

技術實現思路

1、本申請的目的是提供一種時空尺度滑坡災害誘發因子探測方法及相關裝置，可以解決單一los向形變結果在部分場景下無法滿足災害誘發因子探測要求以及難以完成時間尺度下災害誘發因子探測的問題，提高滑坡災害誘發因子探測精度。

2、為實現上述目的，本申請提供了如下方案：

3、第一方面，本申請提供了一種時空尺度滑坡災害誘發因子探測方法，所述時空尺度滑坡災害誘發因子探測方法包括：

4、獲取研究區域的sar影像集、dem數據、多源遙感數據集和gacos大氣數據集；所述sar影像集包括預設時間段中每一第一時刻的sar影像，所述多源遙感數據集包括預設時間段中每一第二時刻的多源遙感數據，所述gacos大氣數據集包括預設時間段中每一第一時刻的gacos大氣數據；所述預設時間段包括衛星升軌過程和衛星降軌過程；

5、基于sar影像集、dem數據、多源遙感數據集和gacos大氣數據集，計算得到研究區域的los向形變數據集；所述los向形變數據集包括預設時間段中每一第一時刻的los向形變數據；

6、以los向形變數據集作為輸入，結合三維解算模型和spfm模型計算得到研究區域的三維形變速率數據集；所述三維形變速率數據集包括預設時間段中每一第一時刻的三維形變速率數據，所述三維形變速率數據包括研究區域中每一位置點的三維形變速率的取值，所述三維形變速率包括垂向形變速率、南北向形變速率和東西向形變速率；

7、基于三維形變速率數據集，利用相關系數法在時間尺度下從第一類潛在影響因子中探測第一影響因子，利用地理探測器在空間尺度下從第二類潛在影響因子中探測第二影響因子，并以第一影響因子和第二影響因子作為災害誘發因子；所述第一類潛在影響因子包括降水和溫度，所述第二類潛在影響因子包括dem數據、坡度、坡向、曲率、ndvi、土地覆蓋類型、地表溫度、水系和斷層中的至少兩項。

8、第二方面，本申請提供了一種計算機設備，包括：存儲器、處理器以及存儲在存儲器上并可在處理器上運行的計算機程序，所述處理器執行所述計算機程序以實現上述的時空尺度滑坡災害誘發因子探測方法。

9、第三方面，本申請提供了一種計算機可讀存儲介質，其上存儲有計算機程序，該計算機程序被處理器執行時實現上述的時空尺度滑坡災害誘發因子探測方法。

10、第四方面，本申請提供了一種計算機程序產品，包括計算機程序，該計算機程序被處理器執行時實現上述的時空尺度滑坡災害誘發因子探測方法。

11、根據本申請提供的具體實施例，本申請具有以下技術效果：

12、本申請提供了一種時空尺度滑坡災害誘發因子探測方法及相關裝置，基于研究區域的sar影像集、dem數據、多源遙感數據集和gacos大氣數據集，計算得到研究區域的los向形變數據集，以los向形變數據集作為輸入，結合三維解算模型和spfm模型計算得到研究區域的三維形變速率數據集，后續利用三維形變速率數據集來完成滑坡災害誘發因子探測，通過引入三維形變速率，能夠解決單一los向形變結果在部分場景下無法滿足災害誘發因子探測要求的問題。且在利用三維形變速率數據集來完成滑坡災害誘發因子探測時，基于三維形變速率數據集，利用相關系數法在時間尺度下從第一類潛在影響因子中探測第一影響因子，利用地理探測器在空間尺度下從第二類潛在影響因子中探測第二影響因子，并以第一影響因子和第二影響因子作為災害誘發因子，通過引入時間尺度，能夠融合時空尺度來完成滑坡災害誘發因子探測，解決難以完成時間尺度下災害誘發因子探測的問題。本申請能夠提高滑坡災害誘發因子探測精度。

技術特征：

1.一種時空尺度滑坡災害誘發因子探測方法，其特征在于，所述時空尺度滑坡災害誘發因子探測方法包括：

2.根據權利要求1所述的時空尺度滑坡災害誘發因子探測方法，其特征在于，基于sar影像集、dem數據、多源遙感數據集和gacos大氣數據集，計算得到研究區域的los向形變數據集，具體包括：

3.根據權利要求1所述的時空尺度滑坡災害誘發因子探測方法，其特征在于，三維解算模型為：

4.根據權利要求1所述的時空尺度滑坡災害誘發因子探測方法，其特征在于，基于三維形變速率數據集，利用相關系數法在時間尺度下從第一類潛在影響因子中探測第一影響因子，利用地理探測器在空間尺度下從第二類潛在影響因子中探測第二影響因子，具體包括：

5.根據權利要求4所述的時空尺度滑坡災害誘發因子探測方法，其特征在于，以每一個第一類潛在影響因子的因子數據集和垂向形變速率數據集作為輸入，利用相關系數法在時間尺度下從第一類潛在影響因子中探測第一影響因子，具體包括：

6.根據權利要求4所述的時空尺度滑坡災害誘發因子探測方法，其特征在于，獲取研究區域的每一個第一類潛在影響因子的因子數據集和每一個第二類潛在影響因子的取值，具體包括：

7.根據權利要求2所述的時空尺度滑坡災害誘發因子探測方法，其特征在于，利用轉換后dem數據計算得到起伏相位和平地相位，并在每一個差分干涉圖中去除起伏相位和平地相位，得到多個初步去除后差分干涉圖，具體包括：

8.一種計算機設備，包括：存儲器、處理器以及存儲在存儲器上并可在處理器上運行的計算機程序，其特征在于，所述處理器執行所述計算機程序以實現權利要求1-7中任一項所述的時空尺度滑坡災害誘發因子探測方法。

9.一種計算機可讀存儲介質，其上存儲有計算機程序，其特征在于，該計算機程序被處理器執行時實現權利要求1-7中任一項所述的時空尺度滑坡災害誘發因子探測方法。

10.一種計算機程序產品，包括計算機程序，其特征在于，該計算機程序被處理器執行時實現權利要求1-7中任一項所述的時空尺度滑坡災害誘發因子探測方法。

技術總結
本申請公開了一種時空尺度滑坡災害誘發因子探測方法及相關裝置，涉及災害誘發因子探測技術領域，該方法包括：基于研究區域的SAR影像集、DEM數據、多源遙感數據集和GACOS大氣數據集，計算研究區域的LOS向形變數據集，以LOS向形變數據集作為輸入，結合三維解算模型和SPFM模型計算研究區域的三維形變速率數據集，基于三維形變速率數據集，利用相關系數法在時間尺度下從第一類潛在影響因子中探測第一影響因子，利用地理探測器在空間尺度下從第二類潛在影響因子中探測第二影響因子，并以第一影響因子和第二影響因子作為災害誘發因子。本申請能夠提高滑坡災害誘發因子探測精度。

技術研發人員：楊紅磊,王澤平,石臻瀚,閆傲,趙松雪
受保護的技術使用者：中國地質大學（北京）
技術研發日：
技術公布日：2025/4/28