近矩形平面状工业产品表面瑕疵检测方法与流程
本发明属于工业产品表面瑕疵检测方法技术领域，特别是一种计算量小、实时性好的近矩形平面状工业产品表面瑕疵检测方法。
背景技术：
流水线生产过程中有可能在产品外壳上形成划痕、擦伤等表面瑕疵，表面完好度检测是很多工业产品出厂合格检测的必检项目之一。
表面瑕疵检测传统上主要依靠人工目测检查的方式。长时间从事这种单调的外观瑕疵检测极易引起质检员视觉疲劳，从而影响检测效率和正确率。工业生产线光照条件稳定，单一流水线产品类型固定，比较适合采用机器视觉的方式在线自动检测产品表面瑕疵。基于机器视觉的工业品外观瑕疵检测方案在纺织、通信等行业都有应用案例，已有应用案例大都在空间域进行数字图像处理，通过与预存的标准产品外观模板对比评估待检产品表面是否存在瑕疵。为了发现一些较为细小的划痕、污点等瑕疵，用于表面瑕疵检测的机器视觉系统一般选择分辨率较高的工业相机采集数字图像。
目前大部分检测系统采用空间域处理的方式比对待检产品和标准图像，在图像分辨率较高的情况下，计算量较大，实时性差。
技术实现要素：
本发明的目的在于提供一种近矩形平面状工业产品表面瑕疵检测方法,计算量小、实时性好。
实现本发明目的的技术解决方案为：
一种近矩形平面状工业产品表面瑕疵检测方法,包括以下步骤：
(10)模板图像生成：在流水线正常照明状态下采集表面完好的标准产品的外观灰度图像，检测该图像边缘，获取产品外轮廓线，依据外轮廓线裁剪出标准产品外观模板图像；
(20)待检测产品图像生成：采集待检测产品外观灰度图像，寻找产品外轮廓矩形，裁剪出待检测产品区域，调整分辨率，得到待检测产品图像，将模板图像与待检测产品图像拼接形成一张整体图像；
(30)疑似瑕疵区域检测：将整体图像进行二维离散傅里叶变换和傅里叶反变换，在经反变换后的图像中剪裁出拼接时的待检测产品图像对应区块，将其二值化得到疑似瑕疵区域；
(40)瑕疵区域检测：将疑似瑕疵区域与模板图像中相同位置附近区域进行对比，如最小差异超出阈值范围，则确定该处存在表面瑕疵。
本发明与现有技术相比，其显著优点为：
1、计算量小：本发明首先通过高效的频域分析检测疑似瑕疵区域，然后仅针对这些疑似瑕疵区域实施计算量较大的瑕疵区域精细比较过程。由于一张待检测产品图像中疑似瑕疵区域数量较少，因此与常用的逐像素对比方法相比明显减小了计算量。
2、实时性好：本发明使用的图像快速傅里叶变换有很高的计算效率，同时得益于提取的疑似瑕疵区域待对比数据量较小，因此比常见的基于空间域分析的方法效率更高。
下面结合附图和具体实施方式对本发明作进一步的详细描述。
附图说明
图1为本发明近矩形平面状工业产品表面瑕疵检测方法的主流程图。
图2为图1中模板图像生成步骤的流程图。
图3为模板图像生成结果示例图。
其中，图3a为标准产品灰度图像，图3b为边缘检测与外轮廓矩形检测结果，图3c为裁剪得到的产品模板图像。
图4为图1中待检测产品图像生成步骤的流程图。
图5为待检测产品图像生成结果示例图。
其中，图5a为待检测产品灰度图像，图5b为边缘检测与外轮廓矩形检测结果，图5c为裁剪得到的待检测产品图像。
图6为图1中疑似瑕疵区域检测步骤的流程图。
图7为模板图像与待检测产品图像拼接结果示例图。
图8为图1中瑕疵区域检测步骤的流程图。
图9为待检测产品外观瑕疵检测结果示例图。
其中，图9a为疑似瑕疵区域，图9b为确定瑕疵区域，图9c为瑕疵检测结果。
具体实施方式
如图1所示，本发明近矩形平面状工业产品表面瑕疵检测方法,包括以下步骤：
(10)模板图像生成：在流水线正常照明状态下采集表面完好的标准产品的外观灰度图像，检测该图像边缘，获取产品外轮廓线，依据外轮廓线裁剪出标准产品外观模板图像；
如图2所示，所述(10)模板图像生成步骤包括：
(11)边缘图像获取：在流水线正常照明状态下，采集表面完好的标准产品的外观灰度图像itemp-o，测量外观灰度图像中产品外轮廓矩形的宽度、高度，使用标准canny边缘检测算子检测外观灰度图像itemp-o边缘，得到二值化边缘图像itemp-edge；
(12)外轮廓线获取：在二值化边缘图像itemp-edge中使用标准hough变换方法检测直线，在检测得到的直线簇中选择长度与产品外轮廓矩形长和宽差异不超过δl的直线段作为疑似外轮廓线，先选择二值化边缘图像itemp-edge中位于最左恻和最上恻的疑似外轮廓线作为实际轮廓线，然后以步骤1-1中测量的产品外轮廓矩形在itemp-o中的宽和高作为先验知识，选择与左侧和上侧轮廓线平行且与其距离最接近产品宽和高的直线段作为右侧和下侧轮廓线；
(13)外观模板图像获取：用标准图像旋转方法旋转四条轮廓线围成的区域，使上侧轮廓线旋转后呈水平，然后将旋转后轮廓线围成的区域保存为产品外观模板图像itemp。
图3所示为模板图像生成的结果示例图。其中，图3a为标准产品灰度图像，图3b为边缘检测与外轮廓矩形检测结果，图3c为裁剪得到的产品模板图像。
(20)待检测产品图像生成：采集待检测产品外观灰度图像，寻找产品外轮廓矩形，裁剪出待检测产品区域，调整分辨率，得到待检测产品图像，将模板图像与待检测产品图像拼接形成一张整体图像；
如图4所示，所述(20)待检测产品图像生成步骤包括：
(21)外观灰度图像采集：采集待检测产品的外观灰度图像；
(22)外轮廓矩形寻找：使用canny边缘检测算子，在外观灰度图像中寻找产品外轮廓矩形；
(23)待检测产品图像获取：旋转外观灰度图像，使产品外轮廓矩形上边缘水平，依据外轮廓裁剪出待检测产品区域，并将该区域缩放至与模板图像itemp同分辨率，得到待检测产品图像itest，
(24)图像拼接：将m*n-1个模板图像itemp与1个待检测产品图像itest整体按照m*n排列拼接形成一张整体图像idet。
所述(24)图像拼接步骤包括：
(241)图像缩放：在拍摄的包含待检测图像的原始图像中裁剪出待检测产品图像itest并缩放至与模板图像itemp4相同分辨率；
(242)图像排列：将m*n-1个模板图像itemp与1个待检测产品图像itest分别作为整体图像块，按照m*n排列拼接形成一张图像idet，其中将itest放置于
图5为待检测产品图像生成结果示例图。
其中，图5a为待检测产品灰度图像，图5b为边缘检测与外轮廓矩形检测结果，图5c为裁剪得到的待检测产品图像。
(30)疑似瑕疵区域检测：将整体图像进行二维离散傅里叶变换和傅里叶反变换，在经反变换后的图像中剪裁出拼接时的待检测产品图像对应区块，将其二值化得到疑似瑕疵区域；
如图6所示，所述(30)疑似瑕疵区域检测步骤包括：
(31)傅里叶变换：将整体图像idet进行二维离散傅里叶变换得到idet-f，如式(1)
其中mdet，ndet为idet的高和宽，udet-f、vdet-f、xtest、ytest为像素坐标，j为虚数单位，r(udet-f,vdet-f)为傅里叶变换后的实数部分，i(udet-f,vdet-f)为虚数部分，以下公式中类似参数含义相似仅以下标区别，式(1)中：
(32)幅度谱计算：使用式(2)计算傅里叶变换后idet-f的幅度谱
(33)傅里叶逆变换：然后使用式(4)对idet-f'进行二维离散傅里叶逆变换，即
在idet'中剪裁出原itest放置位置的图像块保存为itest'；
(34)疑似瑕疵区域获取：使用标准otsu算法将itest'二值化，二值化图像中亮度为255的像素点为疑似瑕疵区域。
图7为模板图像与待检测产品图像拼接结果示例图。
(40)瑕疵区域检测：将疑似瑕疵区域与模板图像中相同位置附近区域进行对比，如最小差异超出阈值范围，则确定该处存在表面瑕疵。
如图8所示，所述(40)瑕疵区域检测步骤包括：
(41)待对比图像块选取：依次选择疑似瑕疵区域中的疑似瑕疵点作为中心像素点，取边长wn像素的正方形窗口内的近邻像素与中心像素点一起作为待对比图像块；
(42)待对比图像块差异计算：记该待对比图像块为imgpatchtest(i)，其中心像素在itest'中的图像坐标为(xi,yi)，在模板图像itemp中，以(xi,yi)为中心，取边长wt的正方形对比窗口sw，依次以窗口sw中的像素为中心像素，在模板图像itemp中取中心像素与边长wn像素的正方形窗口内的近邻像素组成对比图像块imgpatchtemp(j)，0<j<wt2+1，将wt2个imgpatchtemp(j)分别与imgpatchtest(i)对比计算差异；
所述(42)待对比图像块差异计算步骤中，对比计算差异具体为：
假设当前带对比图像块为imgpatchtest(i)，对比模板图像块为imgpatchtemp(j)，每个图像块中包含wn2个像素，使用式(5)计算imgpatchtest(i)与imgpatchtemp(j)的差异值：
式中，gmax为所使用灰度图像格式的可取灰度级最大值，itest'中的每个imgpatchtest(i)与item中的wt2个对比图像块做对比得到wt2个差异值。
(43)瑕疵判定：记取最小的差异值为emin，如果emin大于设定的差异阈值eth，则待检测产品图像itest中(xi,yi)处即确定为瑕疵，否则不为瑕疵。
图9为待检测产品外观瑕疵检测结果示例图。
其中，图9a为疑似瑕疵区域，图9b为确定瑕疵区域，图9c为瑕疵检测结果。
本发明使用快速频域变换技术，利用频域滤波检测疑似瑕疵区域，然后在空间域中仅把疑似瑕疵区域与模板图像对比，通过频率域与空间域混合分析方法显著减少了对比时的计算量；利用单张图像重复排列拼接方法人为构造图像中的整体重复模式，使用简单的频域滤波算法即可发现待检测图像与模板图像不一致的区域，从而实现瑕疵检测；与一般使用sift等特征算子匹配方法匹配模板图像与待检测图像不同，本方法通过检测产品轮廓矩形的方法在图像中裁剪出产品区域并调整角度使产品上边缘平行，然后将裁剪出的待检测产品图像缩放至和裁剪出的模板图像同分辨率实现图像匹配，在确定的生产线环境下，检测摄像机安装位置固定，产品区域在模板图像和待检测图像中实际缩放比例很小，因此这种方法比稠密特征算子匹配方法更加简便、高效。
- 还没有人留言评论。精彩留言会获得点赞！