(武汉大学电子信息学院，湖北武汉430072)
关键词：高频地波雷达；射频干扰抑制；YOLO；高阶奇异值分解；深度学习
本文提出将深度学习中的YOLO(YouOnlyLookOnce)算法，结合高阶张量分析的HOSVD算法，实现RFI自动识别检测与抑制一体化工作。实验结果表示，与之前的RFI抑制技术相比，所提出的YOLO-HOSVD算法具有显著优势，既能够自动识别RFI并定位，又能在不损坏所需信号的情况下快速抑制RFI，能够实现自动实时化HFSWR信号处理工作，可以实际落地到HFSWR系统应用中进行实时精细化处理。
近年来，深度学习算法逐渐用于雷达图像识别和信号处理方面，如基于SAR的目标检测[13]，利用CNN和Nanodet等网络对车辆[14]、飞机[15]、舰船[16]等目标进行检测，深度学习应用在雷达领域上相比传统方法具有快速、准确、效率高的优点[17-18]。目前，深度学习主要应用在SAR目标探测，鲜见应用于高频雷达干扰识别与抑制。HFSWR回波数据经过两次快速傅里叶变换(FastFourierTransform,FFT)后，数据处理结果反映在距离多普勒谱(Range-DopplerSpectrum,RDS)也是以图像形式，因此可以应用图像处理的方法来进行特征识别。
YOLO-HOSVD算法与传统的干扰检测算法相比，具有自动提取深度特征、准确率高等优点。以YOLO为代表的目标检测算法将检测问题转化为回归问题，实现“端到端”处理[19]，直接预测对象类别和位置，速度和精度取得了良好平衡，可进行实时检测。目前，YOLO已成功应用于农业[20]、地质[21]、遥感[22]和医学[23]等领域。该算法基于RD谱中RFI的特征，能够利用小样本数据训练后实时输出检测结果，定位RFI数据，结合HOSVD算法处理多通道张量数据，自动实时抑制RFI。
以YOLOv5s为例，模型由4部分组成：输入端负责数据预处理、主干部分提取特征、检测部分采集特征图和预测端输出结果[24]。YOLOv5的网络体系如图1所示。YOLOv5检测时将RDS平均拆分为S2个网格，若待检测对象的中心在某一网格上，此网格就负责预测该对象。每个网格产生多个候选框，每个候选框包含回归位置信息和置信度值，置信度值代表候选框中包含物体的置信度以及候选框的预测准确度的双重信息。模型训练时会预设大小不一的候选框，微调置信度高的候选框后输出预测框，计算二者差值后反向更新，候选框与预测框如图2所示。
图1YOLOv5s网络体系结构图
图2候选框与预测框示意图
(1)
公式(1)表示候选框和预测框关系，tx和ty表示待检测目标所在网格的中心位置，候选框的长宽分别为pw和ph，预测框的长宽分别为bw和bh，cx和cy表示偏移左上角值，tw和th代表尺度缩放值，最终得到的预测坐标值是(bx,by,bw,bh)。
YOLOv5使用GIoU作为损失函数，精度(Precision,P)、召回率(Recall,R)、平均精度均值(meanAveragePrecision,mAP)和每秒传输帧数(FramesPerSecond,FPS)作为评估指标[25]。
HFSWR采用线性调频中断连续波(FrequencyModulatedInterruptedContinuousWave,FMICW)，通过两次FFT提取RDS信息。RFI可以看作是一系列单频信号的线性组合，在第n个扫频周期的t时刻，某一信道的单频RFI信号Sn(t)表示为
Sn(t)=Ajcos(2πfj(t+(n-1)Ts+φj))
(2)
式中，Aj、fj和φj分别为第j个RFI的振幅、载波频率和初始相位，Ts为扫频周期。经过正交混频和低通滤波之后，RFI分量转换为
(3)
RFI在RDS中通常表现为沿距离轴的明亮条纹，根据RDS特征，结合YOLOv5进行处理并输出模型。构建训练所需的数据集步骤如下：第一步是对回波数据进行两次FFT，得到RDS；第二步是将图像分为两类；第三步是使用标注软件标记对象。
标注软件用矩形框标记RD谱，将对应矩形框的结果整理成PascalVOC格式可直接被YOLOv5读取，如图3所示，Bragg峰和RFI用矩形框标记，框标签信息保存为xml文件，文件中包含沿矩形框对角线的两点坐标和类别名称等信息。
图3标记方法和标签信息
Dopplerin=2·(2·Xcenter-1)
(4)
Rangebin=90·(1-2·Ycenter)
(5)
YOLOv5检测结果输出为图4右侧数据结构形式，经过公式(4)和公式(5)映射到RDS上范围，图4中标注了RFI和一个Bragg峰在RDS上和YOLOv5结果对应数据，结合RFI和Bragg峰位置，可以自适应输出训练张量和待抑制张量数据以供HOSVD处理。
图4YOLOv5检测结果图
1.4.1HOSVD投影
HOSVD充分利用RFI的多维特性，构造训练张量A∈I1×I2×I3，I1、I2和I3分别代表距离元、多普勒单元和通道数。为了准确估计干扰子空间，所构造的训练张量一般包含RFI而不包含海洋回波，训练张量和待抑制张量通过YOLOv5检测结果自动匹配，图5为RFI数据进行HOSVD示意图。
图5HOSVD分解示意图
由文献[27]，A写作模n积形式可得
A=S×1U(1)×2U(2)×3U(3)
(6)
S=A×1U(1)T×2U(2)T×3U(3)T
(7)
HOSVD将A分解为核心张量S∈S1×S2×S3和模n积形式，U(n)(n=1,2,3)分别是沿3个维度的展开矩阵，均为正交矩阵，U(n)U(n)T代表U(n)的投影矩阵。
1.4.2RFI子空间估计
由于在多普勒维，RFI与海洋回波并不相互独立，结合RFI的空、频域特点对其分量进行联合估计，对模态2和模态3的展开矩阵进行SVD可得
(8)
(9)
U(2)和U(3)进行SVD分解后的左奇异矩阵分别是U2=[u21,u22,…,u2I2]和U3=[u31,u32,…,u3I3]，列向量分别为u2i∈I2×1,i=1,2,…,I2和u3i∈I3×1,i=1,2,…,I3，Σ2I2×I1I3和Σ3I3×I1I2是奇异值矩阵，V2∈I1I3×I1I3和V3∈I1I2×I1I2是右奇异矩阵。
RFI在三种模态展开矩阵中对应的投影矩阵为
(10)
(11)
(12)
1.4.3张量投影抑制RFI
图6构造多维结构RDS中的训练张量
YOLOv5输出的训练张量为A∈I1×I2×I3，待处理张量为χ，如图6所示，其中χ中的多普勒单元和通道数与A相同，I1个距离单元仅存在RFI不存在海洋回波。χ和A中的RFI的空、频域子空间相同，将待处理向量χ向估计的多域干扰子空间投影，即可得χ的RFI子空间投影分量χr是：
χr=χ×1Proj1×2Proj2×3Proj3
(13)
RFI抑制后，期望信号的张量χs为
χs=χ-χr
(14)
基于YOLO-HOSVD算法的自动识别检测与抑制RFI算法基本步骤如图7所示，具体步骤如下：
1)对HFSWR回波数据进行两次FFT处理后，输出RDS，标记对象制作数据集，输入YOLOv5模型中进行训练，输出训练模型；
3)计算张量各模态展开矩阵，应用HOSVD算法计算各模态展开矩阵的左奇异矩阵Un(n=1,2,3)，得到多域RFI投影子空间矩阵Projn(n=1,2,3)；
4)将待处理张量往估计的子空间投影，得到χ在RFI子空间投影分量χr；
5)信号张量χ减去投影分量χr，即可得到抑制干扰后的期望信号。
图7自动识别检测抑制RFI流程图
通过仿真对基于YOLO-HOSVD的自动抑制RFI算法性能进行测试，根据中国福建省漳州市东山县(23.6575°N，117.4863°E)的2016年12月27日的HFSWR探测海洋数据为基础，模拟单频信号组成的RFI和目标添加到原始无干扰的数据中，HFSWR的系统参数如表1所示。
表1雷达参数表
图8高频地波雷达距离多普勒谱
原始HFSWR数据的RDS如图8(a)所示，将仿真单频RFI和目标信号依次加入到RDS中，RFI覆盖范围为多普勒频率0.30～0.45Hz范围内所有距离元，仿真目标位于第34距离元，速度为4m/s，多普勒频率为0.36Hz，目标被RFI覆盖，无法从RDS中直接观测，如图8(b)所示，RFI掩盖一阶Bragg峰和目标时无法提取有用信息，对海态反演和目标检测都有较大的影响。将加入RFI和目标的数据输入到训练好的YOLOv5模型，自动输出训练张量和待抑制张量数据，图8(c)为经过YOLOv5s处理后的输出结果，训练张量A的范围是图8(c)中红色实线框内的8个通道的数据，待抑制张量范围为图8(c)中蓝色实线框内所有通道的数据。图8(d)是经过HOSVD处理后的RDS，结果表明RFI显著减轻，被掩盖的Bragg峰出现，海洋回波得以保留，目标也凸显出来。综上所述，YOLOv5结合HOSVD算法能够在不丢失有用信息的情况下实时自动识别并抑制RFI，即使RFI覆盖海洋回波和目标，HOSVD也能有效抑制干扰。
仿真实验中，对RDS进行识别分类用了0.01s，抑制RFI用了0.04s。分析可得，YOLO-HOSVD算法在几乎不损失有用信号的同时可以自动地检测并抑制RFI，即使在RFI覆盖海洋回波和目标的情况下依然适用，增加了HFSWR系统自动识别抑制RFI的能力，为海态反演和目标探测提供有力帮助。
图9东山和龙海雷达站位置
表2双站双频地波组网模式
两个雷达站空间直线距离约90km，分别位于中国福建省的龙海市(24.2674°N，118.1353°E)和漳州市东山县(23.6575°N，117.4863°E)，每个站包含两路高低频发射天线和4个宽带接收天线(6～13.5MHz)，双站通过GPS同步协同工作，构成双基地、单基地、多频和MIMO的混合系统。
本次试验中YOLOv5模型的运行环境：CPU为AMDRyzen75800HwithRadeonGraphics，GPU为NVIDIAGeForceRTX3060，操作系统为window11。将YOLOv5s应用在本次实测数据随机挑选的200张包含RFI的RD谱，进行训练并输出模型。训练集放入YOLO模型训练，以YOLOv5s为例，在一般配置的笔记本电脑上训练300轮仅需0.5h左右，训练好的模型可供后续数据处理持续使用，处理单场数据仅需11ms左右。FPS越大，模型检测率越快。当FPS≥30时，可以实现实时处理速度[28]。
图10YOLOv5s模型的训练结果
利用YOLOv5s模型计算2021年8月8日至9月8日31天内双频双站干扰存在率，结果展示在图12中。
从图12可以看出，频率不同、位置不同的站点受干扰污染的情况也不同，且在一个月连续不间断的实测中，双站双频雷达几乎每天都会受到RFI污染，只是程度不一样。说明抑制RFI对于实测数据处理是不可避免的一项工作，抑制RFI可以保证有用信息的提取，提高HFSWR数据利用率，有利于海态反演和目标信号检测处理工作。
图12双频双站RFI存在率
图13雷达距离多普勒谱图
选取4场不同特征的RFI抑制结果展示，图13所示为8月12日上午11时35分的雷达RDS图，图13(a)显示该场数据存在一较强能量的RFI，部分海洋回波被掩盖。应用YOLOv5s模型后输出结果为多普勒范围在0.33～0.74Hz范围内红色框覆盖的距离元的4通道张量作为训练数据A，待抑制张量为多普勒0.33～0.74Hz范围内所有距离元数据，如图13(b)所示，抑制RFI后结果如图13(c)所示，识别和抑制过程共用了0.41s，经过处理，RFI被消除，之前被RFI掩盖的海洋回波有效保留。
从RDS中可以看出，距离元位置不同，海洋回波的能量也不同，为了验证HOSVD对海洋回波能量的影响，计算RDS中覆盖海洋回波的0～10距离元的RFI抑制前后信号与干扰加噪声比(SignaltoInterferenceplusNoiseRatio,SINR)变化，如图14所示。结果表明，系统的SINR有了明显的提高。增加幅度最高是第6个距离元的18.32dB，SINR平均增加约14.42dB。
图14RFI抑制前后信噪比变化曲线
图15为8月12日凌晨2时10分的RFI抑制前后RDS图，图15(a)中包括4个明显的RFI污染和大范围电离层干扰，图15(b)为YOLOv5输出张量结果图，图15(c)为经过抑制后处理的结果，处理该场数据用了1.78s，结果显示4个RFI都得到了明显抑制，对电离层和期望信号数据没有影响。
图15雷达距离多普勒谱图
图16雷达距离多普勒谱图
图17雷达距离多普勒谱图
从图17(b)可以看出，经过一次YOLO-HOSVD处理后，部分强烈RFI被抑制，再次通过YOLOv5s检测后的结果如图17(c)所示，标记出了第一次未被检测出的RFI，应用HOSVD处理后的结果如图17(d)所示，可以看出剩余RFI已被抑制。
从图18可以看出，该算法抑制RFI后干扰存在率降低了56.60%，抑制效率达到93.68%。YOLO-HOSVD算法没有达到100%完全抑制效果可能存在的原因有2种，一是存在其他强烈干扰与RFI干扰形态相似，导致RFI误判，另一种是部分RFI并没有被完全抑制。经过人工对比后，第二次经过YOLOv5算法检测后的结果的确存在RFI干扰，说明YOLO并未存在误判情况，可能存在原因是部分RFI污染没有完全抑制，例如图17中展示情况，解决方案是将抑制后数据重新送入YOLO模型中检测，如若重复出现干扰即进行二次抑制即可提高RFI抑制效果。
本文将深度学习算法应用于HFSWR干扰抑制中，提出了基于YOLO-HOSVD的HFSWR系统自动识别检测与抑制RFI的算法，充分利用了YOLO和HOSVD算法的优点，利用YOLO模型识别RDS中的RFI和Bragg，自动划分训练张量和待抑制张量，实时输出张量数据，将HOSVD算法应用于多通道HFSWR数据的三维张量来估计干扰子空间，利用了多通道HFSWR数据的固有结构和射频信号在频域和空间域的特性，从而提高了RFI抑制性能。仿真和实测结果表明，YOLO-HOSVD算法实现了对HFSWR中RFI的自动识别与抑制，应用在普通性能的笔记本处理单场数据不超过1.8s，将该方法应用于长期遭受严重RFI影响的HFSWR数据上，可以实时有效准确地统计回波谱遭受RFI污染情况，并自动抑制干扰。此算法打破了传统RFI抑制算法不能在工程应用中大规模处理数据的限制，可作为HFSWR信号处理的系统方法，进行实时精细化处理，也为雷达数据质量评估提供参考，有利于后续海流反演、目标探测等信号处理工作。
参考文献：
[1]CHUX,ZHANGJ,JIY,etal.ExtractionofWindDirectionfromtheHFHybridSky-SurfaceWaveRadarSeaEchoes[J].IEEEAerospaceandElectronicSystemsMagazine,2018,33(3):42-47.
[2]PABLOL,MARCOSS,LOTFIA,etal.ExtremeWaveHeightEventsinNWSpain:ACombinedMulti-SensorandModelApproach[J].RemoteSensing,2018,10(1):1.
[3]WYATTLR.MeasuringtheOceanWaveDirectionalSpectrum‘FirstFive’withHFRadar[J].OceanDynamics,2019,69(1):123-144.
[4]LIC,WUX,YUEX,etal.ExtractionofWindDirectionSpreadingFactorfromBroad-BeamHigh-FrequencySurfaceWaveRadarData[J].IEEETransonGeoscienceandRemoteSensing,2017,55(9):5123-5133.
[5]AIX,LUOY,ZHAOG.TransientInterferenceExcisioninOver-the-HorizonRadarbyRobustPrincipalComponentAnalysiswithaStructuredMatrix[J].IEEEGeoscience&RemoteSensingLetters,2015,13(1):48-52.
[6]GURGELKW,BARBINY,SCHLICKT.RadioFrequencyInterferenceSuppressionTechniquesinFMCWModulatedHFRadars[C]∥OCEANS2007-Eu-rope,Aberdeen,UK:IEEE,2007:1-4.
[7]李敏乐,毕大平.宽带波束域频率聚焦自适应波束形成算法[J].雷达科学与技术,2016,14(6):594-598.
LIMinle,BIDaping.AFrequency-FocusingAdaptiveBeamformingAlgorithminWidebandBeamspace[J].RadarScienceandTechnology,2016,14(6):594-598.(inChinese)
[8]NAZARIME,HUANGW,ZHAOC.RadioFrequencyInterferenceSuppressionforHFSurfaceWaveRadarUsingCEMDandTemporalWindowingMethods[J].IEEEGeoscienceandRemoteSensingLetters,2020,17(2):212-216.
[9]杨俊.高频地波雷达干扰抑制方法研究[D].武汉：武汉大学,2004.
[10]刘爱军.基于极化信息的高频地波雷达干扰抑制方法研究[D].哈尔滨：哈尔滨工业大学,2011.
[11]ZHENT,WENB,JINL,etal.RadioFrequencyInterferenceSuppressionAlgorithminSpatialDomainforCompactHigh-FrequencyRadar[J].IEEEGeoscienceandRemoteSensingLetters,2017,15(1):102-106.
[12]LIY,YUEX,WUX,etal.AHigher-OrderSingularValueDecomposition-BasedRadioFrequencyInterferenceMitigationMethodonHigh-FrequencySurfaceWaveRadar[J].IEEETransonGeoscienceandRemoteSensing,2020,58(4):2770-2781.
[13]夏勇，田西兰，常沛，等.基于深度学习的复杂沙漠背景SAR目标检测[J].雷达科学与技术,2019,17(3):305-309.
XIAYong，TIANXilan，CHANGPei，etal.SARTargetDetectioninComplexDesertBackgroundImagesBasedonDeepLearning[J].RadarScienceandTechnology,2019,17(3):305-309.(inChinese)
[14]常沛，夏勇，李玉景，等.基于CNN的SAR车辆目标检测[J].雷达科学与技术,2019,17(2):220-224.
CHANGPei,XIAYong,LIYujing,etal.SARVehicleTargetDetectionBasedonCNN[J].RadarScienceandTechnology,2019,17(2):220-224.(inChinese)
[15]王思雨，高鑫，孙皓,等，基于卷积神经网络的高分辨率SAR图像飞机目标检测方法[J].雷达学报，2017,6(2):195-203.
[16]刘方坚,李媛.基于视觉显著性的SAR遥感图像NanoDet舰船检测方法[J].雷达学报,2021,10(6):885-894.
[17]张云,穆慧琳,姜义成,等.基于深度学习的雷达成像技术研究进展[J].雷达科学与技术,2021,19(5):467-478.
ZHANGYun,MUHuilin,JIANGYicheng,etal.OverviewofRadarImagingTechniquesBasedonDeepLearning[J].RadarScienceandTechnology,2021,19(5):467-478.(inChinese)
[18]毕海霞，魏志强.深度学习在极化SAR图像分类上的应用综述[J].雷达科学与技术,2021,19(5):539-551.
BIHaixia,WEIZhiqiang.ASurvey:theApplicationofDeepLearninginPolSARImageClassification[J].RadarScienceandTechnology,2021,19(5):539-551.(inChinese)
[19]聂凡杰.基于端到端的深度学习目标检测算法研究[D].北京:北京邮电大学,2018.
[20]陈斌,张漫,徐弘祯,等.基于改进YOLOv3-tiny的全景图像农田障碍物检测[J].农业机械学报,2021,52(S1):58-65.
[21]GAOJ,CHENY,WEIY,etal.DetectionofSpecificBuildinginRemoteSensingImagesUsingaNovelYOLO-S-CIOUModel.Case:GasStationIdentification[J].Sensors,2021,21(4):1375.
[22]PANGL,LIUH,CHENY,etal.Real-TimeConcealedObjectDetectionfromPassiveMillimeterWaveImagesBasedontheYOLOv3Algorithm[J].Sensors,2020,20(6):1678-1693.
[23]YAOSJ,CHENYW,TIANX,etal.AnImprovedAlgorithmforDetectingPneumoniaBasedonYOLOv3[J].AppliedSciences,2020,10(5):1818.
[24]LISW,GUXY,XUXG,etal.DetectionofConcealedCracksfromGroundPenetratingRadarImagesBasedonDeepLearningAlgorithm[J].ConstructionandBuildingMaterials,2021,273:121949.
[25]YANB,FANP,LEIX,etal.AReal-TimeAppleTargetsDetectionMethodforPickingRobotBasedonImprovedYOLOv5[J].RemoteSensing,2021,13(9):1619.
[26]杨绍麟，柯亨玉，侯杰昌，等.OSMAR2000基于MUSIC的超分辨率海洋表面流算法[J].武汉大学学报(理学版),2001(5):601-608.
[27]NGUYENVD,KAIIMAM,LINH-TRUNGN.FastTensorDecompositionsforBigDataProcessing[C]∥2016InternationalConferenceonAdvancedTechnologiesforCommunications(ATC),Hanoi,Vietnam:IEEE,2016:215-221.
[28]MAH,LIUY,RENY,etal.DetectionofCollapsedBuildingsinPost-EarthquakeRemoteSensingImagesBasedontheImprovedYOLOv3[J].RemoteSensing,2020,12(1):44.
DONGSonglin，YUEXianchang，WUXiongbin，GAOYubin
(SchoolofElectronicInformation,WuhanUniversity,Wuhan430072,China)
Abstract:Thedetectionperformanceofhighfrequencysurfacewaveradar(HFSWR)isextremelysusceptibletoradiofrequencyinterference(RFI).Currently,thesuppressionofRFIinHFSWRechoesisconductedframebyframethroughmanualidentification,andtherearefewresearchesonautomaticidentificationandsuppressionofRFI.Withthedevelopmentofdeeplearninginradarimageprocessing,thisarticletriestointroduceitintohighfrequencyradarRFIsuppression,andusestheYOLO(youonlylookonce)modeltoidentifytheRFIintheradarrangeDopplerspectrum,andthensuppressedbythehigherordersingularvaluedecomposition(HOSVD)method.TheresultsofsimulationandexperimentshowthattheYOLO-HOSVDjointalgorithmcanrealizetheautomaticidentificationandsuppressionofRFIinHFSWRechoes,andtheprocessingtimeofasinglefielddatadoesnotexceed1.8s.ThismethodisapplicabletotheoperationalHFSWRtoautomaticallyidentifyandsuppresstheRFIinobservingseastate.
Keywords:HFSWR;RFIsuppression;YOLO;HOSVD;deeplearning
DOI:10.3969/j.issn.1672-2337.2022.03.004
收稿日期：2022-01-05；修回日期：2022-02-07
基金项目：广东省重点领域研发计划(No.2020B1111020005)；湖北省重点研发计划项目(No.2020BCA080)；科技部重点研发计划(No.2016YFC1401101)
中图分类号：TN972
文献标志码：A
作者简介：
东松林女，1997年生，河南周口人，武汉大学电子信息学院在读硕士研究生，主要研究方向为高频地波雷达干扰抑制。
岳显昌男，1975年生，辽宁人，博士，副教授，主要研究方向为高频雷达遥感海洋风浪反演、雷达信号处理。
吴雄斌男，1968年生，湖北人，博士，教授，主要研究方向为高频地波雷达海洋监测以及无线电波传播、电离层层析成像算法。
高玉斌男，1998年生，硕士研究生，主要研究方向为高频地波雷达干扰抑制。