[发明专利]一种顾及地形起伏的InSAR残余运动误差估计方法及装置

权利要求书

1.一种顾及地形起伏的InSAR残余运动误差估计方法，其特征在于，包括以下步骤：

估计从影像对应的地形相位：获取研究区域的N+1景SAR影像，依次进行配准、重采样、去平地相位和干涉处理生成N*(N+1)/2个干涉对，假设公共主影像的相位为0，根据干涉得到的N*(N+1)/2个干涉对估计N个从轨迹中包含InSAR残余运动误差的地形相位；

获取差分相位观测值：利用覆盖研究区域的外部DEM数据和InSAR干涉几何参数，计算模拟地形相位；将包含InSAR残余运动误差的地形相位与模拟地形相位求差，并对差值进行相位解缠处理，再将解缠差分干涉相位作为观测值；

构建顾及地形起伏的差分干涉相位参数化模型：根据包含InSAR残余运动误差的地形相位和模拟地形相位分别关于InSAR干涉几何的计算式，计算得到差分干涉相位关于InSAR干涉几何、DEM误差以及基线误差分量的表达式，并将该表达式作为差分干涉相位参数化模型；

构建多基线观测方程组：针对每个从影像的每一距离行像素，均将差分干涉相位的观测值与差分干涉相位参数化模型进行匹配构建等式，得到观测方程；所有观测方程构成多基线观测方程组；

解算残余运动误差：基于已知的InSAR干涉几何参数和观测值，解算多基线观测方程组中的DEM误差和基线误差分量；其中的基线误差分量即为待估计的残余运动误差。

2.根据权利要求1所述的方法，其特征在于，还包括按预设迭代次数，迭代计算新的残余运动误差和DEM误差，将最后一轮迭代计算得到的残余运动误差作为最终的残余动作误差；其中，每一轮迭代计算方法为：

(1)使用上一轮解算得到的DEM误差，对当前的外部DEM数据的地形高程进行补偿，得到更新的外部DEM数据；

(2)使用更新的外部DEM数据和InSAR干涉几何参数，重新构建观测方程和解算残余运动误差，得到新一轮的残余运动误差和DEM误差。

3.根据权利要求1所述的方法，其特征在于，模拟地形相位的计算式为：

式中，表示第i个模拟地形相位，h代表外部DEM数据中的地形高程，B_⊥,i代表正交于第i个从轨迹的斜距的基线分量，λ表示波长，R_i0代表第i个从轨迹成像几何下的斜距，θ₀为公共主影像的视角；

包含InSAR残余运动误差的地形相位关于InSAR干涉几何参数的计算式为：

式中，表示第i个从轨迹的干涉地形相位，ΔB_y,i和ΔB_z,i代表第i个从轨迹的基线误差水平分量和垂直分量，θ_i代表第i个从轨迹的入射角，代表真实地形高程；

计算得到的差分干涉相位参数化模型为：

式中，表示第i个从轨迹对应的差分干涉相位，Δh为DEM误差。