
Citation: | Tongjun WANG, Feng WU, Wei XU. A Novel Imaging Approach for Improving Azimuth Angular Resolution of Automotive Radars[J]. Journal of Electronics & Information Technology, 2020, 42(8): 2037-2044. doi: 10.11999/JEIT190618 |
As the azimuth angular resolution is limited by the antenna length in automotive radars, a novel imaging approach for improving azimuth angular resolution of automotive radars is proposed based on multi-beam real-aperture radar images combination processing. Firstly, the antenna beam of the phased array antenna is electronically scanned to obtain forward-looking real-aperture radar images. Afterwards, multiple real-aperture radar images are coherent accumulated according to the imaging geometry of automotive radar to improve azimuth angular resolution. Simulation results validate the proposed imaging approach to improve the azimuth angular resolution of automotive radar.
随着现代工业技术的发展,交通安全问题已成为世界各国亟需解决的问题,而交通安全问题主要由于驾驶员们的疏忽大意和无法及时反应所导致。无人驾驶技术通过车载传感系统感知道路环境,自动控制车辆的转向和速度,从而使车辆能够安全、可靠地在道路上行驶[1-3]。前视微波成像雷达是车载传感系统的重要组成部分,与超声波、摄像头、激光雷达相比,微波雷达不仅具有高距离分辨率、高角度分辨率、高速度分辨率,还具有不受恶劣天气及极端光线影响的优势[4-9]。
微波前视成像雷达距离向通过宽带信号脉冲压缩技术来获取高距离分辨率,方位向通过实孔径成像或数字阵列合成处理来获取较高的方位分辨率[10-13]。受雷达角分辨瑞利限的约束,要提高雷达方位角分辨率就需要提高雷达工作频率和雷达方位孔径天线长度,这些给低成本微波成像雷达的实现带来了较大的工程实现难度[2]。数字阵列超分辨处理技术也能够提升雷达方位角分辨率,但阵列超分辨处理技术运算量大一般难以满足汽车雷达实时性的要求。
针对汽车雷达方位角分辨率受方位向天线长度限制的问题,本文提出一种基于多波束实孔径雷达图像融合来提升汽车雷达方位角分辨率的成像方法,本方法通过多波束实孔径雷达图像相参合成来提升雷达方位分辨率。从本质上讲,本文方法利用合成孔径技术提升了雷达方位角分辨率。相比于传统合成孔径处理,本文方法的处理可以在实孔径雷达图像域进行,而且雷达数据合成不受单个雷达指向波束条件约束。因此,本文方法更加有利于提升雷达的方位角分辨率。
本文首先介绍了提升汽车雷达方位角分辨率的成像方式和雷达回波信号模型;然后根据汽车雷达成像方式和信号模型,提出了一种提升汽车雷达方位角分辨率的成像方法;最后通过计算机仿真验证了本文所提方法的准确性。
车载微波前视成像雷达的距离向与汽车平台的运动方向一致,而方位向与汽车平台的运动方向垂直。雷达距离向通过发射大时宽-带宽信号并利用脉冲压缩技术实现距离向的高分辨率,方位向通常通过实孔径成像或数字阵列合成处理来获取较高的方位分辨率。本文所提车载微波前视成像雷达工作方式如图1所示,汽车雷达首先通过方位窄波束匀速扫描实现对前方区域的观测,雷达波束扫描范围为
由于汽车雷达最大作用距离一般要小于300 m[2],因此通常采用收发分置的连续波雷达工作体制。假设车载毫米波前视成像雷达发射连续波可以写成
s(t,τ)=wa(θs(t))⋅rect[ττp]exp(j2πfcτ+jπKrτ2) | (1) |
其中,
如图1所示,假设车载前视成像雷达在t时刻,点目标(x, y)对应的雷达回波可以写成
sr(t,τ)=wa(θs(t)−θ(t))⋅rect[τ−2R(t)/cτp]⋅exp[j2πfc(τ−2R(t)/c)+jπKr(τ−2R(t)/c)2] | (2) |
其中,
R(t)=√x2+(y−vt)2 | (3) |
θ(t)=arctan(xy−vt) | (4) |
其中,c为光速,
sref(t,τ)=rect[τ−2Rref/cτref]⋅exp[j2πfc(τ−2Rref/c)+jπKr(τ−2Rref/c)2] | (5) |
此时,去斜后的雷达回波信号可以写成
sr,1(t,τ)=sr(t,τ)⋅s∗ref(t,τ)=wa(θs(t)−θ(t))⋅rect[τ−2R(t)/cτp]⋅exp[−j4πKrc(τ−2Rref/c)⋅ΔR(t)−j4πfcc⋅ΔR(t)+j4πKrc2ΔR2(t)] | (6) |
其中,
ΔR(t)=R(t)−Rref | (7) |
由式(7)可以看出,连续波信号经过去斜处理后,距离信号为快时间
ρr=0.886c2⋅(1τp⋅1Kr)=0.886c2Br | (8) |
此外,雷达的方位向角分辨率则由方位向双程天线方向图
为了提升垂直于视线方向的几何分辨率,这里采用合成孔径技术来提升方位几何分辨率。如图2所示,汽车雷达观测目标的合成孔径长度为
以目标观测中心时刻为基准,根据合成孔径原理,目标垂直于雷达中心时刻观测视线方向的几何分辨率
ρa=0.886λ2(θN−θ1)=0.443λarctanyx−Ls/2−arctanyx+Ls/2 | (9) |
此时,目标对应的方位角分辨率为
根据上述车载微波前视成像雷达利用多角度合成孔径技术来提升方位角分辨率的工作方式和雷达回波信号模型,图3给出了本文所提成像方法的处理流程图,本成像方法主要包含距离向处理和方位向处理两大处理步骤。其中,通过距离向脉冲压缩处理即可得到实孔径雷达图像,通过方位向多张实孔径雷达图像相参累加即可获得方位角分辨率得到明显提升的高分辨率雷达图像。
(1) 移除RVP相位:对去斜后的雷达回波进行快速傅里叶变换(Fast Fourier Transform, FFT)即可直接得到脉冲压缩的雷达数据。但由于式(6)中的RVP项存在,去斜后的回波脉冲在时间上无法对齐,这给时域加权抑制脉冲压缩旁瓣带来困难。因此,需要通过传递函数
Hc(f)=exp(−jπf2Kr) | (10) |
其中,f表示去斜接收后的距离残差频率。
(2) 距离向时域加窗处理:在移除RVP项后,雷达回波数据通过逆快速傅里叶变换(Inverse Fast Fourier Transform, IFFT)转换到时域,此时利用时域加窗函数来抑制距离脉冲压缩信号的旁瓣,时域余弦窗函数为
Hw(τ)=α+(1−α)⋅cos(ττp⋅π) | (11) |
其中,
sreal(t,f)≈A⋅wa(θs(t)−θ(t))⋅sinc[1k⋅τp(f+2KrcΔR(t))]⋅exp[−j4πfcc⋅ΔR(t)] | (12) |
其中,
(3) 成像区域打网格处理:为了实现通过多角度雷达观测数据相参积累来提升目标的方位几何分辨率的目标,这里采用后向投影(Back Projection, BP)处理算法。BP算法首先需要对目标区域进行打网格处理,成像区域打网格处理如图4所示,其中网格的距离方位为
(4) 多角度雷达回波数据相位校正:从前面提到的车载微波前视成像雷达工作方式来看,雷达通过相控阵波束电扫描来实现对一定角度范围内的目标观测成像。相控阵天线通过移相器和延时线来实现雷达波束电扫描,此时会引入附加的相位,这在一定程度上会影响后续多角度雷达观测数据相参积累。因此,在针对每个网格进行雷达回波数据相参积累之前,需要对不同角度观测的雷达回波数据进行相位校正如图3所示,相位校正函数
(5) 网格内多角度雷达回波相参积累:针对每个成像网格,先利用相位补偿函数
Haz(t)=exp[j4πfcc⋅ΔR(t)] | (13) |
此外,由于距离向雷达回波数据采样率较为有限,在选择指定网格雷达回波数据时,需要进行雷达回波数据的sinc插值或者升采样处理,从而能够精确完成距离徙动校正处理(Range Cell Migration Correction, RCMC)。在网格内多角度雷达回波数据相参积累过程中,各个网格的计算相互独立,这种情况可以采用GPU并行计算提升运算速度[15-17],从而满足汽车雷达实时性的需求。
为了验证本文所提方法的有效性,这里开展点目标仿真试验,雷达仿真参数如表1所示。根据表1中参数,可以计算得到方位向单程法向3 dB波束宽度为0.53°,双程3 dB波束宽度为0.38°,完成整个方位区域扫描共需要扫描100个波位,即获取整个方位区域实孔径雷达图像的时间为2.5 ms,此时对应的方位合成孔径的等效PRF为400 Hz。根据本文所提利用多角度实孔径雷达图像相参合成处理的工作方式和原理,图5(a)给出了在合成孔径长10 m条件下不同成像区域的方位角分辨率结果,图5(b)对比了不同合成孔径长度下可以获得的方位角分辨率。从图5中的仿真结果可以看出,多角度实孔径雷达图像相参积累可以大幅提升车载前视成像雷达的方位角分辨率,但该方法只能对斜前方区域方位角分辨率改善较为明显,而对于正前方目标无法提高其方位角分辨率。
参数 | 数值 |
雷达载频 | 96 GHz |
方位向天线长度 | 0.3 m |
系统PRF | 4000 Hz |
脉冲宽度 | 80 μs |
信号带宽 | 1 GHz |
去斜接收后雷达信号采样率 | 150 MHz |
AD量化位数 | 12位 |
方位波束扫描角度范围 | ±15° |
雷达作用距离 | 20~300 m |
波束跃度 | 0.3° |
合成孔径长度 | 10 m |
汽车速度 | 15 m/s |
为了进一步验证本文方法的有效性,这里设置了一个9个点的点阵成像场景目标,它们之间的几何位置关系如图6所示。图7给出了该点阵目标的实孔径成像结果和多角度相参合成孔径成像的处理结果。从图7中的成像结果看出,位于方位两侧的6个点目标方位角分辨率得到了十分明显提升,而正前方3个点目标的方位角分辨率基本没有变化,这与图5的分析结果完全一致。
图8、图9和图10分别给出了实孔径10 m车载前视雷达图像、15 m合成孔径和5 m合成孔径高分辨率雷达图像中的3个点目标P1, P2和P3的8倍插值后的等高线。从图中的仿真结果可以看出,实孔径雷达图像和合成孔径雷达图像都取得了很好的聚焦效果。此外,实孔径雷达图像中点目标方位向旁瓣较低,这主要由于收发双程天线方向图的缘故。表2给出了实孔径雷达图像、图8和图9点目标的性能测试指标(分辨率、峰值旁瓣比(Peak SideLobe Ratio, PSLR)、积分旁瓣比(Integrated SideLobe Ratio, ISLR))和理论指标,技术指标与理论分析指标基本相一致。
模式 | 目标 | 距离向 | 方位向 | |||||
分辨率(m) | PSLR(dB) | ISLR(dB) | 分辨率(°) | PSLR(dB) | ISLR(dB) | |||
实孔径图像 | P1 | 0.09 | –13.26 | –9.98 | 0.390 | –26.40 | –22.24 | |
P2 | 0.09 | –13.26 | –9.98 | 0.380 | –26.52 | –22.47 | ||
P3 | 0.09 | –13.26 | –9.98 | 0.390 | –26.40 | –22.24 | ||
实孔径理论值 | P1 | 0.09 | –13.26 | –9.80 | 0.380 | –26.60 | –22.30 | |
P2 | 0.09 | –13.26 | –9.80 | 0.380 | –26.60 | –22.30 | ||
P3 | 0.09 | –13.26 | –9.80 | 0.380 | –26.60 | –22.30 | ||
合成处理图像(10 m合成孔径) | P1 | 0.09 | –13.26 | –9.98 | 0.010 | –13.18 | –9.74 | |
P2 | 0.09 | –13.26 | –9.98 | 0.390 | –26.52 | –22.29 | ||
P3 | 0.09 | –13.26 | –9.98 | 0.010 | –13.18 | –9.74 | ||
合成处理理论值(10 m合成孔径) | P1 | 0.09 | –13.26 | –9.80 | 0.010 | –13.26 | –9.80 | |
P2 | 0.09 | –13.26 | –9.80 | 0.390 | –26.60 | –22.30 | ||
P3 | 0.09 | –13.26 | –9.80 | 0.010 | –13.26 | –9.80 | ||
合成处理图像(15 m合成孔径) | P1 | 0.09 | –13.06 | –9.72 | 0.007 | –13.08 | –9.44 | |
P2 | 0.09 | –13.26 | –9.98 | 0.390 | –26.52 | –22.29 | ||
P3 | 0.09 | –13.06 | –9.72 | 0.007 | –13.07 | –9.40 | ||
合成处理理论值(15 m合成孔径) | P1 | 0.09 | –13.26 | –9.80 | 0.007 | –13.26 | –9.80 | |
P2 | 0.09 | –13.26 | –9.80 | 0.390 | –26.60 | –22.30 | ||
P3 | 0.09 | –13.26 | –9.80 | 0.007 | –13.26 | –9.80 |
相比于与超声波、摄像头、激光雷达相比等传感器,微波前视雷达具备在恶劣天气以及极端光线条件下高分辨率成像和运动目标检测的能力,因此其将在未来无人驾驶技术上发挥重要作用。毫米波雷达在方位向通常采用实孔径成像方法,雷达方位分辨受到天线孔径长度限制。针对这个问题,本文利用多角度观测合成孔径处理技术来提升方位角分辨率方法,给出了本成像方法详细的成像处理流程,最后通过计算机仿真验证了方法的有效性。从计算机仿真结果可以看出,本文方法可以提升毫米波前视成像雷达方位角分辨率,在两侧区域方位角分辨率提升效果更佳明显。
XU Zhihuo, SHI Quan, SUN Ling, et al. Novel orthogonal random phase-coded pulsed radar for automotive application[J]. Journal of Radars, 2018, 7(3): 364–375. doi: 10.12000/JR17083
|
PATOLE S M, TORLAK M, WANG Dan, et al. Automotive radars: A review of signal processing techniques[J]. IEEE Signal Processing Magazine, 2017, 34(2): 22–35. doi: 10.1109/MSP.2016.2628914
|
KRONAUGE M and ROHLING H. New chirp sequence radar waveform[J]. IEEE Transactions on Aerospace and Electronic Systems, 2014, 50(4): 2870–2877. doi: 10.1109/TAES.2014.120813
|
TAK J, JEONG E, and CHOI J. Metamaterial absorbers for 24-GHz automotive radar applications[J]. Journal of Electromagnetic Waves and Applications, 2017, 31(6): 577–593. doi: 10.1080/09205071.2017.1297257
|
GUO Kunyi, HOARE E G, JASTEH D, et al. Road edge recognition using the stripe hough transform from millimeter-wave radar images[J]. IEEE Transactions on Intelligent Transportation Systems, 2015, 16(2): 825–833. doi: 10.1109/TITS.2014.2342875
|
MAO Xuesong, INOUE D, MATSUBARA H, et al. Demonstration of in-car doppler laser radar at 1.55 μm for range and speed measurement[J]. IEEE Transactions on Intelligent Transportation Systems, 2013, 14(2): 599–607. doi: 10.1109/TITS.2012.2230325
|
LEE J E, LIM H S, JEONG S H, et al. Enhanced iron-tunnel recognition for automotive radars[J]. IEEE Transactions on Vehicular Technology, 2016, 65(6): 4412–4418. doi: 10.1109/TVT.2015.2460992
|
KELLNER D, BARJENBRUCH M, KLAPPSTEIN J, et al. Tracking of extended objects with high-resolution Doppler radar[J]. IEEE Transactions on Intelligent Transportation Systems, 2016, 17(5): 1341–1353. doi: 10.1109/TITS.2015.2501759
|
WANG Xiao, XU Linhai, SUN Hongbin, et al. On-road vehicle detection and tracking using MMW radar and monovision fusion[J]. IEEE Transactions on Intelligent Transportation Systems, 2016, 17(7): 2075–2084. doi: 10.1109/TITS.2016.2533542
|
WANG H N, HUANG Yingwei, and CHUN S J. Spatial diversity 24-GHz FMCW radar with ground effect compensation for automotive applications[J]. IEEE Transactions on Vehicular Technology, 2017, 66(2): 965–973. doi: 10.1109/TVT.2016.2565608
|
ASKELAND S A and EKMAN T. Tracking with a high-resolution 2D spectral estimation based automotive radar[J]. IEEE Transactions on Intelligent Transportation Systems, 2015, 16(5): 2418–2423. doi: 10.1109/TITS.2015.2407571
|
LEE M S and KIM Y H. Design and performance of a 24-GHz switch-antenna array FMCW radar system for automotive applications[J]. IEEE Transactions on Vehicular Technology, 2010, 59(5): 2290–2297. doi: 10.1109/TVT.2010.2045665
|
HU Chenxi, LIU Yimin, MENG Huadong, et al. Randomized switched antenna array FMCW radar for automotive applications[J]. IEEE Transactions on Vehicular Technology, 2014, 63(8): 3624–3641. doi: 10.1109/TVT.2014.2308895
|
保铮, 邢孟道, 王彤. 雷达成像技术[M]. 北京: 电子工业出版社, 2005: 169–170.
BAO Zheng, XING Mengdao, and WANG Tong. Radar Imaging Technology[M]. Beijing: Publishing House of Electronics Industry, 2005: 169–170.
|
DEVADITHYA S, PEDROSS-ENGEL A, WATTS C M, et al. GPU-accelerated enhanced resolution 3-D SAR imaging with dynamic metamaterial antennas[J]. IEEE Transactions on Microwave Theory and Techniques, 2017, 65(12): 5096–5103. doi: 10.1109/TMTT.2017.2766060
|
CHANG W Y, WU Mengche, CHANG Yanglang, et al. GPU acceleration of adaptive local kriging applied to retrieving slant-range surface motion maps[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2018, 11(11): 4317–4325. doi: 10.1109/JSTARS.2018.2871877
|
WIJAYASIRI A, BANERJEE T, RANKA S, et al. Dynamic data-driven SAR image reconstruction using multiple GPUs[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2018, 11(11): 4326–4338. doi: 10.1109/JSTARS.2018.2873198
|
参数 | 数值 |
雷达载频 | 96 GHz |
方位向天线长度 | 0.3 m |
系统PRF | 4000 Hz |
脉冲宽度 | 80 μs |
信号带宽 | 1 GHz |
去斜接收后雷达信号采样率 | 150 MHz |
AD量化位数 | 12位 |
方位波束扫描角度范围 | ±15° |
雷达作用距离 | 20~300 m |
波束跃度 | 0.3° |
合成孔径长度 | 10 m |
汽车速度 | 15 m/s |
模式 | 目标 | 距离向 | 方位向 | |||||
分辨率(m) | PSLR(dB) | ISLR(dB) | 分辨率(°) | PSLR(dB) | ISLR(dB) | |||
实孔径图像 | P1 | 0.09 | –13.26 | –9.98 | 0.390 | –26.40 | –22.24 | |
P2 | 0.09 | –13.26 | –9.98 | 0.380 | –26.52 | –22.47 | ||
P3 | 0.09 | –13.26 | –9.98 | 0.390 | –26.40 | –22.24 | ||
实孔径理论值 | P1 | 0.09 | –13.26 | –9.80 | 0.380 | –26.60 | –22.30 | |
P2 | 0.09 | –13.26 | –9.80 | 0.380 | –26.60 | –22.30 | ||
P3 | 0.09 | –13.26 | –9.80 | 0.380 | –26.60 | –22.30 | ||
合成处理图像(10 m合成孔径) | P1 | 0.09 | –13.26 | –9.98 | 0.010 | –13.18 | –9.74 | |
P2 | 0.09 | –13.26 | –9.98 | 0.390 | –26.52 | –22.29 | ||
P3 | 0.09 | –13.26 | –9.98 | 0.010 | –13.18 | –9.74 | ||
合成处理理论值(10 m合成孔径) | P1 | 0.09 | –13.26 | –9.80 | 0.010 | –13.26 | –9.80 | |
P2 | 0.09 | –13.26 | –9.80 | 0.390 | –26.60 | –22.30 | ||
P3 | 0.09 | –13.26 | –9.80 | 0.010 | –13.26 | –9.80 | ||
合成处理图像(15 m合成孔径) | P1 | 0.09 | –13.06 | –9.72 | 0.007 | –13.08 | –9.44 | |
P2 | 0.09 | –13.26 | –9.98 | 0.390 | –26.52 | –22.29 | ||
P3 | 0.09 | –13.06 | –9.72 | 0.007 | –13.07 | –9.40 | ||
合成处理理论值(15 m合成孔径) | P1 | 0.09 | –13.26 | –9.80 | 0.007 | –13.26 | –9.80 | |
P2 | 0.09 | –13.26 | –9.80 | 0.390 | –26.60 | –22.30 | ||
P3 | 0.09 | –13.26 | –9.80 | 0.007 | –13.26 | –9.80 |
参数 | 数值 |
雷达载频 | 96 GHz |
方位向天线长度 | 0.3 m |
系统PRF | 4000 Hz |
脉冲宽度 | 80 μs |
信号带宽 | 1 GHz |
去斜接收后雷达信号采样率 | 150 MHz |
AD量化位数 | 12位 |
方位波束扫描角度范围 | ±15° |
雷达作用距离 | 20~300 m |
波束跃度 | 0.3° |
合成孔径长度 | 10 m |
汽车速度 | 15 m/s |
模式 | 目标 | 距离向 | 方位向 | |||||
分辨率(m) | PSLR(dB) | ISLR(dB) | 分辨率(°) | PSLR(dB) | ISLR(dB) | |||
实孔径图像 | P1 | 0.09 | –13.26 | –9.98 | 0.390 | –26.40 | –22.24 | |
P2 | 0.09 | –13.26 | –9.98 | 0.380 | –26.52 | –22.47 | ||
P3 | 0.09 | –13.26 | –9.98 | 0.390 | –26.40 | –22.24 | ||
实孔径理论值 | P1 | 0.09 | –13.26 | –9.80 | 0.380 | –26.60 | –22.30 | |
P2 | 0.09 | –13.26 | –9.80 | 0.380 | –26.60 | –22.30 | ||
P3 | 0.09 | –13.26 | –9.80 | 0.380 | –26.60 | –22.30 | ||
合成处理图像(10 m合成孔径) | P1 | 0.09 | –13.26 | –9.98 | 0.010 | –13.18 | –9.74 | |
P2 | 0.09 | –13.26 | –9.98 | 0.390 | –26.52 | –22.29 | ||
P3 | 0.09 | –13.26 | –9.98 | 0.010 | –13.18 | –9.74 | ||
合成处理理论值(10 m合成孔径) | P1 | 0.09 | –13.26 | –9.80 | 0.010 | –13.26 | –9.80 | |
P2 | 0.09 | –13.26 | –9.80 | 0.390 | –26.60 | –22.30 | ||
P3 | 0.09 | –13.26 | –9.80 | 0.010 | –13.26 | –9.80 | ||
合成处理图像(15 m合成孔径) | P1 | 0.09 | –13.06 | –9.72 | 0.007 | –13.08 | –9.44 | |
P2 | 0.09 | –13.26 | –9.98 | 0.390 | –26.52 | –22.29 | ||
P3 | 0.09 | –13.06 | –9.72 | 0.007 | –13.07 | –9.40 | ||
合成处理理论值(15 m合成孔径) | P1 | 0.09 | –13.26 | –9.80 | 0.007 | –13.26 | –9.80 | |
P2 | 0.09 | –13.26 | –9.80 | 0.390 | –26.60 | –22.30 | ||
P3 | 0.09 | –13.26 | –9.80 | 0.007 | –13.26 | –9.80 |