郑晓玲， 刘斌

(1. 数字化视觉测量厦门市重点实验室， 福建 厦门 361021;

2. 黎明职业大学 机电及自动化学院， 福建 泉州 362000)



采用机器视觉的铝压铸件表面缺陷检测

郑晓玲1,2， 刘斌1

(1. 数字化视觉测量厦门市重点实验室， 福建 厦门 361021;

2. 黎明职业大学 机电及自动化学院， 福建 泉州 362000)

摘要：对铝压铸件常见的缺陷形态进行分析，提出适用于铝压铸件表面缺陷的检测算法.首先，采用阈值分割与形态学相结合的方法分割出可疑区域；然后，根据基于面积、亮度均值、亮度均值差、灰度曲线分析的4个剔除原则，剔除伪缺陷的干扰.试验结果表明：文中的检测方法具有低成本、高精度、可操作性强等优点，能有效提高铝压铸件生产过程中的检测效率.

关键词：铝压铸件; 缺陷检测; 伪缺陷; 机器视觉

在铝压铸件生产过程中，由于种种原因，铝压铸件表面不可避免地产生气孔、裂纹、划痕等缺陷.这些缺陷严重影响产品的表面质量和物理机械性能.因此，生产过程中的缺陷检测具有十分重要的意义.目前，国内铸造生产车间多采用人工目视检测法，主要依赖于人工经验，劳动强度大且工作效率低.随着市场对压铸件质量要求的日益提高，以及企业劳动力成本的提升，依赖于人工检测的传统质量监控体系已无法适应现代生产发展的要求.在铸件的大批量生产中，基于机器视觉的实时在线无损检测成为压铸件质量检测的发展方向.国内外许多学者进行了相关研究[1-4]，但针对铝压铸件表面缺陷检测的研究成果较少.由于铝压铸件形状复杂，最常见的气孔缺陷尺寸细微，难以准确分割出缺陷区域，并且在生产过程中出现的水渍、油污等易被误识为缺陷，对检测结果造成干扰.因此，本文对铝压铸件常见缺陷形态进行分析，提出适用于铝压铸件表面缺陷识别的算法.

图1　机器视觉检测流程Fig.1　Process of machine vision

1检测过程

表面缺陷视觉检测流程，如图1所示.图像处理算法可分为两个步骤.1) 采用较快速的算法检测到感兴趣的可疑区域，将其分割出来,保存至可疑图像缓存区.这一步骤不仅可以初步判断是否有可疑的目标，还能缩小需要后续处理的图片数据.2) 对可疑图像缓存区的对象进行筛选.在铸件生产过程中出现的水渍、油污等不属于铸件缺陷，但由于其外观形貌在图像上与缺陷非常类似，易被检测系统误识为缺陷.对这些伪缺陷的处理不仅增加系统检测的负担，而且影响对严重缺陷的正确识别.因此，在判断出可疑目标的基础上，还要进一步深入研究可疑区域上是否有真实缺陷，剔除伪缺陷的干扰.

1.1可疑区域分割

1.1.1分割算法的选择基于边缘检测的可疑区域分割算法和影像相减法是缺陷检测系统常用的检测算法[5-7]. 但这两种算法各有局限之处， 并不适用于复杂形状铸件细微缺陷的检测. 文中采用阈值分割和形态学相结合的方法.气孔缺陷处亮度较图片平均亮度偏暗，因此，通过阈值分割获得灰度值较低的像素点.阴影和噪声点区域的灰度值和缺陷区域相似，在分割出缺陷区域时，也分割出干扰区域.利用形态学中的膨胀技术将气孔缺陷与阴影、噪声点的特征进行区分，从而准确分割出缺陷区域.

1.1.2阈值分割检测黑点有气孔缺陷的汽车启动机盖图像，如图2所示.为了提高检测结果的精确性，在RGB模型下选取一定阈值，将低于此阈值的赋值为1，其余赋值为0.RGB模型下检测黑点结果，如图3所示.由图3可知：受成像过程中出现的阴影及噪声点的影响，检测出缺陷区域的同时，也出现了大量的干扰区域.因此，需要通过后续途径去除干扰，准确分离出缺陷区域.

图2　有气孔缺陷的汽车启动机盖图像　　　　　　　　图3　RGB模型下检测黑点结果图　Fig.2　Car boot cover image with porosity defects Fig.3　RGB model for the detection of black spot′s result

1.1.3形态学图像处理在图像获取过程中，由于相机噪声、光照不均、被测物表面细小灰尘等因素影响，获得的图像会产生一些噪声点，对正确识别出可疑区域造成干扰.将数学形态学中的灰度膨胀应用于缺陷检测算法中，可有效地去除噪声点的干扰.

数学形态学是用具有一定形态的结构元素提取和度量图像中的对应形状，以达到分析识别图像的目的.因其具有天然并行结构实现的优点，可实现处理算法和形态学的并行，进而快速提高图像分析和处理速度.数学形态学的基本运算有腐蚀、膨胀、开启、闭合等4种.

若输入图像为f(x,y)，结构元素为B(i,j)，则形态学膨胀定义为

(1)

式(1)中:Df为图像大小；DB为结构元素大小.

结构元素的选取除了考虑形状的相似性及尺寸的覆盖性外，还应考虑到腐蚀的结果是否能突显出目标区域.因后续的检测方法基于连通区域的几何特征，故结构元素像素大小的选取尤为重要.若结构元素过小，则无法区分目标区域与噪声点；若结构元素过大，则目标区域与背景区域融合为同一个连通区域.经过多次实验效果比对，最终选用15个像素的方块结构元素对图像进行膨胀处理.

1.1.4连通区域面积计算标记连通区域，参考文献[8]提出的算法，按列对图像进行扫描编码，每一个非零像素团记为一个run，并赋予其初始标签.记录所有run的起始行、结束行及所在列.若某一个run的相邻列存在与之在行上交叠的run，则判定两者属于同一个连通域，赋予相同标签.若某一列的第k个run与相邻前一列的n个run(n>1)都存在行上的交叠，则判定以上n+1个run都是属于同一个连通区域[8].

完成所有连通区域标记后，计算每个标签对应的连通区域面积，保留面积值在阈值范围内的连通块，删除其余连通块，得到缺陷掩模图像，如图4所示.将掩模图像与图2相乘，成功分割出缺陷区域，如图5所示.

图4　缺陷掩模图像 图5　分割出来的缺陷区域　Fig.4　Defect mask image Fig.5　Segmented defect area

通过试验验证，该算法能检测到直径为0.1mm的气孔缺陷，平均每幅1 280px×960px的图像处理时间为0.061s，可获得准确的缺陷焦点位置，完整保存缺陷形状，以利于后续分析.

1.2伪缺陷剔除

对从原图分割出来的可疑区域样本进行研究，可以发现：真缺陷区域大小在一定范围内，且较周围存在一定范围的亮度差，区域内部的灰度变化较均匀.

分别对气孔和油污样本的区域中心提取行灰度曲线，如图6所示.图6中：g为灰度；s为位置.由图6可知：气孔区域中心行灰度曲线存在一个很明显的谷底，而油污的区域中心行灰度曲线则存在大量的尖峰抖动.铝压铸件气孔为一个凹坑区域，其表面的反射光线由于方向的改变，无法进入相机，在成像效果上呈现暗区域效果，相对周围像素点的灰度值有明显变化，区域内部处于一个一定范围的灰度渐变区，理论上沿任何方向扫描，其灰度曲线都能存在一个明显的山峰区域；而油污伪缺陷区域表面其实是平整的，并不存在反射光线方向的改变，其成像效果上的暗区域是由油污本身的物质结构所决定.在生产过程中，油污很少处处均匀地附着在铝压铸件表面，其区域内部的灰度值并不是有规律地均匀渐变，故提取行灰度曲线会发现存在多处尖峰抖动.

针对铝压铸件生产过程中容易出现的油污、水渍等伪缺陷，提出以下4个剔除原则.1) 计算可疑区域面积，剔除在阈值范围外的区域.2) 计算可疑区域亮度均值，剔除在阈值范围外的区域.3) 计算可疑区域与周围背景的均值亮度差，剔除差值在阈值范围外的区域.4) 用数学的方法分析气孔区域和油污区域中心行灰度曲线波形差异，计算能够表征此差异的特征值，剔除特征值在阈值范围外的区域.

(a) 气孔缺陷区域中心行灰度曲线1 (b) 气孔缺陷区域中心行灰度曲线2

(c) 气孔缺陷区域中心行灰度曲线3 (d) 气孔缺陷区域中心行灰度曲线4

(e) 油污区域中心行灰度曲线1 (f) 油污区域中心行灰度曲线2

(g) 油污区域中心行灰度曲线3 (h) 油污区域中心行灰度曲线4图6　真伪缺陷的中心行灰度曲线Fig.6　Center line gray curve of true and false defect

原则4)的波形分析具体有以下5个步骤.

步骤1对可疑区域提取中心行灰度曲线，如图7所示.

(a) 伪缺陷样本 (b) 真缺陷样本

(c) 伪缺陷区域 (d) 真缺陷区域图7　中心行灰度曲线Fig.7　Center line gray curve

步骤2提取灰度曲线所有极大、极小值点.将得到的行灰度曲线离散化，得到曲线的各个点为(x1,y1)，(x2,y2),…,(xn,yn).其中:n>1;x,y分别为各点的像素位置及像素值.记曲线点集为C,有

(2)

式(2)中：yi为C集合中第i个点的像素值.当m>0时，所对应的点(xi,yi)为极值点，将其储存入极值点集合M中.

步骤3最小二乘法拟合曲线.假设集合M中有m个点，记为(ti,yi)，j=1,2,…,m.通过这些数据寻找一个m次近似多项式，有

(3)

式(3)中:φk(t)=tk,k=0,1,…,m.求系数a0,a1,…,am，使误差的偏差平方和

(4)

最小，则P(t)称为最小二乘拟合多项式，或称为t,r的经验公式.可用多元函数极值的方法求其最小点和最小值，得驻点方程为

(5)

即

(6)

利用内积的定义式，式(6)可写为

(7)

式(7)称为正规方程组，通过这个正规方程组解出ak的值，代入多项式P(t)中，即可得到拟合的曲线方程[9].图8(a)，(b)分别为图7(c)，(d)的拟合结果.

(a) 真缺陷 (b) 伪缺陷图8　真伪缺陷拟合结果Fig.8　Fitting result of false and true defect

步骤4计算拟合曲线的最大曲率值为

图9　真伪缺陷特征曲率的差异　　Fig.9　Differences in the characteristics 　　curvature of true and false defect

(8)

根据式(8)计算拟合曲线上每一点的曲率，取最大值作为该样本图片的特征值.

步骤5选取曲率阈值.真伪缺陷区域的特征曲率有着较为明显的区别，如图9所示.图9中：k为曲率；n为样本编号.将特征值在曲率阈值以下的图片判定为伪缺陷图片，从缓存区域删除.

2试验分析

为了验证铝压铸件缺陷检测方法的有效性，从生产车间现场采集137幅含有若干真伪缺陷的样本图像进行试验.其中，气孔缺陷区域109处，油污42处，水渍4处.具体检测结果，如表1所示.表1中：n为实际数目；n1为检出数；η1为检出率；n2为误检个数；η2为误检率.由表1可知：有102处气孔缺陷被成功检测，有4处油污被误识为缺陷，检出率和误检率分别为 93.5%，9.5%.造成漏检的原因是少部分气孔缺陷出现在很靠近边缘的区域，在区域分割时被算法判别为阴影或噪声点.在后续工作中，可通过改进图像采集平台对算法进行辅助.

表1　检测结果

3结束语

气孔是铝压铸件铸造过程中最常见的缺陷，其尺寸通常在毫米以下级别.采用阈值分割和形态学相结合的方法能准确分割出气孔区域，适用于在灰度不均匀,且存在噪声干扰的图像中分割出小目标区域的情况.对于油污水渍这类伪缺陷，再通过增加4个筛选原则从可疑区域中剔除干扰，从而检测出真正的缺陷.相比于传统人工检测方法，文中的检测方法具有低成本、高精度、可操作性强等优点，能有效提高铝压铸件生产过程中的检测效率.

参考文献：

[1]徐科.基于光度立体学的金属板带表面微小缺陷在线检测方法[J].机械工程学报,2013,49(4):25-29.

[2]周新星.基于非下采样Contourlet变换和PCNN的表面缺陷自动识别方法[J].应用基础与工程学科学报,2013,21(1):174-184.

[3]段志娟.基于独立成分分析的冷轧带钢表面缺陷识别[J].钢铁研究学报,2011,23(10):63-66.

[4]ELBEHIERY H.Surface defects detection for ceramic tiles using[J].World Academy of Science,2007,12(10)：158-162.

[5]张鑫.基于图像处理的注塑制品缺陷检测方法研究[D].沈阳:东北大学,2009：34-38.

[6]彭向前.产品表面缺陷在线检测方法研究及系统实现[D].武汉：华中科技大学,2009：25-27.

[7]吴炜峰.基于水平集的灰度不均匀图像分割研究[D].广州:华南理工大学,2013：41-42.

[8]HARALICK,ROBERT M,LINDA G.ShAPIRO, computerand robot vision[J].Addison-Wesley,1992,1(8):28-48.

[9]孙志忠.计算方法与实习[M].南京:东南大学出版社,2011:116-120.

(责任编辑： 钱筠英文审校： 杨建红)

Surface Defect Detection of Aluminum Die Casting Using Machine Vision

ZHENG Xiaoling1,2, LIU Bin1

(1. Xiamen Key Laboratory of Digital Vision Measurement, Xiamen 361021, China;2. Institute of Electrical and Automation， Liming Vocational College, Quanzhou 362000, China)

Abstract：The common defects of aluminum die casting were analyzed in this paper, then a detection algorithm suitable for aluminum die casting surface defect was put forward. First, by using the combined method of threshold segmentation and suspicious area. Then, based on four removing principles: the area, the mean brightness, the mean brightness difference, gray curve analysis, false defects interference was eliminated. The experimental results showed that the detection method had the advantages of low cost, high precision, strong operability, and could effectively improve the detection efficiency in the process of aluminum die castings production.

Keywords：aluminum die casting; defect detection; false defect; machine vision

中图分类号：TP 391

文献标志码：A

基金项目：国家自然科学基金面上资助项目(51175191)

通信作者：刘斌(1970-)，男，教授，博士，主要从事数字化设计的研究.E-mail:mold_bin@hqu.edu.cn.

收稿日期：2014-10-29

doi:10.11830/ISSN.1000-5013.2016.02.0139

文章编号：1000-5013(2016)02-0139-06