[发明专利]基于修正嵌入式容积卡尔曼滤波的相位展开算法

说明书

技术领域

本发明涉及干涉合成孔径雷达(InSAR)技术测量应用中干涉相位展开算法领域，尤其涉及一种基于修正嵌入式容积卡尔曼滤波的相位展开算法。

背景技术

干涉合成孔径雷达(InSAR)可以高精度、高可靠性地获取地表三维信息和高程变化信息，被广泛应用在大地测绘、海洋监测、火山监测和地震检测等领域。干涉相位展开是干涉合成孔径雷达(lnSAR)技术应用中尤为关键的环节，其相位展开精度直接影响着干涉合成孔径雷达(InSAR)系统高程测量精度,故干涉图展开问题受到越来越多的关注。传统相位展开算法，包括枝切法、质量图引导法、网络流法、最小二乘法等，要在相位展开之前使用前置预滤波器抑制干涉图中的相位噪声，以提高展开相位的精度，但在尽可能地抑制干涉相位噪声的同时却易模糊掉干涉条纹的边缘特性。

随后，一类基于数据融合的相位展开算法陆续被提出，包括扩展卡尔曼滤波相位展开算法(EKFPU)、不敏卡尔曼滤波相位展开算法(UKFPU)、粒子滤波相位展开算法等。这类算法利用自身具有抑制噪声的能力，可以在展开缠绕像元的同时抑制相位噪声，在一定程度上弥补了传统相位展开算法使用前置预滤波器的不足，但其相位展开精度有待进一步提高。

发明内容

针对上述现有技术，本发明要解决的技术问题是如何进一步提高相位展开算法的精度。

为解决上述技术问题，本发明提供的技术方案是提供一种基于修正嵌入式容积卡尔曼滤波的相位展开算法(MECKFPU)，把经Levenberg-Marquardt方法修正后的嵌入式容积卡尔曼滤波器、基于修正矩阵束模型(AMPM)的局部相位梯度估计算法以及量化跟踪策略结合起来，利用修正嵌入式容积卡尔曼滤波器(MECKF)同时执行干涉相位的递推估计和相位噪声抑制，利用基于修正矩阵束模型(AMPM)的局部相位梯度估计算法快速和精确地从复干涉图中获取相位梯度信息，以及利用量化跟踪策略从高质量像元到低质量像元快速地搜索展开路径，以保证精确和高效地展开缠绕像元，主要步骤如下：

步骤1利用干涉图相干系数图以及微分偏差图构造路径引导图，并将其进行归一化和量化操作，同时创建附链表的优先队列和一个全局的指针，指针初始值设为零；

步骤2在路径引导图中选择一个质量最好的非边界像元作为初始像元，并默认它为已展开像元，把其缠绕相位作为展开相位；将与初始展开像元相邻的非边界缠绕像元标记为待展开像元，并根据待展开像元的量化路径引导值，把待展开像元插入到优先队列对应链表的顶部，同时指针指向量化路径引导值最大的待展开像元；

步骤3根据指针标记，从优先队列相应链表顶部获取最佳待展开像元x，并利用基于修正嵌入式容积卡尔曼滤波的相位展开算法(MECKFPU)展开此缠绕像元；从优先队列相应链表中删除像元x；把邻接像元x的非边界缠绕像元标记为待展开像元，并根据其量化路径引导值，把它们分别插入到优先队列相应链表的顶部；并使指针始终指向量化路径引导值最大的待展开像元；

步骤4优先队列中是否存在待展开像元？存在，转步骤3，否则转步骤5；

步骤5结束。

利用基于修正嵌入式容积卡尔曼滤波的相位展开算法(MECKFPU)同时执行干涉相位的递推估计和相位噪声抑制，过程如下：

利用干涉图中相邻干涉相位之间的关系，以及把归一化复干涉的同向分量和正交分量分别作为干涉相位的两个观测值，在沿某一确定路径下，相位展开系统模型如下：

其中，x_(m,n)和x_(a,s)分别表示(m,n)像元和(a,s)像元的真实干涉相位；表示干涉图中 (m,n)像元与(a,s)像元之间的相位梯度估计值，本文通过AMPM局部相位梯度估计算法获得；表示干涉图中(m,n)像元与(a,s)像元之间的相位梯度估计误差；v_1,(m,n)和v_2,(m,n)分别表示附加在复干涉正交分量和同向分量上的噪声。

嵌入式容积卡尔曼滤波器采用三阶嵌入式容积准则，使用带有权值的积分点实现对非线性系统后验均值和方差的逼近，具有良好的状态估计精度和数值稳定性，相应的积分点和权值计算如下：

上式中，n_x表示状态矢量维数；[δ]_j表示[δ]的第j列，其[δ]的详细取值为现有技术。针对上述系统方程，利用修正嵌入式容积卡尔曼滤波器(MECKF)执行干涉相位的递推估计，按如下进行：