1、第4期2023年8月Vol.21 No.4August 2023雷达科学与技术Radar Science and TechnologyDOI:10.3969/j.issn.16722337.2023.04.008基于改进PFA的车载大斜视SAR成像算法王佳慧，刘衍琦，王纪平，金秋，毛新华（南京航空航天大学电子信息工程学院雷达成像与微波光子技术教育部重点实验室，江苏南京 211106）摘要：为了应对车载毫米波雷达大斜视成像困难的问题，本文提出一种改进的极坐标格式算法（PolarFormat Algorithm,PFA）对条带车载毫米波雷达斜视回波进行基于子孔径拼接的成像。该算法从条带数据与聚束数
2、据的特点出发，将全孔径回波划分为子孔径，利用PFA处理子孔径数据。由于PFA存在波前弯曲误差，子图像不能直接拼接，因此对每一幅子图像进行几何失真校正。同时，以重叠子孔径的划分方式保证成像结果的高分辨率。最后截取子图像进行拼接得到条带SAR成像结果。所提方法解决了车载毫米波雷达大斜视情况下两维耦合严重的问题。通过对点目标和实测数据的仿真与分析验证了所提方法的有效性。关键词：车载毫米波雷达斜视SAR；子孔径划分；极坐标格式算法；重叠子孔径；子图像拼接中图分类号：TN959文献标志码：A文章编号：16722337（2023）04041109Automotive High Squint SAR Ima
3、ging Algorithm Based on Improved PFAWANG Jiahui,LIU Yanqi,WANG Jiping,JIN Qiu,MAO Xinhua（Key Laboratory of Radar Imaging and Microwave Photonics,Ministry of Education,College of Electronic and Information Engineering,Nanjing University of Aeronautics and Astronautics,Nanjing 211106,China）Abstract:Th
4、is paper proposes an improved polar format algorithm（PFA）,which is used to process the highsquint stripmap automotive millimeter wave radar echoes.Based on the characteristics of stripmap data and spotlight data,the fullaperture is divided into subapertures,and the subapertures are processed by PFA.
5、Due to the wavefrontbending error,subimages cannot be mosaicked directly,so geometric distortion correction is carried out for each subaperture.The high resolution is ensured by overlapped subapertures.Finally,subimages are intercepted for mosaickingto obtain stripmap SAR imaging result.The algorith
6、m solves the severe problem of twodimensional coupling in highsquint imaging of stripmap SAR.Simulation and processing of the point targets and the measured data demonstrate theeffectiveness of the proposed method.Key words:automotive millimeter wave radar squint SAR;subaperture division;polar forma
7、t algorithm（PFA）;overlapped subapertures;subimages mosaicking0引言随着智能驾驶辅助系统（Advanced Driver Assistance System,ADAS）的兴起，车载毫米波雷达在保障驾驶员和乘客的安全性上起到不可替代的作用。车载毫米波雷达通过接收回波信号感知周围环境，是实现自适应巡航控制、自动紧急刹车等功能必不可少的传感器1。与摄像头相比，毫米波雷达不受环境限制，在雨雪等复杂环境下也可发挥作用。低成本、高可靠性的特点使得毫米波雷达更胜一筹。基于以上优势，再加之技术趋于成熟，毫米波雷达已被装载到越来越多的汽车上投入使用2。
8、实际应用中为了保证行车安全，不仅要对正侧方场景进行观察，同时要将天线倾斜以方便对斜视方向的情况进行探测。雷达斜视时，距离向与方位向两维耦合严重，距离多普勒算法（RangeDoppler Algorithm,RDA）存在近似，斜视会使误差变大，成像困难。后向投影算法（BackProjectionAlgorithm,BPA）适应于斜视情况，但是考虑到车载环境下对算法的实时性要求很高，而BPA计算复杂，效率低，也难以应用。与之相比，PFA在时域进行运动补偿，适应于非线性航迹，且成像实时性好，因此得到广泛应用。但其只适应于小场景下收稿日期：20221101；修回日期：20221209基金项目：国家自然
9、科学基金（No.62071225）；航空科学基金（No.2020Z017052001）；江苏省自然科学基金（No.BK20211181）雷达科学与技术第 21 卷 第 4 期的聚束数据的处理，因此在车载环境下需要做出改进。市场的扩大和用户需求量的增加对车载雷达提出了更高的要求，因此车载毫米波雷达成像已经得到了社会各界研究者的广泛关注。将多传感器融合技术用于车载SAR系统是一种比较先进的方式，但是该方法高度依赖于惯导信息，成本较高3。有关正侧视的车载毫米波雷达成像研究较多，也都获得了较好的结果和应用，比如短距情况下的人体和车辆检测等46。BPA适应于斜视，因此关于车载毫米波雷达斜视成像的处理应用
10、BPA较多，但是考虑到实时性则需要做出进一步的研究78。为了解决当前车载环境下斜视成像困难的问题，在综合分析现有各种雷达成像算法特点的前提下，本文提出了一种改进的极坐标格式算法，将全孔径条带SAR回波划分为子孔径，对子孔径数据进行PFA处理和几何失真校正，将子图像进行截取、拼接得到条带SAR成像结果。同时，一定条件下以重叠子孔径的方式划分全孔径保证成像的高分辨率。通过用所提算法对点目标和实测数据进行仿真处理，目标均成像在正确的位置上，并得到了完整的高分辨率图像。1基于子孔径拼接的极坐标格式成像算法1.1子孔径划分车载毫米波雷达以条带模式采集数据，其特点是天线指向不变，波束照射区域与雷达同向运动
11、，能够获取较大场景的回波。聚束模式下天线波束指向时刻变化，波束照射区域固定，获取数据的场景较小。PFA适用于小场景下的聚束模式，原因如下：一是 PFA 存在平面波假设，存在几何失真，且场景越大误差越大，而通常车载环境的场景较近，波前弯曲误差更加严重；二是PFA以场景中心作为参考点进行运动补偿，距离中心点越远误差越大。为了使用PFA处理条带数据，需要先将条带数据转换为聚束数据。转换示意图如图1所示，车辆从Q行驶到Q的过程中，对应长度为Lg的阴影部分始终处于波束照射下，符合聚束数据的特点。条带合成孔径长度为L，等效聚束合成孔径长度为Ls，斜视角为s。当车辆行驶速度也即雷达运动的速度为Vr，子孔径回
12、波个数为Naa时，Ls=VrNaaPRF，要求L Ls。由于毫米波雷达波束较宽，对其进行滤波后得到方位向波束宽度为。YXQQLssLgOLR0图1条带数据转换为聚束数据的几何模型子孔径脉冲个数Naa的选择首先要满足相邻子孔径间无间隔，其次为了使成像结果看起来无拉伸感，应使距离向与方位向分辨率相匹配。1.2子孔径数据的极坐标格式算法处理毫米波雷达发射线性调频连续波（Linear Frequency Modulation Continuous Wave,LFMCW）信号为s()tr,ta=A rect()trTr rect()taTa exp()jkrtr2exp()j2 fctr（1）式中A表示
13、信号幅度，kr为距离向线性调频斜率，fc为信号载频，tr表示距离向快时间变量，ta表示方位向慢时间变量，Tr为脉冲宽度，Ta为目标在波束下的照射时间。忽略信号传输过程中的幅度变化，设回波延迟为0，则经过解调后天线接收到的回波信号为s()tr,ta=A rect()tr-0Tr rect()taTaexp()jkr()tr-02 exp()-j2fc0（2）毫米波雷达带宽较大，达到几个GHz，同时在车辆行驶过程中，探测场景相对较小，为了降低带宽、简化AD采样、降低硬件要求，车载毫米波雷达通过在模拟阶段使用 Dechirp 的方式进行脉冲压缩，脉冲压缩后的表达式为412王佳慧：基于改进PFA的车载
14、大斜视SAR成像算法2023 年第 4 期s()tr,ta=Arect()tr-0Trrect()taTaexp()j2krtr0exp()-jkr02 exp()j2 fc0（3）将式（3）进行距离向FFT，得到s()fr,ta=A Tr sincTr()fr-kr0 rect()taTaexp()-j2fr0 exp()jkr02 exp()j2fc0（4）式中fr为距离频率，sinc 函数的峰值位于fr=kr0处，经过残余视频相位去除（RVP）后，得到s()fr,ta=A Tr sincTr()fr-kr0rect()taTa exp()j2fc0（5）忽略俯仰角，图2为车载毫米波雷达子
15、孔径在斜视情况下采集数据的几何模型，以场景中心O为坐标原点，平行于雷达运动方向为X轴，垂直于雷达运动方向为Y轴建立坐标系。图中P(xt,yt)为目标位置，s为斜视角，Rt和Ra分别为雷达天线相位中心到目标和场景中心的距离，雷达方位角为，雷达航迹到X轴的距离为y0。YQYXXsRtRaP(xt,yt)(xa,ya)OQrty0yc图2子孔径采集数据的几何模型将式（5）作距离向IFFT并忽略幅度得到s()tr,ta=exp()j4()fc+krtrRtc（6）对每一个子孔径的场景中心单独做运动补偿后s()tr,ta=exp()j4()fc+krtrRc（7）式中R=Rt-Ra为雷达天线相位中心到目
16、标与场景中心的距离差。在平面波假设条件下R=Rt-Ra xtsin+ytcos（8）则波数域的信号表示为S()Kx,Ky=exp()jxtKx+jytKy（9）式中，Kr=4()fc+krtrcKx=KrsinKy=Krcos（10）由式（9）及式（10）可知，进行两维FFT后即可得到成像结果，但雷达信号在波数域是以极坐标格式进行的非均匀采样，需要先经过距离向和方位向插值实现两维均匀重采样后才能通过两维FFT实现聚焦。1.3几何失真校正在平面波条件下用式（8）近似，但雷达波前为球面波，所以平面波近似会带来波前弯曲误差。随着对分辨率要求的提高、成像场景的扩大和成像距离的减小，波前弯曲误差愈加不可
17、忽略9。将式（9）进行两维Taylor展开得到S()Kx,Ky=expjKr()Rt-Ra=expja01Kx+a10()Ky-Kyc+（11）式中Kyc=4fcc为Ky的偏置，c表示电磁波在空气中的传播速度。在图2中，建立了XOY平面直角坐标系，对其进行旋转得到新的坐标系XOY，旋转角度为s，雷达位置为(xa,ya)，孔径中心时刻雷达到X轴的距离为yc。坐标旋转如下：()XY=()coss-sinssinscoss()XY（12）为了得到Taylor展开式的系数，用表示所有项，在图2中有Rt=x2a+y2aRa=()xa-xt2+()ya-yt2（13）413雷达科学与技术第 21 卷 第
18、4 期同时在橙色三角形中Vrtasin=yccos(s+)，得到 xa=yctancosscoss-sinstanya=yc+xatans=yc+yctansinscoss-sinstanKxKy=tan（14）将式（13）和式（14）代入式（10）得到相应的系数为a10=rco-rcta01=rcorctxt+yc-rco()yc-ytrcttans（15）式中rco=yc，rct=x2t+(yc-yt)2，分别为孔径中心时刻天线相位中心到场景中心和目标的距离，由此得到目标在PFA图像中的位置(a01,a10)，通过插值找到其“像”，实现几何失真校正1011。至此完成了子孔径数据的PFA处理
19、及其几何失真校正，步骤如图3所示。5C.IFFT4FFTD=C.,Ds(tr,ta)=exp()-j4(fc+kr0)RacS(Kx,Ky)=exp(jxtKx+jytKy)()XY=()coss-sinssinscoss()XYa10=rco-rcta01=rcorctxt+()yc-rco(yc-yt)rcttans图3子孔径成像及几何失真校正过程1.4子孔径图像拼接由上述分析可知，子孔径回波个数的选择取决于波束宽度和分辨率12。将子孔径回波数据利用PFA成像并进行几何失真校正后得到子图像，需要按照顺序将子图像拼接起来得到全部区域的成像。子孔径的划分与拼接详细介绍如下。如图4（a）所示。l
20、2与l3之间为实际成像需要的部分。当LsLe时，直接截取对应的虚线框内的成像场景部分进行拼接，拼接示意图如图4（b）所示，其中W1=Lsx，对应的是雷达也即车辆在子孔径数据时间内运动的距离，x为坐标系内相邻点的横向方位间隔。QQQLLsLel1l2l3s（a）子孔径划分示意图W1123K（b）子孔径拼接示意图图4Ls Le时子孔径划分与拼接示意图当Ls Le时，为了在实现连续拼接的同时保证分辨率，采用重叠子孔径的方法13。重叠子孔径的划分方式如图 5（a）所示，第一个子孔径为Q1Q1，第二个子孔径为Q2Q2，其回波个数为Naa，长度为Ls=VrNaa/PRF。设置两个子孔径间重叠的回波个数为N
21、aa，则对应的长度为Lc=VrNaa/PRF，那么相邻两个子孔径间的重叠率rate=Naa/Naa=Lc/Ls。子孔径拼接如图5（b）所示，其中W2=(Ls-Lc)/x，对应的是未重叠部分雷达运动的距离。设全孔径回波信号个数为Na，由上述分析可知，当子孔径无重叠时，划分的子孔个数为K1=414王佳慧：基于改进PFA的车载大斜视SAR成像算法2023 年第 4 期Na/Naa。当子孔径重叠时，划分的子孔径个数为K2=()Na-Naa/()Naa-Naa+1，式中数学符号代表向上取整。基于子孔径拼接的条带车载毫米波雷达斜视SAR全孔径成像流程如图6所示。LsLcLeLsQ2Q1Q1Q2l1l2l3
22、W2123K（a）重叠子孔径划分示意图（b）重叠子孔径拼接示意图图5Ls Le时重叠子孔径划分与拼接示意图DPFA,D4+F(E45+%(DD(EPFA*LU,图6基于子孔径拼接的条带车载毫米波雷达斜视成像过程成像效率也是车载环境中需要考虑的重要因素。利用文献 14 的方法进行成像复杂度的计算如表1所示，计算复杂度主要来自于成像过程中的复乘运算，其中M为全孔径回波脉冲个数，M为一个子孔径中的脉冲个数，K为划分的子孔径数，N为每个脉冲的采样点数。表1计算复杂度操作DechirpRVP距离向FFT/IFFT方位向FFT运动补偿距离向插值方位向插值几何失真校正计算复杂度MNMNlog2N+MN2K
23、Nlog2NK Mlog2MK MK Nlog2NK Mlog2MK (Mlog2M+Nlog2N)2仿真和实测数据处理2.1点目标仿真为了验证算法的有效性，设置3个点目标位置分别为（-5，8），（0，10.5），（2，13），利用所提算法进行仿真后得到的结果如图7所示，点目标放大结果和距离向及方位向剖面图如图8、图9所示。可以看出3个点目标经过PFA成像及几何失真校正后均在正确位置上，误差很小，可以忽略不计，说明上述基于子孔径拼接的成像方法是有效的。4C./m/m-10252015105XY-4.96781 8.02485XY 2.06386 12.9997XY 0.0548149 10.4
24、3690-505图7点目标仿真结果-6.0-5.5-5.0-4.5-4.0m7.07.58.08.59.09.5C.m（a）（-5,8）放大结果415雷达科学与技术第 21 卷 第 4 期-1.0-0.500.51.09.510.010.511.011.5mC.m（b）（0，10.5）放大结果1.01.52.02.53.012.012.513.013.514.014.5mC.m（c）（2，13）放大结果图8点目标放大结果01002003004005006007002-100-80-60-40-200L0100200300400500600700C.2-100-80-60-40-200dB dB
25、 C.L（a）（-5,8）剖面图0100200300400500600700-100-80-60-40-2000100200300400500600700-100-80-60-40-2002LC.2C.LdB dB（b）（0,10.5）剖面图0100200300400500600700-100-80-60-40-2000100200300400500600700-100-80-60-40-2002LC.2dB dB C.L（c）（2,13）剖面图图9方位向和距离向剖面图2.2面目标仿真设置如图10（a）所示的面阵，相邻点间方位向间隔和距离向间隔均为2 m，距离向长度共6 m，方位向长度14 m
26、，共计28个点目标。其成像结果如图10（b）所示，方位向回波脉冲数共有16 384个，子孔径脉冲设置为1 024个，划分为16个子孔径，其中第1，3，5，7个子孔径几何失真校正前后的成像结果如图 11和图 12所示。从图 11可以看出，PFA成像使用的平面波假设与实际球面波不符，存在较大的误差，使成像结果扭曲、变形、错位，又由于车载情况下一般是对近场景成像，平面波假设带来的影响更大，因此要弥补这些误差。如图12所示，经几何失真校正后误差几乎被补偿掉了，目标成像结果与设置相同。YXO2 m2 m（a）面阵（b）成像结果图10面阵及其成像结果416王佳慧：基于改进PFA的车载大斜视SAR成像算法2
27、023 年第 4 期(a)01(b)03(c)05(d)07图11经PFA后子孔径成像结果(a)01(b)03(c)05(d)07图12几何失真校正后子孔径成像结果2.3实测数据处理本文实测数据使用IWR1843毫米波雷达模块发射和接收毫米波，经DCA1000高速捕获原始数据后通过网线传给计算机。在实际应用场景中，IWR1843毫米波雷达使用的参数如表2所示。表2实测数据参数参数脉冲间隔载频PRF带宽采样频率参数值60 s77 GHz2 000 Hz3 240 MHz10 MHz首先在校内停车坪采集了一组斜视数据，斜视角为 30，车速为 10 km/h，全孔径共划分为 113个子孔径，每个子孔
28、径包含512个脉冲。图13和图14为第9195个子孔径几何失真校正前后的成像结果对比。由图13可以看出，由于子孔径回波在经过PFA处理后存在一定角度的旋转以及几何位置的失真，导致图像混乱，目标不在正确的位置上，不能直接拼接。经几何失真校正后，目标位置正确，子图像可连续拼接。（a）第91个子孔径（b）第92个子孔径（c）第93个子孔径（d）第94个子孔径（e）第95个子孔径图13几何失真校正前子孔径成像结果（a）第91个子孔径（b）第92个子孔径（c）第93个子孔径（d）第94个子孔径（e）第95个子孔径图14几何失真校正后子孔径成像结果图15为条带实测数据成像结果及其与实际成像场景的对应情况。
29、在图15中可以看到，经过子孔径划分、PFA成像、几何失真和子孔径拼接后，得到了停车坪连续完整的成像结果，车辆和草坪均成像在正确的位置，路灯（红色框内）、指示牌（黄色框内）清晰可见，实现了良好的聚焦。由于实现过程中滤波后的波束仍然较宽，所以未重叠子孔径。417雷达科学与技术第 21 卷 第 4 期图15校内停车坪成像结果与实际场景的对应情况图16为在校外公路采集的一组斜视数据的成像结果及其与实际场景的对应情况，将波束滤波较窄，使用重叠子孔径的方式对其进行成像。图16全孔径回波脉冲共有95 342个，成像时子孔径回波脉冲为 512 个，重叠率rate 52%，即每两个相邻子孔径间有约270个脉冲相
30、互重叠。成像之后，路边车辆均停放在其实际所在位置上，说明了算法的有效性。图16校外公路成像结果与实际场景的对应情况3结束语本文提出了一种改进的 PFA算法，依据聚束数据与条带数据间的联系，将全孔径数据划分为子孔径，利用PFA进行子孔径成像，并通过几何失真校正解决了成像扭曲变形的问题。同时，为了保证成像的高分辨率，在一定条件下以重叠子孔径的方式划分全孔径。根据点目标及实测数据的成像结果可以看出，该算法能够解决车载毫米波雷达大斜视时两维耦合严重的问题，在车载毫米波雷达成像领域中具有良好的应用前景。本文重在研究成像算法，虽然上述分析均基于雷达匀速运动进行，但是成像结果满足当前分辨率的要求。随着对分辨
31、率要求的提高，实际行车过程中由于速度不均而带来的目标散焦等问题将在后续工作中进行深入研究。参考文献：1朱玉成.智能驾驶汽车毫米波雷达应用分析 J.中国新通信，2021，23（2）:9899.2韩宝石，王峥.车载毫米波雷达国内外发展现状综述J.数字通信世界，2019（9）:1516.3刘明亮.基于多传感器融合的智能车道路目标检测D.镇江：江苏大学，2021.4TAGLIAFERRI D,RIZZI M,NICOLI M,et al.NavigationAided Automotive SAR for High Resolution Imaging ofDriving Environments J
32、.IEEE Access,2021,S（9）:3559935615.5PANDEY N,RAM S S.Automated Parking Test Using418王佳慧：基于改进PFA的车载大斜视SAR成像算法2023 年第 4 期ISAR Images from Automotive Radar C/2021 IEEE International Conference on Autonomous Systems,Montral,Qubec,Canada:IEEE,2021:15.6TAVANTI E,RIZIK A,FEDELI A,et al.A ShortRangeFMCW Rada
33、rBased Approach for MultiTarget HumanVehicle Detection J.IEEE Trans on Geoscience and Remote Sensing,2022,60:116.7IQBAL H,LFFLER A,MEJDOUB M N,et al.RealisticSAR Implementation for Automotive ApplicationsC/2020 17th European Radar Conference,Utrecht,Netherlands:IEEE,2020:306309.8GAO Xiangyu,ROY S,XI
34、NG Guanbin.MIMOSAR:AHierarchical HighResolution Imaging Algorithm for MMWave FMCW Radar in Autonomous DrivingJ.IEEETrans on Vehicular Technology,2021,70（8）:73227334.9薛鹏凤，郭苹.基于改进极坐标格式算法的大斜视SAR成像方法 J.空军工程大学学报（自然科学版），2022，23（3）:8288.10毛新华，朱岱寅，叶少华，等.一种基于图像后处理的极坐标格式算法波前弯曲补偿方法 J.电子学报，2010，38（1）:244250.11陈
35、旭，张大炜.大斜视角SAR的方位多级重叠子孔径算法 J.无线电通信技术，2013，39（4）:6063.12叶铮，朱岱寅，吴迪.一种基于重叠子孔径回波信息的SAR图像配准算法 J.航空学报，2023，44（8）:196207.13TANG Yu,ZHANG Bo,XING Mengdao,et al.AzimuthOverlapped Subaperture Algorithm in Frequency Domain for Highly Squinted Synthetic Aperture Radar J.IEEE Geoscience and Remote Sensing Letters
36、,2013,10（4）:692696.14WALTER G,RONALD M,RON S.Spotlight SyntheticAperture RadarM.Boston:Artech House Publishers,1995:479481.作者简介：王佳慧女，1999年生，江苏徐州人，南京航空航天大学硕士研究生，主要研究方向为车载毫米波雷达信号处理。刘衍琦男，1997年生，湖南岳阳人，南京航空航天大学博士研究生，主要研究方向为合成孔径雷达成像技术。王纪平男，1999年生，浙江杭州人，南京航空航天大学硕士研究生，主要研究方向为合成孔径雷达成像技术。金秋男，1998年生，江苏镇江人，南京航空航天大学硕士研究生，主要研究方向为高空高速平台雷达成像。毛新华男，1979 年生，湖南涟源人，博士，教授、博士生导师，主要研究方向为合成孔径雷达成像技术。陶海红女，1976年生，陕西宝鸡人，西安电子科技大学雷达信号处理全国重点实验室教授，主要研究方向为雷达信号处理、先进阵列处理、网络化雷达系统。畅鑫男，1990年生，河北石家庄人，西安电子科技大学博士研究生，主要研究方向为雷达信号对抗分析、电磁信号安全研究、多智能体强化学习。（上接第410页）419
