[发明专利]一种基于半参数模型的InSAR干涉相位解缠方法

权利要求书

1.一种基于半参数模型的InSAR干涉相位解缠方法，用于在InSAR或DInSAR数据处理过程中从干涉相位获取解缠相位，其特征在于，包括以下步骤：

1)半参数相位解缠模型的建立；

2)半参数相位解缠模型的解算及解缠相位的估计；

所述步骤1)中，半参数相位解缠模型具体为：

其中，V为观测值的残差，B为观测值的系数矩阵，为缠绕相位和解缠相位之间的关系，为解缠相位的估值，L为缠绕相位的相位差，为缠绕相位相位差中的系统误差的估值，m为缠绕相位的行数，n为缠绕相位的列数；

所述步骤2)具体为：

201)建立基于补偿最小二乘基本原则的半参数模型解算模型：

其中，α为一个标量，称作正则参数，R为正则矩阵，P为观测值的权阵，V^T为观测值的残差V的转置，为系统误差的估值的转置，所述正则矩阵R采用时间序列的方式构建；

202)基于拉格朗日极值基本原则对所述解算模型进行求解，获取半参数相位解缠模型的误差方程的法方程：

203)获取关于解缠相位的方程：

其中：

204)采用预共轭梯度迭代算法对所述关于解缠相位的方程进行求解，获取解缠相位的估值。

2.根据权利要求1所述的基于半参数模型的InSAR干涉相位解缠方法，其特征在于，所述缠绕相位的相位差利用InSAR或者DInSAR干涉相位计算获得。

3.根据权利要求1所述的基于半参数模型的InSAR干涉相位解缠方法，其特征在于，所述正则参数α采用L曲线的方法进行最优化选取，具体为：

A1)基于预共轭梯度迭代算法获取解缠相位的估值则系统误差的估值为：

A2)计算其中i为迭代次数；

A3)如果α所对应的L_i最小，则α为最优的正则化参数。

4.根据权利要求1所述的基于半参数模型的InSAR干涉相位解缠方法，其特征在于，所述步骤204)具体为：

B1)对进行初始化令同时获取预共轭矩阵

B2)计算z_i＝H^-1·r_i，判断当前迭代是否为第一次迭代，若是，则令p_i＝z_i，若否，则令w_i＝ρ_i/ρ_i-1，p_i＝z_i-1+w_ip_i-1，i为迭代次数；

B3)计算r_i＝r_i-1-v_iAp_i；

B4)判断是否存在||r_i||＞ε||r_i-1||，若是，则重复B2)、B3)进行迭代计算，若否，则执行步骤B5)；

B5)迭代终止获取解缠相位的估计值。