[发明专利]基于混合整数优化模型的多基线InSAR相位解缠算法

权利要求书

1.基于混合整数优化模型的多基线InSAR相位解缠算法，包括如下步骤：

(1)输入SAR复图像数据；

(1a)输入干涉合成孔径雷达InSAR主天线获取的主图像数据；

(1b)输入干涉合成孔径雷达InSAR辅天线获取的已完成与主图像完全配准的辅图像数据；

(1c)输入干涉合成孔径雷达InSAR成像处理的处理参数和系统参数；

(2)复干涉相位图：

(2a)对输入的主、辅SAR图像进行干涉处理，得到多幅干涉相位图；

(3)获取局部最优窗口；

局部窗口的大小是混合整数优化模型的关键参数，根据干涉包相与假定线性地形的偏差，计算局部最优窗口的大小；所述局部窗口的大小即为方位和距离方向的像素数；步骤(3)所述的获取局部最优窗口是利用下式实现：

其中，P表示局部窗口方位向的像素点数量，S1和S2表示配准后的复SAR图像对，(m，n)表示取中心像素点的坐标，(·)*表示取共轭操作；

通过类似的方法可获得局部窗口距离向的像素点数量；

(4)构建混合整数优化模型；

利用InSAR系统参数和步骤(3)获得的局部最优窗口参数，将局部窗口的地形近似为斜面，构建混合整数优化模型；

(5)计算复干涉相位图的模糊整数；

对于每一个局部窗口，利用步骤(4)建立的混合整数优化模型求解两个干涉相位图的模糊参数；

(6)计算各个干涉图的绝对相位；

(7)输出整个场景的DEM。

2.根据权利要求1所述的基于混合整数优化模型的多基线InSAR相位解缠算法，其特征在于：步骤(4)所述构建混合整数优化模型利用下式实现：

minimising(|ε₁|+|ε₂|+|ε₃|)从属于

其中，|ε₁|和|ε₂|表示真实高度与斜面之间的最大距离，|ε₃|表示各个干涉图之间的高度差，a，b和c是构建一个斜平面方程的固有参数，X_i，j和R_i，j分别是相对像素中心的方位向和距离向，α₁和α₂是各个干涉图的绝对相位与高度的比例因子，和表示各个干涉图中通过精准SAR图像的共轭相乘法直接测量到的干涉相位，和表示有效的跨轨道的基线长度，和表示每个相位图的模糊数。

3.根据权利要求1所述的基于混合整数优化模型的多基线InSAR相位解缠算法，其特征在于：步骤(5)所述的计算复干涉相位图的模糊整数：

利用权利要求1所确定的最优窗口和权利要求2所述混合整数优化模型计算所得。