衍射层析模型下穿墙雷达三维学习成像方法

卞 粱， 晋良念,2， 刘庆华

(1.桂林电子科技大学信息与通信学院, 广西桂林 541004；2.广西无线宽带通信与信号处理重点实验室， 广西桂林 541004)

0 引 言

超宽带穿墙雷达三维成像不仅能提供距离向和方位向上的目标位置信息，还能提供俯仰等更多维度的信息，满足了建筑内部的结构特征、人体目标的姿势等信息的需求。其现有的三维成像有后向投影算法[1]和衍射层析(Diffraction Tomographic，DT)算法[2-3]等。但是这些成像算法不仅需要进行大量的数据采样、存储和处理，而且所成像还存在主瓣较宽、旁瓣较高的问题。

为此，国内外学者已将压缩感知(Compressive Sensing，CS)理论应用于穿墙雷达成像(Through- the-Wall Radar Imaging，TWRI)中，不仅减少了数据的采集时间、采集量和存储代价，还可以获得高分辨成像。文献[4]首次将CS应用于穿墙场景中，利用稀疏随机矩阵的采样方法可在随机的采样天线中随机选择频点以减少数据采集量。在此基础上，文献[5]通过正交匹配追踪(OMP)算法对墙后目标成像，提高了成像质量和减少了旁瓣干扰。由于OMP重构算法易受到信号信噪比影响导致成像质量较差，文献[6]在贝叶斯算法的基础上通过噪声先验信息，同时考虑扩展目标像素间的结构信息，对目标进行高分辨成像。考虑到贝叶斯算法复杂度高，计算量大，文献[7]在低秩稀疏约束下提出了一种基于奇异值和l1范数的快速软阈值迭代(FISTA)的凸优化算法重构，该算法具有计算量小、收敛速度快的特点。但是，这种稀疏重构方法的字典矩阵是通过目标到成像网格的时延构建的，对于三维成像，字典矩阵所需内存太大，难以提供硬件需求；而且它往往需要预先确定阀值参数，而恰当地选择该参数对重建图像质量具有非常重要的影响。

目前，深度神经网络(Deep Neural Networks， DNNs)已经在图像去噪、图像复原等领域展现出优势。而模型驱动[8-9]作为DNNs中特殊的一类，将凸优化算法中的单次迭代过程映射为网络结构的每一层，进而利用少量的样本即可学习出最优参数。此外，因网络结构的显示表达，所以也保证了高质量的图像重建。鉴于此，本文基于这种模型驱动的思路，通过衍射层析成像算法构造快速傅里叶变换算子，避免了字典矩阵的构建，减少了所需内存；通过修正近似消息传递(Approximate Message Passing，AMP)算法进行重构，并将其迭代过程映射成多层神经网络，最后通过数据驱动学习近似消息传递算法中可调参数，实现高分辨、高效率的三维图像的重建。

1 信号模型

如图1所示，使用收发天线共置的雷达系统平行于前墙的方向对墙后目标进行探测。系统发射步进频率信号，扫频范围从f0开始，步进频率为Δf，频率点数为K；采用M行N列二维收发共置天线阵列，天线空间坐标向量记为(xR,yR,zR)。空间区域被分成3个区域，区域0和区域2为自由空间，区域1为墙体，自由空间的介电常数和磁导率记为ε0和μ0，墙体介电常数和磁导率记为εb和μb，所述隐蔽目标位于区域2中；设x方向为方位向，y方向为高度向，z方向为距离向，墙体的厚度为d。

图1 穿墙雷达探测场景

根据文献[2]和[3]，图1所示场景的衍射层析三维成像表达式为

式中，O(x,y,z)代表目标像函数,k=2πf/c代表自由空间波数(这里的c为光速)，kx代表方位向波数，ky代表高度向波数，kz代表距离向波数，ξ0=120π代表自由空间中的波阻抗,F(kx,ky,k)代表并矢函数，即

F(kx,ky,k)=

exp(-jkxxR-jkyyR)dxRdyR

(2)

J(kx,ky,k,z)=

则式(1)可以简化为

O(x,y,z)=

J(kx,ky,k,z)exp(jkxx+jkyy)

这里，等式右边的第2和第3个积分项也是一个二维空域傅里叶逆变换，则式(2)就继续简化为

J(kx,ky,k,z)}

(3)

由于使用步进频率雷达信号，自由空间波数k积分可以转化为求和形式，则式(3)可进一步变为

J(kx,ky,k,z)}

(4)

可以看出，式(4)是将回波散射数据变换为目标图像。因此，我们可以借鉴文献[10]的方法，将其利用傅里叶变换反演，得到衍射层析成像正向变换，即

J-1(kx,ky,k,z)}

(5)

E=Ψ(O)

(6)

式中，假设O与E的维度相同，均为M×N×K；当O和E的维度在不等的情况可以通过补零来实现。

在压缩感知理论框架下，对频点和阵元位置进行随机降采样，但对阵元位置进行随机选取后要对未被选取的阵元作补零处理，则式(6)改写为

Y=ΦE=ΦΨ(O)+N

(7)

(8)

式中，λ为正则化参数。

2 算法描述

2.1 LAMP网络构建

对于式(8)，可以采用AMP算法[11]求解。但是，其迭代过程还必须考虑衍射层析模型中的快速傅里叶变换算子Ψ，所以得到修正的AMP算法为

(9)

(10)

(11)

(12)

式中，ΨH为Ψ的共轭转置算子，由式(4)可得

根据文献[11]的模型驱动学习方法，我们将式(9)～(12)的AMP迭代过程映射成多层神经网络，然后通过数据驱动进行学习，即所谓的学习近似消息传递算法(Learning Approximate Message Passing，LAMP)。LAMP的网络结构如图2所示，设定有T层。根据式(9)～(12)的迭代过程，它的每一层应设计成如图3所示的3个隐层，分别是残差层、算子更新层和非线性变换层，具体如下：

1) 残差层：对上一层残差Dt-1进行更新，输出本层的残差Dt。为简化网络结构，将式(9)代入式(10)，可得

(13)

(14)

(15)

η(Rt;λt)=sign(Rt)ReLU{|Rt|-λt}

(16)

式中， ReLU为整流线性单元。

2.2 LAMP网络训练

(17)

图2 LAMP网络结构

图3 第t层LAMP的内部结构

(18)

3 算法复杂度分析

从前文可知，衍射层析算法构建的压缩感知模型与文献[4]构建的压缩感知模型区别在于字典矩阵的构建。对于模型空间复杂度，一般字典矩阵的大小为天线数×频点数×成像网格数[4]，即矩阵维度为MNK×MNK(假设目标像O与回波数据具有相同大小)；对于本文的衍射层析模型而言，快速傅里叶变换算子的J决定了它的大小，J的维度降为MNK×K，明显小于字典矩阵所占内存。对于模型时间复杂度，字典矩阵构建时间采用大O记法可以表示为O(M2N2K)，衍射层析模型构建时间为O(MNK2)；对于凸优化算法运行时间主要体现在计算Ψ(O)上，文献[4]计算所需M2N2K2+MNK次，本文如表1所示，总体小于文献[4]所需时间。

表1 衍射层析模型运行次数

4 仿真与实验数据结果分析

4.1 仿真数据测试结果分析

图4 均方误差

下面，设验证场景为放置在空间坐标为(0.3,0.5,0.5)，(0.7, 0.5,0.5)的点目标，其对应反射系数均为0.5，信噪比为0 dB，对数据50%降采样。LAMP网络成像结果如图5所示，左上图为沿距离向z=0.5 m处方位-高度切片图，右上图为沿高度向y=0.5 m 处方位-距离切片图，左下图为沿方位向x=0.3 m处距离-高度切片图，右下图为三维立体图，可以看出在没有手工调整阈值参数的情况下，成像网络得到的重建图像去除了旁瓣和杂波，准确对点目标成像。

图5 LAMP网络成像结果

图6给出了不同算法的成像结果对比。其中图6(a)为DT算法成像结果，图6(b)为迭代软阈值算法(ISTA)算法成像结果，图6(c)～图6(e)给出了α=1,5,10时AMP算法的成像结果。与图5的成像结果对比，衍射层析算法成像中旁瓣的存在导致成像结果杂波较多；ISTA算法图像与图5像基本相同，但是ISTA算法阈值参数需要手动调整；AMP算法在不同的α下对旁瓣滤除效果明显不同，表明人工设定阈值参数将很大程度上决定最终的成像结果。

(a) DT算法成像结果

为了更直观地对算法性能进行对比，表2给出其在TCR、ENT、MFOD、PRT方面的对比，分别反映了目标杂波比、图像熵、目标的边缘轮廓特性和算法运行时间。从表中可以看出，DT算法成像的PRT值最小，程序运行时间较短，但其TCR、ENT和MFOD较大，图像聚焦度不够高，杂波较多，轮廓特性不明显；ISTA算法TCR、ENT和MFOD较小，成像质量较高，但其PRT值为691 s，反映出程序迭代次数过多，运行时间较长；AMP算法在不同α值下TCR、ENT和MFOD不同，表明了阈值参数对重建图像质量的影响；LAMP网络TCR值比手动调参下小3 dB，ENT与MFOD也小于其他算法，所成图像最清晰，此外，AMP单次迭代时间约为40 s，而LAMP单层网络使用GPU运行时间约为4 s，PRT值小于ISTA和AMP算法，表明了本文方法运行时间小于其他方法。

表2 TCR、ENT、MFOD、PRT的比较

4.2 实测数据结果分析

采用维拉诺瓦(Villanova)大学高级通信中心(CAC)与空军研究实验室(AFRL)合作收集得到的穿墙雷达成像实验数据，探测场景如图7所示[14]。选取6557平面阵列所得到的数据进行三维穿墙成像，发射信号是载波频率为2 GHz，带宽为1 GHz的步进频率波形，步进间隔为5 MHz。对数据进行预处理，使用背景相消法消除墙体杂波和噪声。使用Populated Scene场景回波数据生成训练样本，加入信噪比为-5～5 dB范围内噪声，一共生成50组样本。使用Calibrated Scene场景进行实验验证，DT算法成像结果如图8所示，LAMP网络下回波数据50%降采样成像结果如图9所示。在图8、图9中，图(a)为高度为1.2 m处切片图，图(b)为三维成像立体图。通过对比可知，DT算法存在目标周围旁瓣较多、受噪声干扰严重、数据量大等缺点；LAMP网络的成像目标周围基本不存在旁瓣干扰，随机噪声也得到了很好的抑制，同时没有损坏目标的边缘特性。表3给出了DT算法和LAMP网络所成像的TCR、ENT、MFOD和PRT值。

(a) 探测场景图

(a) 方位-距离向切片图

(a) 方位-距离向切片图

表3 TCR、ENT、MFOD、PRT的比较

5 结束语

本文提出了一种衍射层析稀疏模型下的学习成像方法，通过构建衍射层析稀疏模型降低了存储空间，并通过设计LAMP网络自动学习可调参数，提高了成像质量，适用于三维穿墙雷达成像。还需注意的是：在衍射层析稀疏模型中没有使用深度学习卷积层替换其他参数，也没有考虑墙体杂波与目标回波的分离，下一步拟对网络结构和训练样本作一定的更改，实现它们的分离。