[发明专利]基于后向序贯最小二乘的地表形变监测方法及系统

权利要求书

1.一种基于后向序贯最小二乘的地表形变监测方法，其特征在于，包括：

确定监测区域在当前时间段内每一个时间单元的第一形变速率及对应的协因数矩阵；

确定监测区域在历史时间段内的SAR影像对应的多幅历史解缠干涉图；

根据多个所述第一形变速率、协因数矩阵和所述多幅历史解缠干涉图，利用后向序贯最小二乘方法，确定所述历史时间段内的每一个时间单元的第二形变速率并更新第一形变速率；

根据多个更新后的所述第一形变速率和多个所述第二形变速率，确定所述监测区域在所述当前时间段和所述历史时间段之间的累积形变。

2.根据权利要求1所述的基于后向序贯最小二乘的地表形变监测方法，其特征在于，确定监测区域在当前时间段内每一个时间单元的第一形变速率及对应的协因数矩阵，具体包括：

确定监测区域在当前时间段内的SAR影像对应的多幅解缠干涉图；

根据所述多幅解缠干涉图，利用SBAS-InSAR技术确定监测区域在当前时间段内每一个时间单元的第一形变速率及对应的协因数矩阵。

3.根据权利要求2所述的基于后向序贯最小二乘的地表形变监测方法，其特征在于，根据第一公式确定监测区域在当前时间段内每一个时间单元的第一形变速率，所述第一公式为：

V₁＝A₁X₁-L₁,P₁

式中，V₁为第一残差向量，A₁为M₁×N₁的设计矩阵，该矩阵由当前时间段内的多个时间单元构成，M₁为所述解缠干涉图的数量，N₁为当前时间段内的所述SAR影像的数量，L₁为所述解缠干涉图的差分干涉相位组成的向量，P₁为L₁对应的权矩阵，X₁为多个所述第一形变速率组成的向量，X₁＝(A₁^TP₁A₁)^-1A₁^TP₁L₁。

4.根据权利要求3所述的基于后向序贯最小二乘的地表形变监测方法，其特征在于，根据第二公式确定监测区域在当前时间段内每一个时间单元的第一形变速率的协因数矩阵，所述第二公式为：

式中，为X₁的协因数矩阵，为A₁的转置矩阵。

5.根据权利要求1-4任一项所述的基于后向序贯最小二乘的地表形变监测方法，其特征在于，根据第三公式确定所述历史时间段内的每一个时间单元的第二形变速率，所述第三公式为：

式中，V₂为第二残差向量，Y为多个所述第二形变速率组成的向量，X₁为多个所述第一形变速率组成的向量，X₂为X₁重新平差后的向量，即为更新后的第一形变速率，B为Y的系数矩阵，A₂为X₂的系数矩阵，L₂为所述历史解缠干涉图的差分干涉相位组成的向量，P₂为L₂对应的权矩阵。

6.根据权利要求5所述的基于后向序贯最小二乘的地表形变监测方法，其特征在于，所述第三公式对应的平差准则如第四公式，所述第四公式为：

式中，V₂为第二残差向量，X₂为X₁重新平差后的向量，即为更新后的第一形变速率，P₂为L₂对应的权矩阵，为X₁的协因数矩阵。