② (北京环境特性研究所 北京 100854)
② (Beijing Institute of Environmental Features, Beijing 100854, China)
近年来,随着无人机产业快速发展,无人机在航拍、农业、监测、测绘等领域的作用日益显现。然而,无人机的日益增多引发了一些安全担忧,如其对空中交通管理的潜在威胁,以及其滥用于犯罪行为甚至恐怖袭击等[1–3]。因此,对不同类型的无人机进行准确的检测和分类具有重要的意义。
由于雷达具有全天候全天时工作的能力,并且具有实现远距离识别的优势,近几年来利用雷达对无人机进行分类识别开始受到重视。文献[4]提出了从窄带雷达跟踪目标的回波信号中提取飞行轨迹模式与速度特征,以区分无人机与其他鸟类、载人飞机。文献[5,6]通过数值仿真和实验测量分析研究了无人机雷达散射截面(Radar Cross Section, RCS)的变化从而实现对无人机进行分类。文献[7]提出了采用仿真计算与室内测量研究了无人机极化散射特性,揭示了无人机与金属化仿真模型的直接的全极化散射特性差异。文献[8]提出了采用多极化雷达对目标进行观测的方案,实现对鸟类与无人机的区分。此外,基于微多普勒效应的无人机检测与识别是一个重要的研究方向[9–13]。
所谓“微多普勒效应”,指的是雷达目标在进行平动的同时,其整体或部件伴随有振动、旋转、摇摆或其他形式的微运动,使得雷达回波在平动多普勒频率之外还包含附加的时变频率调制的物理现象[9]。微多普勒特征是目标内部部件精细运动在雷达回波中的反映,近年来被广泛应用于人类行为分类[10,11]、手势识别[12,13],目标分类[14–16]等各类雷达目标的分类及识别中。文献[17]通过提取无人机的微多普勒信号特征实现对无人机的分类识别。文献[18]利用卷积神经网络(Convolutional Neural Networks, CNN)与无人机的微多普勒时频图进行无人机分类识别,准确率达到了94.7%。文献[19]通过实验表明时频图的相位信息也可以被利用来提高无人机的识别率。文献[20,21]提出利用多角度或者多站点的雷达观测的微多普勒特征相对于单站点雷达信号而言,有助于提高无人机的识别鲁棒性。
针对提高无人机识别的准确率,本文重点研究了基于多角度雷达传感器观测和微动特征的无人机分类方法。首先采用时频分析方法获得无人机的微动特性时频图,然后利用主成分分析 (Principal Component Analysis, PCA)从时频图中提取特征,并对多角度雷达传感器观测的数据进行特性融合。相对于单个角度雷达传感器观测而言,多角度雷达传感器观测的融合特征的维数增加了一倍。最后利用支持向量机(Support Vector Machine, SVM)进行训练与识别不同类型的无人机。基于实测数据的实验结果表明,该方法对3种典型无人机的分类效果良好,多角度雷达传感器观测融合得到的分类精度优于单角度雷达传感器。
2 实验数据采集实验采用双波段连续波雷达系统同时从不同角度对无人机进行数据采集。采用两个不同波段雷达是为了避免两个传感器之间的相互干扰。其中K波段雷达的载波频率
![]() |
图 1 实验场景几何示意图 Fig.1 The experimental geometry |
实验采集了3种不同类型的无人机包括四轴无人机、直升机、六轴无人机,如图2所示。由于微多普勒信号对相位角[15]高度敏感,因此采用了6个不同的垂直和水平位置来固定无人机,从而实现对无人机不同角度的观测。每个固定装置的重复测量10次,每次6 s,总共有(2部雷达)×(6 s)×(3类无人机)×(6个固定位置)×(10次测量)=2160 s的测量数据。
![]() |
图 2 3类无人机示意图 Fig.2 Appearance of the drones |
典型的基于微多普勒特征识别的过程一般包括回波信号预处理、目标特征提取与选择、特征数据集构建、分类器训练和识别等如图3所示,其中图4为图3红框部分的流程图,具体处理论述见本节。
![]() |
图 3 典型识别流程 Fig.3 Typical flow chart of radar recognition |
![]() |
图 4 特征提取与融合流程示意图 Fig.4 The flow chart of feature extraction and fusion with dual radar sensors |
首先将实测回波数据分割成若干个分段测量回波样本。本文设置每个分段回波样本的时间长度为0.5 s,这个时间长度足够包含若干个无人机的回波旋转周期,从而确保获得精确的微多普勒调制。然后利用高斯窗的短时傅里叶变换(Short-Time Fourier Transform, STFT)对样本进行时频分析。
![]() |
其中,
![]() |
图 5 3种无人机的时频图 Fig.5 Spectrograms of three types of drones |
从图5可以看出,无人机目标的微多普勒调制显现周期特性,其中直升机的回波信号结构较为简单,微动分量周期特性分明,而多轴无人机的回波信号较为复杂,微动分量周期边缘起伏程度较为明显。这是因为直升机的旋翼结构相对比较简单,主旋翼共用一个转子,而多轴无人机包含多个转子和对应的旋翼桨,尽管转子转动的角速度一般是一致的,但每两个转子之间初始相位差并不固定,因此不同类型的无人机的雷达回波中的多普勒频率调制有着较大的差别,可依据这些特性实现对无人机目标的分类。
3.2 特征提取PCA是一种基于目标统计特性的最佳正交变换。由于它能够有效地消除冗余数据并且本身的计算量小,能很好地用于实时处理,本文采用PCA进行特性提取。首先将2维时频图进行向量化操作,排列为1维向量
![]() |
其中,
![]() |
图6(a)和图6(b)分别展示了仅使用K波段雷达传感器和仅使用X波段雷达传感器的实验数据的前3个特征分布(即Y(1:3))。图6(c)显示了利用t分布邻域嵌入算法(t-distributed Stochastic Neighbor Embedding, t-SNE)[22]将两个角度观测的雷达融合特征分布在3维空间的可视化分布。从图6中可以看出,融合了两个角度两个波段观测的雷达传感器特征相对于单角度单波段观测的雷达传感器而言,进一步增强了特征空间的聚类。
![]() |
图 6 实验数据的特征分布 Fig.6 Feature distribution of the experimental data |
通过上述预处理与特征提取得到训练特征数据集,如何提高目标识别的正确率是分类器需要解决的问题。由于微多普勒信号对相位角较为敏感,而在实际中常常不能得到目标所有相位角度的回波结果,因此导致数据集的不完备性。传统的线性分类方法常常不能得到良好的目标识别性能。基于结构风险最小原理(Structure Risk Minimization, SRM)和核技巧,SVM[23]通过映射原始的特征矢量到高维空间,算法仅依赖于训练库的子集,具有良好的推广能力,计算复杂度低。本文采用SVM分类器来评估性能,其判决函数
![]() |
式中:q为支持向量个数
为了评估该方法的分类性能,每一类目标数据随机抽取30%即216个样本作为训练样本,剩下70%即504个样本作为测试样本进行分类器的训练与测试。对于不同波段雷达数据采用100次蒙特卡洛实验得到的平均结果如表1–表3所示。
![]() |
表 1 仅利用K波段雷达传感器识别混合矩阵 Tab.1 Confusion matrix using only the K-band radar sensor |
![]() |
表 2 仅利用X波段雷达传感器识别混合矩阵 Tab.2 Confusion matrix using only the X-band radar sensor |
![]() |
表 3 融合双观测角度的雷达传感器识别混合矩阵 Tab.3 Confusion matrix fusing dual-angle radar sensors |
如表1–表3所示,直升机更容易与其他两类无人机区分,这是因为直升机的主旋转叶片共用一级转子,以固定相位旋转,而其他两类无人机的旋转叶片之间的旋转相位差并不固定导致每个叶片的回波相位并不相关。此外由于K波段雷达传感器的波长较短,单独使用K波段雷达传感器比单独使用X波段雷达传感器具有更高的分类准确率。当无人机由双波段雷达传感器同时从不同角度测量时,采用双角度双波段雷达观测的融合特征量是单独使用单雷达传感器的两倍,因此融合后提高了目标分类的准确率。
4.1 不同训练样本比例下的识别率分析本实验考察了在不同的训练样本比例下的3类无人机识别率。设置训练样本比例从10%增加到90%,间隔为10%,采用100次蒙特卡洛实验得到的平均识别率如图7所示。
![]() |
图 7 识别率随训练样本比例变化情况 Fig.7 Recognition accuracy under different sizes of training set |
由图7可知,随着训练样本比例增大,基于双雷达传感器融合的方法、基于K波段雷达传感器和基于X波段雷达传感器的方法的识别率均提高。基于双雷达传感器融合方法能够进一步提高基于单雷达传感器方法的识别率。
4.2 不同噪声水平下的识别率分析本实验考察了在不同噪声水平下的3类无人机识别率。实验中,将无人机的实测雷达信号与仿真加性高斯白噪声(AWGN)按照如下方式混合[13]:
![]() |
其中,
![]() |
图 8 识别率随加性噪声强度变化情况 Fig.8 Recognition accuracies under different levels of additive Gaussian white noise |
图8表明,随着加性噪声的强度变大,基于单雷达传感器的方法的识别率下降较快,基于双雷达传感器融合方法的识别率超过了基于单雷达传感器的方法。此外,在训练样本较少的情况下(30%数据作为训练样本),基于双雷达传感器融合方法的识别率在低噪声水平下相对于基于单雷达传感器的方法具有优势,这一结果与图7的结果保持一致。
5 结论本文研究了一种基于多角度双波段雷达观测特征融合提高无人机分类精度的方法。基于观测数据的实验结果表明,采用双角度双波段雷达传感器进行特征融合相对于单个雷达传感器而言可以获得更高的目标分类精度。本文的研究成果可拓展到双/多基地雷达中。在下一步工作中,将进一步研究多角度的相干通道雷达的微多普勒特征融合的提高无人机分类性能的方法。
[1] |
De Wit J J M, Harmanny R I A, and Prémel-Cabic G. Micro-Doppler analysis of small UAVs[C]. Proceedings of the 9th European Radar Conference, Amsterdam, Netherlands, 2012: 210–213.
(![]() |
[2] |
De Wit J J M, Harmanny R I A, and Molchanov P. Radar micro-Doppler feature extraction using the Singular Value Decomposition[C]. Proceedings of 2014 International Radar Conference, Lille, France, 2014: 1–6. DOI: 10.1109/RADAR.2014.7060268.
(![]() |
[3] |
夏铭禹, 赵凯, 倪威. 要地防控反无人机系统及其关键技术[J].
指挥控制与仿真, 2018, 40(2): 53-60, 71. Xia Ming-yu, Zhao Kai, and Ni Wei. Anti-UAV system and key technology for key point defense[J]. Command Control & Simulation, 2018, 40(2): 53-60, 71. DOI:10.3969/j.issn.1673-3819.2018.02.010 ( ![]() |
[4] |
Mohajerin N, Histon J, Dizaji R, et al.. Feature extraction and radar track classification for detecting UAVs in civillian airspace[C]. Proceedings of 2014 IEEE Radar Conference, Cincinnati, USA, 2014: 674–679. DOI: 10.1109/RADAR.2014.6875676.
(![]() |
[5] |
Ritchie M, Fioranelli F, Griffiths H, et al.. Micro-drone RCS analysis[C]. Proceedings of 2015 IEEE Radar Conference, Johannesburg, South Africa, 2015: 452–456. DOI: 10.1109/RadarConf.2015.7411926.
(![]() |
[6] |
Schroder A, Renker M, Aulenbacher U, et al.. Numerical and experimental radar cross section analysis of the quadrocopter DJI phantom 2[C]. Proceedings of 2015 IEEE Radar Conference, Johannesburg, South Africa, 2015: 463–468. DOI: 10.1109/RadarConf.2015.7411928.
(![]() |
[7] |
白杨, 吴洋, 殷红成, 等. 无人机极化散射特性室内测量研究[J].
雷达学报, 2016, 5(6): 647-657. Bai Yang, Wu Yang, Yin Hong-cheng, et al. Indoor measurement research on polarimetric scattering characteristics of UAV[J]. Journal of Radars, 2016, 5(6): 647-657. DOI:10.12000/JR16032 ( ![]() |
[8] |
Torvik B, Olsen K E, and Griffiths H. Classification of birds and UAVs based on radar polarimetry[J].
IEEE Geoscience and Remote Sensing Letters, 2016, 13(9): 1305-1309. DOI:10.1109/LGRS.2016.2582538 (![]() |
[9] |
Chen V C. The Micro-Doppler Effect in Radar[M]. Boston: Artech House, 2011: 35–60.
(![]() |
[10] |
Kim Y and Ling H. Human activity classification based on micro-Doppler signatures using a support vector machine[J].
IEEE Transactions on Geoscience and Remote Sensing, 2009, 47(5): 1328-1337. DOI:10.1109/TGRS.2009.2012849 (![]() |
[11] |
Fioranelli F, Ritchie M, and Griffiths H. Classification of unarmed/armed personnel using the NetRAD multistatic radar for micro-Doppler and singular value decomposition features[J].
IEEE Geoscience and Remote Sensing Letters, 2015, 12(9): 1933-1937. DOI:10.1109/LGRS.2015.2439393 (![]() |
[12] |
Kim Y and Toomajian B. Hand gesture recognition using micro-Doppler signatures with convolutional neural network[J].
IEEE Access, 2016, 4: 7125-7130. DOI:10.1109/ACCESS.2016.2617282 (![]() |
[13] |
Li G, Zhang R, Ritchie M, et al. Sparsity-driven micro-Doppler feature extraction for dynamic hand gesture recognition[J].
IEEE Transactions on Aerospace and Electronic Systems, 2018, 54(2): 655-665. DOI:10.1109/TAES.2017.2761229 (![]() |
[14] |
Stove A G and Sykes S R. A Doppler-based target classifier using linear discriminants and principal components[C]. Proceedings of 2003 International Conference on Radar, Adelaide, SA, Australia, 2003: 107–125. DOI: 10.1109/RADAR.2003.1278734.
(![]() |
[15] |
熊丁丁, 崔国龙, 孔令讲, 等. 基于互相关熵的非高斯背景下微动参数估计方法[J].
雷达学报, 2017, 6(3): 300-308. Xiong Dingding, Cui Guolong, Kong Lingjiang, et al. Micro-motion parameter estimation in non-Gaussian noise via mutual correntropy[J]. Journal of Radars, 2017, 6(3): 300-308. DOI:10.12000/JR17007 ( ![]() |
[16] |
王宝帅. 基于微多普勒效应的空中飞机目标分类研究[D]. [博士论文], 西安电子科技大学, 2015: 4–22.
Wang Bao-shuai. Study on classification of airplane targets based on micro-Doppler effect[D]. [Ph.D. dissertation], Xidian University, 2015: 4–22. ( ![]() |
[17] |
Molchanov P, Egiazarian K, Astola J, et al. . Classification of small UAVs and birds by micro-Doppler signatures[C]. Proceedings of 2013 European Radar Conference, Nuremberg, Germany, 2013: 172–175.
(![]() |
[18] |
Kim B K, Kang H S, and Park S O. Drone classification using convolutional neural networks with merged Doppler images[J].
IEEE Geoscience and Remote Sensing Letters, 2017, 14(1): 38-42. DOI:10.1109/LGRS.2016.2624820 (![]() |
[19] |
Ren J F and Jiang X D. Regularized 2-D complex-log spectral analysis and subspace reliability analysis of micro-Doppler signature for UAV detection[J].
Pattern Recognition, 2017, 69: 225-237. DOI:10.1016/j.patcog.2017.04.024 (![]() |
[20] |
Ritchie M, Fioranelli F, Borrion H, et al. Multistatic micro-Doppler radar feature extraction for classification of unloaded/loaded micro-drones[J].
IET Radar, Sonar & Navigation, 2017, 11(1): 116-124. DOI:10.1049/iet-rsn.2016.0063 (![]() |
[21] |
Zhang R, Li G, Clemente C, et al. Multi-aspect micro-Doppler signatures for attitude-independent L/N quotient estimation and its application to helicopter classification[J].
IET Radar, Sonar & Navigation, 2017, 11(4): 701-708. DOI:10.1049/iet-rsn.2016.0271 (![]() |
[22] |
Van Der Maaten L and Hinton G. Visualizing data using t-SNE[J].
Journal of Machine Learning Research, 2008, 9(11): 2579-2605. (![]() |
[23] |
Vapnik V N. The Nature of Statistical Learning Theory[M]. New York: Springer-Verlag, 1995: 112–150.
(![]() |