InSAR通道联合稀疏贝叶斯特征化成像

侯育星 徐 刚(陕西黄河集团有限公司 西安 710043)

②(东南大学信息科学与工程学院毫米波国家重点实验室 南京 210096)

1 引言

基于传统合成孔径雷达(Synthetic Aperture Radar， SAR) 2维成像，干涉合成孔径雷达(Interferometric SAR， InSAR)通过在高度向配置基线，利用多个天线相位中心在高度上的差异能够获取目标场景多视角观测，具备对观测场景3维测绘的能力[1-3]。相应地，InSAR在数字高程测绘(Digital Elevation Model， DEM)生成、场景高度估计和地形变化检测等方面已经得到了广泛应用[4，5]。在InSAR处理中，获取高精度的干涉相位是其中一项关键技术。然而，干涉相位不可避免的存在噪声是InSAR固有问题，这是因为通道之间视角差异引入SAR复图像之间的去相干问题，对应干涉相位噪声[6]。在幅度方面，相干斑噪声是SAR图像典型性特征，这是由SAR成像的相干积累特性所决定的。直观地讲，相干斑噪声造成相同散射介质的物体在SAR图像幅度表现出“椒盐”噪声特性，从而影响目标的分类和识别。因此，InSAR干涉图像幅度和相位同时存在噪声问题。为了满足后续图像应用的要求，InSAR图像的相干斑抑制和相位噪声滤波处理是必要过程。

近十几年来，许多学者较为系统的研究了干涉相位噪声滤波和相干斑抑制算法[7-10]。在现有大多数研究工作中，干涉相位噪声建模为加性噪声模型[11，12]，而相干斑则基于乘性噪声模型假设。绝大多数算法都是基于场景目标局部平稳假设，即同一散射介质目标相邻散射系数具有平稳特性，并且目标高度上连续变化(除了边缘区域)。近几年来，越来越多的学者研究利用结合幅度和相位提高两者噪声抑制的性能。文献[13，14]提出了一种联合幅度和干涉相位的稀疏正则化算法，通过联合构造幅相差分(Total Variation， TV)字典进行联合稀疏约束，以提升在城区场景的降噪能力。在文献[15]中，作者的前期工作研究了利用自适应表征字典进行InSAR图像干涉相位和幅度的联合稀疏表征，在稀疏成像中能够实现幅相滤波的特征化增强。文献[16]利用InSAR干涉处理，构造稀疏字典进行InSAR图像稀疏化表征，实现在降数据率采样下的图像重构。

在稀疏信号处理中，贝叶斯算法作为一种典型性的算法理论，相比其他算法，具有更高的稳健性和精度。相应的，稀疏贝叶斯算法在雷达成像中得到了较为广泛的研究和应用。进一步，结构化稀疏贝叶斯成像被提出用于提升雷达成像性能。文献[17]研究了联合稀疏模型在稀疏ISAR(Inverse SAR，ISAR)成像中的应用。文献[18]研究了结构化统计模型在ISAR成像中的应用。在文献[19]中，作者的前期工作研究了结构化稀疏贝叶斯模型在ISAR成像中的应用。文献[17-19]都是利用ISAR图像中散射目标相邻结构相关性以提升成像性能，其区别在于统计建模和算法求解的差异。对于InSAR而言，通道之间具有相干性，利用其相干性，通过结构化贝叶斯方法联合通道图像处理，能够进一步提升干涉相位滤波和相干斑抑制的性能，对应本文的主要研究内容。

本文提出了一种基于贝叶斯联合通道稀疏的InSAR成像算法，通过贝叶斯结构化稀疏化表征，提升干涉相位噪声滤波和相干斑抑制性能。首先，基于贝叶斯准则建立InSAR成像模型，通过多层级统计建模，将InSAR复图像联合幅度和干涉相位，建立为结构化稀疏统计模型。然后，利用最大期望(EM， Expectation Maximization)算法进行InSAR图像重构和多层级统计参数估计。相比较文献[15]的方法，本文方法由于利用贝叶斯结构化稀疏模型，能够进一步提升InSAR图像幅度和相位噪声滤波的能力。最后，通过实验分析进一步验证本文算法的有效性。

2 InSAR信号模型

图1所示InSAR成像几何，在SAR成像基础上，InSAR通过在高度向增加通道数，即通道(本文主要考虑两个通道情况，此时)，可实现对目标场景不同高度视角的观测。InSAR与SAR 2维成像几何相同，平台飞行方向为方位向，对应图1中X轴，与X轴垂直的波束指向为距离向，对应图1中的Y轴。如图1所示，在方位时间时刻，通道与目标P的距离可以表示为。假设已经利用SAR算法进行成像处理，此时数据在距离时域和方位多普勒域，第()个通道的信号可以表示为

3 所提方法

3.1 成像模型

在第2节中，式(2)实质是一个求逆以及降噪问题。在文献[15]中，作者的前期工作已经研究了利用自适应表征字典实现InSAR图像稀疏化表征，其稀疏成像模型为

基于式(5)和式(8)，由贝叶斯准则求得最大后验概率密度函数

至此，基于贝叶斯建模，本文已经建立了联合通道稀疏成像模型。假设给定一组估计参数，的最大后验概率估计为。在实际中，由于InSAR数据量通常较大，因此式(10)中的矩阵求逆运算量巨大，为了避免矩阵求逆，可以利用求解线性方程求得

式(11)的问题可以通过共轭梯度算法进行求解，参考文献[20]，利用句柄操作直接对2维图像或者进行距离和方位分维度处理，而不拉成向量操作，能有效提高算法的实际可行性。贝叶斯方法利用EM算法将问题求解分解为图像重构和参数的估计两个部分，两个部分交替迭代，构成一个整体。

3.2 算法求解

运算量分析：通过分析可知，本文算法的主要运算量在于式(10)和式(11)中和的计算，即矩阵的求逆操作。为了避免矩阵求逆，前面已经分析，可以利用式(11)通过共轭梯度算法求解。对于式(11)而言，其运算复杂度主要依赖于和操作，分别对应于SAR投影算子以及其逆过程操作，在本文中为方位傅立叶变换矩阵及其逆变换矩阵，其运算复杂度为。假设EM算法迭代次数为，在每次迭代中，求解式(11)的共轭梯度迭代次数为，那么本文算法的运算量为。其中，运算量与通道数成正比。

表1 本文算法流程框图Tab. 1 Algorithm flow chart in this paper

4 实验分析

下面通过实验分析以验证本文算法的有效性，并通过与现有方法比较，揭示本文算法的优势。在本节中，将实验分为2个部分。第1个实验部分为仿真数据实验，第2个部分为实测数据实验，包含山区数据和城区数据。下面所有实验均在个人电脑进行运行，电脑配置为CPU 3.50 GHz双核处理器，编程软件为R2015b版本MATLAB。

在仿真数据实验中，理想情况下图像幅度和干涉相位值如图2(a)和图2(b)所示。其中，幅度值对应平地和山区两个区域，每个区域对应相同的散射系数。图2(c)和图2(d)对应存在相干斑和相位噪声的图像。为了抑制噪声，利用本文算法进行稀疏化特征成像处理，其成像结果如图2(e)和图2(f)所示。可见，幅度的相干斑噪声和干涉相位噪声都得到了很好的抑制。为了更好地评估本文算法的性能，此处本文利用小波变换相位滤波[21]算法进行比较，其结果如图2(g)所示。通过比较图2(e)和图2(g)可知，本文算法在地形平坦区域和高度陡变区域具有更好的相位滤波效果。图2(h)和图2(i)进一步给出两种算法结果与真实相位的差异图。通过比较可知，本文算法的相位估计误差更小。同时，计算了两种算法相位估计的均方误差，本文算法和文献[21]算法结果分别为0.0082 rad2和0.0405 rad2，与图2(h)和图2(i)一致。在运算量方面，EM算法的迭代次数设置为20，算法可近似达到收敛，此时两个通道图像处理时间约为118 s。

在实测数据实验中，我们利用RADARSAT-2重复航过干涉数据进行实验，每个通道为单视复(Single Look Complex， SLC)图像。两个通道的录取回波数据间隔时间为24 h，并被应用于DEM生成，可以用来高程估计。在第一个实验中，选取山区场景的部分数据。通道1幅度的图像如图3(a)所示，两个通道的原始干涉相位如图3(b)所示。明显可见，干涉相位存在噪声，SAR幅度图像存在相干斑。然后，利用稀疏特征化增强的方法，在稀疏成像中进行噪声滤波处理。其中，本文利用文献[15]中作者前期工作的算法作为对比。在稀疏化表征中，都选择双树小波变换作为稀疏表征字典。图3(d)-图3(i)分别为文献[15]和本文方法的成像结果。相比传统SAR成像，两种稀疏方法都能够较为有效的进行干涉相位噪声滤波和相干斑抑制。同时，我们给出了传统SAR成像、文献[15]和本文方法结果的相位“残点”，如图3(c)、图3(f)和图3(i)所示。可见，文献[15]和本文方法在山区地形陡变部分，能够保留一定的细节信息。相比于文献[15]，本文方法在地形平坦局域具有更好的滤波效果。需要指出的是，文献[15]和本文方法由于稀疏字典构造等局限性因素，其滤波效果具有一定的低通特性，这是后续工作需要改进的地方。

第2个实验部分，本文选取城区场景进行实验。文献[15]和本文的成像结果如图4所示。由图4可知，本文的方法在在建筑物区域具有更好的噪声抑制能力，同时在能够较好地避免相位的不连续性，相应地，干涉相位图像的突变区域对应幅度图像的突变部分，表现为观测场景建筑物的边缘区域。为了定性分析，本文通过计算图像的等效视数(Equivalent Number of Looks， ENL)评估相干斑抑制的均匀场景区域强度图像的均值和方差。本文选取2个不同场景区域，以目标像素为中心，通过加窗(大小)选取相邻窗口像素计算ENL，其结果如表2所示。由表2可知，传统SAR成像的ENL结果近似为1，与使用的单视复图像一致。可见，本文的算法能够更加有效实现在同一散射介质区域的相干斑噪声抑制。相比于文献[15]的方法，本文算法可以取得更大的ENL值，对应较好的相干斑抑制性能，这得益于算法的结构化稀疏贝叶斯模型。为了评估算法对于独立散射点的分辨率保持能力，利用孤立散射点进行分辨率测量。传统SAR成像2维分辨率为2.67 m×4.13 m(理论值为1.92 m×3.27 m)。文献[15]和本文算法分辨率与传统成像比较接近，图5给出了方位脉冲响应形式，图中蓝色、红色和黑色分别对应传统成像、文献[15]和本文方法结果。

表2 ENL估计Tab. 2 ENL estimates

5 结论

本文针对InSAR稀疏成像，提出了一种通道联合的结构化稀疏贝叶斯成像算法。基于贝叶斯准则建立了联合通道结构化稀疏模型，并通过EM算法进行了图像重构和参数估计，相比于传统稀疏算法，本文算法利用了通道之间联合稀疏的相关特性，能够更加有效地提高InSAR成像算法中的噪声滤波特性，降低幅度噪声和相位噪声对成像结果的影响，实验结果证明了本文算法的有效性。