[发明专利]基于整数规划模型的MIMO InSAR相位解缠方法

权利要求书

1.本发明提供了一种基于整数规划模型的MIMO InSAR相位解缠的方法，其特征在于，主要步骤如下：

步骤R1，根据基线比、波长比和高程相位之间的对应关系构造整数规划模型；

步骤R2，采用分枝定界法，求解满足整数规划模型的最小解k_i(i＝1,2,3,···,N²)；

步骤R3，计算相应干涉图对应像元处的干涉相位微分的解模糊值dΦ_i；

步骤R4，分别以给定的相位解缠种子对各干涉图中所有像元处的干涉相位微分解模糊值dΦ_i，按照一定的路径积分，得到各干涉图的解缠结果Φ_i(x,y)；

在所述步骤R1中，将MIMO InSAR系统获得的N²幅SAR图像精确配准，且已去除了平地效应，根据基线比B₁:B₂:···:B_N和波长比λ₁:λ₂:···:λ_N，可得高程相位对应关系

其中，基线的最小公倍数为波长的最小公倍数为为干涉相位微分；k_i为的模糊数，且：

根据式(1)和式(2)，可将MIMO InSAR相位解缠问题转化为整数规划模型：

2.根据权利要求1所述的基于整数规划模型的MIMO InSAR相位解缠的方法，其特征在于，在所述步骤R2，采用分枝定界法对式(3)所描述的整数规划问题求解，可得满足S最小的解：k_i(i＝1,2,3,···,N²)，其具体步骤如下：

步骤1，求解与整数规划相对应的线性规划S，若S无可行解，则整数规划也没有可行解，计算停；若S的最优解是整数解，则该解即为整数规划的最优解，计算停；若S的最优解不是整数解，则转步骤2；

步骤2，分枝：在S的最优解中任选一个不符合整数条件的变量其值为(B^-1b)_i，[(B^-1b)]为小于(B^-1b)_i的最大整数，构造两个约束条件和将这两个约束条件分别加在问题S的约束条件上，形成两个子问题S₁和S₂，并求解S₁和S₂；

步骤3，定界：取整数解中最大目标值为界限值T下界，如果计算中尚无整数解，则取S＝-∞；检查分枝S_i，若它的最优解不是整数解，且S_i＞T，则重复步骤2，若S_i≤T，则S_i不再分枝；

步骤4，重复步骤2-3，直至所有分枝都不能再分解为止，这时界限值T对应的整数解即为原问题的最优解。

3.根据权利要求2所述的基于整数规划模型的MIMO InSAR相位解缠的方法，其特征在于,在所述步骤R3中，根据k_i(i＝1,2,3,···,N²)，可得各干涉图相应像元处的干涉相位微分的解模糊值dΦ_i；

4.根据权利要求3所述的基于整数规划模型的MIMO InSAR相位解缠的方法，其特征在于,在所述步骤R4中，分别以给定的相位解缠种子对各干涉图中所有像元处的干涉相位微分解模糊值dΦ_i，按照一定的路径积分，得到各干涉图的解缠结果Φ_i(x,y)；