买专利卖专利找龙图腾，真高效！ 查专利查商标用IPTOP,全免费！专利年费监控用IP管家,真方便！

申请/专利权人：东南大学

摘要：本发明公开了一种基于ELM与DSmT的飞机姿态识别方法。将飞机的360度飞行姿态分为7大类。对于待识别的飞机图像，经过去噪、灰度化、二值化处理后，提取具备平移、旋转、缩放不变性的SIFT特征与轮廓局部奇异值特征，可以有效解决单一特征提供的目标信息有限而造成算法识别率差这一问题；然后，采用学习速度快的ELM算法训练泛化能力强的SLFNs，将训练好的SLFNs作为算法的分类器；最后，采用DSmT融合理论中的PCR5公式，实现多特征的决策级融合，给出最终的识别结果。基于ELM与DSmT的飞机姿态识别方法可以解决飞机姿态识别中由于飞机姿态种类繁多、外界干扰等造成的识别率低下这一问题，同时在实时性方面，本方法也有一定的优越性。

主权项：一种基于ELM与DSmT的飞机姿态识别方法，其特征在于，具体步骤如下：第一步：对于飞机图像进行去噪、灰度化和二值化预处理；第二步：对经过预处理的飞机图像提取二值化图像的SIFT特征描述子和飞机轮廓特征；第三步：对SIFT特征描述子采用经典的BoW模型进行描述并最终得到BoWSIFT向量；第四步：对第二步中得出的飞机轮廓特征进行处理得到飞机目标轮廓局部奇异值向量；第五步：根据BoWSIFT向量，构建SLFNs网络，采用实时性好的ELM算法训练SLFNs网络，得到训练好的SLFNs1；对于轮廓局部奇异值向量，构建SLFNs，然后采用实时性好的ELM算法训练SLFNs网络，得到训练好的SLFNs2；第六步：将两个训练好的SLFNs1和SLFNs1作为算法的分类器，采用DSmT融合理论中的PCR5公式，实现多特征的决策级融合，得到最终的识别结果。

全文数据：一种基于ELM与DSmT的飞机姿态识别方法技术领域[0001]本发明属于自动目标识别技术领域，特别涉及一种基于ELM与DSmT的飞机姿态识别方法。背景技术[0002]自动目标识别技术AutomaticTargetRecognition，ATR是取得战场控制信息权的关键因素之一。从二十世纪八十年代开始，国内外很多的专家、学者积极参与了ATR技术的研究，使得ATR理论方面得到了完善，应用领域也得到很大的扩展。飞机目标识别作为ATR的重要领域之一，在现代战争和未来战争中都将占有着极为重要的地位。现代军事活动中，飞机的姿态识别越来越发挥着重要的作用，通过定位飞机的姿态，能够更方便准确的定位飞机的局部关键部位如螺旋浆、机舱、油箱）。[0003]早期的研究主要集中于利用单特征结合模式分类方法对飞机目标做分类和识别。刘刚等刘刚，梁晓庚，张京国.基于红外图像的飞机目标关键攻击部位识别[J].计算机工程与应用，2011,4724:174-178基于红外图像，将飞机图像分割成背景、机身、尾焰三部分，通过计算一种组合矩，同时结合面积比特征来实现飞机识别。涂建平等（涂建平，彭应宁，庄志洪.弹道终端飞机目标红外图像瞄准点识别方法[J].光学技术，2003，293:sei¬ses通过对飞机目标进行骨架提取、霍夫变换等过程，检测飞机机轴，识别飞机。[0004]基于单一特征的飞机目标识别在面对飞机种类繁多时，其算法性能很差，同时基于单一特征的算法在面对图像存在噪声如遮挡、缩放等情况下，误识别率很低。为了解决这个问题，一部分学者提出了基于多特征融合的方法来提高飞机目标识别算法的鲁棒性。曾接贤等（曾接贤，季康.基于多特征融合的多视点飞机目标识别算法[J].南昌航空大学学报：自然科学版，2016，302:9-15提出一种基于DSmT理论和SVM分类器相结合的多特征融合额多视点飞机目标识别算法，该方法首先对处理后的图像提取仿射不变矩与傅里叶描述子，利用SVM分类器对这两个证据源构建信度赋值，最后利用DSmT在决策级进行融合，从而完成对飞机目标的识别。[0005]基于决策级的多特征融合的目标识别算法流程为：[0006]1图像特征提取与表示；[0007]2分类器的构建与训练；[0008]3多特征融合；[0009]多特征融合识别算法在特征方面弥补了单一特征在识别方面的不足，但是面对飞机在飞行过程中姿态多变、外界环境复杂等情况，会给识别结果带来较大的影响，导致识别错误。同时在实际军事应用中，算法的实时性是很关键的指标之一。综上，如何选择鲁棒性好的多个特征，同时降低算法的时间复杂度是很有意义的工作。发明内容[0010]发明目的:针对现有技术中存在的由于飞机种类繁多、姿态多变导致识别率较低，同时算法的实时性普遍很差等问题，本发明提供一种将SIFT算法引入飞机姿态识别方法中，提取飞机图像的SIFT特征与轮廓局部奇异值特征;将ELM算法训练后的SLFNs作为分类器，可以有效的提高算法的实时性;最终在决策级，采用DSmT融合理论中的PCR5公式实现多特征融合，可以在有效的提高算法的鲁棒性的基于ELM与DSmT的飞机姿态识别方法。[0011]技术方案:为解决上述技术问题，本发明提供一种基于ELM与DSmT的飞机姿态识别方法，具体步骤如下：[0012]第一步:对于飞机图像进行去噪、灰度化和二值化预处理；[0013]第二步:对经过预处理的飞机图像提取二值化图像的SIFT特征描述子和飞机轮廓特征；[00M]第三步:对SIFT特征描述子采用经典的BoW模型进行描述并最终得到BoWsift向量；[0015]第四步:对第二步中得出的飞机轮廓特征进行处理得到飞机目标轮廓局部奇异值向量；[0016]第五步:根据BoWsift向量，构建SLFNs网络，采用实时性好的ELM算法训练SLFNs网络，得到训练好的SLFNs1;对于轮廓局部奇异值向量，构建SLFNs，然后采用实时性好的ELM算法训练SLFNs网络，得到训练好的SLFNs2;[0017]第六步：将两个训练好的SLFNsl和SLFNsl作为算法的分类器，采用DSmT融合理论中的PCR5公式，实现多特征的决策级融合，得到最终的识别结果。[0018]进一步的，所述第二步中对SIFT特征进行提取的具体步骤如下：[0019]2.1:SIFT特征点检测[0020]1，7表示一幅图像，〇是尺度空间因子，则其尺度空间1^，7,〇定义为：[0021]Lx，y，〇=Gx，y，〇*Ix，y[0022]其中Gx，y，〇为高斯核函数，二维高斯核函数定义为：[0023][0024]Lowe用计算简单、运算速度快的DoG算子替代LoG算子，DoG算子计算公式为：[0025]Dx,y,〇=Lx,y,k〇-Lx,y,〇[0026]其中，k为相邻两个高斯尺度空间的比例因子。[0027]通过DoG算子构建图像差分金字塔后，对每个采样点与同尺度的8个邻点以及上下两个尺度的9个邻点共计26个点进行比较，如果该采样点是最大值或者最小值，则认为其是图像在该尺度下的一个特征点，然后对筛选出来的特征点中的不稳定的边界点进行剔除；[0028]2.2:SIFT特征点描述：[0029]关键点检测完后，利用关键点邻域像素的梯度方向分布特性为每个关键点确定主方向和辅助方向，关键点x，y处的梯度幅值mx，y和方向Θx，y计算公式为：[0032]利用直方图统计关键点邻域内方向分布，将0〜360度的角度范围分为36个柱，每柱10度，直方图的峰值方向为关键点的主方向，将大于主方向峰值80%的方向作为该关键点处的辅助方向；[0033]将每个关键点的邻域点分配到4X4的子区域中，计算影响子区域的采样点的梯度方向和幅值，将每个子区域上计算得到的梯度分配到8个方向上，利用直方图统计每个方向的邻域点对关键点的贡献;对于每个关键点，形成了一个4X4X8=128维的特征向量，该向量为每个特征点的特征描述子，最后对特征描述子进行归一化处理。[0034]进一步的，所述第二步中提取飞机轮廓特征的具体步骤如下：首先提取出飞机图像的轮廓，同时计算图像的质心，然后选定固定起始点，从起始点将整条轮廓进行η等分，采样点记为匕^2，"_411]，计算每个采样点到质心的距离，并进行归一化，得到归一化距离向量[di,cb，…，dn]O[0035]进一步的，所述第三步中对SIFT特征采用经典的BoW模型进行描述并最终得到BoWsift向量的具体步骤如下：[0036]3.1:提取训练样本库中所有图像的SIFT特征描述子；[0037]3.2:使用KMeans聚类算法，将所有SIFT特征描述子聚为K类，生成K个视觉单词；[0038]3.3:计算样本中每个SIFT特征描述子与K个视觉单词的欧氏距离，通过欧式距离最小的原则统计所有SIFT特征描述子，然后得至IjK维向量，即为样本的BoWsift向量。[0039]进一步的，所述第四步中对轮廓特征进行处理得到轮廓局部奇异值的具体步骤如下:对于得到的归一化距离向量[cU，^，···，^]，借鉴数字信号处理中的采样思想，利用一个长度为w的采样窗口对归一化距离向量进行采样，两个采样窗口的重叠部分的长度为m，可以得1个采样后维度为w的距离向量，将这1个w维的归一化距离向量构造出1个wXw的循环矩阵，分别计算奇异值，最终按顺序可以组合为一个IXw的轮廓局部奇异值向量，完成图像轮廓局部奇异值的表不。[0040]与现有技术相比，本发明的优点在于：[0041]本发明将对噪声有不错稳定性的SIFT特征引入飞机姿态识别方法中，可以提高算法的识别准确率，同时在发生外界干扰如光照、遮挡等情况下，算法也具有一定的准确率；由于ELM算法相比于基于梯度法的训练算法有收敛速度快、不会陷入局部最优解等优点，将ELM算法训练的SLFNs作为算法的分类器，可以在提升算法识别率的同时，大大降低算法的训练时间开销。附图说明[0042]图1为本发明的总体流程图；[0043]图2为具体实施例中飞机7个姿态说明图；[0044]图3为具体实施例中SIFT特征128维向量生成过程图；[0045]图4为具体实施例中飞机目标轮廓离散采样图；[0046]图5为具体实施例中物体的BoW模型描述过程图；[0047]图6为具体实施例中基于BoW模型的物体识别过程图；[0048]图7为具体实施例中SLFNs网络结构图；[0049]图8为具体实施例中经典DSmT组合规则在yfΘ的描述图。具体实施方式[0050]下面结合附图和具体实施方式，进一步阐明本发明。[0051]如图1所示一种基于ELM与DSmT的飞机姿态识别方法，首先，对于待识别的飞机姿态图像，提取具备平移、旋转、缩放不变性的SIFT特征与轮廓局部奇异值特征，可以有效解决单一特征提供的目标信息有限而造成算法鲁棒性差这一问题;然后，采用学习速度很快的ELM算法训练泛化能力强的SLFNs，将训练好的SLFNs作为算法的分类器;最终，采用DSmT融合理论中的PCR5公式，实现多特征的决策级融合。下面对本发明加以具体说明。[0052]1、图像特征提取[0053]如图2所示，将飞机飞行的360度姿态分为具有代表性的7类。飞机在飞行过程中姿态多变，同时由于外界环境的影响，容易造成遮挡、缩放的发生，这就要求提取的特征在飞机图像发生平移、旋转、缩放时，具备不变性。对于特定物体识别，局部特征相比于全局特征有更好的鲁棒性，因为全局特征易受外界噪声如光照、遮挡等因素的影响。本发明使用了在目标跟踪与识别方法中应用广泛的SIFT特征，同时提取了轮廓局部奇异值特征，以这两种特征来描述飞机目标。[0054]I.ISIFT特征提取[0055]应用于飞机姿态识别方法中的SIFT特征提取可以分为两个步骤[0056]1.1.1SIFT特征点检测[0057]设Ix，y表示一幅图像，〇是尺度空间因子，则其尺度空间Lx，y，〇定义为：[0058]Lx，y，〇=Gx，y，〇*Ix，y[0059]其中Gx，y，〇为高斯核函数，二维高斯核函数定义为：[0060][0061]Lowe用计算简单、运算速度快的DoG算子替代LoG算子，DoG算子计算公式为：[0062]Dx,y,〇=Lx,y,k〇-Lx,y,〇[0063]其中，k为相邻两个高斯尺度空间的比例因子。[0064]通过DoG算子构建图像差分金字塔后，通过比较每个采样点与同尺度的8个邻点以及上下两个尺度的9个邻点共计26个点进行比较，如果此采样点是最大值或者最小值，则认为其是图像在该尺度下的一个特征点，之后还需要对筛选出来的特征点进行剔除一些不稳定的边界点。[0065]1·1·2SIFT特征点描述[0066]关键点检测完后，为了让算子具备旋转不变形，需要为每个特征点计算一个方向，然后利用关键点邻域像素的梯度方向分布特性为每个关键点确定主方向和辅助方向。关键点X，y处的梯度幅值mX，y和方向θX，y计算公式为：[0069]利用直方图统计关键点邻域内方向分布，可以将0〜360度的角度范围分为36个柱，每柱10度，直方图的峰值方向为关键点的主方向，为增加鲁棒性，将大于主方向峰值80%的方向作为该关键点处的辅助方向。[0070]接下来就是通过一组向量将每个关键点描述出来，将每个关键点的邻域点分配到4X4的子区域中，计算影响子区域的采样点的梯度方向和幅值，将每个子区域上计算得到的梯度分配到8个方向上，利用直方图统计每个方向的邻域点对关键点的贡献。对于每个关键点，形成了一个4X4X8=128维的特征向量，如图3所示。最后还需要对向量进行归一化处理，以降低光照对特征描述子的影响。[0071]1.2轮廓奇异值[0072]对于飞机的二值化图像，首先提取出轮廓，同时计算图像的质心，选定固定起始点，从起始点将整条轮廓进行η等分，采样点记为[pi，P2，…，pn]，计算每个采样点到质心的距离，并进行归一化，得到归一化距离向量[cU，d2,···，!„]，飞机轮廓的离散采样过程如图4所示。[0073]2图像特征描述[0074]2.ISIFT特征表示[0075]每张图像经过SIFT特征点提取与表示后，均由多个IX128的SIFT特征描述子表示，不同图像的特征点个数可能不同，难以直接采用分类器进行分类。可以采用经典的BoW模型进行描述，可以有效的解决这个问题。基于BoW模型的物体描述如图5所示，具体步骤如下[0076]2.1.1:首先提取训练样本库中所有图像的SIFT特征描述子；[0077]2.1.2:接着使用KMeans聚类算法，将所有SIFT特征描述子聚为K类，生成K个视觉单词；[0078]2.1.3:最后计算样本中每个SIFT特征描述子与K个视觉单词的欧氏距离，通过欧式距离最小的原则统计所有SIFT特征描述子，如此可以得到K维向量，即为样本的BoWsift向量。[0079]基于BoW的物体识别过程如图6所示。[0080]2.2轮廓局部奇异值[0081]对于得到的归一化距离向量[cU，^，···，^]，借鉴数字信号处理中的采样思想，利用一个长度为w的采样窗口对归一化距离向量进行采样，两个采样窗口的重叠部分的长度为m，可以得1个采样后维度为w的距离向量，将这1个w维的归一化距离向量构造出1个wXw的循环矩阵，分别计算奇异值，最终按顺序可以组合为一个IXw的轮廓局部奇异值向量，完成图像轮廓局部奇异值的表示。[0082]3分类器训练环节[0083]单隐层前馈神经网络（SingleHidden-layerFeedforwardNeuralNetworks具有较强的学习能力，能够逼近复杂非线性函数，能够解决传统方法无法解决的问题等优点，其网络结构如图7所示。训练SLFNs通常采用的算法是基于梯度下降法，这类学习算法普遍具有收敛速度慢、易陷入局部最优点等缺点，黄广斌等（HuangGB，ZhuQY，SiewCK.ExtremeLearningMachine:ANewLearningSchemeofFeedforwardNeuralNetworks[C]·In:ProceedingsoftheIJCNN·2004，2:985-990提出极限学习机算法，可以改善单隐层神经网络的学习算法以及泛化过程存在的问题。[0084]设HX为隐层输出矩阵，β为输出权重，T为给定的监督矩阵，SLFNs—般模型可以表示为：[0085]Ηβ=Τ[0086]其中[0089]极限学习机算法可以总结为如下流程：[0090]给定训练集山er，激活函数g，隐层节点个数为^[0091]1随机生成输入权重Wi和隐层偏置.¾¾…，爲；[0092]2计算隐层输出矩阵H;[0093]3计算输出权重β:[0094]β=Η+Τ[0095]H+是隐藏层输出矩阵H的Moore-Penrose增广逆，简称为伪逆。[0096]当β计算完成后，一个单隐层反馈神经网络的训练环节就结束了。对于一个测试样本X，它的标签可以通过训练好的单隐层反馈神经网络推测：[0097]fLx=hx0[0098]其中hX是神经网络隐层关于此样本的响应，初始化参数后，hX即保持不变。[0099]4多特征融合识别环节[0100]实际中进行飞机目标识别时，由于飞机种类繁多，可能会出现同类飞机差距很大，不同类飞机在某个姿态下很相似等情况，飞机目标识别在本质上具有高冲突和模糊性的特点，DSmT证据融合理论针对这一特性，为基于决策级融合的飞机目标识别提供了一个非常好的思路。经典DSmT组合规则在VΘ的描述如图8所示。本文针对基于决策级融合的飞机目标识别，在实验过程中做了一定的简化，在超幂集D0的基础上，采用完全排他性的约束条件下的DSmT模型，即仅Q1，θ2，…θη具有信度赋值。本文中只有轮廓局部奇异值和SIFT两种特征，采用PCR5组合规实现两特征对应的决策级融合。[0101]以上所述仅为本发明的实施例子而已，并不用于限制本发明。凡在本发明的原则之内，所作的等同替换，均应包含在本发明的保护范围之内。本发明未作详细阐述的内容属于本专业领域技术人员公知的已有技术。

权利要求：1.一种基于ELM与DSmT的飞机姿态识别方法，其特征在于，具体步骤如下：第一步:对于飞机图像进行去噪、灰度化和二值化预处理；第二步：对经过预处理的飞机图像提取二值化图像的SIFT特征描述子和飞机轮廓特征；第三步:对SIFT特征描述子采用经典的BoW模型进行描述并最终得到BoWsift向量；第四步：对第二步中得出的飞机轮廓特征进行处理得到飞机目标轮廓局部奇异值向量；第五步:根据BoWsift向量，构建SLFNs网络，采用实时性好的ELM算法训练SLFNs网络，得到训练好的SLFNsl;对于轮廓局部奇异值向量，构建SLFNs，然后采用实时性好的ELM算法训练SLFNs网络，得到训练好的SLFNs2;第六步：将两个训练好的SLFNsl和SLFNsl作为算法的分类器，采用DSmT融合理论中的PCR5公式，实现多特征的决策级融合，得到最终的识别结果。2.所述第二步中对SIFT特征进行提取的具体步骤如下：2.1:SIFT特征点检测1，7表示一幅图像，〇是尺度空间因子，则其尺度空间1^，7,〇定义为：Lx，y，〇=Gx，y，〇*Ix，y其中Gx，y，σ为高斯核函数，二维高斯核函数定义为：Gx,y\a=-^-e-{x2+y2l2cx22πσιLowe用计算简单、运算速度快的DoG算子替代LoG算子，DoG算子计算公式为：Dx,y,〇=Lx,y,k〇-Lx,y,〇其中，k为相邻两个高斯尺度空间的比例因子。通过DoG算子构建图像差分金字塔后，对每个采样点与同尺度的8个邻点以及上下两个尺度的9个邻点共计26个点进行比较，如果该采样点是最大值或者最小值，则认为其是图像在该尺度下的一个特征点，然后对筛选出来的特征点中的不稳定的边界点进行剔除；2.2:SIFT特征点描述：关键点检测完后，利用关键点邻域像素的梯度方向分布特性为每个关键点确定主方向和辅助方向，关键点X，y处的梯度幅值mX，y和方向θX，y计算公式为：m{x,y—+1,ν·~L{x—I,y+L{x,v+1—L{x,v—I^x,v=tan丨产-+DU-1'Lx+l,.v-LX-Ky利用直方图统计关键点邻域内方向分布，将0〜360度的角度范围分为36个柱，每柱10度，直方图的峰值方向为关键点的主方向，将大于主方向峰值80%的方向作为该关键点处的辅助方向；将每个关键点的邻域点分配到4X4的子区域中，计算影响子区域的采样点的梯度方向和幅值，将每个子区域上计算得到的梯度分配到8个方向上，利用直方图统计每个方向的邻域点对关键点的贡献;对于每个关键点，形成了一个4X4X8=128维的特征向量，该向量为每个特征点的特征描述子，最后对特征描述子进行归一化处理。3.所述第二步中提取飞机轮廓特征的具体步骤如下：首先提取出飞机图像的轮廓，同时计算图像的质心，然后选定固定起始点，从起始点将整条轮廓进行η等分，采样点记为匕1^2广_411]，计算每个采样点到质心的距离，并进行归一化，得到归一化距离向量[11，d2,***,dn]〇4.所述第三步中对SIFT特征采用经典的BoW模型进行描述并最终得到BoWsift向量的具体步骤如下：3.1:提取训练样本库中所有图像的SIFT特征描述子；3.2:使用KMeans聚类算法，将所有SIFT特征描述子聚为K类，生成K个视觉单词；3.3:计算样本中每个SIFT特征描述子与K个视觉单词的欧氏距离，通过欧式距离最小的原则统计所有SIFT特征描述子，然后得至IjK维向量，即为样本的BoWsift向量。5.所述第四步中对轮廓特征进行处理得到轮廓局部奇异值的具体步骤如下:对于得到的归一化距离向量[cU，^，···，1η]，借鉴数字信号处理中的采样思想，利用一个长度为w的采样窗口对归一化距离向量进行采样，两个采样窗口的重叠部分的长度为m，可以得1个采样后维度为w的距离向量，将这1个w维的归一化距离向量构造出1个wXw的循环矩阵，分别计算奇异值，最终按顺序可以组合为一个IXw的轮廓局部奇异值向量，完成图像轮廓局部奇异值的表不。

百度查询： 东南大学 一种基于ELM与DSmT的飞机姿态识别方法