[发明专利]一种快速更新InSAR变形序列结果的方法

权利要求书

1.一种快速更新InSAR变形序列结果的方法，其特征在于，包括如下步骤：

A、根据获取的相干目标点信息及时序相位信息，只对新增SAR影像的差分干涉图进行相位提取；

所述步骤A具体包括：

A1、新增加的SAR影像与历史SAR影像集合的主影像进行高精度配准；

A2、将配准后的新增SAR影像与历史SAR影像集合进行干涉处理生成干涉图，并根据数字高程模型的数据去除干涉图中的平地相位和地形相位；

A3、利用永久散射体的位置对差分干涉图进行相位提取，通过时空同质滤波和极大似然估计对同分布散射体的差分干涉相位进行提取；

所述步骤A3具体包括：

A31、假设历史SAR影像集合的差分干涉图共有N景，所述同分布散射体的历史时序相位为其中k的取值范围是1到N的正整数；

A32、基于同分布散射体对新增SAR影像与历史SAR影像集合生成的差分干涉图θ_k进行时空同质滤波处理，获取滤波相位值其中k的取值范围是1到N的正整数；

A33、采用傅立叶变换法或小波变换法去除差分干涉图中的残余系统相位，计算同分布散射体的相干性γ_k；

A34、获得新增SAR影像所对应的第N+1景差分干涉图的时序相位的极大似然估计值

其中，i为复数单位，arg代表计算相位，*代表共轭转置，采用非线性优化算法来解算

B、利用解算参数和解缠相位对新增SAR影像的差分干涉图进行相位解缠；

C、利用新增SAR影像的差分干涉图的相位解缠结果，更新解算参数和变形序列。

2.根据权利要求1所述的快速更新InSAR变形序列结果的方法，其特征在于，所述相干目标点包括永久散射体的像素点和同分布散射体的像素点。

3.根据权利要求2所述的快速更新InSAR变形序列结果的方法，其特征在于，所述非线性优化算法采用拟牛顿法。

4.根据权利要求3所述的快速更新InSAR变形序列结果的方法，其特征在于，所述解算参数包括变形速率和高程。

5.根据权利要求4所述的快速更新InSAR变形序列结果的方法，其特征在于，所述差分干涉图的相位信息包括大气相位、地形相位、变形相位以及相位噪声，所述地形相位利用高程解算参数进行去除。

6.根据权利要求5所述的快速更新InSAR变形序列结果的方法，其特征在于，所述步骤B具体包括：

B1、假设历史SAR影像的差分干涉图共有N景，所对应的时序相位解缠结果表示为Ψ_k，k∈[1,N]；

B2、所述第N+1景差分干涉图与第N景差分干涉图的时序相位差值

B3、采用空间相位解缠方法对时序相位差值进行处理，获取解缠时序相位值Ψ_r；

B4、获得第N+1景差分干涉图的时序相位解缠结果：Ψ_N+1＝Ψ_N+Ψ_r。

7.根据权利要求6所述的快速更新InSAR变形序列结果的方法，其特征在于，所述空间相位解缠方法包括枝切线法和最小费用流法的其中之一。

8.根据权利要求7所述的快速更新InSAR变形序列结果的方法，其特征在于，所述步骤C具体包括：

C1、利用时间域高通滤波和空间域低通滤波方法对大气相位进行估计并将其从差分干涉图中去除，然后利用时序滤波方法对残余大气相位和相位噪声进行削弱，得到监测点的变形序列；

C2、根据观测时间和变形序列对监测点的变形速率v进行更新：

Δt_k·v＝d_k，k∈[1,N+1]

其中，Δt_k＝t_k-t₁，t_k为第k景SAR影像的观测时间，d_k为第k景SAR影像的累计变形量，通过最小二乘法来解算d_k。