CN108596156A - 一种智能sar雷达空中飞行目标识别系统 - Google Patents
一种智能sar雷达空中飞行目标识别系统 Download PDFInfo
- Publication number
- CN108596156A CN108596156A CN201810457998.6A CN201810457998A CN108596156A CN 108596156 A CN108596156 A CN 108596156A CN 201810457998 A CN201810457998 A CN 201810457998A CN 108596156 A CN108596156 A CN 108596156A
- Authority
- CN
- China
- Prior art keywords
- feature
- target
- module
- sar
- indicate
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V20/00—Scenes; Scene-specific elements
- G06V20/10—Terrestrial scenes
- G06V20/13—Satellite images
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
- G06T7/62—Analysis of geometric attributes of area, perimeter, diameter or volume
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/20—Image preprocessing
- G06V10/26—Segmentation of patterns in the image field; Cutting or merging of image elements to establish the pattern region, e.g. clustering-based techniques; Detection of occlusion
- G06V10/267—Segmentation of patterns in the image field; Cutting or merging of image elements to establish the pattern region, e.g. clustering-based techniques; Detection of occlusion by performing operations on regions, e.g. growing, shrinking or watersheds
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Multimedia (AREA)
- Geometry (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Astronomy & Astrophysics (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种智能SAR雷达空中飞行目标识别系统,包括SAR雷达、数据库以及上位机;SAR雷达、数据库、上位机依次相连,所述SAR雷达对空中进行实时监测,并将SAR雷达获得的图像数据存储到所述的数据库中,所述的上位机包括图像预处理模块、特征提取模块、特征选择模块、分类器训练模块、智能寻优模块以及结果显示模块。本发明提供一种实现在线识别、精度高的空中飞行目标识别系统。
Description
技术领域
本发明涉及雷达数据处理领域,特别地,涉及一种智能SAR雷达空中飞行目标识别系统。
背景技术
现代雷达的作用是通过探测获取目标位置及运动方面的数据,而目标的形状、尺寸等物理特征数据,只有通过特征识别、成像识别或者分类识别才可以得到。无源雷达可以获得目标的电磁辐射信号用来实现对目标位置的确定,跟踪和识别,可以弥补常规雷达识别机型和飞机架数中出现的错误。国外在目标识别方面的研究已经形成了一套成熟的理论,包括理论和软硬件多方面的研究,并产生了一批非常经典的识别体系。目前各国主要的研究对象是无源雷达识别、雷达识别、红外识别、反侦察识别系统、侦察识别和综合识别系统。另外,国外正在研究电子接收机、红外、激光传感器、视频摄像等光电设备,主要对目标成像进行识别。美国在一次战争中使用一种近红外装置,配合地面作战运载工具,因此较好地获得了制空权,正确地得到了目标的属性。其中基于SAR图像的飞行目标检测已经有了广泛的研究,而飞行目标分类识别由于SAR图像分辨率的限制才刚刚起步,已有的一些研究成果也由于研究不够透彻,模型的效果并不是很好。因此,积极开展基于高分辨率SAR图像的飞行目标分类识别研究具有极其重大的意义。
发明内容
为了克服目前基于SAR图像的空中飞行目标识别准确率不高的不足,本发明的目的在于提供一种实现实时分析的SAR雷达空中飞行目标识别系统。
本发明解决其技术问题所采用的技术方案是:一种智能SAR雷达空中飞行目标识别系统,包括SAR雷达、数据库以及上位机,SAR雷达、数据库和上位机依次相连,所述SAR雷达对空中进行实时监测,并将SAR雷达获得的图像数据存储到所述的数据库中,所述的上位机包括:
图像预处理模块,用以进行SAR雷达图像数据预处理,采用如下过程完成:
1)从数据库中传来的SAR图像灰度级为L,f(x0,y0)为像素点(x0,y0)处的灰度值,g(x0,y0)为像素点(x0,y0)的N×N邻域内像素的平均值,其中x0,y0分别表示像素点的横坐标和纵坐标;
2)通过计算满足f=m和g=n的像素数目h(m,n),得到二维联合概率密度pmn:
pmn=p(m,n)=h(m,n)/M
其中,M表示图像像素的总数目;
3)计算二维直方图的均值向量μ:
4)分别计算图像中目标和背景出现的概率P0,1和均值向量μ0,1:
其中,t、s、下标0、下标1分别表示f分割阈值、g分割阈值、目标区域、背景区域;
5)计算类间方差BCV:
BCV=P0(μ0-μ)(μ0-μ)′+P1(μ1-μ)(μ1-μ)′;
其中,μ表示均值向量,上标’表示矩阵的转置。
6)最佳阈值即为使得BCV为最大值时的二维阈值向量[s0,t0]:
特征提取模块,用以进行飞行目标典型特征的提取,采用如下过程完成:
1)从图像预处理模块传来的只包含一个飞行目标的SAR图像切片I(m,n),其中只包含目标区域的二值图为B(m,n),则只包含目标的图像T(m,n):
T(m,n)=I(m,n)×B(m,n)
其中,×表示对应像素相乘;
2)在B(m,n)中根据飞行目标个体的主轴方向求得飞行目标主体区域的最小外接矩形,则该矩形的长边长度Length即为飞行目标个体的长度,矩形的短边长度Width即为飞行目标个体的宽度;
3)计算得到几何结构特征,其中包括周长、面积、长宽比、形状复杂度、目标质心位置以及转动惯量:
周长面积长宽比R=Length/Width;形状复杂度C=Length2/4πS;目标区域的质心位置
转动惯量式中,r代表了目标像素点与质心之间的距离,
4)计算得到灰度统计特征,其中包括质量、均值、方差系数、标准差、分形维数、加权填充比:
质量均值方差系数标准差式中分别表示灰度对数和、灰度对数平方和;分形维数H=(log10N1-log10N2)/(log10d1-log10d2),该特征的计算方法是:用分割后的SAR图像切片构建一个保留了目标区域的K(这里取K=50)个最亮像素点的二值图B2(m,n),首先将一个大小为d1×d2的窗口在这个二值图中连续滑动,记下窗口中包含亮点的窗口总数记为N1,接着再用一个大小为d2×d2的窗口在这个二值图中连续滑动,记下在窗口中包含亮点的窗口总数记为N2;加权填充比
特征选择模块,用以选择出最优特征子集,采用如下过程完成:
1)计算类内距离类间距离以及类内类间距Ji:
其中,i表示特征标号,ω表示飞行目标类别的标号,表示特征向量的2范数,表示训练集样本的总体均值,Nω表示第ω类飞行目标的数量,N表示训练集中飞行目标总数,E表示期望,下标W、下标B分别表示类内、类间。
2)计算得到归一化方差系数
其中,i表示特征标号,ω表示飞行目标类别的标号,表示特征向量的2范数,和分别表示特征的平方均值以及均值的平方。特征的方差系数越小,表明该目标特征的稳定性越好;
3)计算得到相关系数ri,j:
其中,i,j表示特征标号,||Fi||2表示特征Fi的2范数,和分别表示Fi和Fj的均值,σi,i和σj,j分别表示Fi和Fj的标准差。由相关系数的性质可知,0≤ri,j≤1;如果两个特征完全不相关,ri,j=0;如果两个特征完全相关,ri,j=1;如果两个特征之间的相关性很低,即特征之间的信息冗余非常少,那么ri,j就会越接近0;反之,如果两个特征之间的相关性很高,即特征之间的信息冗余非常多,那么ri,j就会越接近1;
4)通过上述得到的类内类间距、归一化方差系数、相关系数筛选出最优特征子集,构造最优输入特征向量;
分类器训练模块,用以进行分类器训练,采用如下过程完成:
1)从特征选择模块中采集N个SAR雷达图像xi作为训练样本,i=1,2,…,N;
2)对训练样本进行归一化处理,得到归一化样本
3)将归一化后的训练样本重构,分别得到输入矩阵X和对应的输出矩阵Y:
其中,D表示重构维数,D为自然数,且D<N,D的取值范围为50-70;
4)得到分类器训练集合其中xi代表n×1的输入特征向量,ti代表m×1的目标向量。给定一个激活函数g(x)以及隐含层的节点个数那么ELM分类器为:
其中,ωi代表第i个隐含层节点和输入层之间的权重向量,bi代表第i个隐含层节点的偏置,βi代表第i个隐含层节点和输出层之间的权重向量,oj代表第j个输入数据的目标输出。另外,ωi·xj代表了ωi和xj的内积。
该网络的输出能够无限地接近于输入的N个样本,即:
可得:
上式可以表达为矩阵形式:Hβ=T
其中,H表示隐含层的输出矩阵,H的第i列分别表示隐含层的第i个节点对应于N个输入x1,x2,…,xN的输出值。单隐层前馈神经网络(SLFNs)的输入权值和隐含层的偏差在网络训练的过程中不需要调整,可以任意给定。基于上述理论,输出权重可以通过计算Hβ=T的最小二乘解求得:
可以利用线性方法快速求得方程的解,如式所示:
其中,代表H的Moore-Penrose广义逆矩阵,代表最小范数最小二乘解,它正好是最小二乘解中范数最小的解。相比很多已有的分类器系统,极限学习机通过这种Moore-Penrose广义逆的求解能够以非常快的速度达到很好的训练效果。
智能寻优模块,用以采用遗传算法对分类器的核参数θ和惩罚因子γ进行优化,采用如下过程完成:
1)采用实数编码方法对θ和γ进行编码;
2)随机产生初始种群;
3)采用下式计算种群每个个体的适应度,并判断是否符合算法终止条件,若符合,输出最佳个体及其代表的最优解,并结束计算,否则继续迭代;
其中,为待估计函数,Cmax为f(x)的最大估计;
4)采用正态分布概率选择个体;
5)通过单点线性交叉产生新个体;
6)通过均匀变异方式产生新个体;
7)产生新一代种群,返回3)进行迭代;
其中,初始种群大小为50-200,最大代数50-300,最佳个体选择概率为0.05-0.1,交叉概率为0.5-0.9,变异概率为0.001-0.01,个体适应度选择模型的均方根误差,终止条件为达到最大代数或者连续五次适应度不变。
结果显示模块,用以进行识别结果的显示,即将输入SAR图像中飞行目标的类型显示在屏幕当中。
本发明的技术构思为:本发明针对SAR雷达全天时、全天候工作以及强穿透的特性,对SAR雷达监测到的空中图像进行图像预处理,接着进行特征的提取以及特征的选择,最后通过分类器的训练过程建立空中飞行目标识别模型,从而实现SAR雷达空中飞行目标的识别。
本发明的有益效果主要表现在:1、可实时识别空中飞行目标;2、所用的识别方法只需较少的训练样本;3、智能化、受人为因素干扰小。
附图说明
图1是本发明所提出的系统的整体结构图;
图2是本发明所提出的上位机的功能模块图。
具体实施方式
下面根据附图具体说明本发明。上述实施例用来解释说明本发明,而不是对本发明进行限制,在本发明的精神和权利要求的保护范围内,对本发明作出的任何修改和改变,都落入本发明的保护范围。
实施例
参照图1、图2,一种智能SAR雷达空中飞行目标识别系统,包括SAR雷达1、数据库2及上位机3,SAR雷达1、数据库2和上位机3依次相连,所述SAR雷达1对所监测空中进行照射,并将SAR雷达图像存储到所述的数据库2,所述的上位机3包括:
图像预处理模块4,用以进行SAR雷达图像数据预处理,采用如下过程完成:
1)从数据库中传来的SAR图像灰度级为L,f(x0,y0)为像素点(x0,y0)处的灰度值,g(x0,y0)为像素点(x0,y0)的N×N邻域内像素的平均值,其中x0,y0分别表示像素点的横坐标和纵坐标;
2)通过计算满足f=m和g=n的像素数目h(m,n),得到二维联合概率密度pmn:
pmn=p(m,n)=h(m,n)/M
其中,M表示图像像素的总数目;
3)计算二维直方图的均值向量μ:
4)分别计算图像中目标和背景出现的概率P0,1和均值向量μ0,1:
其中,t、s、下标0、下标1分别表示f分割阈值、g分割阈值、目标区域、背景区域;
5)计算类间方差BCV:
BCV=P0(μ0-μ)(μ0-μ)′+P1(μ1-μ)(μ1-μ)′;
其中,μ表示均值向量,上标’表示矩阵的转置。
6)最佳阈值即为使得BCV为最大值时的二维阈值向量[s0,t0]:
特征提取模块5,用以进行飞行目标典型特征的提取,采用如下过程完成:
1)从图像预处理模块传来的只包含一个飞行目标的SAR图像切片I(m,n),其中只包含目标区域的二值图为B(m,n),则只包含目标的图像T(m,n):
T(m,n)=I(m,n)×B(m,n)
其中,×表示对应像素相乘;
2)在B(m,n)中根据飞行目标个体的主轴方向求得飞行目标主体区域的最小外接矩形,则该矩形的长边长度Length即为飞行目标个体的长度,矩形的短边长度Width即为飞行目标个体的宽度;
3)计算得到几何结构特征,其中包括周长、面积、长宽比、形状复杂度、目标质心位置以及转动惯量:
周长面积长宽比R=Length/Width;形状复杂度C=Length2/4πS;目标区域的质心位置
转动惯量式中,r代表了目标像素点与质心之间的距离,
4)计算得到灰度统计特征,其中包括质量、均值、方差系数、标准差、分形维数、加权填充比:
质量均值方差系数标准差式中分别表示灰度对数和、灰度对数平方和;分形维数H=(log10N1-log10N2)/(log10d1-log10d2),该特征的计算方法是:用分割后的SAR图像切片构建一个保留了目标区域的K(这里取K=50)个最亮像素点的二值图B2(m,n),首先将一个大小为d1×d2的窗口在这个二值图中连续滑动,记下窗口中包含亮点的窗口总数记为N1,接着再用一个大小为d2×d2的窗口在这个二值图中连续滑动,记下在窗口中包含亮点的窗口总数记为N2;
加权填充比
特征选择模块6,用以选择出最优特征子集,采用如下过程完成:
1)计算类内距离类间距离以及类内类间距Ji:
其中,i表示特征标号,ω表示飞行目标类别的标号,表示特征向量的2范数,表示训练集样本的总体均值,Nω表示第ω类飞行目标的数量,N表示训练集中飞行目标总数,E表示期望,下标W、下标B分别表示类内、类间。
2)计算得到归一化方差系数
其中,i表示特征标号,ω表示飞行目标类别的标号,表示特征向量的2范数,和分别表示特征的平方均值以及均值的平方。特征的方差系数越小,表明该目标特征的稳定性越好;
3)计算得到相关系数ri,j:
其中,i,j表示特征标号,||Fi||2表示特征Fi的2范数,和分别表示Fi和Fj的均值,σi,i和σj,j分别表示Fi和Fj的标准差。由相关系数的性质可知,0≤ri,j≤1;如果两个特征完全不相关,ri,j=0;如果两个特征完全相关,ri,j=1;如果两个特征之间的相关性很低,即特征之间的信息冗余非常少,那么ri,j就会越接近0;反之,如果两个特征之间的相关性很高,即特征之间的信息冗余非常多,那么ri,j就会越接近1;
4)通过上述得到的类内类间距、归一化方差系数、相关系数筛选出最优特征子集,构造最优输入特征向量;
分类器训练模块7,用以进行分类器训练,采用如下过程完成:
1)从特征选择模块中采集N个SAR雷达图像xi作为训练样本,i=1,2,…,N;
2)对训练样本进行归一化处理,得到归一化样本
3)将归一化后的训练样本重构,分别得到输入矩阵X和对应的输出矩阵Y:
其中,D表示重构维数,D为自然数,且D<N,D的取值范围为50-70;
4)得到分类器训练集合其中xi代表n×1的输入特征向量,ti代表m×1的目标向量。给定一个激活函数g(x)以及隐含层的节点个数那么ELM分类器为:
其中,ωi代表第i个隐含层节点和输入层之间的权重向量,bi代表第i个隐含层节点的偏置,βi代表第i个隐含层节点和输出层之间的权重向量,oj代表第j个输入数据的目标输出。另外,ωi·xj代表了ωi和xj的内积。
该网络的输出能够无限地接近于输入的N个样本,即:
可得:
上式可以表达为矩阵形式:Hβ=T
其中,H表示隐含层的输出矩阵,H的第i列分别表示隐含层的第i个节点对应于N个输入x1,x2,…,xN的输出值。单隐层前馈神经网络(SLFNs)的输入权值和隐含层的偏差在网络训练的过程中不需要调整,可以任意给定。基于上述理论,输出权重可以通过计算Hβ=T的最小二乘解求得:
可以利用线性方法快速求得方程的解,如式所示:
其中,代表H的Moore-Penrose广义逆矩阵,代表最小范数最小二乘解,它正好是最小二乘解中范数最小的解。相比很多已有的分类器系统,极限学习机通过这种Moore-Penrose广义逆的求解能够以非常快的速度达到很好的训练效果。
智能寻优模块9,用以采用遗传算法对分类器的核参数θ和惩罚因子γ进行优化,采用如下过程完成:
1)采用实数编码方法对θ和γ进行编码;
2)随机产生初始种群;
3)采用下式计算种群每个个体的适应度,并判断是否符合算法终止条件,若符合,输出最佳个体及其代表的最优解,并结束计算,否则继续迭代;
其中,为待估计函数,Cmax为f(x)的最大估计;
4)采用正态分布概率选择个体;
5)通过单点线性交叉产生新个体;
6)通过均匀变异方式产生新个体;
7)产生新一代种群,返回3)进行迭代;
其中,初始种群大小为50-200,最大代数50-300,最佳个体选择概率为0.05-0.1,交叉概率为0.5-0.9,变异概率为0.001-0.01,个体适应度选择模型的均方根误差,终止条件为达到最大代数或者连续五次适应度不变。
结果显示模块8,用以进行识别结果的显示,即将输入SAR图像中飞行目标的类型显示在上位机当中。
所述上位机3的硬件部分包括:I/O元件,用于数据的采集和信息的传递;数据存储器,存储运行所需要的数据样本和运行参数等;程序存储器,存储实现功能模块的软件程序;运算器,执行程序,实现指定的功能;显示模块,显示设置的参数和识别结果。
上述实施例用来解释说明本发明,而不是对本发明进行限制,在本发明的精神和权利要求的保护范围内,对本发明作出的任何修改和改变,都落入本发明的保护范围。
Claims (5)
1.一种智能SAR雷达空中飞行目标识别系统,其特征在于:包括SAR雷达、数据库以及上位机,SAR雷达、数据库和上位机依次相连,所述SAR雷达对空中进行实时监测,并将SAR雷达获得的图像数据存储到所述的数据库中,所述的上位机包括图像预处理模块、特征提取模块、分类器训练模块、分类器训练模块、群智优化模块和结果显示模块,所述图像预处理模块、特征提取模块、分类器训练模块、分类器训练模块和结果显示模块依次相连,分类器训练模块与群智优化模块相连。
2.根据权利要求1所述智能SAR雷达空中飞行目标识别系统,其特征在于:所述图像预处理模块用以进行SAR雷达图像数据预处理,采用如下过程完成:
1)从数据库中传来的SAR图像灰度级为L,f(x0,y0)为像素点(x0,y0)处的灰度值,g(x0,y0)为像素点(x0,y0)的N×N邻域内像素的平均值,其中x0,y0分别表示像素点的横坐标和纵坐标;
2)通过计算满足f=m和g=n的像素数目h(m,n),得到二维联合概率密度pmn:
pmn=p(m,n)=h(m,n)/M
其中,M表示图像像素的总数目;
3)计算二维直方图的均值向量μ:
4)分别计算图像中目标和背景出现的概率P0,1和均值向量μ0,1:
其中,t、s、下标0、下标1分别表示f分割阈值、g分割阈值、目标区域、背景区域;
5)计算类间方差BCV:
BCV=P0(μ0-μ)(μ0-μ)′+P1(μ1-μ)(μ1-μ)′;
其中,μ表示均值向量,上标’表示矩阵的转置。
6)最佳阈值即为使得BCV为最大值时的二维阈值向量[s0,t0]:
3.根据权利要求1所述智能SAR雷达空中飞行目标识别系统,其特征在于:所述特征提取模块,用以进行飞行目标典型特征的提取,采用如下过程完成:
1)从图像预处理模块传来的只包含一个飞行目标的SAR图像切片I(m,n),其中只包含目标区域的二值图为B(m,n),则只包含目标的图像T(m,n):
T(m,n)=I(m,n)×B(m,n)
其中,×表示对应像素相乘;
2)在B(m,n)中根据飞行目标个体的主轴方向求得飞行目标主体区域的最小外接矩形,则该矩形的长边长度Length即为飞行目标个体的长度,矩形的短边长度Width即为飞行目标个体的宽度;
3)计算得到几何结构特征,其中包括周长、面积、长宽比、形状复杂度、目标质心位置以及转动惯量:
周长面积长宽比R=Length/Width;形状复杂度C=Length2/4πS;目标区域的质心位置
转动惯量式中,r代表了目标像素点与质心之间的距离,
4)计算得到灰度统计特征,其中包括质量、均值、方差系数、标准差、分形维数、加权填充比:
质量均值方差系数标准差式中分别表示灰度对数和、灰度对数平方和;分形维数H=(log10N1-log10N2)/(log10d1-log10d2),该特征的计算方法是:用分割后的SAR图像切片构建一个保留了目标区域的K(这里取K=50)个最亮像素点的二值图B2(m,n),首先将一个大小为d1×d2的窗口在这个二值图中连续滑动,记下窗口中包含亮点的窗口总数记为N1,接着再用一个大小为d2×d2的窗口在这个二值图中连续滑动,记下在窗口中包含亮点的窗口总数记为N2;加权填充比
4.根据权利要求1所述智能SAR雷达空中飞行目标识别系统,其特征在于:所述特征选择模块用以选择出最优特征子集,采用如下过程完成:
1)计算类内距离类间距离以及类内类间距Ji:
其中,i表示特征标号,ω表示飞行目标类别的标号,||Fi (ω)||2表示特征向量Fi (ω)的2范数,表示训练集样本的总体均值,Nω表示第ω类飞行目标的数量,N表示训练集中飞行目标总数,E表示期望,下标W、下标B分别表示类内、类间。
2)计算得到归一化方差系数ρi (ω):
ρi (ω)=E[||Fi (ω)||2 2]-E2[||Fi (ω)||2]/E[||Fi (ω)||2 2]
其中,i表示特征标号,ω表示飞行目标类别的标号,||Fi (ω)||2表示特征向量Fi (ω)的2范数,E[||Fi (ω)||2 2]和E2[||Fi (ω)||2]分别表示特征的平方均值以及均值的平方。特征的方差系数ρi (ω)越小,表明该目标特征的稳定性越好;
3)计算得到相关系数ri,j:
其中,i,j表示特征标号,||Fi||2表示特征Fi的2范数,和分别表示Fi和Fj的均值,σi,i和σj,j分别表示Fi和Fj的标准差。由相关系数的性质可知,0≤ri,j≤1;如果两个特征完全不相关,ri,j=0;如果两个特征完全相关,ri,j=1;如果两个特征之间的相关性很低,即特征之间的信息冗余非常少,那么ri,j就会越接近0;反之,如果两个特征之间的相关性很高,即特征之间的信息冗余非常多,那么ri,j就会越接近1;
4)通过上述得到的类内类间距、归一化方差系数、相关系数筛选出最优特征子集,构造最优输入特征向量。
5.根据权利要求1所述智能SAR雷达空中飞行目标识别系统,其特征在于:所述分类器训练模块用以进行分类器训练,采用如下过程完成:
1)从特征选择模块中采集N个SAR雷达图像xi作为训练样本,i=1,2,…,N;
2)对训练样本进行归一化处理,得到归一化样本
3)将归一化后的训练样本重构,分别得到输入矩阵X和对应的输出矩阵Y:
其中,D表示重构维数,D为自然数,且D<N,D的取值范围为50-70;
4)得到分类器训练集合其中xi代表n×1的输入特征向量,ti代表m×1的目标向量。给定一个激活函数g(x)以及隐含层的节点个数那么ELM分类器为:
其中,ωi代表第i个隐含层节点和输入层之间的权重向量,bi代表第i个隐含层节点的偏置,βi代表第i个隐含层节点和输出层之间的权重向量,oj代表第j个输入数据的目标输出。另外,ωi·xj代表了ωi和xj的内积。
该网络的输出能够无限地接近于输入的N个样本,即:
可得:
上式可以表达为矩阵形式:Hβ=T
其中,H表示隐含层的输出矩阵,H的第i列分别表示隐含层的第i个节点对应于N个输入x1,x2,…,xN的输出值。单隐层前馈神经网络(SLFNs)的输入权值和隐含层的偏差在网络训练的过程中不需要调整,可以任意给定。基于上述理论,输出权重可以通过计算Hβ=T的最小二乘解求得:
可以利用线性方法快速求得方程的解,如式所示:
其中,代表H的Moore-Penrose广义逆矩阵,代表最小范数最小二乘解,它正好是最小二乘解中范数最小的解。相比很多已有的分类器系统,极限学习机通过这种Moore-Penrose广义逆的求解能够以非常快的速度达到很好的训练效果。
所述智能寻优模块用以采用遗传算法对分类器的核参数θ和惩罚因子γ进行优化,采用如下过程完成:
1)采用实数编码方法对θ和γ进行编码;
2)随机产生初始种群;
3)采用下式计算种群每个个体的适应度,并判断是否符合算法终止条件,若符合,输出最佳个体及其代表的最优解,并结束计算,否则继续迭代;
其中,为待估计函数,Cmax为f(x)的最大估计;
4)采用正态分布概率选择个体;
5)通过单点线性交叉产生新个体;
6)通过均匀变异方式产生新个体;
7)产生新一代种群,返回3)进行迭代;
其中,初始种群大小为50-200,最大代数50-300,最佳个体选择概率为0.05-0.1,交叉概率为0.5-0.9,变异概率为0.001-0.01,个体适应度选择模型的均方根误差,终止条件为达到最大代数或者连续五次适应度不变。
所述结果显示模块用以进行识别结果的显示,即将输入SAR图像中飞行目标的类型显示在屏幕当中。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810457998.6A CN108596156A (zh) | 2018-05-14 | 2018-05-14 | 一种智能sar雷达空中飞行目标识别系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810457998.6A CN108596156A (zh) | 2018-05-14 | 2018-05-14 | 一种智能sar雷达空中飞行目标识别系统 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108596156A true CN108596156A (zh) | 2018-09-28 |
Family
ID=63637471
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810457998.6A Pending CN108596156A (zh) | 2018-05-14 | 2018-05-14 | 一种智能sar雷达空中飞行目标识别系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108596156A (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110288001A (zh) * | 2019-05-28 | 2019-09-27 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | 基于目标数据特征训练学习的目标识别方法 |
CN112926381A (zh) * | 2021-01-08 | 2021-06-08 | 浙江大学 | 一种基于激光检测技术的水下目标识别系统 |
CN112926383A (zh) * | 2021-01-08 | 2021-06-08 | 浙江大学 | 一种基于水下激光图像的目标自动识别系统 |
CN112925932A (zh) * | 2021-01-08 | 2021-06-08 | 浙江大学 | 一种高清晰的水下激光图像处理系统 |
CN112946682A (zh) * | 2021-01-08 | 2021-06-11 | 浙江大学 | 一种高准确率高敏捷性的水下激光目标识别仪 |
CN113534093A (zh) * | 2021-08-13 | 2021-10-22 | 北京环境特性研究所 | 飞机目标的螺旋桨叶片数量反演方法及目标识别方法 |
CN113705415A (zh) * | 2021-08-23 | 2021-11-26 | 中国电子科技集团公司第十五研究所 | 基于雷达情报的空情目标特征提取方法及装置 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1916792A (zh) * | 2006-09-12 | 2007-02-21 | 浙江大学 | 一种工业生产过程小样本条件下的软测量方法 |
CN102147463A (zh) * | 2011-03-03 | 2011-08-10 | 浙江大学 | 一种群智雷达海杂波预报系统及方法 |
CN102147464A (zh) * | 2011-03-03 | 2011-08-10 | 浙江大学 | 一种鲁棒智能雷达海杂波预报系统及方法 |
CN102183745A (zh) * | 2011-03-03 | 2011-09-14 | 浙江大学 | 一种智能雷达海杂波预报系统及方法 |
CN102183751A (zh) * | 2011-03-03 | 2011-09-14 | 浙江大学 | 一种智能雷达海上目标检测系统及方法 |
CN105260786A (zh) * | 2015-09-08 | 2016-01-20 | 哈尔滨工程大学 | 一种电力推进系统仿真可信度评估模型综合优化方法 |
CN105354363A (zh) * | 2015-10-09 | 2016-02-24 | 上海大学 | 基于核极限学习机的脉动风速预测方法 |
CN106295686A (zh) * | 2016-08-01 | 2017-01-04 | 安徽农业大学 | 参数优化支持向量机的茶叶香气分类方法 |
-
2018
- 2018-05-14 CN CN201810457998.6A patent/CN108596156A/zh active Pending
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1916792A (zh) * | 2006-09-12 | 2007-02-21 | 浙江大学 | 一种工业生产过程小样本条件下的软测量方法 |
CN102147463A (zh) * | 2011-03-03 | 2011-08-10 | 浙江大学 | 一种群智雷达海杂波预报系统及方法 |
CN102147464A (zh) * | 2011-03-03 | 2011-08-10 | 浙江大学 | 一种鲁棒智能雷达海杂波预报系统及方法 |
CN102183745A (zh) * | 2011-03-03 | 2011-09-14 | 浙江大学 | 一种智能雷达海杂波预报系统及方法 |
CN102183751A (zh) * | 2011-03-03 | 2011-09-14 | 浙江大学 | 一种智能雷达海上目标检测系统及方法 |
CN105260786A (zh) * | 2015-09-08 | 2016-01-20 | 哈尔滨工程大学 | 一种电力推进系统仿真可信度评估模型综合优化方法 |
CN105354363A (zh) * | 2015-10-09 | 2016-02-24 | 上海大学 | 基于核极限学习机的脉动风速预测方法 |
CN106295686A (zh) * | 2016-08-01 | 2017-01-04 | 安徽农业大学 | 参数优化支持向量机的茶叶香气分类方法 |
Non-Patent Citations (2)
Title |
---|
JUN WU 等: "A novel ship classification approach for high resolution SAR images based on the BDA-KELM classification model", 《INTERNATIONAL JOURNAL OF REMOTE SENSING》 * |
陈文婷: "SAR图像舰船目标特征提取与分类识别方法研究", 《中国优秀硕士学位论文全文数据库(电子期刊)基础科学辑》 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110288001A (zh) * | 2019-05-28 | 2019-09-27 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | 基于目标数据特征训练学习的目标识别方法 |
CN110288001B (zh) * | 2019-05-28 | 2023-09-05 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | 基于目标数据特征训练学习的目标识别方法 |
CN112926381A (zh) * | 2021-01-08 | 2021-06-08 | 浙江大学 | 一种基于激光检测技术的水下目标识别系统 |
CN112926383A (zh) * | 2021-01-08 | 2021-06-08 | 浙江大学 | 一种基于水下激光图像的目标自动识别系统 |
CN112925932A (zh) * | 2021-01-08 | 2021-06-08 | 浙江大学 | 一种高清晰的水下激光图像处理系统 |
CN112946682A (zh) * | 2021-01-08 | 2021-06-11 | 浙江大学 | 一种高准确率高敏捷性的水下激光目标识别仪 |
CN112926381B (zh) * | 2021-01-08 | 2022-06-21 | 浙江大学 | 一种基于激光检测技术的水下目标识别系统 |
CN112926383B (zh) * | 2021-01-08 | 2023-03-03 | 浙江大学 | 一种基于水下激光图像的目标自动识别系统 |
CN113534093A (zh) * | 2021-08-13 | 2021-10-22 | 北京环境特性研究所 | 飞机目标的螺旋桨叶片数量反演方法及目标识别方法 |
CN113534093B (zh) * | 2021-08-13 | 2023-06-27 | 北京环境特性研究所 | 飞机目标的螺旋桨叶片数量反演方法及目标识别方法 |
CN113705415A (zh) * | 2021-08-23 | 2021-11-26 | 中国电子科技集团公司第十五研究所 | 基于雷达情报的空情目标特征提取方法及装置 |
CN113705415B (zh) * | 2021-08-23 | 2023-10-27 | 中国电子科技集团公司第十五研究所 | 基于雷达情报的空情目标特征提取方法及装置 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108596156A (zh) | 一种智能sar雷达空中飞行目标识别系统 | |
CN108717568B (zh) | 一种基于三维卷积神经网络的图像特征提取与训练方法 | |
CN106999107B (zh) | 身体成像 | |
Li et al. | An intelligent optimization method of motion management system based on BP neural network | |
CN107229904A (zh) | 一种基于深度学习的目标检测与识别方法 | |
CN110472627A (zh) | 一种端到端的sar图像识别方法、装置及存储介质 | |
CN107169435A (zh) | 一种基于雷达仿真图像的卷积神经网络人体动作分类方法 | |
CN108490415A (zh) | 一种智能sar雷达海上舰船目标识别系统 | |
CN109978882A (zh) | 一种基于多模态融合的医疗影像目标检测方法 | |
Salehi et al. | Generative adversarial networks (GANs): An overview of theoretical model, evaluation metrics, and recent developments | |
CN109214366A (zh) | 局部目标重识别方法、装置及系统 | |
Jiang et al. | An eight-layer convolutional neural network with stochastic pooling, batch normalization and dropout for fingerspelling recognition of Chinese sign language | |
CN109815814B (zh) | 一种基于卷积神经网络的人脸检测方法 | |
CN115661943B (zh) | 一种基于轻量级姿态评估网络的跌倒检测方法 | |
CN106326857A (zh) | 基于人脸图像的性别识别方法及装置 | |
CN110060233A (zh) | 一种玉米果穗破损检测方法 | |
CN106778474A (zh) | 3d人体识别方法及设备 | |
CN104298974A (zh) | 一种基于深度视频序列的人体行为识别方法 | |
CN108363959A (zh) | 一种针对乒乓球或羽毛球运动动作识别方法 | |
CN110084293A (zh) | 一种全明格局房屋的确定方法和装置 | |
CN107016371A (zh) | 基于改进的深度置信网络的无人机着陆地貌分类方法 | |
CN108596155A (zh) | 一种sar雷达空中飞行目标识别系统 | |
CN110111426A (zh) | 一种动静分离格局房屋的确定方法和装置 | |
CN108629382A (zh) | 一种群智优化sar雷达海上舰船目标识别系统 | |
Arif et al. | Human pose estimation and object interaction for sports behaviour |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180928 |
|
RJ01 | Rejection of invention patent application after publication |