幅度和相位估计的ISAR超分辨成像方法

， ， ，

(1.南京航空航天大学雷达成像与微波光子技术教育部重点实验室， 江苏南京 211100；2.中国人民解放军96764部队， 河南洛阳 471000)

0 引言

逆合成孔径雷达可全天候、全天时对远距离运动目标进行高分辨率成像，是一种有效的目标识别途径。RD算法是ISAR成像的常用算法，该方法首先对脉冲压缩后的雷达回波数据进行运动补偿，接着在方位向进行快速傅里叶变换(Fast Fourier Transform, FFT)实现成像。RD算法假设相干积累时间内目标相对雷达转动引起的多普勒频率恒定，因此在实际中限制了成像积累时间和方位向分辨率。在这种情况下，采用超分辨方法[1-13]，在不增加成像角度的情况下，可以获得高于理论分辨率的图像分辨率。

超分辨方法可以分为参数化法和非参数化法。参数化法要求获得散射点的个数、位置和幅度，对数据模型误差非常敏感，该类方法包括旋转不变子空间算法(Estimation of Signal Parameters via Rotational Invariance Techniques, ESPRIT)[1]、多重信号分类法(Multiple Signal Classification, MUSIC)[2-4]、基于带宽外推的线性预测法[5]、基于自适应模型的方法[6]等。非参数化法假设场景为连续的函数，不需要获取散射点的个数，包括最小方差法(Minimum Variance Method, MVM)[7-9]、空间变迹法(Spatially Variant Apodization, SVA)[10-12]、自适应旁瓣抑制法(Adaptive Sidelobe Reduction, ASR)[13]等。实际中，精确数据模型的获取相当困难，所以非参数化法的适用性更强，但MVM法频谱估计不够精确，降低了算法的准确性；SVA法对信噪比和频谱的平坦程度有较高要求；ASR法的性能较好，但是运算量较大。

本文将非参数化的APES法用于ISAR方位向和距离向成像，给出具体的成像算法，并通过仿真和实测数据处理分析该成像方法的性能。仿真结果和实测ISAR数据成像结果验证了该算法可以获得超分辨，优于RD成像结果，并可降低旁瓣，提高对比度和信噪比，成像效果与MVM法相当。

1 基于RD的ISAR成像方法

成像雷达多采用线性调频信号作为发射波形。在“停-走-停”假设下，经过解调和距离压缩后，目标上某一散射点回波信号可写为

(1)

式中，t为距离向快时间，τ为方位向慢时间，kr为发射波形的调频率，Ap为包含散射点散射率和天线方向图调制的回波幅度因子，rp(τ)表示τ时刻散射点和雷达之间的距离，Tcp表示压缩后的脉冲宽度，λ为雷达波长。

ISAR成像中，目标与雷达间的相对运动通常被分解为目标上某一参考点相对于雷达的平动和目标上散射点绕该参考点的转动。因此，散射点和雷达之间的距离[14]可表示为

rp(τ)=r0(τ)+(Ω(τ)×r(τ))·R(τ)

(2)

式中，r0(τ)为τ时刻雷达和参考中心点的距离，R为雷达视线方向的单位矢量，r为由参考中心点指向散射点的矢量，Ω为目标相对于参考中心点的转动矢量。平动分量r0(τ)对ISAR成像没有贡献，相反会影响成像，需要对其补偿，即所谓的ISAR运动补偿。

将式(2)代入式(1)中，并假设运动补偿已完成，可得

s(t,τ)=

(3)

在小转角假设下，散射点所在距离单元内的方位向回波信号sc(τ)[15]可以近似表示为

(4)

式中，Ω0为有效转动角速率，将Ω分解为与R正交的分量即可得到有效转动矢量Ω0，y为散射点在ISAR图像投影平面(Image Projection Plane, IPP)内的方位向坐标。

由式(4)可见，在小转角假设下，散射点回波方位向为单频信号，频率为fd,y=2Ω0y/λ，与散射点在IPP内的方位坐标一一对应。因此，对每个距离单元的回波信号作傅里叶变换，即频谱分析，即可在方位向对目标进行成像。这就是传统ISAR成像广泛使用的距离-多普勒成像方法的原理。

2 基于APES的ISAR成像方法

由上节可知，ISAR方位向成像的本质就是频谱分析。APES是一种自适应有限冲击响应(Finite Impulse Response, FIR)方法[16]，它依据最小二乘估计的准则对输入的数据进行处理，构造出相应的有限冲击响应滤波器来完成谱估计，可提供更准确、副瓣更低的超分辨性能。在本节给出基于APES频谱估计的ISAR成像方法来提升ISAR成像质量。

假设目标由Nd个散射点构成，且考虑噪声，由上节给出的ISAR成像模型可知，距离压缩后，每一距离单元内第n个脉冲信号可表示为

(5)

式中，Sc(n)=sc(nTr)为sc(τ)的离散形式(注意这里已扩展为多个散射点的情形)，Tr为脉冲重复周期，ωi为散射点对应的数字多普勒角频率：

(6)

式中，yi为第i个散射点在图像投影平面方位向的坐标。式(5)中，Ap,i是与第i个散射点散射率相关的回波幅度，e(n)为噪声和干扰。

APES成像方法可归纳为如下步骤：

m=0,1,…,N-M

(7)

式中，zm表示Sc(n)的第m个采样值，(·)T表示转置。

m=0,1,…,N-M

(8)

式中，(·)*表示共轭。

(9)

(10)

式中，(·)H表示共轭转置。为了提高估计精度，由式(9)和式(10)可以进一步获得前后向采样协方差矩阵[17]：

(11)

则噪声和干扰的协方差矩阵[18]为

(12)

式中，

(13)

(14)

式中,L=N-M+1。

4) 构造长度为M的FIR滤波器。

根据近似最大似然法，FIR滤波器的脉冲响应函数[14]选为

(15)

5) 获得方位向回波频谱的最小二乘估计[18]：

(16)

对距离向也可以应用上述步骤进行超分辨成像，处理过程与方位向类似。首先将脉压后的一维距离像进行逆傅里叶变换，返回数据域。然后构造自适应滤波器来完成距离向幅度相位估计，具体过程这里不再赘述。APES法对雷达回波的处理流程如图1所示。

3 RD，MVM，APES对比

因此，RD成像方法是MVM方法和APES方法的一种特殊情况。MVM方法和APES方法均依赖于输入数据，属于自适应滤波方法。构造滤波器使用的协方差矩阵的不同导致MVM和APES二者性能的差异。

4 实验验证

本节采用RD法、MVM法和APES法对点目标仿真数据、卫星目标仿真数据，以及飞机和舰船的实测ISAR数据进行处理，验证APES成像算法有效性和优越性。为了更好地评价成像方法的性能，采用两种定量评估图像的指标。

1) 图像对比度为图像强度的标准差和均值之间的定量关系，定义如下[20]：

(17)

式中，A{·}为空间平均算子，I为图像的幅度。对比度可以反映图像的聚焦程度，对比度越大表示图像聚焦效果越好。

2) 图像信噪比的计算方法如下[21]：

(18)

式中，α表示目标区域，β表示背景区域，Nα表示目标区域的像素个数，Nβ表示背景区域的像素个数。

4.1 点目标仿真

为了验证算法的正确性，采用RD法、MVM法和APES法分别对点目标的仿真数据进行处理，并对点目标的成像质量进行比较。

为了便于对成像结果进行分析，ISAR点仿真采用经典的转台模型。ISAR点仿真的参数为脉冲宽度10 μs，信号带宽200 MHz，采样频率480 MHz，工作波长0.03 m，径向距离5 000 m，脉冲重复频率200 Hz，点目标的转动角速度为0.02 rad/s。采用RD法、MVM法和APES法对点仿真进行成像的结果如图2所示。

表1给出了点目标距离向和方位向的积分旁瓣比(Integrated Side Lobe Ratio, ISLR)和峰值旁瓣比(Peak Side Lobe Ratio, PSLR)。

表1 点目标成像质量比较

表2给出了点目标的对比度和信噪比。

表2 点目标的定量分析结果

经过计算得到的距离向和方位向的理论分辨率分别为0.75 m和0.3 m，而APES法的距离向和方位向的分辨率分别为0.492 m和0.178 m，因此APES法的分辨率超越了理论分辨率，算法的正确性得到验证。表1说明了相比于RD法和MVM法，APES法能更好地抑制旁瓣。由图像性能指标可以看出，APES法成像的图像质量最好。

4.2 卫星仿真

仿真所用“长曲棍球”卫星模型共923个散射点，反射系数为单位值，采用卫星工具包(Satellite ToolKit, STK)软件导出卫星轨道，卫星模型如图3所示。

卫星仿真时发射信号带宽为1 GHz，雷达信号波长为0.03 m，距离向采样率为1.21 GHz，脉冲重复频率为400 Hz。采用RD法、MVM法和APES法对卫星进行成像的结果如图4所示。

表3给出了卫星数据的对比度和信噪比结果。

表3 卫星数据的定量分析结果

经过计算，卫星距离向和方位向的理论分辨率分别为0.15 m和0.159 m，而APES法的距离向和方位向的分辨率分别为0.134 m和0.115 m。说明本文的算法可以实现超分辨成像。从卫星的成像图中可以看出，与RD法相比，MVM法和APES法得到的卫星的主体与天线罩成像效果较好，天线的聚焦效果也更好，且APES法得到的图像旁瓣比MVM法得到的更低。由图像对比度和信噪比可以看出，APES法成像的图像质量与MVM法相当。

4.3 实测数据处理

为了验证APES方法的有效性和优越性，本节分别对飞机和舰船的ISAR实测数据进行处理，并对成像结果进行分析。

使用工作在C波段的地面雷达进行飞机数据采集，发射波的带宽为400 MHz，成像积累时间为1.28 s。使用工作在X波段的岸基雷达进行舰船数据采集，发射波的带宽为170 MHz，成像积累时间为0.4 s。采用全局最小熵法进行距离对准，相位梯度自聚焦法进行相位补偿。

采用RD法、MVM法和APES法对飞机和舰船数据的成像结果分别如图5和图6所示。

从图中可以看出，相比传统的RD法，MVM法和APES法聚焦更好，旁瓣明显降低，APES法与MVM法的成像效果相当。

经过计算得到，飞机和舰船理论的距离向分辨率分别为0.375 m和0.882 m，APES法得到的距离向分辨率分别为0.342 m和0.735 m，超过了距离向的理论分辨率。

由于实测数据缺少方位向的转角信息，无法定量比较各种方法所得图像的方位向分辨率，因此图7和图8分别给出了飞机和舰船实测数据的强散射点方位向的3 dB主瓣宽度图，从而定性比较分辨率。从图中可以看出，相比于RD法，APES法的分辨率提高。

表4和表5给出了飞机和舰船成像结果的对比度和信噪比分析。

表4 飞机实测数据的定量分析结果

表5 舰船实测数据的定量分析结果

从表中可以看出， APES法重建图像的对比度和信噪比最大，说明APES法成像质量在三者中较优。

5 结束语

本文给出了一种基于APES的ISAR超分辨成像算法。该算法对雷达回波数据脉压、运动补偿后进行处理，采用空间平滑法得到前后向采样协方差矩阵，再通过构造脉冲响应函数重建目标散射点。仿真结果、飞机和舰船的实测数据成像结果表明，APES法可以抑制旁瓣，改善成像效果，实现ISAR超分辨成像，成像效果与常用的MVM方法相当。