[發明專利]一種用于大壩及滑坡變形GB-SAR監測的數據后處理方法及系統

在審

• 申請號: 202011596817.1
• 申請日: 2020-12-29
• 公開（公告）號: CN112685819A
• 公開（公告）日: 2021-04-20
• 發明（設計）人: 項霞;陳辰;劉超;聶銳化;王輝;魯恒;劉鐵剛;陳建康;楊正麗
• 申請（專利權）人: 四川大學
• 主分類號: G06F30/13
• 分類號: G06F30/13;G06T17/05
• 代理公司: 成都玖和知識產權代理事務所(普通合伙) 51238
• 代理人: 胡琳梅;王海權

1.一種用于大壩及滑坡變形GB-SAR監測的數據后處理方法，其特征在于，包括以下步驟：

步驟1：獲取N個監測影像子集，對監測影像子集內部數據進行處理分析，得到N幅平均影像；

步驟2：以其中一幅平均影像作為主影像，其余N-1幅平均影像作為從影像進行干涉處理，得到N-1幅干涉影像；

步驟3：選取PS點，并建立Delaunay三角網；

步驟4：PS點差分相位建模及參數估計，包括：

建立相鄰兩PS點基線差分干涉相位的最佳模型；

根據所述差分干涉相位的最佳模型估算Delaunay三角網中每條基線的線性變形速率增量和殘余相位增量；

步驟5：將估算得到的所述線性變形速率增量在空間上采用區域網最小二乘法間接平差處理，得到各PS點的絕對線性變形速率和線性變形量；

步驟6：將估算得到的所述殘余相位增量在空間上按照最小二乘法間接平差處理，得到各PS點的絕對殘余相位，對所述絕對殘余相位進行分離和提取，得到非線性變形量；

步驟7：將各個PS點的線性變形量與非線性變形量相加，獲得各個PS點上的總體變形量，并采用插值方式獲取整個監測物體全域范圍的密集變形場。

2.根據權利要求1所述的變形GB-SAR監測的數據后處理方法，其特征在于，所述建立相鄰兩PS點基線差分干涉相位的最佳模型方法為：

在GB-SAR大壩或滑坡監測的第i幅干涉影像中，任意相鄰的兩PS點，P、Q點，連接基線的差分干涉相位模型表示為：

式中，Δv(x_p,y_p；x_q,y_q；T_i)是兩PS點P(x_p,y_p)、Q(x_q,y_q)間的線性變形速率增量，Δw_i(x_p,y_p；x_q,y_q；T_i)是兩PS點P(x_p,y_p)、Q(x_q,y_q)間的殘余相位增量，i＝1,2,...N，T_i為第i幅從影像與主影像之間的時間間隔；

將線性變形速率增量Δv進行優化處理估算，此時，差分干涉相位的模型即為最佳模型。

3.根據權利要求2所述的變形GB-SAR監測的數據后處理方法，其特征在于，所述優化處理估算包括：

使整體復相關系數γ_k最大化，所述整體復相關系數γ_k最大化模型表示如下：

式中，Δω_i為觀測值與估計值之差，N為每對PS點對應的N幅差分干涉相位圖；

根據取得的整體復相關系數γ_k最大值，估算出線性變形速率增量Δv。