[发明专利]一种迭代最大似然估计多基线InSAR相位解缠方法

说明书

技术领域

本发明属于雷达技术领域，它特别涉及干涉合成孔径雷达(InSAR)相位解缠技术领域。

背景技术

干涉合成孔径雷达(InSAR)是合成孔径雷达(SAR)成像技术的扩展，可以测量得到目标的高度信息，从而实现对目标场景的三维测绘，该技术主要利用两个或两个以上的天线对同一观测场景进行不同角度观测得到的两幅以上SAR复图像，通过两幅SAR复图像信息可以得到干涉相位信息，经过干涉处理后的相位信息通过高程反演后得到场景的地形高程信息。InSAR技术不但可以全天候、全天时成像，并且成像区域大、成像精度高，已经成为当今世界获取地形三维图像和地形高程变化信息的一项重要遥感技术。其中多基线InSAR有着测量精度更高、能更有效地处理复杂地形等优点，多基线InSAR也得到了越来越多的关注，也是当今InSAR研究的重要方向之一。

在多基线InSAR中，每个副天线的复图像和主天线复图像做干涉处理可以得到一组干涉相位，由于每个副天线的位置不同，所以导致了每个副天线和主天线得到的干涉相位存在一定的差异，而基线长度与干涉相位的条纹个数有关，并且长基线对应的副天线得到的干涉相位具有更多的干涉条纹，所以能够反演得到精度更高的地形信息。

相位解缠是InSAR处理的关键步骤，InSAR处理过程中相位提取得到的干涉相位都是缠绕的，进行地形高程重构前必须对缠绕的干涉相位进行解缠，相位解缠精度的高低，直接影响后面高程反演的精度。为了对复杂地形的缠绕相位和低信噪比条件下的缠绕相位进行相位解缠，高精度的多基线InSAR相位解缠方法有了深入的研究。

最大似然估计多基线相位解缠方法利用多个通道的缠绕相位进行相位解缠，根据最大似然函数值求得的解缠相位，解缠精度较高，然而，在有噪声存在时，最大似然函数值对应的解缠相位一定是准确的，有可能是另外一个较大的似然函数值对应的解缠相位，这时得到的解缠相位就有较大误差。详见文献“G.Fornaro.Maximum likelihood multi-baseline SAR interferometry.2005.”。

发明内容

本发明提出了一种迭代最大似然估计多基线InSAR相位解缠方法，该方法首先采用最大似然函数估计方法，得到多基线InSAR初始的解缠相位，再求出几组较大似然函数值对应的解缠相位，然后通过多组较大似然函数对应的解缠相位迭代优化初始解缠相位，减小相位解缠过程中噪声带来的误差，为后续的高精度InSAR高程反演提供保证。为了方便描述本发明的内容，首先作以下术语定义：

定义1、干涉合成孔径雷达(InSAR)

干涉合成孔径雷达(InSAR)指利用同一观测场景不同观测角度的两组或者两组以上SAR数据进行干涉成像处理，然后结合雷达系统参数和雷达平台几何位置信息反演地形高度及高程变化信息的遥感技术，详见文献“合成孔径雷达成像原理”，皮亦鸣等编著，电子科技大学出版社出版。

定义2、相位解缠

一切将相位由主值或相位差恢复为真实值的过程统称为相位解缠。详见文献“星载合成孔径雷达干涉测量”，王超等编著，科学出版社。

定义3、似然函数：

设总体X服从分布P(x；θ)(当X是连续型随机变量时为概率密度，当X为离散型随机变量时为概率分布)，θ为待估参数，X₁,X₂,….X_n是来自于总体X的样本，x₁,x₂,…,x_n为样本X₁,X₂,….X_n的一个观察值，则样本的联合分布(当X是连续型随机变量时为概率密度，当X为离散型随机变量时为概率分布)称为似然函数。详见文献“概率论与数理统计”，徐全智等编著，电子科技大学出版社。

定义4、缠绕相位：

缠绕相位是指在实际干涉处理中，经过三角运算，得到的干涉相位的主值。详见文献“星载合成孔径雷达干涉测量”，王超等编著，科学出版社。

定义5、解缠相位：

解缠相位是指干涉处理中从缠绕相位经过相位解缠处理恢复的真实相位。详见文献“星载合成孔径雷达干涉测量”，王超等编著，科学出版社。

定义6、高斯分布函数randn(m,n)

高斯分布函数randn(m,n)是MATLAB中产生均值为0，标准差为1的高斯分布，并且为m×n的随机矩阵。