[发明专利]一种利用史赖伯规则的InSAR干涉参数区域网平差方法

权利要求书

1.一种利用史赖伯规则的InSAR干涉参数区域网平差方法，其特征在于，该平差方法包括以下步骤：

1)获取图像覆盖范围内的高程控制点和各个干涉图的解缠结果；

2)以每个干涉图像对为基本平差单元，利用高程控制点和不同像对重叠处的同名连接点，根据InSAR干涉参数定标原理得到控制点误差方程和连接点误差方程；

3)建立连接点误差方程的等效误差方程式，以该等效误差方程进行InSAR干涉参数区域网平差的平差求解，以得到相应的干涉参数；

若整个区域参与平差的干涉像对数为t，对某个控制点j，将设它处于m个干涉像对重叠处，m≤t，则可以列误差方程式，误差方程中的未知数即为m个干涉像对的干涉参数改正数Δ_O1,Δ_O2,…,Δ_Om，误差方程式为：

式(14)简化为：

V_GCPj＝A_GCPjΔ_GCPj-L_GCPj P_GCPj (15)

式(15)中，A_GCPj为干涉参数改正数系数矩阵，P_GCPj为权矩阵；V_GCPj＝[v₁ v₂ … v_m]^T_m×1为残差向量；Δ_GCPj＝[Δ_O1 Δ_O2 … Δ_Om]^T_3m×1为干涉参数改正数向量；L_GCPj＝[l_1j l_2j … l_mj]^T_m×1为常数项向量，且有：

其中，Δ_Oi＝[ΔB_i Δα_i Δφ_0i]^T为第i个干涉像对的干涉参数改正数向量，A_ij＝[b_0ijb_1ij b_2ij]第i个干涉像对上控制点j的干涉参数改正数系数向量；

对于某个连接点k，假设它处于n个干涉像对重叠处，n≤t，由于高程测量值未知，误差方程中未知数除了n个像对的干涉参数改正数Δ_O1,Δ_O2,…,Δ_On外，还包含该连接点的高程改正数Δ_Gk，原始误差方程式为：

式(16)简化为：

V_TPk＝A_TPkΔ_TPk+C_TPkΔ_Gk-L_TPk P_TPk (17)

式(17)中，V_TPk＝[v₁ v₂ … v_n]^T_n×1为残差向量；Δ_TPk＝[Δ_O1 Δ_O2 … Δ_On]^T_3n×1为干涉参数改正数向量；Δ_Gk＝[Δh_k]为连接点k的高程改正数；L_TPk＝[l_1k l_2k … l_nk]^T_n×1为常数项向量；A_TPk为干涉参数改正数系数矩阵，C_TPk为连接点高程改正数系数矩阵，P_TPk为权矩阵，且有：

其中，A_ik＝[b_0ik b_1ik b_2ik]为第i个干涉像对上连接点k的干涉参数改正数系数向量；C_ik＝[b_3ik]为第i个干涉像对上连接点k的高程改正数系数；

由于区域网平差需要答解干涉参数和连接点高程值两类未知数，通过连接点数量较大，需求解干涉参数，再解算连接点的高程；依据史赖伯规则，对连接点的原始误差方程式构建等效误差方程，从而求解干涉参数；式(17)的等效公式如式(18)和式(19)所示：

式(18)中，

式(18)和式(19)可以简化为：

V_VTPk＝A_VTPkΔ_VTPk-L_VTPk P_VTPk (20)

其中，A_ik、C_ik的含义与式(17)相同；

按照公式(15)和(20)构建成误差方程式，并简化为：

V_EE＝A_EEΔ_O-L_EE P_EE (24)

公式(24)即为等效误差方程式，其中Δ_O为干涉参数改正数向量，V_EE为残差向量，A_EE为干涉参数改正数系数矩阵，L_EE为常数项向量，P_EE为权矩阵，且有

对等效误差方程式法化，得到的法方程N_EEΔ_O＝U_EE的系数矩阵和常数项矩阵为：