[发明专利]一种InSAR干涉相位真值确定及差分干涉测量方法

权利要求书

1.一种InSAR干涉相位真值确定方法，其特征在于，用于获得InSAR差分干涉图中一个像素点观测值的真值，所述的方法包括：

步骤1、通过建立观测方程求取当前像素点的m个观测值各自的闭合残差，根据闭合残差的大小从所有观测值中选择n+k个初选观测值，n为观测方程中设计矩阵的秩，m1，nm,k＝[0,1,…,m-n]；

步骤2、将m个观测值按照各自真误差绝对值的大小从小到大进行排序，获得观测值序列；

若|dΔ|_max＝||Δ_n+k|-|Δ_n+k-1||，将步骤1中选择的n+k个初选观测值作为复选观测值，执行步骤3；否则，k＝k+1，从所述观测值序列中选择前n+k个观测值作为初选观测值，重新执行步骤2；

其中，dΔ为所述观测值序列中相邻观测值的真误差绝对值之间的差值，Δ_n+k为观测值序列中第n+k个观测值的真误差，Δ_n+k-1为观测值序列中第n+k-1个观测值的真误差；

步骤3、根据每一个观测值的评价参数绝对值的大小从m个观测值中选取y个观测值，获得一个y×1维的观测值向量L_y；

其中，第i个观测值的评价参数i＝1，2，…，m，Δ_i为第i个观测值的真误差，Δ_z为第z个观测值的真误差，z＝1，2，…，n+k；

步骤4、利用式VII计算当前像素点观测值的解缠误差向量

式VII中，L_r为r×1维的观测值向量，r＝m-y，R为m×1维的观测值向量L的平差因子矩阵，其中R_y×y与R_y×r为根据识别出来的y个解缠误差位置对R所作出的分析，R表达式为，

步骤5、利用式X获得当前像素点观测值的真值L'：

其中，为对填充后的m×1维解缠误差向量，填充方式为将其余的不是含有解缠误差的m-y个观测值的解缠误差向量用0填充。

2.如权利要求1所述的InSAR干涉相位真值确定方法，其特征在于，所述的根据闭合残差的大小从所有观测值中选择n+k个初选观测值时，将m个观测值按照其闭合残差的大小从小到大进行排序，选取前n+k个观测值作为初选观测值。

3.如权利要求2所述的InSAR干涉相位真值确定方法，其特征在于，将m个观测值按照其闭合残差的大小从小到大进行排序，选取前n+k个观测值作为初选观测值时，采用单位权的方法进行选取。

4.如权利要求1所述的InSAR干涉相位真值确定方法，其特征在于，所述的根据每一个观测值的评价参数绝对值的大小从m个观测值中选取y个观测值时，从m个观测值中选取评价参数绝对值大于阈值的y个观测值。

5.一种差分干涉测量方法，其特征在于，所述的方法包括：

步骤A、在SBAS-InSAR框架下，设置时空基线阈值，将SAR影像组网，并分别进行干涉测量，获得干涉组合的InSAR差分干涉图；

步骤B、采用权利要求1-4中任一项权利要求所述的InSAR干涉相位真值确定方法对所述的InSAR差分干涉图中的一个像素点进行处理，获得该像素点观测值的真值；

步骤C、重复步骤B，直至获得InSAR差分干涉图中所有像素点观测值的真值。