合成孔径雷达(SAR)[1]因其具备全天候、全天时的工作能力以及分辨率优势,被广泛应用于战场侦察和地理测绘。星载[2]和机载[3]SAR 凭借远程探测和大带宽的特点,成为高效的战略侦察和地球测绘工具。然而,星载和机载SAR 的调度过程较为复杂,无法应对快速变化的战场环境或灾害救援。
近年来,随着国家对低空经济的积极推动,无人机(Unmanned Aerial Vehicle,UAV)行业得到了迅速发展。如图1所示,将微型SAR(NanoSAR)[4-5]系统搭载于中小型无人机上能够充分发挥它们各自的优势。无人机载NanoSAR 通常用于低空遥感作业,能够快速便捷地获取高分辨率SAR 图像,其便利性和灵活性远超传统的星载和机载SAR系统[6]。
图1 无人机载SAR成像几何示意图
传统SAR 系统基于高精度ADC 采样,具有数据量大、功耗高的缺点。我们利用频移单比特成像算法进行SAR 成像,在保持成像质量的同时,大大降低了数据量,简化了SAR 系统的软硬件资源,使得本文设计NanoSAR 系统的重量更小、功耗更低。
对于距离雷达为Rt 的目标,调频连续波(FMCW)的接收回波可建模为
式中:Tp, fc和Kr分别是所发射信号的脉冲宽度、载波频率和调频斜率;σ 和τ 分别表示散射体的散射系数和时间延迟;tr 表示快时间;目标的回波延迟τ = 2Rt/c。接收机中的参考信号为
通过对s(tr )与sref( tr)混频处理可以得到中频信号:
其中φ表示相位项:
根据文献[7],对中频信号sif(tr)进行单比特量化:
式中,sign(·)表示取数据符号的符号函数。单比特量化导致数据的实部在1 和-1 之间变化,虚部在j和-j 之间变化,从而产生高阶谐波分量。式(5)可以进一步表示为
式中,表示m 阶谐波的幅值,Jm(·) 为m 阶贝塞尔函数,其中m ≠0 时ε0 = 1,εm = 2。当m为偶数时Am = 0,式(6)可重建模为
单比特量化信号由原信号的奇次谐波组成。
一般来说,在许多应用中,三次谐波和五次谐波的能量较大,是不可忽略的。为此,可使用频谱搬移的方法将三次和五次谐波推到更高的频段,保留有效频率信息。定义基带范围为0 ∼fmax,采样率为fs。从上述分析中可以看出中频信号的基带频率范围为0~fmax,五次谐波的频率范围为0~5fmax,三次谐波的频率范围为0~-3fmax。五次谐波的频带与原始频带重叠,会对SAR 成像造成影响。因此,我们考虑添加一个固定的频偏fΔ,使用该频偏对中频信号进行调制,最后再进行单比特采样。这个过程可以表示为
增加采样率以确保采样带宽足够宽就能将三次和五次谐波的频谱移动到远离基带的位置。考虑到混叠的情况,频率偏移和采样率应满足以下条件:
在随后的信号处理中,丢弃0~fΔ 的带宽,保留fΔ~fΔ + fmax 的带宽。如上所述,适当的参数设置使我们能够保留单比特量化的优势的同时,抑制单比特量化带来的谐波对成像的影响。
无人机载SAR成像的几何关系如图1所示,无人机平台在距离地面高度为H 的空中,以速度v 沿规划的航线飞行,进行侧视成像。
单比特NanoSAR 是一款专门为中小型无人机设计的雷达系统,其核心技术参数详见表1。该雷达采用FMCW 体制,工作频段位于Ka 波段,中心频率为35 GHz。选择该频段是为了减小射频相关硬件的尺寸和重量,同时更容易提升雷达的带宽,从而增强雷达的距离分辨能力。这款雷达样机支持最高达2 GHz的射频带宽,距离分辨率能够达到7.5 cm。此外,该雷达通过FPGA 和DAC 生成基带线性调频信号,使其能够根据使用场景和预设无人机的飞行速度,灵活调整调频斜率和脉冲重复频率等雷达参数。
表1 雷达样机主要设计指标
名称雷达体制中心频率带宽分辨率探测距离重量无人机飞行高度数值LFMCW 35 GHz 2 GHz(可调)7.5 cm 500~1 000 m 1.85 kg 120~500 m
NanoSAR 系统的硬件结构如图2 和图3 所示,主要包括数字板、射频板、电源板、收发天线以及惯性导航设备(惯导)。整个系统的结构支撑材料为碳纤维和铝合金,兼顾了轻量化和强度。
图2 雷达系统实物照片
图3 雷达系统的主要硬件构成
数字板是系统的核心控制单元,负责与地面控制设备通信,接收地面发出的指令并回传雷达的工作状态信息。同时,数字板还负责生成基带信号、采集中频回波信号、数据打包、存储以及传输等功能。
射频板的功能包括为雷达系统提供工作时钟,将数字板输出的基带信号上变频至Ka 波段,经过滤波和放大后传输至发射天线。此外,射频板负责将天线接收的回波信号与发射端的信号进行混频和滤波,生成中频信号并输出给数字板进行处理。
SAR 成像理想情况下要求平台进行匀速直线运动,但在实际应用中,无人机平台的运动往往存在不稳定性,主要表现为运动速度和欧拉角的变化。这种不稳定性会对SAR 成像产生干扰,导致图像失真。为此,该雷达样机配备了具备实时动态测量(Real Time Kinematic,RTK)定位功能的惯导模组,可以提供厘米级的定位精度和精确的角度变化信息,为成像算法中的运动补偿步骤提供了数据支持。
雷达样机的重量和功耗分配如表2所示,雷达样机包括机械结构在内的总体重量约1 850 g,雷达开启探测时功耗达到峰值,峰值功耗约60 W,雷达系统使用无人机提供电源,不需要额外的储能单元。
表2 雷达样机的重量和功耗
名称数字电路板(含散热外壳)射频收发机(含散热外壳)收发天线电压转换板雷达吊舱框架描述生成、采集、存储、传输基带信号发射和接收射频信号两个喇叭天线输入12~48 V,输出12 V和5 V固定其他组件的碳纤维板和铝合金重量750 g 450 g 180 g 100 g 370 g功耗46 W 12 W 2 W
雷达样机的硬件接口方案如图4 所示。数字板作为雷达系统的控制和数据处理中枢,通过多种接口与其他组件进行通信。具体来说,数字板与射频板之间通过SMA 接口传输时钟和雷达模拟信号,通过J30J-25接口控制收发增益。此外,射频板通过2.92 mm接口与收发天线相连。
图4 硬件接口方案示意图
数字板通过LoRA 与地面控制设备传输控制命令和工作状态信息,通过千兆以太网或无线网桥传输雷达数据。此外,数字板通过RS232/TTL接口与惯导通信,通过M2 接口与板载固态硬盘通信。
在数字板内部,时钟芯片、DAC 芯片以及ADC芯片均通过SPI接口进行寄存器配置。FPGA 通过JESD204B 接口向DAC 传输雷达基带信号,通过LVDS接口接收ADC采集的中频信号。
数字板的设计示意图如图5所示,数字板的主控芯片是XC7Z100 系列SoC,该SoC 由FPGA(称为PL)以及ARM 子系统(称为PS)组成,该ARM 子系统的核心是Cortex A9 处理器。FPGA 部分控制时钟分配芯片、DAC、ADC、DDR3、固态硬盘以及IMU等外设并进行数据传输。Cortex A9 处理器负责解析来自上位机的控制指令、传输雷达数据以及控制雷达系统的工作参数与功能开关。
图5 数字板的设计示意图
在雷达开始采集数据之前,地面端通过LoRA向数字板发送一系列配置参数和初始化指令。数字板上的处理器(PS)接收到这些指令后,会将其解析并传递给FPGA 执行。为了优化能耗管理,数字板中一些功耗较高的芯片默认处于断电状态,FPGA 根据使用需求动态控制电源管理芯片,对这些芯片进行上电和断电,并保证上电时序。
雷达开启探测时,其发射端的工作流程如下:首先,数字板上的FPGA 根据预设的调频斜率、起始频率和脉冲重复频率等参数,利用多路并行DDS 生成带宽为500 MHz、等效采样率为2 GS/s 的LFMCW数字信号。
FPGA 产生的基带数字信号通过JESD204B 接口传输至DAC。在DAC 内部,基带信号进行上采样,并根据设定的参数进行数字IQ 调制,将信号调制至2.5~3 GHz 频段。经过滤波和DA 转换后,信号被输出至射频板。
雷达接收端的工作流程如下:ADC 负责采集来自射频板的中频信号,并将采样数据传输至FPGA。FPGA 将这些采样数据与惯导数据及雷达工作参数结合,构建完整的雷达数据包,并将其临时缓存在DDR3 内存中。同时,FPGA 构建NVMe 的提交队列(Submission Queue,SQ)请求,并通过门铃寄存器(Doorbell Register,DB)通知固态硬盘从DDR3 中读取雷达数据包进行存储。在存储过程中,雷达会将起始扇区和结束扇区的信息回传至地面端,以便后续读取。
当进行雷达数据回传时,地面端首先依据记录中的扇区地址范围,通知雷达回传指定扇区范围的数据。数字板中的PS 接收到数据回传指令后,会尝试通过有线或无线链路与地面端建立TCP连接。TCP 连接建立成功后,PS指示FPGA 从固态硬盘中读取雷达数据包并缓存在DDR3中,并响应来自FPGA 的中断请求将这些数据包回传给地面端。地面端的数据处理软件则基于雷达数据包中记录的雷达参数、采样数据以及无人机平台的运动信息进行成像处理,生成SAR图像。
射频板的设计示意图如图6所示,该射频板支持输入雷达基带信号的频率范围为2.5~3 GHz。其发射功率大于18 dBm,接收增益约50 dB。
图6 射频板的设计示意图
与传统的FMCW架构的收发机不同,该射频板内部使用两路本振信号同时与数字板中DAC 输出的信号进行混频,混频的结果其中一路用于发射,另一路用于与接收端进行混频,从而实现将中频信号的起始频率从零频移动至5 MHz或15 MHz,这有利于在进行单比特采样时减少三次和五次谐波分量的干扰。
使用KEYSIGHT 公司N9040B 频谱分析仪对雷达样机的收发性能进行测试。图7 展示了射频板发射端的测试结果,测试中使用的发射带宽是1.8 GHz,测量结果表明发射信号的频率范围是34.1~35.9 GHz,发射功率大于18 dBm,通带的频谱平坦度较好,峰值旁瓣比大于23 dB。图8 展示了射频板输出的中频信号频谱,该频谱的通带范围是5~16 MHz,带宽约11 MHz。
图7 发射功率和带宽测试
图8 中频信号频谱
雷达天线的水平辐射方向图和俯仰辐射方向图的测试结果分别如图9和图10所示,该天线的频率范围是34~38 GHz,测试中使用的频率是35 GHz。测试结果表明,雷达天线水平方向的3 dB 波束宽度是15°,俯仰方向的3 dB 波束宽度是50°,水平方向的旁瓣低于-12 dB,俯仰方向的旁瓣低于-20 dB。
图9 水平辐射方向图
图10 俯仰辐射方向图
无人机载SAR 的运动补偿相较于传统星载和机载SAR 存在两个主要的挑战:1) 由于无人机平台对载荷体积和功耗的限制,无人机载SAR 通常选用较高的雷达波段,更短的波长会导致成像结果对平台非规则运动产生的相位误差更敏感;2)无人机平台飞行的稳定性较差,易受气流等环境因素的干扰。
本文采用能同时提供卫星导航与RTK 定位功能的惯导模组采集无人机平台的运动信息。在成像过程中,首先使用卡尔曼滤波对惯导模组采集的位置信息进行平滑处理,滤波后的位置信息作为粗略的运动补偿用于初步成像,随后使用相位梯度自聚焦方法对初步的成像结果做进一步聚焦处理,对于子图像之间的配准和拼接,本文结合使用了平台的运动信息与图像特征点的匹配信息来估计子图像间的几何变化关系,从而提高图像配准的精度。
为了测试雷达样机的成像效果,我们将其安装在Matrice 600 PRO 无人机上进行成像测试。挂飞实验的照片如图11 所示,测试中使用的无人机运动参数和雷达系统参数详见表3。图12 展示了测试场景的卫星光学图像,图13(a)和图13(b)分别展示了高精度(16-bit)和单比特(1-bit)数据的SAR 成像结果,图14 展示了两种采样数据成像结果的局部对比。SAR 成像结果表明,本文设计的雷达样机具备对较大场景进行SAR 成像的能力,成像场景中的草坪、树木、人行天桥和篮球场等目标均具有较高的辨识度。
表3 成像测试参数
参数飞行高度/m飞行速度/(m·s-1)擦地角/(°)带宽/GHz中心频率/GHz脉冲宽度/μs采样率/MHz脉冲重复频率/Hz参数值120 7.5 45 1.5 35 450 20 2 000
图11 雷达样机挂飞照片
图12 卫星光学图像
图13 SAR成像测试结果
图14 SAR成像测试结果(局部)
对于图13 中展示的成像结果,以高精度数据得到的SAR 图像作为评价基准,使用结构相似度(Structural Similarity,SSIM)作为评价指标,单比特数据成像结果的SSIM 为88.57%,与高精度SAR 图像质量十分接近。使用不同倍数的过采样频率,单比特成像所需的数据量可降低到高精度的1/3甚至1/5,可以有效降低雷达系统对数据存储和传输的性能要求,降低雷达系统的功耗和成本等方面存在优势。
此外,对于需要实时SAR 成像的应用场景,高精度量化数据由于数据率大,无法实时无线传输雷达数据,需要在雷达系统中部署成像单元,将成像结果传回地面端。而单比特量化后的数据由于数据率的降低,雷达系统能无线实时回传采样数据,在地面端进行成像处理,避免在雷达系统中部署成像单元,有利于简化雷达系统的结构,并降低系统的成本和功耗。
将本文设计的NanoSAR 挂载在无人机上单比特成像实验,取得与预期相符的成像效果,验证了本系统设计的可行性和稳定性。相较于传统SAR系统,单比特NanoSAR 大幅减少采样数据量,降低了传输和处理的负担,降低了对无人机载重量的要求,在低空地质测绘,资源勘探等领域有着十分广泛的前景。
[1] HOVANESSIAN SEVAK A. Introduction to Synthetic Array and Imaging Radars[M]. USA:Artech House, 1980.
[2] STOFAN E R, EVANS D L, SCHMULLIUS C, et al. Overview of Results of Spaceborne Imaging Radar-C, X-Band Synthetic Aperture Radar(SIR -C/X -SAR)[J]. IEEE Trans on Geoscience and Remote Sensing, 1995, 33(4):817-828.
[3] GRAY A L, FARRIS-MANNING P J. Repeat-Pass Interferometry with Airborne Synthetic Aperture Radar[J].IEEE Trans on Geoscience and Remote Sensing, 1993, 31(1):180-191.
[4] 王岩飞,刘畅,詹学丽,等.无人机载合成孔径雷达系统技术与应用[J].雷达学报,2016,5(4):333-349.
[5] 吕泽鑫,仇晓兰,丁赤飚.无人机载PolInSAR 城区散射机制初步分析[J].雷达科学与技术,2022,20(5):475-479.
[6] 朱岱寅,张营,俞翔,等.微型合成孔径雷达成像信号处理技术[J].雷达学报,2019,8(6):793-803.
[7] PAN Tianyi, ZHAO Bo, HUANG Ling, et al. One-Bit SAR Imaging and Moving Target Detection Based on Frequency Shift of Dechirp Echo[J].IET International Radar Conference, Online, 2021:233-238.
UAV-Mounted 1-bit Nano Synthetic Aperture Radar System
王 伟 男,博士研究生,主要研究方向为SAR 系统设计和FPGA工程实现。
龙天尧 男,博士研究生,主要研究方向为雷达信号处理和稀疏阵列设计。
黄 磊 男,博士,教授、博士生导师,主要研究方向为阵列信号处理、稀疏信号处理、统计信号处理、压缩感知、频谱感知、机器学习、雷达信号处理、雷达通信一体化。