CN114254141A - 一种基于深度分割的端到端雷达信号分选方法 - Google Patents
一种基于深度分割的端到端雷达信号分选方法 Download PDFInfo
- Publication number
- CN114254141A CN114254141A CN202111563815.7A CN202111563815A CN114254141A CN 114254141 A CN114254141 A CN 114254141A CN 202111563815 A CN202111563815 A CN 202111563815A CN 114254141 A CN114254141 A CN 114254141A
- Authority
- CN
- China
- Prior art keywords
- pdw
- pulse
- pulses
- image
- pixel point
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 52
- 230000011218 segmentation Effects 0.000 title claims abstract description 47
- 238000012545 processing Methods 0.000 claims abstract description 8
- 230000005855 radiation Effects 0.000 claims description 31
- 230000006870 function Effects 0.000 claims description 21
- 239000011159 matrix material Substances 0.000 claims description 21
- 238000012549 training Methods 0.000 claims description 21
- 230000008569 process Effects 0.000 claims description 19
- 238000005070 sampling Methods 0.000 claims description 15
- 230000009466 transformation Effects 0.000 claims description 10
- 230000004913 activation Effects 0.000 claims description 8
- 238000001208 nuclear magnetic resonance pulse sequence Methods 0.000 claims description 8
- 238000006243 chemical reaction Methods 0.000 claims description 6
- ORILYTVJVMAKLC-UHFFFAOYSA-N Adamantane Natural products C1C(C2)CC3CC1CC2C3 ORILYTVJVMAKLC-UHFFFAOYSA-N 0.000 claims description 5
- 238000010606 normalization Methods 0.000 claims description 5
- 238000005259 measurement Methods 0.000 claims description 3
- 238000012216 screening Methods 0.000 claims description 3
- 239000013598 vector Substances 0.000 claims description 3
- 238000004422 calculation algorithm Methods 0.000 description 14
- 238000013135 deep learning Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000018109 developmental process Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000003064 k means clustering Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000003909 pattern recognition Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/50—Information retrieval; Database structures therefor; File system structures therefor of still image data
- G06F16/51—Indexing; Data structures therefor; Storage structures
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/50—Information retrieval; Database structures therefor; File system structures therefor of still image data
- G06F16/55—Clustering; Classification
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/214—Generating training patterns; Bootstrap methods, e.g. bagging or boosting
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/22—Matching criteria, e.g. proximity measures
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/241—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
- G06F18/2415—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches based on parametric or probabilistic models, e.g. based on likelihood ratio or false acceptance rate versus a false rejection rate
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10032—Satellite or aerial image; Remote sensing
- G06T2207/10044—Radar image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20081—Training; Learning
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Evolutionary Computation (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Computational Biology (AREA)
- Artificial Intelligence (AREA)
- Life Sciences & Earth Sciences (AREA)
- Databases & Information Systems (AREA)
- Probability & Statistics with Applications (AREA)
- Software Systems (AREA)
- Image Analysis (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明属于雷达信号分选技术领域,具体涉及一种基于深度分割的端到端雷达信号分选方法。本发明通过对PDW的处理和深度分割网络的结合实现未知复杂环境下的端到端雷达信号分选,操作简单,无需按照传统的预分选和主分选每一步都预设参数、阈值等信息。本发明能够分选脉间调制类型多变的雷达信号、时域频域交叠严重的PDW,适应复杂环境和未知信号。此外本发明还能够分选出从未训练过且没有先验信息的信号,对于脉冲丢失20%的情况也适用良好。
Description
技术领域
本发明属于雷达信号分选技术领域,具体涉及一种基于深度分割的端到端雷达信号分选方法。
背景技术
雷达信号分选技术是关键技术之一,它主要是将接收机接收到的未知的、交叠的脉冲信号按照不同的辐射源个体区分开,如图1所示,为接下来的进一步分析和采取对抗策略提供基础。学者们根据传统的统计模式识别算法研制了许多雷达信号分选方法。
实际应用中比较广泛的信号分选算法是利用脉冲描述字(pulse descriptionword,PDW)进行分析。PDW包含了每个脉冲的不同属性特征,如到达方向(direction ofarrival,DOA)、载频(carrier frequency,CF)、脉冲宽度(pulse width,PW)、脉冲幅度(pulse amplitude,PA)、脉冲到达时间(time of arrival,TOA)等。常见的信号分选算法大多都是利用这些信息对脉冲流进行去交错,采用预分选和主分选相结合的方式,将交叠混乱的PDW分选出不同辐射源脉冲描述字(emitter description word,EDW)。对于预分选,有较多的学者采用聚类算法将PDW进行初步的聚类,如k-means聚类算法及其改进算法、基于密度的聚类算法DBSCAN等,这些预分选方法需要事先设定聚类中心、阈值、容差等参数;在预分选的基础上,有学者采用统计直方图算法进行主分选,如一些学者利用CDIF、SDIF等算法统计聚类后每组脉冲的脉冲重复间隔,然后配合脉冲搜索完成主分选;还有很多学者在预分选的基础上,采用脉冲重复间隔(pulse repetition interval,PRI)变换法进行主分选,同样利用PRI变换算法对聚类后的脉冲分组估计PRI值,然后配合脉冲搜索完成主分选,算法中采用了可变的移动因子,并选用了固定的箱分辨率以及合理的门限,能够适应PRI值分布在较大范围内雷达信号。此外,由于深度学习的良好的处理大数据能力,有些学者将深度学习相关算法引入到信号分选中来,但是仅仅是针对雷达信号的PRI调制方式进行识别。
这些传统的基于模板和固定顺序流程的算法,针对传统的雷达信号取得了不错的效果。但是随着现代电子对抗和雷达发展的发展,电磁环境时域、频域和空域上也更复杂。表现为雷达信号脉冲流密度已经达到数百万级别,信号交叠严重;雷达信号参数调制复杂多变,多样的调制方式使得从复杂的交叠脉冲序列分选辐射源描述字EDW更为艰难。这些传统算法的固定参数和人工经验设定的阈值无法自适应很多复杂调制类型的信号,对没有先验信息的信号也无法分选成功。预分选加主分选的前后结构互相制约,缺乏灵活性,无法实现端到端的处理过程。利用深度学习方法得到脉冲组的PRI调制方式只能作为分选过程中的一个参考,目前还未见有人将脉冲聚类,脉冲搜索等与深度学习相结合,具体实现由交叠的PDW到有序的EDW分选过程。
近年来,学者们将深度学习应用到图像处理领域,并取得了非常好成果,其中深度分割相关模型可以实现像素级分类。本发明受此启发,利用PDW构成PDW序列图像,单个PDW表现为像素点,EDW的PDW序列表现为像素分类结果,从而可以考虑将深度分割方法引入到信号分选中,比如U-Net网络。避免了预分选、主分选的传统分选中的前后制约的影响,构建端到端的实现方案,从而能灵活应对未知复杂电磁环境。
发明内容
本发明的目的在于提供一种基于深度分割的端到端雷达信号分选方法。
一种基于深度分割的端到端雷达信号分选方法,包括以下步骤:
步骤1:PDW采集,截取时间长度为T、带宽为B的PDW并保存为PDW脉冲序列;
步骤2:PDW序列图像生成;
选择PDW中CF、PW、TOA这三个参数进行图像变换处理,TOA、CF值用于描述像素点位置,PW值用于描述像素点颜色深度,从而构成tl(TOA)×fw(CF)×pd(PW)大小的三维脉冲矩阵P;同时为适应深度分割网络的图像的分辨率大小输入需求,需完成图像比例变换,构建图像矩阵Q,并根据对应的比例变换关系,构建像素点与脉冲之间的索引表;
步骤3:利用训练好的U-Net深度分割模型进行像素点分类,得到像素点分类结果图像矩阵Q2;
步骤4:PDW序列搜索;根据像素点的分类结果图像矩阵Q2和步骤2中得到的像素点与脉冲之间的索引表,按照像素点的分类情况对脉冲进行搜索归类;对于一个像素点索引多个脉冲的情况,采用参数匹配法,实现PDW的归类,获得最终不同辐射源信号的EDW。
进一步地,所述步骤2中构建像素点与脉冲之间的索引表的方法具体为:
假设每个维度上脉冲与像素点对应的比例变换关系为β1,β2,β3,根据三维脉冲矩阵P以及脉冲序列图像矩阵Q,可以得到不同维度上的变换关系分别为β1,β2,β3:
若第n个脉冲的TOA和CF为{TOAn,CFn},它对应像素点在PDW序列图像的二维位置坐标则为:
(x,y)=(β1×TOAn,β2×CFn)
这样脉冲流中每个脉冲信号都对应PDW序列图像上一个像素点,进而构建像素点索引表,如表1所示:
表1像素点与脉冲之间的索引表
表中K表示K个脉冲像素点,每个像素点最多对应M个脉冲;脉冲参数{TOA,CF,PW}k,m的下角标k,m表示由第k个像素点可以搜索到m个脉冲,同时得到该脉冲的TOA、CF、PW参数值。
进一步地,所述步骤3中U-Net深度分割模型具体为:
U-Net深度分割模型对输入的PDW序列图像的像素点进行深度分割,得到像素点分类结果图像Q2,U-Net深度分割模型的详细结构如表2所示:
表2U-Net网络参数
U-Net深度分割过程分为下采样和上采样,各自包含5个网络层,每次下采样包括两个卷积核为3x3的卷积,每个卷积后面跟着一个激活函数ReLU和一个归一化函数BatchNormlization,两次卷积能够提取出图像的深层特征,2x2的最大池操作使图像尺寸缩小;在每个下采样步骤中,特征通道的数量增加一倍,图像尺寸减半;上采样过程采用反卷积结构,每个反卷积的卷积核为2x2,后面继续进行两个卷积核为3x3的卷积,同样每个卷积后面跟着一个激活函数ReLU和一个归一化函数BatchNormlization;每个上采样特征通道减半,图像尺寸放大一倍;同时为了优化输出和提高特征的利用率,需要将反卷积的所得深层特征和前面下采样对应提取的特征向量进行拼接融合;网络的最后一层是一个3x3卷积,接着是一个激活函数Sigmoid。
进一步地,所述步骤3中U-Net深度分割模型的训练方法具体为:
在训练过程中,训练集中输入为PDW序列图像,目标输出为像素点分类好的图像;训练的目的是让U-Net更注重区分不同像素点间的差距,像素点不会丢失,PDW序列中像素点的位置也不会改变,即U-Net的输入Q和输出Q2对应;训练在基于PyTorch框架的基础上,采用了衡量模型输出与目标输出之间的交叉损失函数熵:
其中,oi代表网络模型的输出,对应的是中PDW图进入网络后的预测结果;ti代表目标输出,对应PDW图像的像素点分类结果标签;I是网络计算一次损失所需要的样本数;σi是网络的需要学习的权重参数;模型参数优化器选择Adam函数,它利用梯度的一阶矩和二阶矩自适应调整每个参数的学习速率;设定学习率初始值、对数据集的迭代次数、损失阈值,训练时当损失收敛达到损失阈值时停止训练保存训练好的模型。
进一步地,所述步骤4中进行PDW序列搜索的方法具体为:
步骤4.1:初筛;先将一个像素点索引一个脉冲的情况优先挑选归类出来,形成该辐射源信号的EDW;
步骤4.2:统计该辐射源信号的EDW范围;遍历EDW所有脉冲的参数信息,统计好该辐射源信号CF、PW、TOA的范围[CFmin,CFmax],[PWmin,PWmax],[TOAmin,TOAmax],若不存在一个像素点对应多个脉冲的情况则完成该辐射源信号的分选,可继续对下一类像素点进行搜索归类;若存在一个像素点对应多个脉冲的情况,则需要进参数匹配;
步骤4.3:参数匹配:将一个像素点对应的多个脉冲依次与上述统计辐射源信号的参数范围进行匹配,若某个脉冲参数均匹配成功,则将该脉冲挑选出来归类到这个辐射源信号的EDW中,完成对这类像素点的脉冲搜索归类操作;匹配不成功则存放至剩余脉冲组中,当所有类像素点完成对应脉冲搜索归类后,将个别剩余的脉冲与不同的辐射源EDW参数范围再次匹配归类,最终PDW搜索结束完成分选。
本发明的有益效果在于:
本发明通过对PDW的处理和深度分割网络的结合实现未知复杂环境下的端到端雷达信号分选,操作简单,无需按照传统的预分选和主分选每一步都预设参数、阈值等信息。本发明能够分选脉间调制类型多变的雷达信号、时域频域交叠严重的PDW,适应复杂环境和未知信号。此外本发明还能够分选出从未训练过且没有先验信息的信号,对于脉冲丢失20%的情况也适用良好。
附图说明
图1是雷达信号分选过程的示意图。
图2是本发明的流程框图。
图3是PDW序列图像生成过程示意图。
图4是深度分割模型U-Net结构图。
图5是不同PDW输入的分选结果图。
图6是不同方法分选结果对比图。
图7是脉冲丢失下分选结果图。
图8是不同方法脉冲丢失情况分选结果对比图。
具体实施方式
下面结合附图对本发明做进一步描述。
为减少信号分选对预置的参数和流程以及人工经验的依赖,加强信号分选方法对参数捷变,如频率捷变雷达、脉组捷变雷达、重频参差、抖动雷达等雷达,对时频域混叠严重、脉冲丢失严重,对未知复杂电磁环境的分选能力,同时加深深度学习相关网络模型与实际分选应用之间的结合,本发明设计了一种基于深度分割的端到端雷达信号分选方法,通过对PDW的处理和深度分割网络的结合实现未知复杂环境下的端到端雷达信号分选。
本发明的实现方案步骤如图2所示:
第一步:PDW采集。可截取时间长度为T、带宽为B的PDW并保存为PDW脉冲序列。
第二步:PDW序列图像生成。选择PDW中CF、PW、TOA这三个参数进行图像变换处理:TOA、CF值用于描述像素点位置,PW值用于描述像素点颜色深度。从而构成tl(TOA)×fw(CF)×pd(PW)大小的三维脉冲矩阵P。同时为适应深度分割网络的图像的分辨率大小输入需求,如l(长)×w(宽)×d(深度),需完成图像比例变换,构建图像矩阵Q。并根据对应的比例变换关系,构建像素点与脉冲之间的索引表。
第三步:U-Net深度分割。利用训练好的U-Net深度分割模型像素点分类,同时得到输出像素点分类结果图像矩阵Q2。
第四步:PDW序列搜索。根据像素点的分类结果图像矩阵Q2和第二步的像素点索引表,按照像素点的分类情况对脉冲进行搜索归类。对于一个像素点索引多个脉冲的情况,可采用参数匹配法,实现PDW的归类,获得最终不同辐射源信号的EDW。
在网络训练过程中,输入可采用多部雷达辐射源交叠的PDW序列。现有的深度分割模型有许多,本发明选取10层U-Net网络对PDW序列图像进行像素点分类。U-Net网络最初提出用于细胞图像的分割,通过下采样和上采样的结构产生非常好的分割效果。
由以上步骤可看出,该分选方案实现了端到端的过程,且由于引入了深度网络,可适应未知环境。
本发明适用的条件:
(1)足够的样本。需要采样固定时长、固定带宽的大量PDW数据对U-Net网络进行训练。
(2)PDW序列图像的大小需要适应具体深度网络的输入需求。需要利用PDW矩阵和PDW序列图像的对应比例变换关系,适应具体的深度网络输入。
本发明实现了端到端分选,操作简单,无需按照传统的预分选和主分选每一步都预设参数、阈值等信息。本发明能够分选脉间调制类型多变的雷达信号,时域频域交叠严重的PDW,适应复杂环境和未知信号。此外还能够分选出从未训练过且没有先验信息的信号,对于脉冲丢失20%的情况也适用良好,如图5至图8所示,
实施例1:
(1)PDW采集
在空间中存在许许多多信号,接收机可以采集到空间中的一些信号。假设在固定时长T内,固定带宽B条件下,接收机截获N个雷达脉冲信号形成PDW,如图3所示(a)到(b)过程所示。PDW之间都是相互独立,每个PDW包含这个脉冲的许多特征参数,比如载频CF、脉宽PW、到达时间TOA等,若根据这三个参数的PDW三维示意图如图(c)所示。接下来就是将采集到的PDW保存起来为PDW序列图像的生成提供数据。
(2)PDW序列图像生成
PDW序列图像生成示意图如图3中(c)到(d)映射过程,将每一个脉冲可视化为二维平面图像上的一个像素点。选择PDW中CF、PW、TOA这三个参数进行图像变换处理:TOA、CF值用于描述像素点位置,PW值用于描述像素点颜色深度。从而构成tl(TOA)×fw(CF)×pd(PW)大小的三维脉冲矩阵P。同时为适应深度分割网络的图像的分辨率大小输入需求,如l(长)×w(宽)×d(深度),需完成图像比例变换,构建图像矩阵Q。
假设每个维度上脉冲与像素点对应的比例变换关系为β1,β2,β3,根据三维脉冲矩阵P以及脉冲序列图像矩阵Q则,可以得到不同维度上的变换关系分别为β1,β2,β3:
对于一个脉冲而言,它的TOA描述其在PDW序列图像横轴方向上的位置,CF描述其纵轴方向上的位置,若第n个脉冲的TOA和CF为{TOAn,CFn},它对应像素点在PDW序列图像的二维位置坐标则为:
(x,y)=(β1×TOAn,β2×CFn) (2-4)
这样脉冲流中每个脉冲信号都对应PDW序列图像上一个像素点,可以构建像素点索引表,如表1所示。
表中K表示K个脉冲像素点,每个像素点最多对应M个脉冲(M<<N)。表中脉冲参数{TOA,CF,PW}k,m的下角标k,m表示由第k个像素点可以搜索到m个脉冲,同时得到该脉冲的TOA、CF、PW参数值。
表1像素点索引表
像素点索引包含了PDW脉冲序列和PDW序列图像之间的连接关系,为后续利用U-Net得到像素点分类结果索引出同属于一类信号的辐射源脉冲描述字EDW提供条件。
(3)U-Net深度分割模块
U-Net深度分割模块是利用U-Net对输入的PDW序列图像的像素点进行深度分割,得到像素点分类结果图像Q2。U-Net网络为详细结构如表2所示,U-Net深度分割过程如图4所示,主要分为下采样和上采样,各自包含5个网络层。
表2U-Net网络参数
每次下采样包括两个卷积核为3x3的卷积(无填充卷积),每个卷积后面跟着一个激活函数(ReLU)和一个归一化函数(BatchNormlization)。两次卷积能够提取出图像的深层特征,2x2的最大池操作(步幅为2),使图像尺寸缩小。在每个下采样步骤中,特征通道的数量增加一倍,图像尺寸减半。上采样过程采用反卷积结构,每个反卷积的卷积核为2x2,后面继续进行两个卷积核为3x3的卷积,同样每个卷积后面跟着ReLU和BatchNormlization。每个上采样特征通道减半,图像尺寸放大一倍。同时为了优化输出和提高特征的利用率,需要将反卷积的所得深层特征和前面下采样对应提取的特征向量进行拼接融合。网络的最后一层是一个3x3卷积,接着是一个激活函数Sigmoid。
在训练过程中,训练集中输入为PDW序列图像,目标输出为像素点分类好的图像。训练的目的是让U-Net更注重区分不同像素点间的差距,像素点不会丢失,PDW序列中像素点的位置也不会改变,即U-Net的输入Q和输出Q2对应。训练在基于PyTorch框架的基础上,采用了衡量模型输出与目标输出之间的交叉损失函数熵:
oi代表网络模型的输出,对应的是中PDW图进入网络后的预测结果,ti代表目标输出,对应PDW图像的像素点分类结果标签,I是网络计算一次损失所需要的样本数,σi是网络的需要学习的权重参数。
训练时模型的超参设置如表3所示:
表3模型超参设置
模型参数优化器选择适应性矩估计(adaptive moment estimation,Adam)函数,它利用梯度的一阶矩和二阶矩自适应调整每个参数的学习速率。默认的学习率初始值为0.001,对数据集的迭代次数为10,训练时当损失收敛达到0.01时停止训练保存训练好的模型。
(4)PDW搜索模块
经过深度分割模型对PDW序列图像的每个像素点的深度分割,可以得到PDW序列图像像素点分类结果。知道了像素点的分类结果,也就知道了脉冲的分类结果。最后需要做的就是将同属于一个辐射源信号的所有脉冲从PDW脉冲流中挑选出来,即选出每一类信号的辐射源描述字EDW,这就是PDW搜索。PDW搜索依靠的就是前面第二个步骤中构建的像素点索引表。
在搜索时遇到的大部分都是一个像素点对应一个脉冲,但也存在一个像素点索引多个脉冲情况,比如表1中像素点3对应2个脉冲。出现这种现象的原因可能是这两个脉冲的到达时间TOA和载频CF的数值相同或者相近,导致这两脉冲对应的像素点重叠,在这种情况下,利用像素点挑选脉冲的时候需要进行另外的参数匹配。
下面具体介绍一下对每一类像素点的PDW搜索方法:
(1)初筛:先将一个像素点索引一个脉冲的情况优先挑选归类出来,形成该辐射源信号的EDW;
(2)统计该辐射源信号的EDW范围:遍历EDW所有脉冲的参数信息,统计好该辐射源信号CF、PW、TOA的范围[CFmin,CFmax],[PWmin,PWmax],[TOAmin,TOAmax],若不存在一个像素点对应多个脉冲的情况则完成该辐射源信号的分选,可继续对下一类像素点进行搜索归类。若存在一个像素点对应多个脉冲的情况,则需要进参数匹配;
(3)参数匹配:将一个像素点对应的多个脉冲依次与上述统计辐射源信号的参数范围进行匹配,若某个脉冲参数均匹配成功,则将该脉冲挑选出来归类到这个辐射源信号的EDW中,完成对这类像素点的脉冲搜索归类操作。匹配不成功则存放至剩余脉冲组中,当所有类像素点完成对应脉冲搜索归类后,将个别剩余的脉冲与不同的辐射源EDW参数范围再次匹配归类,最终PDW搜索结束完成分选。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (5)
1.一种基于深度分割的端到端雷达信号分选方法,其特征在于,包括以下步骤:
步骤1:PDW采集,截取时间长度为T、带宽为B的PDW并保存为PDW脉冲序列;
步骤2:PDW序列图像生成;
选择PDW中CF、PW、TOA这三个参数进行图像变换处理,TOA、CF值用于描述像素点位置,PW值用于描述像素点颜色深度,从而构成tl(TOA)×fw(CF)×pd(PW)大小的三维脉冲矩阵P;同时为适应深度分割网络的图像的分辨率大小输入需求,需完成图像比例变换,构建图像矩阵Q,并根据对应的比例变换关系,构建像素点与脉冲之间的索引表;
步骤3:利用训练好的U-Net深度分割模型进行像素点分类,得到像素点分类结果图像矩阵Q2;
步骤4:PDW序列搜索;根据像素点的分类结果图像矩阵Q2和步骤2中得到的像素点与脉冲之间的索引表,按照像素点的分类情况对脉冲进行搜索归类;对于一个像素点索引多个脉冲的情况,采用参数匹配法,实现PDW的归类,获得最终不同辐射源信号的EDW。
2.根据权利要求1所述的一种基于深度分割的端到端雷达信号分选方法,其特征在于:所述步骤2中构建像素点与脉冲之间的索引表的方法具体为:
假设每个维度上脉冲与像素点对应的比例变换关系为β1,β2,β3,根据三维脉冲矩阵P以及脉冲序列图像矩阵Q,可以得到不同维度上的变换关系分别为β1,β2,β3:
若第n个脉冲的TOA和CF为{TOAn,CFn},它对应像素点在PDW序列图像的二维位置坐标则为:
(x,y)=(β1×TOAn,β2×CFn)
这样脉冲流中每个脉冲信号都对应PDW序列图像上一个像素点,进而构建像素点索引表,如表1所示:
表1像素点与脉冲之间的索引表
表中K表示K个脉冲像素点,每个像素点最多对应M个脉冲;脉冲参数{TOA,CF,PW}k,m的下角标k,m表示由第k个像素点可以搜索到m个脉冲,同时得到该脉冲的TOA、CF、PW参数值。
3.根据权利要求1所述的一种基于深度分割的端到端雷达信号分选方法,其特征在于:所述步骤3中U-Net深度分割模型具体为:
U-Net深度分割模型对输入的PDW序列图像的像素点进行深度分割,得到像素点分类结果图像Q2,U-Net深度分割模型的详细结构如表2所示:
表2 U-Net网络参数
U-Net深度分割过程分为下采样和上采样,各自包含5个网络层,每次下采样包括两个卷积核为3x3的卷积,每个卷积后面跟着一个激活函数ReLU和一个归一化函数BatchNormlization,两次卷积能够提取出图像的深层特征,2x2的最大池操作使图像尺寸缩小;在每个下采样步骤中,特征通道的数量增加一倍,图像尺寸减半;上采样过程采用反卷积结构,每个反卷积的卷积核为2x2,后面继续进行两个卷积核为3x3的卷积,同样每个卷积后面跟着一个激活函数ReLU和一个归一化函数BatchNormlization;每个上采样特征通道减半,图像尺寸放大一倍;同时为了优化输出和提高特征的利用率,需要将反卷积的所得深层特征和前面下采样对应提取的特征向量进行拼接融合;网络的最后一层是一个3x3卷积,接着是一个激活函数Sigmoid。
4.根据权利要求1所述的一种基于深度分割的端到端雷达信号分选方法,其特征在于:所述步骤3中U-Net深度分割模型的训练方法具体为:
在训练过程中,训练集中输入为PDW序列图像,目标输出为像素点分类好的图像;训练的目的是让U-Net更注重区分不同像素点间的差距,像素点不会丢失,PDW序列中像素点的位置也不会改变,即U-Net的输入Q和输出Q2对应;训练在基于PyTorch框架的基础上,采用了衡量模型输出与目标输出之间的交叉损失函数熵:
其中,oi代表网络模型的输出,对应的是中PDW图进入网络后的预测结果;ti代表目标输出,对应PDW图像的像素点分类结果标签;I是网络计算一次损失所需要的样本数;σi是网络的需要学习的权重参数;模型参数优化器选择Adam函数,它利用梯度的一阶矩和二阶矩自适应调整每个参数的学习速率;设定学习率初始值、对数据集的迭代次数、损失阈值,训练时当损失收敛达到损失阈值时停止训练保存训练好的模型。
5.根据权利要求1所述的一种基于深度分割的端到端雷达信号分选方法,其特征在于:所述步骤4中进行PDW序列搜索的方法具体为:
步骤4.1:初筛;先将一个像素点索引一个脉冲的情况优先挑选归类出来,形成该辐射源信号的EDW;
步骤4.2:统计该辐射源信号的EDW范围;遍历EDW所有脉冲的参数信息,统计好该辐射源信号CF、PW、TOA的范围[CFmin,CFmax],[PWmin,PWmax],[TOAmin,TOAmax],若不存在一个像素点对应多个脉冲的情况则完成该辐射源信号的分选,可继续对下一类像素点进行搜索归类;若存在一个像素点对应多个脉冲的情况,则需要进参数匹配;
步骤4.3:参数匹配:将一个像素点对应的多个脉冲依次与上述统计辐射源信号的参数范围进行匹配,若某个脉冲参数均匹配成功,则将该脉冲挑选出来归类到这个辐射源信号的EDW中,完成对这类像素点的脉冲搜索归类操作;匹配不成功则存放至剩余脉冲组中,当所有类像素点完成对应脉冲搜索归类后,将个别剩余的脉冲与不同的辐射源EDW参数范围再次匹配归类,最终PDW搜索结束完成分选。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111563815.7A CN114254141B (zh) | 2021-12-20 | 2021-12-20 | 一种基于深度分割的端到端雷达信号分选方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111563815.7A CN114254141B (zh) | 2021-12-20 | 2021-12-20 | 一种基于深度分割的端到端雷达信号分选方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114254141A true CN114254141A (zh) | 2022-03-29 |
CN114254141B CN114254141B (zh) | 2024-08-23 |
Family
ID=80793210
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111563815.7A Active CN114254141B (zh) | 2021-12-20 | 2021-12-20 | 一种基于深度分割的端到端雷达信号分选方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114254141B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114925721A (zh) * | 2022-04-27 | 2022-08-19 | 中国电子科技集团公司第二十九研究所 | 一种基于U-Net分割网络的信号分选方法 |
CN116821097A (zh) * | 2023-07-05 | 2023-09-29 | 中国人民解放军92728部队 | 全脉冲数据处理方法 |
CN117272086A (zh) * | 2023-11-22 | 2023-12-22 | 中国电子科技集团公司第二十九研究所 | 一种基于dbscan的雷达信号扫描包络分割方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110110738A (zh) * | 2019-03-20 | 2019-08-09 | 西安电子科技大学 | 一种基于多特征融合的雷达辐射源识别方法 |
CN111190146A (zh) * | 2020-01-13 | 2020-05-22 | 中国船舶重工集团公司第七二四研究所 | 一种基于视觉图形特征的复杂信号分选方法 |
EP3674741A1 (fr) * | 2018-12-27 | 2020-07-01 | Avantix | Système et procédé d'identification de source radar |
CN113156391A (zh) * | 2021-04-25 | 2021-07-23 | 电子科技大学 | 一种雷达信号多维特征智能分选方法 |
CN113759323A (zh) * | 2021-07-19 | 2021-12-07 | 西安电子科技大学 | 基于改进K-Means联合卷积自编码器的信号分选方法及装置 |
-
2021
- 2021-12-20 CN CN202111563815.7A patent/CN114254141B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP3674741A1 (fr) * | 2018-12-27 | 2020-07-01 | Avantix | Système et procédé d'identification de source radar |
CN110110738A (zh) * | 2019-03-20 | 2019-08-09 | 西安电子科技大学 | 一种基于多特征融合的雷达辐射源识别方法 |
CN111190146A (zh) * | 2020-01-13 | 2020-05-22 | 中国船舶重工集团公司第七二四研究所 | 一种基于视觉图形特征的复杂信号分选方法 |
CN113156391A (zh) * | 2021-04-25 | 2021-07-23 | 电子科技大学 | 一种雷达信号多维特征智能分选方法 |
CN113759323A (zh) * | 2021-07-19 | 2021-12-07 | 西安电子科技大学 | 基于改进K-Means联合卷积自编码器的信号分选方法及装置 |
Non-Patent Citations (1)
Title |
---|
陈涛;王天航;郭立民;: "基于时差多参分选的多层感知器网络脉间识别", 电子与信息学报, no. 07, 9 April 2018 (2018-04-09) * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114925721A (zh) * | 2022-04-27 | 2022-08-19 | 中国电子科技集团公司第二十九研究所 | 一种基于U-Net分割网络的信号分选方法 |
CN116821097A (zh) * | 2023-07-05 | 2023-09-29 | 中国人民解放军92728部队 | 全脉冲数据处理方法 |
CN116821097B (zh) * | 2023-07-05 | 2024-01-26 | 中国人民解放军92728部队 | 全脉冲数据处理方法 |
CN117272086A (zh) * | 2023-11-22 | 2023-12-22 | 中国电子科技集团公司第二十九研究所 | 一种基于dbscan的雷达信号扫描包络分割方法 |
CN117272086B (zh) * | 2023-11-22 | 2024-02-13 | 中国电子科技集团公司第二十九研究所 | 一种基于dbscan的雷达信号扫描包络分割方法 |
Also Published As
Publication number | Publication date |
---|---|
CN114254141B (zh) | 2024-08-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110443143B (zh) | 多分支卷积神经网络融合的遥感图像场景分类方法 | |
CN114254141B (zh) | 一种基于深度分割的端到端雷达信号分选方法 | |
CN108647602B (zh) | 一种基于图像复杂度判定的航空遥感图像场景分类方法 | |
CN109446922B (zh) | 一种实时鲁棒的人脸检测方法 | |
CN103927531A (zh) | 一种基于局部二值和粒子群优化bp神经网络的人脸识别方法 | |
CN112990282B (zh) | 一种细粒度小样本图像的分类方法及装置 | |
CN104657980A (zh) | 一种改进的基于Meanshift的多通道图像分割算法 | |
CN110610717A (zh) | 复杂频谱环境混合信号的分离方法 | |
CN109948722B (zh) | 一种识别空间目标的方法 | |
Bi et al. | Multiple instance dense connected convolution neural network for aerial image scene classification | |
CN112084897A (zh) | 一种gs-ssd的交通大场景车辆目标快速检测方法 | |
CN109508639A (zh) | 基于多尺度带孔卷积神经网络的道路场景语义分割方法 | |
Xu et al. | TransMIN: Transformer-guided multi-interaction network for remote sensing object detection | |
Sun et al. | Deep learning based pedestrian detection | |
CN107451617A (zh) | 一种图转导半监督分类方法 | |
CN112966544B (zh) | 一种采用ICGAN与ResNet网络的雷达辐射源信号分类识别方法 | |
Nguyen et al. | FoodMask: Real-time food instance counting, segmentation and recognition | |
Tang et al. | A multi-view SAR target recognition method using feature fusion and joint classification | |
Zheng et al. | Boundary adjusted network based on cosine similarity for temporal action proposal generation | |
CN116500546A (zh) | 一种基于点云分割网络的雷达信号分选方法 | |
Wang et al. | Hierarchical Kernel Interaction Network for Remote Sensing Object Counting | |
Egbueze et al. | Weather Recognition Based on Still Images Using Deep Learning Neural Network with Resnet-15 | |
Qu et al. | Frequency-Hopping Signal Sorting Based on Deep Leaming | |
Basrur et al. | Food Recognition using Transfer Learning | |
CN114241380A (zh) | 一种基于类别标签和属性注释的多任务属性场景识别方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |