一种活性毁伤元爆裂毁伤快速精确图像识别方法

1.本发明涉及弹药毁伤及图像识别技术领域，具体涉及一种活性毁伤元爆裂毁伤快速精确图像识别方法。


背景技术：

2.活性材料是近十年来快速发展的一种新型含能材料，是武器弹药与装备研究的热点之一，其兼具类金属材料的力学性能与含能材料的爆炸性能，具有独特的性能优势。杀爆战斗部是毁伤目标运用最为广泛的战斗部，其利用爆炸后产生的破片飞散对目标进行毁伤。活性材料通过冷压与高温烧结制备而成的活性破片毁伤元以一定的初速碰撞目标后，活性破片在强冲击作用下会发生爆燃反应，释放大量化学能，从而实现对目标的侵爆联合高效毁伤。由于活性破片毁伤元是通过动能与爆炸产生的化学能时序联合毁伤目标，其在后效靶上爆裂产生的穿孔往往数量较多且呈现不规则的形状，导致破片爆裂毁伤的面积难以测量计算，穿孔面积计算一旦误差较大，将极大影响毁伤的后效评估，因此精确识别计算活性破片作用靶板的穿孔面积对于毁伤评估具有重要意义。
3.在图像处理与模式识别中，面积特征是十分重要的一个特征，不规则形状面积的识别与测定在现代工业生产具有广泛的应用领域，例如叶面积作为农作物栽培育种中常用的指标，被用于评价农作物产量与品质；医学中患者受伤皮肤与鼓膜穿孔面积的计算反馈准确而快速，给医生带来方便。相比于前两者，由于活性破片对后效靶作用穿孔轮廓的数量与不规则程度大大增加，虽然随着计算机相关技术的不断发展，实现了对不规则穿孔图像面积的快速识别与计算，但都存在着对环境要求高、机动性差、步骤复杂以及识别效率不高的缺点。


技术实现要素：

4.有鉴于此，本发明提供了一种活性毁伤元爆裂毁伤快速精确图像识别方法，通过图像处理技术，对图像中标尺进行识别和测定，对穿孔进行识别与检测，实现对后效靶爆裂穿孔面积的精确计算。
5.本发明的活性毁伤元爆裂毁伤快速精确图像识别方法，包括：
6.s1，将两个相互垂直的标尺放置在活性毁伤元毁伤后的后效靶上，爆裂产生的穿孔位于两个标尺的量程范围内，对后效靶进行成像；
7.s2，利用标尺刻度的颜色，对标尺进行识别，获得成像图片的像素比例尺；
8.s3，以成像图片中标尺量程范围内的区域作为工作区域，进行轮廓识别，并剔除轮廓内部的噪点以及狭长轮廓，获得穿孔轮廓；
9.s4，由穿孔轮廓内像素点数和s2确定的像素比例尺，获得穿孔面积。
10.较优的，所述s2具体包括：
11.s21，将s1获得的成像图片转到hsv空间，利用标尺刻度对应颜色的掩膜提取刻度所在区域，并利用设定的距离阈值对提取的区域进行合并，获得刻度区域；
12.s22，对所有刻度区域进行共线判断，获得两组标尺；
13.s23，针对各标尺，根据该标尺上相邻刻度区域间的像素长度与实际的标尺刻度间距，获得该标尺方向的像素比例尺。
14.较优的，所述s21中，首先对提取的区域进行面积判断，剔除面积小于或等于设定面积阈值的区域后，再进行合并。
15.较优的，利用opencv中轮廓提取函数findcontours()直接提取面积大于设定面积阈值的区域。
16.较优的，所述s22中，选取任意三个刻度区域，计算其中任意一个刻度区域的中心点到另外两个刻度区域中心点连线的垂直距离，若垂直距离小于设定的共线误差阈值d2，则认为该三个刻度区域共线；遍历所有刻度区域，得到各刻度区域的标尺归属。
17.较优的，所述s23中，提取同一标尺上的任意三个刻度区域，计算该三个刻度区域的最长间距，判断所述最长间距是否在设定的标尺长度范围[d
min
,d
max
]内，若是，则该三个刻度区域为相邻的刻度区域，利用相邻的刻度区域间的像素长度和现实标尺刻度间实际长度，确定该标尺方向的像素比例尺。
[0018]
较优的，所述s3具体包括：
[0019]
s31，提取工作区域的rgb分量，找到工作区域中穿孔子区域与靶板子区域差别最大的分量，并放大所述差别最大的分量，减小其他两个分量；转化为灰度图；
[0020]
s32，调整所述灰度图的亮度阈值，以穿孔轮廓最完整清晰时的亮度阈值为划分阈值，对所述灰度图进行二值分割，得到二值化图像；
[0021]
s33，对二值化图像进行轮廓识别，并剔除轮廓内部的噪点以及狭长轮廓，获得穿孔轮廓。
[0022]
较优的，所述s3中，对所有识别到的轮廓边界上的像素点数量进行判断，小于设定的轮廓周长阈值则认为是噪点，予以剔除；
[0023]
计算识别到的轮廓的面积周长比，若小于设定的轮廓圆度阈值cs，则进一步判断是否是狭长轮廓，若是，则予以剔除。
[0024]
较优的，狭长轮廓的判断方法如下：
[0025]
针对识别到的轮廓a，以轮廓a上任一轮廓点a为起点，利用函数findcontours()寻找同属于轮廓a的其他轮廓点，判断是否存在与轮廓点a序号之差大于参数d
l
的轮廓点存在，若存在则认为轮廓a为狭长轮廓。
[0026]
较优的，所述s3还包括s30：
[0027]
s30，对工作区域进行平滑滤波处理，去除图像噪声。
[0028]
有益效果：
[0029]
(1)本发明将标尺放置在后效靶上，首先基于图像处理技术识别测量标尺，获得图像与现实尺寸的比例尺关系；然后基于图像处理技术识别图像中的轮廓，并基于活性毁伤元爆裂产生的穿孔的特点，剔除非穿孔轮廓；最后利用穿孔轮廓面积及比例尺即可获得不规则形状的穿孔面积。本发明可对不规则形状穿孔轮廓的面积进行高效批量识别与计算，大幅提高了识别精度以及处理效率。本发明方法对仪器设备要求较低，具有可操作性强、方法简单、应用场景多样等特点；对于活性破片毁伤穿孔图像识别效率高，测量计算穿孔面积准确率高，与实际结果对比表明，相对误差在3％内，可短时高效、精确测量穿孔面积，便于
对活性毁伤元侵爆耦合毁伤模式进行毁伤评估。
[0030]
(2)本发明利用标尺刻度的颜色，对标尺刻度区域进行识别，同时利用同一刻度区域距离相近、不同刻度的刻度区域距离较远，对识别的区域进行合并，提高了对标尺红色区域的识别成功率；利用同一标尺刻度的共线性，将所有刻度区域归拢至对应标尺，具有较好的可实施性以及较高的准确率；根据同一标尺相邻刻度区域间的像素长度和现实标尺相邻刻度的真实距离，即可得到像素比例尺，利用该比例尺，只需求取穿孔区域的像素数，即可获得实际穿孔面积，可有效解决战斗部靶场试验中后效靶上穿孔数量多、形状复杂、毁伤面积难以准确测量的问题，只需计算经过图像处理后穿孔轮廓内部的的的像素点数量，且计算量低，可实现快速准确计算。
[0031]
(3)提取刻度区域时，首先剔除面积过小的区域，避免环境土壤等具有与标尺刻度相似颜色造成的误差，提高刻度区域识别率。
[0032]
(4)利用设定的d
min
和d
max
即可确认三个刻度区域是否是属于同一个标尺的相邻刻度区域，计算量小，方便快捷。
[0033]
(5)穿孔区域的轮廓提取过程中，首先利用穿孔子区域与靶板子区域的rgb差别最大的分量，通过放大该差别最大分量，减小另外两个分量，使得穿孔子区域和靶板子区域的灰度差别最大化，使得穿孔部分的灰度值趋于黑色，而靶板部分的灰度值趋于白色，去除不同分量差距过小，导致颜色分布不明显所带来的影响；同时，利用调整灰度图的亮度阈值，基于穿孔轮廓最完整清晰时对应的亮度阈值，对灰度图进行二值分割，进一步减少了噪点对于穿孔轮廓识别的影响，同时穿孔区域内部的噪点不出现连接成片的现象，穿孔轮廓又能够被比较准确地识别出来，大大提高了识别的准确性。
[0034]
(6)利用轮廓边界上的像素点数量剔除噪点；基于活性毁伤元爆裂毁伤穿孔的轮廓特点，利用轮廓圆度和狭长度，剔除非活性毁伤元爆裂产生的穿孔，提高穿孔轮廓的识别率。
[0035]
(7)利用函数findcontours()即可实现关于轮廓及其面积的筛选，如狭长孔等，方便快捷。
[0036]
(8)轮廓识别前，先对工作区域进行平滑滤波处理，去除图像噪声，提高轮廓识别准确率。
附图说明
[0037]
图1为本发明图像快速识别方法流程图。
[0038]
图2为本发明中红色区域识别过程示意图。
[0039]
图3为本发明中识别标尺流程图。
[0040]
图4为本发明中靶板图像沿基准线的rgb分量图。
[0041]
图5为本发明中靶板灰度分布直方图。
[0042]
图6为本发明中b不同取值条件下的图像分割结果示意图。
[0043]
图7为本发明中检测狭长孔示意图。
[0044]
图8为本发明中靶板图像处理过程示意图。
具体实施方式
[0045]
下面结合附图并举实施例，对本发明进行详细描述。
[0046]
本发明提供了一种活性毁伤元爆裂毁伤快速精确图像识别方法，通过python等语言，实现对活性破片毁伤后效穿孔面积的快速图像识别与精确计算。包括以下步骤：
[0047]
s1，将2个相互垂直的标尺放置在活性毁伤元毁伤后的后效靶上，后效靶上爆裂产生的穿孔位于两个标尺的量程范围内；对后效靶进行成像。
[0048]
s2，识别标尺，确定待测图像的像素比例尺：
[0049]
由于标尺上的整十刻度均为红色字体，可据此对待测图像中的标尺刻度进行识别，进而获取该标尺方向的像素比例尺。
[0050]
s21，首先，对待测图像中标尺的红色刻度进行识别。
[0051]
本实施例中，将原本为rgb空间的待测图像转到hsv空间，在hsv空间利用掩膜提取红色区域。利用opencv中轮廓提取函数findcontours()提取面积大于设定阈值a的红色区域；由于标尺上10、20、30等刻度由多个数字组成，利用轮廓提取红色区域时会将同属于一个刻度的红色区域分割为多个小区域；同时，单独数字的轮廓也可能在提取时被识别为多个小区域，因此，还需要将这个同属于一个刻度的小区域进行合并。可利用刻度间的间距，通过距离判断参数d1对提取的红色区域进行甄别，判断这些红色区域是否属于同一个刻度，若是，则进行区域合并，即将任意两个中心距离小于d1的红色区域进行合并，直到各区域间的距离均大于或等于d1，得到最终的刻度区域；
[0052]
s22，由于同一标尺上的红色刻度是共线的，因此通过判断刻度区域的中心点是否共线，即可确定刻度区域的标尺归属。本实施例中，选取任意的三个刻度区域，计算任意一个刻度区域的中点到另外两个刻度区域的中点连线的垂直距离，若垂直距离小于设定的阈值d2，则认为该三个刻度区域共线。遍历所有刻度区域，即可得到2个标尺上的刻度区域。若未能找到两组标尺，则向用户发出警告，提示人工处理。
[0053]
s23，获得所有刻度区域的标尺归属后，利用标尺上的相邻刻度区间的像素长度，比对现实标尺的刻度间距，即可得到该标尺方向的像素比例尺。
[0054]
本实施例中，利用设定的标尺长度范围[d
min
,d
max
]，判断标尺上任意三个刻度区间的最长间距是否在[d
min
,d
max
]范围内，若不在此范围内，则表示该三个刻度区间不相邻；由此，获得该标尺上相邻的刻度区间的像素长度。
[0055]
s3，提取后效靶上活性毁伤元造成的穿孔的轮廓。
[0056]
s30，活性破片穿孔面积都在标尺测量范围内，裁剪标尺量程范围内的图像作为工作区域，去除非工作区域图像带来的不必要影响及操作。
[0057]
s31，对裁剪的工作区域进行平滑滤波法处理，去掉所有可能影响图像中轮廓识别的噪声。工作区域分为穿孔子区域和靶板子区域；在工作区域设定一条基准线对比rgb三种分量的分布情况，找到穿孔子区域与靶板子区域差别最大的分量值，通过放大差别最大参量增大分量差距，缩小其他两种参量减小其他分量影响，然后通过灰度变换使得原图像的颜色分量r＝g＝b，经过变换后的穿孔部分灰度值趋于黑色，靶板部分灰度值趋于白色；
[0058]
s32，对工作区域亮度阈值的取值直接决定了穿孔轮廓识别和面积计算的准确性。为进一步扩大穿孔子区域和靶板子区域的差别，对灰度图进行二值分割。二值分割的具体方法为：调节灰度图的亮度阈值，通过对比不同亮度阈值下的穿孔子区域与靶板子区域的
灰度分布直方图，找到穿孔子区域最完整清晰地在图像上显现出来时对应的亮度阈值，以该亮度阈值为分割阈值，对灰度图进行分割处理，得到理想的二值化图像。
[0059]
s33，识别二值化图像中的轮廓，并剔除非活性毁伤元爆裂毁伤造成的穿孔的轮廓。
[0060]
调用opencv类库下的findcontours()识别二值化图像中的轮廓。在对二值化图像进行轮廓识别的过程中，由于客观原因，图像中的穿孔部分内部有很多噪点存在，这些噪点会被识别成许多细小的轮廓；同时，一些噪点连在一起的不规则轮廓，以及人工标记的线条也会被识别成狭长的线条型轮廓，这些会导致识别的不准确。观察发现，活性毁伤元爆裂毁伤造成的穿孔，虽然其轮廓不规则，但其大多是类圆形，且其面积远大于噪点，因此，本发明通过预先设定轮廓周长阈值cn和轮廓圆度阈值cs，来过滤这些噪点和狭长的线条型轮廓。
[0061]
具体的，对所有识别到的轮廓边界上的像素点数量进行判断，小于cn的轮廓即是噪点；当几个噪点轮廓连在一起或图像中有狭长的线条型轮廓时，通过判断其面积周长比来分析该轮廓的圆度，小于cs的轮廓即是不规则轮廓，对于这些轮廓，还需进一步判断是否为狭长型线条。因此需要使用狭长轮廓判断参数d
l
，筛选除去狭长线条。使用findcontours()方法寻找此点周围小范围区域内同属这一轮廓的所有其他轮廓点，判断是否有与此点序号之差大于参数d
l
的轮廓点存在，序号差距过大则说明此轮廓为狭长孔，不能作为识别到的穿孔轮廓。
[0062]
s4、除去以上干扰因素后，可以绘制出剩余的穿孔轮廓，利用contourarea()方法分别计算穿孔轮廓内部的像素点数n，结合s2确定的像素比例尺，即可求出穿孔面积并输出结果图像。
[0063]
最后对于标尺数量不为2，穿孔数量过多或穿孔轮廓与图像边缘相接等问题发出警告，可通过人工处理进一步确定轮廓识别的准确性。
[0064]
实施例1
[0065]
对图像中标尺的红色刻度进行识别，将待测图像从rgb空间转成hsv空间，由于hsv空间颜色由色温决定，在hsv空间中建立掩膜以屏蔽红色区域，如图2所示，通过findcontours()提取面积大于阈值a的红色区域，，借由距离判断参数d1判别识别的红色区域能否合并，对任意两个中心距离小于d1的小区域进行合并，直至红色区域的间隔均大于或等于d1，合并后的红色区域即为刻度区域。
[0066]
如图3所示，从若干个红色区域中任选三个区域，通过公式(1)判断区域是否共线，若三个区域中心点中，其中任意一点到其他两点连线的垂直距离小于共线误差d2即可得出三点共线的结论，然后通过标尺误差参数d
min
与d
max
判断所有共线三点连线的距离是否处于这两个参数之间，则可被确认为是一组相邻的刻度区域，分别找到两组标尺的相邻刻度区域，并根据其对应的实际长度即可计算求出像素比例尺。若未能找到两组标尺则向用户发出警告，提示需人工处理。
[0067][0068]
然后，在图像中工作区域建立一条基准线对比rgb三种分量沿基准线的分布，如图4所示，可知b分量上穿孔与靶板的差别更大，设定b参量为1.2，r、g参量定为-0.2以放大b分量增强分量差距体现区别性，缩小r、g分量减小分量影响，避免干扰。然后通过灰度变换使
得图像工作区域内的所有像素点的颜色分量r＝g＝b(公式(2))。灰度图中，穿孔部分灰度值趋于黑色，而靶板部分趋于白色。
[0069][0070]
经过灰度变换后的图像像素点除亮度外无其余差别，调整灰度图的亮度，进一步扩大穿孔部分与靶板部分的差异，以穿孔部分最完整且清晰地显现时的亮度为分割阈值，对灰度图进行二值分割，二值分割公式见公式(3)，其中f(x,y)为灰度变换图像，g(x,y)为二值化图像，b为分割阈值，分割阈值直接决定了穿孔面积识别与计算的准确性，如图5所示，待测图像的灰度分布直方图中穿孔与靶板部分的亮度为两个波峰处且有明显区别，但穿孔部分与靶板交界的轮廓部分较为模糊，根据图6可以看出随着b逐渐增大，穿孔内部的噪点也越来越少，但轮廓部分面积越来越大，逐渐失去了识别准确性，当b为80的时候，穿孔趋于轮廓较为完整清晰，内部噪点也不连成片状，各自独立，因此b取80时可以得到理想的状态，得到二值化图像。
[0071][0072]
通过opencv类库下的findcontours()方法识别穿孔轮廓部分后，如图7所示，通过对轮廓边界上的像素点数量判断，其大于轮廓周长参数cn时，可筛选掉轮廓内部尺寸较小的噪点；几个噪点轮廓连在一起时，通过判断其面积周长比是否大于轮廓圆度参数cs来分析该轮廓的圆度，从而筛选掉形状不规则的小轮廓；从图6中可以看出图像中还有标记的线条与可能存在的被微小破片细微切向作用到的小轮廓，它们共同组成了狭长型轮廓，会影响穿孔轮廓识别与面积的计算，因而需要借由狭长轮廓判断参数d
l
，筛选除去狭长轮廓，使用findcontours()方法寻找此点周围小区域内同属这一轮廓的所有其他轮廓点，判断是否有与此点序号之差大于参数d
l
的轮廓点存在，序号差距过大则证明此轮廓为狭长孔，不计入识别的穿孔轮廓之中。
[0073]
绘制出穿孔轮廓后，通过contourarea()方法分别计算穿孔轮廓内部的像素点数量，计算穿孔轮廓内部的像素点数n，然后基于由标尺识别确定的像素比例尺，由公式(4)即可得出穿孔面积并输出结果图像，公式(4)中，s为求得的穿孔面积，n为穿孔轮廓内部的像素点数量，l为识别的标尺像素长度，200为识别的标尺对应的实际长度，单位为mm。
[0074]
s＝n/(l/200)2ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(4)
[0075]
最后对于标尺数量不为2，穿孔数量过多或弹孔轮廓与图像边缘相接等问题发出警告，可通过人工处理进一步确定轮廓识别的准确性。靶板图像处理的过程如图8所示。
[0076]
实施例1中所用预设参数见表1。
[0077]
表1预设参数值
[0078]
d1d2d
mindmaxcncsdl
7030160020005010200
[0079]
为验证此方法的准确性，另选五组不同的图像进行测量实验，如表2所示。首先通过绘图工具找到标尺上整十刻度的坐标，计算得出标尺的像素长度。然后为靶板着色，将弹孔以外的部分全部绘制为白色，并通过matlab软件计算弹孔区域的像素点数量，根据标尺
进行比例换算，得到的结果可近似视为实际面积。结果表明，本发明计算得到的靶板穿孔面积与实际结果的相对误差在3％以内，能够达到工程应用的精度要求。
[0080]
表2本发明计算结果与实际结果对比
[0081][0082]
综上所述，以上仅为本发明的较佳实施例而已，并非用于限定本发明的保护范围。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。