张斌，胡庆荣，李爽，韦立登

(1.北京无线电测量研究所，北京 100854;2.中国航天科工集团有限公司 第二研究院，北京 100854)

0 引言

InSAR技术通过利用2幅或多幅SAR图像的相位信息可以获取大范围、高精度的地表三维信息。其中相位解缠是InSAR处理的关键步骤，其处理精度对DEM(digital elevation model)产生重要影响[1-2]。传统的单基线相位解缠，大都假设Itoh条件满足的前提下进行的，其在相位噪声大和陡峭地形情况下很难满足高精度解缠要求，甚至失效，无法准确还原地形[3]。而多基线相位解缠算法能打破Itoh条件限制，获取更多的观测区的信息，融合多幅干涉图，减小干涉相位的模糊区间，降低解缠难度[4]。

近年来，多基线InSAR相位解缠得到了长足发展，Ghiglia[5]等提出多基线相位解缠绕方法的最小二乘法和最大似然函数法，Ferraiuolo[6]等提出基于马尔可夫随机场的最大后验概率估计方法，夏香根[7]等提出一种基于中国余数定理方法相位解缠，于瀚雯[8]提出基于聚类分析的多基线相位解缠绕算法。而在多基线最大后验概率估计方法中，Ferraioli[9]等提出了基于图割的多基线相位解缠，其直接由多幅干涉条纹图直接计算高程值，从而避开了相位解缠过程。但该类方法由于采用Ishikawa图割[10]，其所建立的高度量化值越多，其相位解缠所用的内存和处理所花的时间越大，因此A.shabou[11]等提出了一种基于二分法的多基线相位解缠方法，其所占用的内存减少，但其每一个循环都要访问所有的高度量化值，相应的计算时间却增加。相比直接估计……