[发明专利]一种基于序贯平差的InSAR时序地表形变监测方法

权利要求书

1.一种基于序贯平差的InSAR时序地表形变监测方法，其特征在于，包括以下步骤：

步骤1：获取待监测的地表区域N+1幅时序SAR影像，并获得对应的M幅历史解缠差分干涉相位图；

步骤2：基于历史解缠差分干涉相位结合和W⁰，求解待监测区域的历史时序形变速率v^p,0、历史地形残差

其中，表示地形残差求解系数矩阵，Δt⁰为历史干涉图时间间隔向量，H⁰为历史干涉图对应的高程转换系数向量，W⁰为M×M的单位矩阵；

步骤3：利用去除地形残差相位的M幅历史解缠差分干涉相位结合和W⁰求解历史时序形变相位从而获得历史时序形变

其中，λ表示雷达波长，表示时序形变相位求解系数矩阵，大小为M*N，M是历史解缠差分干涉相位图的数量，矩阵的第k行表示第k幅干涉图，第k行中的第i个元素值为-1，第j个元素值为1，其余元素值均为0，第k幅干涉图是由第i，j两幅SAR影像干涉生成；

步骤4：获取新增SAR影像，并从新增的SAR影像中选取满足时空基线阈值的SAR影像对，得到新增的M₁幅解缠差分干涉相位图；

步骤5：利用新增的M₁幅解缠差分干涉相位和待监测区域的W⁰，基于序贯平差求得历史地形残差相关未知参数向量的改正数向量进而得到更新后的地形残差相关未知参数向量即获得了待监测区域更新后的时序形变速率v^p,a和地形残差

其中，表示历史地形残差相关未知参数向量，表示更新后的地形残差相关未知参数向量，

步骤6：利用W⁰和对历史时序形变相位进行初步更新，得到时序形变相位初步更新值表示待监测区域历史地形残差的改正数向量；

步骤7：基于去除地形残差相位的M₁幅新增解缠差分干涉相位和初步更新的时序形变相位利用序贯平差对地表时序形变相位进行更新，得到待监测区域的实时时序形变相位从而获得对应的时序形变

2.根据权利要求1所述的方法，其特征在于，所述历史地形残差相关未知参数向量的改正数向量按以下公式计算获得：

其中，Δv^p,a和分别表示待监测区域历史时序形变速率和地形残差的改正数向量；J₂、表示中间变量；

表示新增M₁幅解缠差分干涉相位图对应的地形残差求解系数矩阵，Δt^a为新增M₁干涉图对应的时间间隔向量，H^a为新增M₁干涉图对应的高程转换系数向量；

W^a表示新增M₁幅解缠差分干涉相位图的权重，取值为M₁×M₁的单位矩阵。

3.根据权利要求2所述的方法，其特征在于，所述时序形变相位初步更新值的计算过程如下：

其中，表示表示p点处第k幅干涉图的垂直基线，R^p表示p点到卫星的斜距，θ^p表示p点处的雷达入射角。