一种基于isar图像的舰船目标模糊识别方法

文档序号:10512764阅读:460来源:国知局
一种基于isar图像的舰船目标模糊识别方法
【专利摘要】一种基于ISAR图像的舰船目标模糊识别方法,本发明涉及基于ISAR图像的舰船目标模糊识别方法。本发明的目的是为了解决现有技术采用均匀分段,对上层曲线的分段过于简单,同一结构内的曲线被分到不同的段中,造成段内均值和编码的误差,并且结构编码对上层结构曲线的描述也不够充分,造成识别准确率不高的问题。具体过程为:一、提取上层结构曲线;二、计算分段后的分段相对均值特征;三、使用训练样本估计构造目标类别模板;四、计算并得到待识别目标和目标类别模板的模糊集合形式及隶属度函数;五、计算待识别目标与不同类别模板的模糊集的贴近度,将待识别目标判定为贴近度最大的类别。本发明应用于舰船目标识别领域。
【专利说明】
一种基于ISAR图像的舰船目标模糊识别方法
技术领域
[0001] 本发明涉及基于ISAR图像的舰船目标模糊识别方法。
【背景技术】
[0002] ISAR是利用静止的雷达利用运动目标产生的相对运动形成的合成孔径来成像的 一种雷达成像技术,雷达固定在地面或者装载于舰船上(分别对应岸基ISAR与舰载ISAR), 目标多为非合作目标。目标的形态特征在高分辨ISAR图像中得到了较好的保留,因此基于 ISAR图像的舰船目标识别技术在军事与民用范围内都具有广阔的前景。
[0003] 现有的基于ISAR图像的舰船目标研究,主要使用目标的形态与轮廓特征,对于舰 船目标,尤其是长度体型较为近似的目标,最明显的差异体现在舰船上层结构一一舰船甲 板以上的部分。有人提出方法,在提取舰船目标中心线后,得到上层结构曲线,并提取桅杆 和上层结构编码等特征,使用模板匹配方法对目标进行识别。然而,该方法使用均匀分段, 对上层曲线的分段过于简单,同一结构内的曲线被分到不同的段中,造成段内均值和编码 的误差。并且结构编码对上层结构曲线的描述也不够充分,造成识别准确率不高。

【发明内容】

[0004] 本发明的目的是为了解决现有技术采用均匀分段,对上层曲线的分段过于简单, 同一结构内的曲线被分到不同的段中,造成段内均值和编码的误差,并且结构编码对上层 结构曲线的描述也不够充分,造成识别准确率不高的问题,而提出一种基于ISAR图像的舰 船目标模糊识别方法。
[0005] 上述的发明目的是通过以下技术方案实现的:
[0006] 步骤一、采用Radon变换,对目标ISAR图像提取Centre line,并提取上层结构曲线;
[0007] 步骤二、利用均值梯度法对提取的上层结构曲线进行自适应分段,并计算分段后 的分段相对均值特征;
[0008] 步骤三、使用训练样本估计构造目标类别模板;
[0009] 步骤四、对待识别目标和目标类别模板,采用隶属度原则,计算并得到待识别目标 和目标类别模板的模糊集合形式及隶属度函数;
[0010] 步骤五、计算待识别目标与不同类别模板的模糊集的贴近度,贴近度计算采用海 明贴近度,欧几里得贴近度或格贴近度,应用最大贴近度原则,将待识别目标判定为贴近度 最大的类别。
[0011] 发明效果
[0012] 采用自适应分段技术,对舰船上层结构进行自动化地分段,不同目标之间根据结 构不同在不同的位置进行分段,使得分段结果更加合理,接近人眼的直观结果,并且保证各 段内都包含完整的上层结构,有利于特征的提取,解决了均匀分段法造成段内均值和编码 的误差的问题。可以实现上层结构曲线分段的自动化;同时,自适应分段所采用的原则也使 得段内均值更加平稳,有利于特征的提取。此外,将原有的舰船编码变换为分段相对均值特 征,将特征数值化,一方面可以与模糊集理论结合,采用模糊识别算法,提升上层结构曲线 的识别效果,一方面也增大了特征空间,有利于许多分类器的应用。可以应用更多的识别方 法和学习手段。采用模糊集合,取代原有的单一舰船编码,表示为多种舰船编码全集上的模 糊集合,极大地提升了原有舰船编码的表达能力,并采用最大隶属度原则,进行模糊识别, 可以有效地提高现有基于ISAR图像上层结构的舰船目标识别技术的正确率。并改进了模板 匹配方法,采用最大贴进度原则,有效提升了现有识别算法的准确率。解决了现有结构编码 对上层结构曲线的描述也不够充分,造成识别准确率不高的问题。
[0013] 本实验中采用的自适应分段算法一一均值梯度法,相比于原始的均匀分段,有效 地提升了识别的正确率。随后比较本文算法与特征匹配法。很明显,使用单一编码来描述目 标的特征匹配法,其平均识别正确率只有0.592,远远低于本文算法的0.781。特别的,目标C 的识别率只有0.433,在实际应用中这种水平识别率很可能无法应用。因此,本文的模糊识 别算法大大提升了基于舰船上层结构曲线编码的识别效果。最后比较本文算法与不变矩 法。不变矩法的平均正确率稍低于本文算法,然而,不变矩法对目标A的识别正确率只有 0.423,远低于其对目标B和目标C的识别率。同样地,这种水平的识别率在实际应用中很可 能无法使用,使得该算法在多目标情况下的稳定性存在缺陷。因此,本文提出的模糊识别算 法,相比于现有的其他算法,不仅具有最高的识别正确率,同时也具有最好的识别稳定性。 实验结果证明了该算法的有效性,以及算法对现有技术的提升。
【附图说明】
[0014] 图1是本发明的识别流程图;
[0015] 图2为Radon变换结果示意图;
[0016]图3为舰船Centreline提取结果示意图;
[0017]图4a为舰船上层结构曲线a简单分段结果示意图;
[0018]图4b为舰船上层结构曲线b简单分段结果示意图;
[0019]图5a为均值梯度法曲线a分段结果示意图;
[0020]图5b为均值梯度法曲线b分段结果示意图;
[0021 ]图6a为曲线a分段相对均值特征值示意图;
[0022]图6b为曲线b分段相对均值特征值示意图;
[0023]图7为隶属度函数分布形式图;
[0024]图8a为实验渔船数据库示例图;
[0025]图8b为实验货船数据库示例图;
[0026]图8c为实验军舰船数据库示例图;
[0027] 图9a为渔船示例图的上层结构曲线与自适应分段结果示意图;
[0028] 图9b为货船示例图的上层结构曲线与自适应分段结果示意图;
[0029] 图9c为军舰船示例图的上层结构曲线与自适应分段结果示意图;
[0030]图10a为渔船示例图的上层结构曲线与均匀分段结果示意图;
[0031] 图10b为货船示例图的上层结构曲线与均匀分段结果示意图;
[0032] 图10c为军舰船示例图的上层结构曲线与均匀分段结果示意图;
[0033]图11为三种贴近度的模糊识别正确率示意图,A为渔船,B为货船,C为军舰;
[0034] 图12为四种用于比较的识别技术的正确率示意图,A为渔船,B为货船,C为军舰。
【具体实施方式】
【具体实施方式】 [0035] 一:结合图1说明本实施方式,本实施方式的一种基于ISAR图像的舰 船目标模糊识别方法,具体是按照以下步骤制备的:
[0036]步骤一、采用Radon变换(拉登变换),对目标ISAR图像(逆合成孔径雷达)提取 Centreline(中心线),并提取上层结构曲线;
[0037]步骤二、利用均值梯度法对提取的上层结构曲线进行自适应分段,并计算分段后 的分段相对均值特征;
[0038]步骤三、使用训练样本估计构造目标类别模板;
[0039] 训练样本和目标类别的已知情报信息为用户设置的;
[0040] 步骤四、对待识别目标和目标类别模板采用隶属度原则,计算并得到待识别目标 和目标类别模板的模糊集合形式及隶属度函数;
[0041] 步骤五、计算待识别目标与不同目标类别模板的模糊集的贴近度,其中贴近度的 计算采用三种不同的定义,分别为海明贴近度,欧几里得贴近度和格贴近度,应用最大贴近 度原则,将待识别目标判定为贴近度最大的目标类别。
【具体实施方式】 [0042] 二:本实施方式与一不同的是:所述步骤一中采用 Radon变换(拉登变换),对目标ISAR图像(逆合成孔径雷达)提取Centre 1 ine (中心线),并提 取上层结构曲线;具体过程为:
[0043] 步骤--、提取 Centre line;
[0044] 对目标ISAR图像进行全角度、全截距的Radon变换,得到变换域,在变换域内通过 峰值搜索,找到Centreline的倾角和投影点;图2展示了这样的过程。其中Θ轴表示旋转角 度,投影方向单元数表示投影点距离y'的像素数;提取并绘制出Centreline;图3展示了对 图2中实测数据提取绘制Centrel ine的结果;
[0045] 这种方法利用了舰船长度一般远大于高度的特点,因为上层结构的存在,可以实 现对Centre line的无偏估计。
[0046] Radon变换是将图像矩阵对指定方向进行投影,它可以很好地检测出二维平面图 像沿某条直线的能量分布。
[0047]对于一幅二维平面图像f(x,y),其Radon变换为
[0048] - I fix' cos θ(?) - y' sin 0(t), x' sin 9(t·) +y' cos 9(t)] dy' (1)
[0049] 该式表示坐标轴V、/绕x、y旋转角度0(t)后,沿/轴方向进行线积分的结果;
[0050] 步骤一二、在步骤--的基础上提取上层结构曲线;
[0051] 提取出Centreline后,目标ISAR图像被Centreline分为两个部分,如图3所示,分 别为包含桅杆或船舱的上层结构部分和包含甲板与船底的船身部分;首先应当区分出这两 个部分。沿着Centreline,对目标ISAR图像(目标ISAR图像包括目标和背景)所在的距离单 元内,提取每个距离单元内Centrel ine两侧目标ISAR图像与Centrel ine之间的最大偏移 量,可以得到两条曲线,分别为船身曲线和上层结构曲线;根据舰船目标的结构特性,我们 认为,由于甲板普遍较低,上层结构普遍较高,偏移量大的曲线为上层结构曲线,并由此区 分两个部分。
[0052]其它步骤及参数与【具体实施方式】一相同。
[0053]【具体实施方式】三:本实施方式与【具体实施方式】一或二不同的是:所述步骤二中利 用均值梯度法对提取的上层结构曲线进行自适应分段,并计算自适应分段后的分段相对均 值特征;具体过程为:
[0054]二、上层结构曲线的自适应分段与改进编码特征 [0055] 1、自适应分段的意义
[0056] 上层结构曲线自适应分段是对传统上层结构编码特征的改进。传统编码特征是将 舰船的上层结构曲线分为三段,并利用每一段内曲线的均值对各段进行二进制编码,从而 获取固定长度的特征,用以表示舰船的上层结构。相比于直接使用上层结构曲线作为特征, 舰船编码具有特征维数固定的优势,避免了不同视角和分辨率下舰船目标上层结构曲线长 度不同需要进行的归一化或插值处理,也为识别部分的分类器选择提供了更大的灵活性。 此外,曲线均值的平滑作用也可以降低噪声干扰和视角变化的影响,增大特征的稳健性。
[0057] 根据编码特征的原理,其特征值依赖于曲线的分段结果。现有的舰船编码技术采 用简单的均匀分段,每段的长度相同。这在很多情况下是不合理的,比如目标船只具有很大 很长的船舱,占据了上层结构的大部分长度,相同长度的分段会把完整的船舱曲线分割到 不同的段中,不仅破坏了上层结构的整体性,同时也使得段内曲线高度变化,进而影响到编 码结果,降低了识别的准确性。如图4所示,为两种舰船目标上层结构直接平均分段的结果。 图4a中第二分段点略微偏左,使得分段点右侧一段较低的平坦区域与明显较高的船舱结构 分在一段中,可能会影响第三段内的平均曲线高度。而图4b中两个分段点都位于曲线的峰 值点附近,由于上层结构曲线峰值往往对应着桅杆和船舱顶部等特殊结构,这样的分段结 果会严重影响上层结构的完整性和特征的提取。
[0058] 因此我们将从直观与数学的角度提出对舰船上层结构进行自适应分段的算法,来 提高该特征的合理性和稳健性。
[0059] 步骤二一、采用均值梯度法对上层结构曲线进行自适应分段;
[0060]均值梯度法是我们提出的自适应分段算法。该算法的思想是在不同位置对上层结 构曲线进行分段,根据段内曲线均值的变化来确定最佳分段点。直观上,曲线的谷点一般是 分段比较合理的位置。当曲线长度变化时,若增加的点的值比较小,曲线的均值会变小,当 增加的点是谷点或明显与原曲线各点的值相差很大时,曲线均值就会急剧下降。因此,合适 的分段点处段内曲线均值的一阶梯度会是一个极小值。我们就在这样的点处对舰船上层曲 线进行分段,该算法适合从曲线两端分别进行。
[0061 ]具体算法如下:
[0062] 步骤二一一、从提取的上层结构曲线一端C开始,以1/3曲线长度点为中心,取总曲 线长度的2/9或1/3为范围,计算范围内所有的点到端点C之间的上层结构曲线的均值,绘制 出均值曲线;
[0063] 步骤二一二、对步骤二一一绘制出的均值曲线求一阶梯度,步长为1,找到一阶梯 度的最小值;
[0064] 步骤二一三、以一阶梯度的最小值对应的位置作为端点C的分段点位置;
[0065] 步骤二一四、从提取的上层结构曲线另一端点D开始,以1/3曲线长度点为中心,取 总曲线长度的2/9或1/3为范围,计算范围内所有的点到端点D之间的上层结构曲线的均值, 绘制出均值曲线;对均值曲线求一阶梯度,步长为1,找到一阶梯度的最小值;以一阶梯度的 最小值对应的位置作为端点D的分段点位置;完成自适应分段;
[0066]对图4a、图4b中相同的上层结构曲线采用均值梯度法进行分段,分段的结果如图 5a、图5b所示。
[0067]从图5a、图5b可以看出,均值梯度法所得到的分段点的位置基本都在曲线的拐点 或谷点附近,自适应能力较强,分段结果比较符合人的直观分类,各段内上层建筑的完整性 得到保留。
[0068] 步骤二二、在步骤二一的基础上计算自适应分段后的分段相对均值特征;
[0069] 传统上层结构编码特征,分段完成后将采用二进制数字对各段进行编码,用1表示 较高的段内结构,〇表示较低的段内结构。这种标称型的特征过于粗糙,也会限制分类器的 选择。因此我们将其改进为分段相对均值(Segment Comparative Mean,SCM)特征。
[0070] 分段相对均值特征用数值来表示每一段,每一段的数值由段内曲线高度的均值与 总体曲线均值的差除以总体曲线均值得到。
[0071] 因此,该数值为正时,等效于编码1,该数值为负时则等效于编码0,数值的绝对值 大小表现了段内结构'高'与'低'的程度。由于该数值是相对量,ISAR图像的横向分辨率在 计算中被抵消,特征对分辨率不变,因此该特征可以对不同分辨率的目标图像同时使用而 不需要进行归一化。同时,因为特征为数值型,更多的识别方法和学习手段可以被应用,从 而提升了该特征的应用范围和识别效果。对图5a、图5b的分段结果计算改进的特征值,得到 的结果如图6a、图6b所示。
[0072] 其它步骤及参数与【具体实施方式】一或二相同。
【具体实施方式】 [0073] 四:本实施方式与一至三之一不同的是:所述步骤四 中对待识别目标和目标类别模板采用隶属度原则,计算并得到待识别目标和目标类别模板 的模糊集合形式及隶属度函数;具体过程为:
[0074] 1、模糊集理论简介
[0075] 本小节将简要介绍本文所用模糊识别方法所需要的模糊集合基本概念。
[0076] 隶属度函数:隶属度函数表征一个元素 u隶属于一个集合A的程度,记作AUhAU) 的取值范围为[0,1 ],A(U) = 1表示元素 u完全属于集合A,A(u) =0则表示元素 u完全不属于 A,A(u)介于0和1之间时,表示元素 u以A(u)的程度属于集合A。
[0077] 模糊集:对论域U上所有的元素 u,定义其隶属度函数A(u),则A就是定义在论域U上 的模糊集。若论域U包含有限个元素 m,u2,···,un,则模糊集A可以表示为:
[0078] A = T^lA^ui)/ui (2)
[0079] 其中,符号不是传统的分号,只是表示论域上的元素 m与其隶属度A(m)的对应 关系;符号' Σ '也不是传统意义的求和,而是表示论域U上所有的元素与其隶属度关系对应 关系的总括。
[0080] -、待识别目标和目标类别模板的上层结构编码的模糊集;
[0081] 将待识别目标和目标类别模板的上层结构表示为论域上的模糊集,取代原有的单 个编码;模糊集的隶属度函数计算依赖于分段相对均值(SCM)特征,分段相对均值特征的提 取已在上一节中介绍过。此处将使用SCM指代该特征向量,用SCMi指代SCM向量的第i个值, 即第i段内的相对均值。
[0082]由于论域是八个三位二进制编码(000到111)的全集,三位二进制编码为000到 111;因此,对论域U上的模糊集A表示为:
[0085]其中,表示j对应数字的二进制编码,如u2表示010;A为论域U上的待识别目标模 糊集,B为论域U上的类别模板的模糊集;
[0086]例如,一个典型的被编码为010的舰船目标B可以用模糊集表示为:
[0088] 根据模糊集相关定义,这种表示方式可以反应不同编码隶属于该目标的程度,也 就是不同编码是该目标正确编码的可能性;因而该模糊集描述了目标被编码成这八种编码 的可能性。也就是说,该模糊集可以描述目标上层结构与各种编码形状之间的似然度,而不 是给出一个绝对的编码。这样可以提高编码的精确性。
[0089] 这里的关键在于如何得到各种编码对目标的隶属度。我们通过一个两步的方案来 获取待识别目标和目标类别模板的上层结构编码的模糊集的隶属度函数:一、得到各段内 编码0和编码1的隶属度;二、计算各个编码Uj的总隶属度;
[0090] 通过SCM特征计算各段内编码0和编码1的隶属度;
[0091 ] SCM为分段相对均值,SCMi为SCM向量的第i个值,即第i段内的相对均值;
[0092]用menuo和memu分别表示第i段的编码0和编码1的隶属度;这两个隶属度的函数形 式为:
[0095] 式中,X为 SCMi;
[0096] 两个隶属度函数的分布形式如图7所示。
[0097] 因此,对某一待识别目标和目标类别模板的上层结构编码,第i段的编码0隶属度 为memio(SCMi),第i段的编码1隶属度为menuKSCMi);例如,图6b中第3段SCM!的值为-0.266, 因此该段的编码0隶属度m em3Q为0 · 871,编码' Γ隶属度menm为0 · 129。这个结果是很合理 的,因为从图6b中可以看到,该段非常可能被编码为0。
[0098]得到各段内的menuo与memu后,计算编码叫的总隶属度;若叩是编码的第i段的码 元,贝1Juj的总隶属度为:
[0099] A(uj) ^ (7)
[0100] B(uj) = Wimemiu:i
[0101] 其中^为权值;权值的和应当等于1,以保证隶属度的范围为[0,1],从而满足隶属 度的定义。缺省情况下可以使用[1/3 1/3 1/3]的均匀权值,若对某一段的结果更加重视, 可以适当增加该段的权值大小,并相应减少其他段的权值。
[0102] 例如,图6a中目标的SCM为[0.445 -0.603 0.281 ],则该目标的编码0隶属度为 [0.0074 1 0.1137],编码1隶属度为[0.9926 0 0.8863],使用均匀权值时,该目标的编码 115 = 4017的隶属度为:
[0104]因此,该目标可以被表示为如下的模糊集:
[0106]该表示形式的合理性是显而易见的。由图6a可以看到,目标最可能被编码为 '101',对应编码的隶属度为0.960,为所有编码中隶属度最大的。而最不可能的编码'101' 的隶属度只有〇.04。同时,该模糊集还保留了目标对其他编码的可能性。
[0107] 其它步骤及参数与【具体实施方式】一至三之一相同。
【具体实施方式】 [0108] 五:本实施方式与一至四之一不同的是:所述步骤五 中计算待识别目标与不同类别模板的模糊集的贴近度,其中贴近度的计算采用三种不同的 定义,分别为海明贴近度,欧几里得贴近度和格贴近度,应用最大贴近度原则,将待识别目 标判定为贴近度最大的类别;具体过程为:
[0109] 最大贴近度原则
[0110] 模糊识别中的最大贴近度原则,就是将目标判定为与其最贴近的模糊集的识别方 法。贴近度是两个模糊集之间相似性的度量。两个模糊集A与B之间的贴近度记作P(A,B)。假 设B k,k = l,…,c是有限论域U上的类别模板的模糊集,A为有限论域U上的待识别目标模糊 集,若:
[0111] p(A,Bi) =maxk[p(A,Bk) ] (10)
[0112] 则将待识别目标模糊集A的种类判定为出所属的类,即第1类;BiSBk中与待识别目 标模糊集A贴近度最大的值,B为某个类别模板的模糊集,1为B的编号,c为1的最大值,1 < 1 <c。在本节叙述的识别方法中,A是目标模糊集,使用上面介绍的方法就可以得到。而Bk需 要从训练样本中估计或者预先构造。使用训练样本估计时,可以简单地计算所有训练样本 的平均SCM特征,以该特征作为模板,计算该类的模糊集。而实际应用中若没有训练样本,或 者有情报或光学信息,也可以根据目标形态模型计算或构造目标的SCM特征,从而获取目标 类的模糊集。
[0113] 该识别方法的关键在于贴近度的计算。有三种常见的贴近度定义:贴近度的计算 采用海明贴近度PHamming,欧几里得贴近度PEimlid或格贴近度Plattice ;对一个有限论域U = {U1, 1!2,"_,1!11}和其上的两个模糊集4,8,这三种贴近度的定义分别为 :
[0116] piattice(A,B) = (A〇B)A(Ac〇Bc) (13)
[0117] 其中,A(uJ与B(uJ分别是待识别目标模糊集A和类别模板的模糊集B的隶属度函 数;上式中的#和把是待识别目标模糊集A和类别模板的模糊集B的补集,符号θ和Λ分别表 示模糊集的内积和取最小运算;η取值为论域中元素的个数,本发明中为8;
[0118] 对一个待识别目标模糊集Α,其补集的隶属度函数为:
[0119] Ac(x) = l-A(x) (14)
[0120] 式中,X为 SCMi;
[0121] 对于两个实数a和b,aAb表示:
[0122] aAb=min{a,b} (15)
[0123] 式中,a为A〇B;b 为Ac〇Bc;
[0124] 模糊集A和B的内积为:
[0125] A0B = Vxeu[A(x)AB(x)] (16)
[0126] 其中符号V为取最大值运算:
[0127] Vx[A(x) ] =maxx[A(x) ] (17)
[0128] 这三种贴近度的定义在不同情形下具有不同的效果,在实际应用中应根据具体情 况选择最优的贴近度定义。
[0129] 其它步骤及参数与【具体实施方式】一至四之一相同。
[0130] 采用以下实施例验证本发明的有益效果:
[0131] 实施例一:
[0132] 下面通过对实测ISAR数据进行识别来验证本发明的有效性。
[0133] 1、数据库简介
[0134] 实验的ISAR数据库中包含三种具有不同上层结构的舰船,分别为:A.渔船,B.货船 和C.军舰。三种舰船的ISAR示例图片如图8a、8b、8c所示,而相应的上层结构曲线如图9a、 9b、9c所示。在数据库中,渔船的图片共有73张,纵向分辨率均为0.375m;货船的图片共有84 张,其中46张纵向分辨率为0.75m,剩余38张的纵向分辨率为0.375m;军舰图片共有227张, 分辨率均为0.75m。数据录入时雷达与渔船、货船和军舰的距离分别为5km,13km和约40km。
[0135] 2、实验方案
[0136] 根据本节的研究目的,我们设计了如下的实验方案。从数据库中对每类舰船随机 选取1/3的图片作为训练样本,剩余图片将作为测试样本,用来评价识别技术的效果。为了 降低随机抽取样本带来的样本偏差与偶然性,我们将重复1〇〇〇次随机抽取样本和识别测 试。图片预处理与上层结构曲线的提取采用之前介绍的方法。在自适应分段步骤中,使用均 值梯度法进行目标曲线的自适应分段,代表图片的分段结果已展示在图12中。由训练样本 估计出类别模板后,采用上述的识别算法对三种舰船进行识别,并计算识别正确率。实验 中,为了测试贴近度定义在本应用中的相对优劣,我们将分别采用三种贴近度进行实验,每 种贴近度应用时都使用完全相同的训练样本和测试样本,经过1000次实验后,得到平均的 识别正确率。
[0137] 同时,为了验证自适应分段算法的有效性,我们设计了一个参照组,该参照组使用 完全相同的实验过程和实验样本,只是在自适应分段步骤中,采用简单的均匀分段,并由此 计算相应的SCM特征和目标模糊集。图8a、8b、8c的代表图片采用简单均匀分段的结果如图 10a、10b、10c所示。经过1000次实验后,同样得到平均的识别正确率。
[0138] 最后,为了验证该模糊识别技术对现有技术的提升,我们还设置了两个参照组,分 别使用现有的两种基于舰船目标ISAR图像的上层结构提出的识别技术。其中,技术[1] (D.Pastina,C.Spina.Multi-feature based automatic recognition of ship targets in ISAR[J] · IET radar,sonar&navigation.2009,3(4) :406-423.)使用特征匹配方法,将 目标上层结构曲线编码与数据库模板相比较,从而选出候选类别。而技术[2](D.Pastina, C.Spina.Multi-feature based automatic recognition of ship targets in ISAR[J] .IET radar,sonar&navi gat ion · 2009,3(4): 406-423 ·)计算目标上层结构的 Hu 不变矩特 征,并使用分类器进行识别。在这里我们使用kNN作为分类器。实验时,两种技术都使用与模 糊识别技术实验时完全相同的训练样本和测试样本,重复1000次后得到平均的识别正确 率。
[0139] 3、实验结果
[0140] 本小节将展示实验的结果。实验按照方案进行,首先得到了三种贴近度,即海明贴 近度,欧几里得贴近度和格贴近度定义下的识别结果,列于表1的混淆矩阵中。
[0141]三种贴近度定义的识别正确率如图11所示。从图11中可以清楚地看到,采用欧几 里得贴近度定义的模糊识别技术获得了最好的识别结果。从表1的结果中可以看出,如果目 标没有被正确识别,则目标A(渔船)更有可能被识别为目标C(军舰),目标B(货船)和目标C 则更有可能被识别为目标A。目标B和目标C在欧几里得贴近度定义下的正确率要高于海明 贴近度和格贴近度,同时目标A的正确率只是稍稍低于海明贴近度的结果。因此,欧几里得 贴近度的平均识别率在二种定义中是最尚的。
[0144] 由实验结果可以看出,在该应用背景下,欧几里得贴近度是最合适的贴近度定义。 因此我们将用欧几里得贴近度的实验结果与其他参照组进行比较。
[0145] 简单分段参照组中同样使用欧几里得贴近度。得到该参照组、特征匹配以及不变 矩方法的实验结果,与模糊识别的结果相比较,列于表2的混淆矩阵中。
[0146] 表2
[0149] 四种算法的识别正确率如图12所示。首先比较本文算法与均匀分段模糊识别算法 的结果。可以看到,均匀分段模糊识别对目标A的识别率高于本文算法,但目标B与目标C的 识别率均低于本文算法,总的平均识别率也低于本文提出的算法。因此,本实验中采用的自 适应分段算法一一均值梯度法,相比于原始的均匀分段,有效地提升了识别的正确率。随后 比较本文算法与特征匹配法。很明显,使用单一编码来描述目标的特征匹配法,其平均识别 正确率只有0.592,远远低于本文算法的0.781。特别的,目标C的识别率只有0.433,在实际 应用中这种水平识别率很可能无法应用。因此,本文的模糊识别算法大大提升了基于舰船 上层结构曲线编码的识别效果。最后比较本文算法与不变矩法。不变矩法的平均正确率稍 低于本文算法,然而,不变矩法对目标A的识别正确率只有0.423,远低于其对目标B和目标C 的识别率。同样地,这种水平的识别率在实际应用中很可能无法使用,使得该算法在多目标 情况下的稳定性存在缺陷。因此,本文提出的模糊识别算法,相比于现有的其他算法,不仅 具有最高的识别正确率,同时也具有最好的识别稳定性。实验结果证明了该算法的有效性, 以及算法对现有技术的提升。
[0150] 本发明还可有其它多种实施例,在不背离本发明精神及其实质的情况下,本领域 技术人员当可根据本发明作出各种相应的改变和变形,但这些相应的改变和变形都应属于 本发明所附的权利要求的保护范围。
【主权项】
1. 一种基于ISAR图像的舰船目标模糊识别方法,其特征在于一种基于ISAR图像的舰船 目标模糊识别方法具体是按照以下步骤进行的: 步骤一、采用Radon变换,对目标ISAR图像提取Centrel ine,并提取上层结构曲线; 步骤二、利用均值梯度法对提取的上层结构曲线进行自适应分段,并计算分段后的分 段相对均值特征; 步骤三、使用训练样本估计构造目标类别模板; 步骤四、对待识别目标和目标类别模板采用隶属度原则,计算并得到待识别目标和目 标类别模板的模糊集合形式及隶属度函数; 步骤五、计算待识别目标与不同类别模板的模糊集的贴近度,贴近度计算采用海明贴 近度,欧几里得贴近度或格贴近度,应用最大贴近度原则,将待识别目标判定为贴近度最大 的类别。2. 根据权利要求1所述一种基于ISAR图像的舰船目标模糊识别方法,其特征在于:所述 步骤一中采用Radon变换,对目标ISAR图像提取Centrel ine,并提取上层结构曲线;具体过 程为: 步骤--、提取 Centre 1 ine; 对目标ISAR图像进行全角度、全截距的Radon变换,得到变换域,在变换域内通过峰值 搜索,找到Centreline的倾角和投影点;提取并绘制出Centreline; 步骤一二、在步骤一一的基础上提取上层结构曲线; 提取出Centre line后,目标ISAR图像被Centre line分为两个部分,分别为包含桅杆或 船舱的上层结构部分和包含甲板与船底的船身部分;沿着Centreline,对目标ISAR图像所 在的距离单元内,提取每个距离单元内Centreline两侧目标ISAR图像与Centreline之间的 最大偏移量,得到两条曲线,分别为船身曲线和上层结构曲线;则偏移量大的曲线为上层结 构曲线。3. 根据权利要求2所述一种基于ISAR图像的舰船目标模糊识别方法,其特征在于:所述 步骤二中利用均值梯度法对提取的上层结构曲线进行自适应分段,并计算自适应分段后的 分段相对均值特征;具体过程为: 步骤二一、采用均值梯度法对提取的上层结构曲线进行自适应分段;具体如下: 步骤二一一、从提取的上层结构曲线一端C开始,以1/3曲线长度点为中心,取总曲线长 度的2/9或1/3为范围,计算范围内所有的点到端点C之间的上层结构曲线的均值,绘制出均 值曲线; 步骤二一二、对步骤二一一绘制出的均值曲线求一阶梯度,步长为1,找到一阶梯度的 最小值; 步骤二一三、以一阶梯度的最小值对应的位置作为端点C的分段点位置; 步骤二一四、从提取的上层结构曲线另一端点D开始,以1/3曲线长度点为中心,取总曲 线长度的2/9或1/3为范围,计算范围内所有的点到端点D之间的上层结构曲线的均值,绘制 出均值曲线;对均值曲线求一阶梯度,步长为1,找到一阶梯度的最小值;以一阶梯度的最小 值对应的位置作为端点D的分段点位置;完成自适应分段; 步骤二二、在步骤二一的基础上计算自适应分段后的分段相对均值特征; 分段相对均值特征用数值来表示每一段,每一段的数值由段内曲线高度的均值与总体 曲线均值的差除以总体曲线均值得到。4. 根据权利要求3所述一种基于ISAR图像的舰船目标模糊识别方法,其特征在于:所述 步骤四中对待识别目标和目标类别模板采用隶属度原则,计算并得到待识别目标和目标类 别模板的模糊集合形式及隶属度函数;具体过程为: 一、 待识别目标和目标类别模板的上层结构编码的模糊集; 将待识别目标和目标类别模板的上层结构表示为论域上的模糊集,模糊集的隶属度函 数计算依赖于分段相对均值特征,论域U上的待识别目标模糊集A,类别模板的模糊集B表示 为:其中,Uj表示j对应数字的二进制编码;A为论域U上的待识别目标模糊集,B为论域U上的 类别模板的模糊集; 二、 获取待识别目标和目标类别模板的上层结构编码的模糊集的隶属度函数的过程 为: (一) 、得到各段内编码〇和编码1的隶属度; 通过SCM特征计算各段内编码0和编码1的隶属度; SCM为分段相对均值,SCMi为SCM向量的第i个值,即第i段内的相对均值;用memiQ和memu分别表示第i段的编码0和编码1的隶属度;memiQ和memu隶属度的函数形 式为: (5) (6) 式中,X为SCMi; 因此,对某一待识别目标和目标类别模板的上层结构编码,第i段的编码0隶属度为 memio(SCMi),第i段的编码1隶属度为menui(SCMi); (二) 、计算各个编码Uj的总隶属度; 得到各段内的memiQ与memu后,计算编码Uj的总隶属度;若Uji是编码Uj的第i段的码元, 则Uj的总隶属度为: (7) 其中,Wi为权值。5. 根据权利要求4所述一种基于ISAR图像的舰船目标模糊识别方法,其特征在于:所述 步骤五中计算待识别目标与不同类别模板的模糊集的贴近度,贴近度计算采用海明贴近 度,欧几里得贴近度或格贴近度,应用最大贴近度原则,将待识别目标判定为贴近度最大的 类别;具体过程为: 最大贴近度原则为:假设Bk,k=l,…,c是论域U上的类别模板的模糊集,A为论域U上的 待识别目标模糊集,若: p(A,Bi) =maxk[p(A,Bk) ] (10) 则将待识别目标模糊集A的种类判定为出所属的类,即第1类;BiSBk中与待识别目标模 糊集A贴近度最大的值,B为某个类别模板的模糊集,1为B的编号,1 < 1 < c; 贴近度的计算采用海明贴近度Pffamming,欧几里得贴近度PEuclid或格贴近度Plattice ;对一 个有限论域1]={111,112,一,1!11}和有限论域1]上的两个模糊集4,8三种贴近度的定义分别为 :其中,A(Uj)与B(Uj)分别为待识别目标模糊集A和类别模板的模糊集B的隶属度函数;f 和Β°分别为待识别目标_吴糊集A和类别_吴板的_吴糊集B的补集,符号?和/\分别表吴糊集 的内积和取最小运算;η为论域中元素的个数,η为8; 对一个待识别目标模糊集A,其补集的隶属度函数为: Ac(x)=l-A(x) (14) 式中,x为SCMi; 对于两个实数a和b,a Λ b表示: aAb=min{a,b} (15) 式中,a为 A〇B;b 为 Ac〇Bc; 模糊集A和B的内积为: Α θ B= Vxeu[A(x)AB(x)] (16) 其中符号V为取最大值运算: Vx[A(x) ] =maxx[A(x) ] (17)〇
【文档编号】G06K9/62GK105868794SQ201610244344
【公开日】2016年8月17日
【申请日】2016年4月19日
【发明人】王勇, 朱鹏凯, 谢俊好, 李绍滨
【申请人】哈尔滨工业大学
网友询问留言 已有0条留言
  • 还没有人留言评论。精彩留言会获得点赞!
1