CN115908408B - 基于多邻域预测模型的磁芯缺陷检测方法 - Google Patents

基于多邻域预测模型的磁芯缺陷检测方法 Download PDF

Info

Publication number
CN115908408B
CN115908408B CN202310011002.XA CN202310011002A CN115908408B CN 115908408 B CN115908408 B CN 115908408B CN 202310011002 A CN202310011002 A CN 202310011002A CN 115908408 B CN115908408 B CN 115908408B
Authority
CN
China
Prior art keywords
pixel point
pixel
neighborhood
pixel points
prediction model
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202310011002.XA
Other languages
English (en)
Other versions
CN115908408A (zh
Inventor
王宪保
雷雅彧
项圣
翁扬凯
吴帅杰
方卫杰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhejiang University of Technology ZJUT
Original Assignee
Zhejiang University of Technology ZJUT
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Zhejiang University of Technology ZJUT filed Critical Zhejiang University of Technology ZJUT
Priority to CN202310011002.XA priority Critical patent/CN115908408B/zh
Publication of CN115908408A publication Critical patent/CN115908408A/zh
Application granted granted Critical
Publication of CN115908408B publication Critical patent/CN115908408B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02PCLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P90/00Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
    • Y02P90/30Computing systems specially adapted for manufacturing

Landscapes

  • Image Processing (AREA)
  • Investigating Materials By The Use Of Optical Means Adapted For Particular Applications (AREA)

Abstract

本发明公开了一种基于多邻域预测模型的磁芯缺陷检测方法,包括:对待学习无缺陷样本图片中的每个像素点设置特征表示区;确定每个像素点的感受野区域,并将该像素点作为感受野区域的中心像素点;对每个中心像素点,在其感受野范围内,找出k个与中心像素点相关度最高的邻域像素点,作为最相关邻域像素点;建立每个中心像素点与其最相关邻域像素点之间的多邻域映射关系,得到多邻域预测模型;调用多邻域预测模型数据对待测图片进行缺陷检测,得到待测图片缺陷检测结果的二值图像;对得到的缺陷检测结果的二值图像进行去噪优化,得到最终缺陷检测图像。本发明只需要对一幅正样本图片进行学习,具有样本需求量少,训练时间快,检测精度高的优点。

Description

基于多邻域预测模型的磁芯缺陷检测方法
技术领域
本发明涉及机器视觉领域,尤其涉及一种基于多邻域预测模型的磁芯缺陷检测方法。
背景技术
在产品生产过程中,由于工艺和环境的复杂性,许多类型的产品不可避免地存在如凹坑、划痕、异物以及掉角等表面缺陷,这些缺陷虽然面积较小,但严重影响了产品质量和性能,甚至会降低产品的使用寿命。因此,产品表面缺陷检测是工业生产的重要环节之一。
传统的通过人工肉眼检测的方式识别产品外观缺陷,并人工剔除有缺陷产品,易受工人主观因素影响,并且检测花费时间长效率低,已经无法满足人们对检测效率及准确度的要求。
近年来,随着机器视觉,人工智能的发展,许多产品的缺陷检测由传统的人工目测转变为使用自动化机器视觉缺陷检测,其中许多基于深度学习的磁芯缺陷检测方法表现优异,然而在实际的生产中,缺陷种类多,实际缺陷样本少,无法满足深度学习需要大量带标签样本数据的需求,不足以利用深度学习的方法进行训练。
发明内容
针对现有技术的不足,本发明提出一种基于多邻域预测模型的磁芯缺陷检测方法,具体技术方案如下:
一种基于多邻域预测模型的磁芯缺陷检测方法,包括以下步骤:
S1:对待学习无缺陷样本图片中的每个像素点设置特征表示区,每个像素点的特征表示区均为以该像素点为中心,包含n个像素点的区域,n为正整数;
S2:确定每个像素点的感受野区域,并将该像素点作为所述感受野区域的中心像素点,所述感受野区域包含m个像素点,m为正整数,且m>n;
S3:对每个所述中心像素点,在其感受野范围内,找出k个与所述中心像素点相关度最高的邻域像素点,作为最相关邻域像素点,所述k小于m;
S4:建立每个所述中心像素点与其k个所述最相关邻域像素点之间的多邻域映射关系,得到多邻域预测模型;
所述S4通过以下子步骤实现:依次遍历所述待学习无缺陷样本图片的每一个中心像素点,计算该中心像素点和每一个与其对应的k个最相关领域像素点之间的像素差,并将所述像素差进行拟合,得到其高斯分布的均值和方差;将所述中心像素点的坐标、所述中心像素点对应的k个最相关领域像素点的坐标、高斯分布的均值和标准差分别保存;当检测完所有中心像素点后,保存所有中心像素点的坐标、该中心像素点对应的k个最相关领域像素点的坐标、高斯分布的均值和标准差,得到多领域预测模型;
S5:调用所述多邻域预测模型数据对待测图片进行缺陷检测,得到所述待测图片缺陷检测结果的二值图像;
S6:对S5得到的缺陷检测结果的二值图像进行去噪优化,得到最终缺陷检测图像。
进一步地,所述步骤S3具体通过以下步骤实现:
对于某一所述中心像素点p,其特征表示区为U p ,遍历该中心像素点p的感受野范围内的所有m个像素点,将除p点外的所有像素点令为邻域像素点q a a=1,2,…,m-1,每个所述邻域像素点q a 的特征表示区为
Figure 676235DEST_PATH_IMAGE001
,每个
Figure 71313DEST_PATH_IMAGE001
内包含n个像素点,U p
Figure 181351DEST_PATH_IMAGE002
的具体表达式如下:
Figure 257761DEST_PATH_IMAGE003
式中,p 1 p 2 、p 3 、...、p n 分别表示U p 中的某一个像素点,q a1 、q a2 、q a3 、...、q an 分别表示
Figure 199172DEST_PATH_IMAGE001
中的某一个像素点;
采用皮尔逊相关系数判定相关度,计算所述U p
Figure 30731DEST_PATH_IMAGE004
之间的相关性,表达式如下:
Figure 424803DEST_PATH_IMAGE005
式中,
Figure 117952DEST_PATH_IMAGE006
代表U p 的期望值,
Figure 648291DEST_PATH_IMAGE007
代表
Figure 668329DEST_PATH_IMAGE004
的期望值,
Figure 549697DEST_PATH_IMAGE008
代表U p 的标准差,
Figure 46538DEST_PATH_IMAGE009
代表
Figure 696962DEST_PATH_IMAGE001
的标准差,
Figure 604744DEST_PATH_IMAGE010
代表U p
Figure 707829DEST_PATH_IMAGE001
之间的协方差,
Figure 8360DEST_PATH_IMAGE011
代表U p
Figure 778870DEST_PATH_IMAGE004
之间的皮尔逊相关系数;
计算m-1个所述皮尔逊相关系数中较大的k个值,得到与所述中心像素点p最相关的前k个所述邻域像素点q i i=1,2,3,…,k。
进一步地,所述步骤S4具体通过以下子步骤实现:
(4.1)依次遍历所述待学习无缺陷样本图片的每一个中心像素点p,其坐标为(x,y),根据所述步骤S3求出每个所述中心像素点p对应的k个所述最相关邻域像素点q i ,通过以下公式分别求出k对pq i 像素对之间的像素差Δ(p,q i ):
Figure 857554DEST_PATH_IMAGE012
(4.2)对每一个所述中心像素点p,拟合出其与k个所述最相关邻域像素点q i 的像素差值的高斯分布为:
Figure 713514DEST_PATH_IMAGE013
式中,
Figure 552157DEST_PATH_IMAGE014
表示均值为μ i 、方差为
Figure 911594DEST_PATH_IMAGE015
的高斯分布;
(4.3)将所述中心像素点p的位置(x,y)、对应的k个所述最相关邻域像素点q i 的位置(x i ,y i )、两个统计参数μ i σ i 作为模型参数保存;当检测完所有中心像素点后,保存所有中心像素点的坐标、该中心像素点对应的k个最相关领域像素点的坐标、高斯分布的均值和标准差,得到多领域预测模型。
进一步地,所述步骤S5具体通过以下操作实现:
输入待测图片,设置误差范围参数C,阈值T,计数器初始值count=0;依次遍历所述待测图片中的每个像素点,依次作为所述中心像素点p,根据p的位置在所述步骤S4得到的多邻域预测模型中,找到所述对应的k个最相关邻域像素点q i 的位置、高斯分布的均值μ i 、标准差σ i ;然后依次求出所述待测图片中这k个pq i 像素对的像素差,并判断该中心像素点p是否在对应的高斯分布之外,即判断是否有如下公式成立:
Figure 426758DEST_PATH_IMAGE016
若上式(11)成立,则count=count+1,i=i+1,循环计算式(11);若不成立,则count保持不变,i=i+1,循环计算式(11);直到i>k,结束循环;
若循环结束后有
Figure 504436DEST_PATH_IMAGE017
成立,则该中心像素点p判定为缺陷像素,将其像素值置1;若不成立,则该中心像素点p判定为非缺陷像素,将其像素值置0;最终得到整幅所述待测图片缺陷检测结果的二值图像。
进一步地,所述步骤S6具体通过以下子步骤实现:
(6.1)定义一个结构元素B,对所述步骤S5得到的缺陷检测结果的二值图像A进行腐蚀和膨胀操作,消除分散孤立的白色噪点和黑色空洞:
A⊙B={z|(B)z⊆A} (12)
A⊕B={z|(B)z⋂A≠∅} (13)
式中,A⊙B代表使用结构元素B对二值图像A进行腐蚀处理,A⊕B代表使用结构元素B对二值图像A进行膨胀处理;
(6.2)遍历所述二值图像中具有相同像素值并且位置相邻的像素组成的连通域,通过观察真实缺陷区域的面积分布特点,设定是否为噪声的判断阈值T1,判定包含像素个数S小于阈值T1的连通域为小面积噪声,将其内包含的像素值置0;反之,则判定所述连通域为真实缺陷,其内包含的像素值保留,具体表达式如下:
Figure 146770DEST_PATH_IMAGE018
式中,p S 为连通域中的像素值。
本发明的有益效果是:
(1)本发明仅需要一张无缺陷图片进行训练建模,解决了深度学习中大量缺陷样本获取困难、样本标注费时费力的问题,提高了缺陷检测效率,降低了检测成本。
(2)本发明采用像素级的缺陷检测,保证了较高的缺陷检测精度。
附图说明
图1是本发明基于多邻域预测模型的磁芯缺陷检测方法的流程图。
图2是本发明步骤S5的流程图。
图3是本发明实施例1中待学习和待检测的图片,其中(a)为无缺陷的待学习图片,(b)为有掉皮缺陷的磁芯图片,(c)为有粘底缺陷的磁芯图片,(d)为有杂质缺陷的磁芯图片。
图4是本发明实施例1中对无缺陷的待学习图片中的其中一个像素点设置特征表示区和感受野的示意图。
图5是本发明实施例1中对图3中对应图片的缺陷检测结果,其中(a)为无缺陷磁芯图片的缺陷检测结果,(b)为有掉皮缺陷的磁芯图片的缺陷检测结果,(c)为有粘底缺陷的磁芯图片的缺陷检测结果,(d)为有杂质缺陷的磁芯图片的缺陷检测结果。
图6是本发明实施例2中对布匹、太阳能电池板、皮革的缺陷检测结果,其中(a1)为布匹缺陷检测的待学习图片,(a2)为布匹缺陷检测的待测图片,(a3)为布匹缺陷检测的检测结果;(b1)为太阳能电池板缺陷检测的待学习图片,(b2)为太阳能电池板缺陷检测的待测图片,(b3)为太阳能电池板缺陷检测的检测结果;(c1)为皮革缺陷检测的待学习图片,(c2)为皮革缺陷检测的待测图片,(c3)为皮革缺陷检测的检测结果。
具体实施方式
下面根据附图和优选实施例详细描述本发明,本发明的目的和效果将变得更加明白,以下结合附图和实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
如图1所示,本发明为一种基于多邻域预测模型的磁芯缺陷检测方法,具体包括以下步骤:
S1:对待学习无缺陷样本图片中的每个像素点p设置特征表示区U p 。特征表示区U p 为以某一个像素点p为中心像素点,包含n个像素点的微小区域,n为正整数,表达式为:
Figure 626292DEST_PATH_IMAGE019
S2:确定每个像素点p的感受野区域。将该像素点p作为该感受野区域的中心像素点,感受野为以中心像素点p为中心,包含m个像素点且包含诸多相似特征的较大区域,其中m为正整数,且m>n
S3:对每个中心像素点p,在其感受野范围内,找出k个与其相关度最高的邻域像素点q
对于某一中心像素点p,其特征表示区为U p ,遍历该像素点p的感受野范围内的所有m个像素点,将除p点外的所有像素点记为邻域像素点q a a=1,2,…,m-1每个邻域像素点q a 的特征表示区为
Figure 312358DEST_PATH_IMAGE001
,将其作为候选邻域,其中包含n个像素点,表达式为:
Figure 877331DEST_PATH_IMAGE020
式中,p 1 p 2 、p 3 、...、p n 分别表示U p 中的某一个像素点,q a1 、q a2 、q a3 、...、q an 分别表示
Figure 57777DEST_PATH_IMAGE001
中的某一个像素点。
本发明采用皮尔逊相关系数判定相关度,皮尔逊相关系数是以两变量与各自平均值的离差为基础,通过两个离差相乘来反映两变量之间相关程度,皮尔逊相关系数越大,两变量的相关程度越高;皮尔逊相关系数越小,两变量的相关程度越低。皮尔逊相关系数ω(X,Y)的表达式如下:
Figure 391806DEST_PATH_IMAGE021
式中,X为变量一,Y为变量二,Cov(X,Y)为X与Y协方差,Var[X]为X的方差,Var[Y]为Y的方差。
现将U p 看作式(3)中的X,m-1个候选邻域
Figure 514352DEST_PATH_IMAGE022
分别看作式(3)中的Y,通过以下公式求得皮尔逊相关系数来度量U p
Figure 566621DEST_PATH_IMAGE001
之间的相关性:
Figure 19599DEST_PATH_IMAGE023
Figure 739294DEST_PATH_IMAGE024
式中,
Figure 767161DEST_PATH_IMAGE006
代表U p 的期望值,
Figure 306727DEST_PATH_IMAGE007
代表
Figure 828975DEST_PATH_IMAGE004
的期望值,
Figure 137597DEST_PATH_IMAGE008
代表U p 的标准差,
Figure 336366DEST_PATH_IMAGE009
代表
Figure 363228DEST_PATH_IMAGE001
的标准差,
Figure 689167DEST_PATH_IMAGE010
代表U p
Figure 852295DEST_PATH_IMAGE001
之间的协方差,
Figure 972698DEST_PATH_IMAGE011
代表U p
Figure 741983DEST_PATH_IMAGE004
之间的皮尔逊相关系数;
计算m-1个皮尔逊相关系数中较大的k个值,得到与中心像素点p最相关的前k个邻域像素点q i i=1,2,3,…,k:将计算得到的m-1个皮尔逊相关系数按降序排列,再通过如下表达式即可得到与所述中心像素点p最相关的前k个所述邻域像素点:
Figure 606034DEST_PATH_IMAGE025
式中,Uq为降序排列的皮尔逊相关系数最大的k个值对应的Uqi的集合,
S4:建立每个中心像素点p与其k个最相关邻域像素点qi之间的多邻域映射关系,得到多邻域预测模型。
具体通过以下子步骤实现:
(4.1)依次遍历待学习无缺陷样本图片的每一个中心像素点p,其坐标为(x,y),根据步骤S3求出每个中心像素点p对应的k个最相关邻域像素点q i ,通过以下公式分别求出k对pq i 像素对之间的像素差Δ(p,q i ):
Figure 623668DEST_PATH_IMAGE026
(4.2)对每一个中心像素点p,拟合出其与k个最相关邻域像素点q i 的像素差值高斯分布:
Figure 180551DEST_PATH_IMAGE027
式中,N(μ i i )表示均值为μ i 、方差为σ i 2 的高斯分布;
(4.3)将中心像素点p的位置(x,y)、对应的k个最相关邻域像素点q i 的位置(x i , y i )、两个统计参数μ i σ i 作为模型参数存入存储表格中的一行。当待学习无缺陷样本图片的每一个中心像素点p对应的数据都存入存储表格后,即完成多邻域预测模型的构建。
S5:调用多邻域预测模型对待测图片进行缺陷检测,得到待测图片缺陷检测结果的二值图像。
如图2所示,输入待测图片,设置误差范围参数C,阈值T,计数器初始值count=0,依次遍历待测图片中的每个像素点,依次作为中心像素点p,根据p的位置在步骤S4得到的多邻域预测模型中,找到对应的k个最相关邻域像素点q i 的位置、高斯分布的均值μ i 、方差
σi 2;然后依次求出待测图片中这k个pq i 像素对的像素差,并判断该中心像素点p是否在对应的高斯分布之外,即判断是否有如下公式成立:
Figure 165694DEST_PATH_IMAGE028
若上式(13)成立,则count=count+1,i=i+1,循环计算式(13);若不成立,则count保持不变,i=i+1,循环计算式(13);直到i>k,结束循环。
若循环结束后有
Figure 567856DEST_PATH_IMAGE029
成立,则该中心像素点p判定为缺陷像素;若不成立,则该中心像素点p判定为非缺陷像素。对于缺陷像素,将其像素值置1;对于非缺陷像素,将其像素值置0。最终得到整幅待测图片缺陷检测结果的二值图像。
S6:对检测结果进行去噪优化,得到最终缺陷检测图像。
去噪优化的目的是消除步骤S5得到的缺陷检测结果的二值图像中存在的噪声,这些噪声是分散在二值图像各处的、由于非缺陷像素被误检为缺陷像素,而产生的白色噪点,以及少量缺陷像素被误检为非缺陷像素,使得检测出的缺陷区域出现的黑色空洞,根据这些噪声区别于真实缺陷的分散稀疏的特点,采用去除小连通域以及图像形态学操作中的腐蚀膨胀操作,将它们全部去除,具体操作如下:
(6.1)定义一个结构元素B,然后根据以下公式对步骤S5得到的缺陷检测结果的二值图像A进行腐蚀和膨胀操作,消除分散孤立的白色噪点和黑色空洞:
A⊙B={z|(B)z⊆A} (14)
A⊕B={z|(B)z⋂A≠∅} (15)
式中,A⊙B代表使用结构元素B对二值图像A进行腐蚀处理,A⊕B代表使用结构元素B对二值图像A进行膨胀处理。
(6.2)在(6.1)的基础上遍历二值图像中具有相同像素值并且位置相邻的像素组成的连通域,通过观察真实缺陷区域的面积分布特点,设定阈值T1,判定包含像素个数S小于阈值T1的连通域为小面积噪声,将其内包含的像素值置0;包含像素个数S大于阈值T1的连通域为真实缺陷,其内包含的像素值保留,具体表达式如下:
Figure 439997DEST_PATH_IMAGE030
式中,S为连通域所包含的像素个数,T1为判断是否为噪声的设定阈值,ps为连通域中的像素值。
下面通过几个具体实施例验证本发明基于多邻域预测模型的磁芯缺陷检测方法的检测效果。
实施例1
基于多邻域预测模型的磁芯缺陷检测方法的具体步骤如下:
S1:采用如图3中的(a)所示的磁芯图片作为待学习无缺陷样本,该图片的大小为296×296,对图片上的每个像素点p,其特征表示区U p 设置为以该像素点p为中心的包含9个像素点的矩形区域。
S2:对图片上的每个像素点p,其感受野为以该像素点p为中心的101×101的范围区域。以某一像素点p为例,得到其感受野和特征表示区U p 之间的关系如图4所示。
S3:对每个中心像素点p,选取其感受野范围内的10个相关度最高的邻域像素点q i i=1,2,3,…,10,得到该中心像素点p的最相关邻域。为了避免所求得的最相关邻域过于聚集在某一区域,本实例设置各邻域之间的距离为10,此时每一个中心像素点p在其感受野范围内的候选邻域个数为100,在这100个候选邻域中找到10个相关度最高的候选领域,这10个相关度最高的候选领域的中心像素点即为最相关邻域像素点。
S4:建立每个中心像素点p与其10个最相关邻域像素点q i 之间的多邻域映射关系,并将所有中心像素点p的坐标、该中心像素点p对应的10个最相关领域像素点的坐标、高斯分布的均值和标准差保存,得到多邻域预测模型。
S5:输入如图3中的(b)、(c)、(d)所示的待测图片,设置误差范围参数C=2,阈值T=0.5,调用多邻域预测模型数据对待测图片进行缺陷检测,得到待测图片缺陷检测结果的二值图像。
S6:对步骤S5得到的二值图像进行去噪优化,腐蚀膨胀处理中,所用结构元素B为3×3的十字形矩形;去除小连通域的处理过程中,分析真实缺陷的分布特点,设置阈值T1=30,得到最终的缺陷检测结果如图5所示。
实施例2
为进一步验证本发明对图案纹理比较规则、相对简单的产品表面图片的缺陷检测具有通用性、有效性,本实施例分别在布匹、太阳能电池板、皮革这几个数据集上,使用与实施例1相同的缺陷检测方法进行了测试,测试结果如图6所示。结果表明,除对磁芯进行缺陷检测外,本发明同样适用于产品表面图案纹理比较规则、相对简单的领域,能有效检测出产品表面缺陷。
综上所述,本发明通过检测产品表面图片非缺陷像素对之间、非缺陷与缺陷像素对之间的差异,分析学习无缺陷样本中每个像素点与其邻域的像素点的差异特征,并建立稳定的多邻域映射关系,得到多邻域预测模型;在检测时,只需要判断目标像素点及其邻域是否符合这种多领域映射关系,就可以实现样本的缺陷检测与识别。本发明可以快速有效地检测出产品表面缺陷,且对于表面图案纹理比较规则、相对简单的产品具有较高的适用性。
本领域普通技术人员可以理解,以上所述仅为发明的优选实例而已,并不用于限制发明,尽管参照前述实例对发明进行了详细的说明,对于本领域的技术人员来说,其依然可以对前述各实例记载的技术方案进行修改,或者对其中部分技术特征进行等同替换。凡在发明的精神和原则之内,所做的修改、等同替换等均应包含在发明的保护范围之内。

Claims (5)

1.一种基于多邻域预测模型的磁芯缺陷检测方法,其特征在于,包括以下步骤:
S1:对待学习无缺陷样本图片中的每个像素点设置特征表示区,每个像素点的特征表示区均为以该像素点为中心,包含n个像素点的区域,n为正整数;
S2:确定每个像素点的感受野区域,并将该像素点作为所述感受野区域的中心像素点,所述感受野区域包含m个像素点,m为正整数,且m>n;
S3:对于某一所述中心像素点p,其特征表示区为Up,遍历该中心像素点p的感受野范围内的所有m个像素点,将除p点外的所有像素点令为邻域像素点qa,a=1,2,...,m-1,每个所述邻域像素点qa的特征表示区为
Figure FDA0004181136390000013
每个
Figure FDA0004181136390000012
内包含n个像素点;计算所述Up
Figure FDA0004181136390000011
之间的相关度,从而找出k个与所述中心像素点相关度最高的邻域像素点,作为最相关邻域像素点,所述k小于m;
S4:建立每个所述中心像素点与其k个所述最相关邻域像素点之间的多邻域映射关系,得到多邻域预测模型;
所述S4通过以下子步骤实现:依次遍历所述待学习无缺陷样本图片的每一个中心像素点,计算该中心像素点和每一个与其对应的k个最相关领域像素点之间的像素差,并将所述像素差进行拟合,得到其高斯分布的均值和方差;将所述中心像素点的坐标、所述中心像素点对应的k个最相关领域像素点的坐标、高斯分布的均值和标准差分别保存;当检测完所有中心像素点后,保存所有中心像素点的坐标、该中心像素点对应的k个最相关领域像素点的坐标、高斯分布的均值和标准差,得到多领域预测模型;
S5:调用所述多邻域预测模型数据对待测图片进行缺陷检测,得到所述待测图片缺陷检测结果的二值图像;
S6:对S5得到的缺陷检测结果的二值图像进行去噪优化,得到最终缺陷检测图像。
2.根据权利要求1所述的基于多邻域预测模型的磁芯缺陷检测方法,其特征在于,所述步骤S3具体通过以下步骤实现:
Up
Figure FDA0004181136390000014
的具体表达式如下:
Up={p1,p2,p3,…,pn} (1)
Figure FDA0004181136390000015
式中,p1、p2、p3、...、pn分别表示Up中的某一个像素点,qa1、qa2、qa3、...、qan分别表示
Figure FDA00041811363900000218
中的某一个像素点;
采用皮尔逊相关系数判定相关度,计算所述Up
Figure FDA00041811363900000219
之间的相关度,表达式如下:
Figure FDA0004181136390000021
Figure FDA0004181136390000022
Figure FDA0004181136390000023
Figure FDA0004181136390000024
Figure FDA0004181136390000025
Figure FDA0004181136390000026
式中,
Figure FDA0004181136390000027
代表Up的期望值,
Figure FDA0004181136390000028
代表
Figure FDA0004181136390000029
的期望值,
Figure FDA00041811363900000210
代表Up的标准差,
Figure FDA00041811363900000211
代表
Figure FDA00041811363900000212
的标准差,
Figure FDA00041811363900000213
代表Up
Figure FDA00041811363900000214
之间的协方差,
Figure FDA00041811363900000215
代表Up
Figure FDA00041811363900000216
之间的皮尔逊相关系数;
计算m-1个所述皮尔逊相关系数中较大的k个值,得到与所述中心像素点p最相关的前k个所述邻域像素点qi,i=1,2,3,...,k。
3.根据权利要求2所述的基于多邻域预测模型的磁芯缺陷检测方法,其特征在于,所述步骤S4具体通过以下子步骤实现:
(4.1)依次遍历所述待学习无缺陷样本图片的每一个中心像素点p,其坐标为(x,y),根据所述步骤S3求出每个所述中心像素点p对应的k个所述最相关邻域像素点qi,通过以下公式分别求出k对p、qi像素对之间的像素差Δ(p,qi):
Figure FDA00041811363900000217
(4.2)对每一个所述中心像素点p,拟合出其与k个所述最相关邻域像素点qi的像素差值的高斯分布为:
Δ(p,qi)~N(μi,σi) (10)
式中,N(μi,σi)表示均值为μi、方差为σi 2的高斯分布;
(4.3)将所述中心像素点p的位置(x,y)、对应的k个所述最相关邻域像素点qi的位置(xi,yi)、两个统计参数μi、σi作为模型参数保存;当检测完所有中心像素点后,保存所有中心像素点的坐标、该中心像素点对应的k个最相关领域像素点的坐标、高斯分布的均值和标准差,得到多领域预测模型。
4.根据权利要求3所述的基于多邻域预测模型的磁芯缺陷检测方法,其特征在于,所述步骤S5具体通过以下操作实现:
输入待测图片,设置误差范围参数C,阈值T,计数器初始值count=0;依次遍历所述待测图片中的每个像素点,依次作为所述中心像素点p,根据p的位置在所述步骤S4得到的多邻域预测模型中,找到所述对应的k个最相关邻域像素点qi的位置、高斯分布的均值μi、标准差σi;然后依次求出所述待测图片中这k个p、qi像素对的像素差,并判断该中心像素点p是否在对应的高斯分布之外,即判断是否有如下公式成立:
|Δ(p,qi)-μi|>(C*σi)i=1,2,…,k (11)
若上式(11)成立,则count=count+1,i=i+1,循环计算式(11);若不成立,则count保持不变,i=i+1,循环计算式(11);直到i>k,结束循环;
若循环结束后有
Figure FDA0004181136390000031
成立,则该中心像素点p判定为缺陷像素,将其像素值置1;若不成立,则该中心像素点p判定为非缺陷像素,将其像素值置0;最终得到整幅所述待测图片缺陷检测结果的二值图像。
5.根据权利要求4所述的基于多邻域预测模型的磁芯缺陷检测方法,其特征在于,所述步骤S6具体通过以下子步骤实现:
(6.1)定义一个结构元素B,对所述步骤S5得到的缺陷检测结果的二值图像A进行腐蚀和膨胀操作,消除分散孤立的白色噪点和黑色空洞:
Figure FDA0004181136390000032
Figure FDA0004181136390000033
式中,A⊙B代表使用结构元素B对二值图像A进行腐蚀处理,
Figure FDA0004181136390000034
代表使用结构元素B对二值图像A进行膨胀处理;
(6.2)遍历所述二值图像中具有相同像素值并且位置相邻的像素组成的连通域,通过观察真实缺陷区域的面积分布特点,设定是否为噪声的判断阈值T1,判定包含像素个数S小于阈值T1的连通域为小面积噪声,将其内包含的像素值置0;反之,则判定所述连通域为真实缺陷,其内包含的像素值保留,具体表达式如下:
Figure FDA0004181136390000041
式中,pS为连通域中的像素值。
CN202310011002.XA 2023-01-05 2023-01-05 基于多邻域预测模型的磁芯缺陷检测方法 Active CN115908408B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310011002.XA CN115908408B (zh) 2023-01-05 2023-01-05 基于多邻域预测模型的磁芯缺陷检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310011002.XA CN115908408B (zh) 2023-01-05 2023-01-05 基于多邻域预测模型的磁芯缺陷检测方法

Publications (2)

Publication Number Publication Date
CN115908408A CN115908408A (zh) 2023-04-04
CN115908408B true CN115908408B (zh) 2023-06-02

Family

ID=86478972

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310011002.XA Active CN115908408B (zh) 2023-01-05 2023-01-05 基于多邻域预测模型的磁芯缺陷检测方法

Country Status (1)

Country Link
CN (1) CN115908408B (zh)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2017090444A (ja) * 2015-11-09 2017-05-25 株式会社リコー 検査装置、検査方法及びプログラム
CN109934154A (zh) * 2019-03-08 2019-06-25 北京科技大学 一种遥感影像变化检测方法及检测装置

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9661299B2 (en) * 2011-06-30 2017-05-23 Thomson Licensing Outlier detection for colour mapping
US8885077B2 (en) * 2013-01-07 2014-11-11 Altasens, Inc. Dynamic, local edge preserving defect pixel correction for image sensors
CN104990925B (zh) * 2015-06-23 2017-12-19 泉州装备制造研究所 一种基于梯度多阈值优化缺陷检测方法
CN110020676A (zh) * 2019-03-18 2019-07-16 华南理工大学 基于多感受野深度特征的文本检测方法、系统、设备及介质
US11494584B2 (en) * 2021-01-12 2022-11-08 Disney Enterprises, Inc. Automated prediction of pixel error noticeability

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2017090444A (ja) * 2015-11-09 2017-05-25 株式会社リコー 検査装置、検査方法及びプログラム
CN109934154A (zh) * 2019-03-08 2019-06-25 北京科技大学 一种遥感影像变化检测方法及检测装置

Also Published As

Publication number Publication date
CN115908408A (zh) 2023-04-04

Similar Documents

Publication Publication Date Title
WO2021143343A1 (zh) 一种产品质量检测方法及装置
CN111383209A (zh) 一种基于全卷积自编码器网络的无监督瑕疵检测方法
CN110781913B (zh) 一种拉链布带缺陷检测方法
CN111047655A (zh) 基于卷积神经网络的高清摄像机布料疵点检测方法
CN111598801B (zh) 一种用于微弱Mura缺陷的识别方法
CN114723708A (zh) 一种基于无监督图像分割的工艺品外观缺陷检测方法
CN113313704B (zh) 一种集成电路引脚检测方法、装置、终端及存储介质
CN116152242B (zh) 一种篮球用天然皮革缺陷可视化检测系统
CN116883408B (zh) 基于人工智能的积算仪壳体缺陷检测方法
CN115170529A (zh) 基于注意力机制的多尺度微小瑕疵检测方法
CN115082477B (zh) 一种基于去反光效果的半导体晶圆加工质量检测方法
CN116342583A (zh) 一种用于纺纱生产加工的抗起球性能检测方法
CN116778269A (zh) 一种基于自编码器重构产品表面缺陷检测模型构建方法
CN115829984A (zh) 一种棉织品丝光处理缺陷识别方法
CN117611551B (zh) 基于物联网的显示屏异常检测方法及系统
KR101782364B1 (ko) 케이평균군집화를 통한 학습기반의 비전검사 방법
CN113421223B (zh) 基于深度学习和高斯混合的工业产品表面缺陷检测方法
CN115908408B (zh) 基于多邻域预测模型的磁芯缺陷检测方法
CN114820597B (zh) 一种基于人工智能的冶炼品缺陷检测方法、装置及系统
CN116503340A (zh) 一种micro oled面板缺陷检测方法、装置以及设备
CN116363584A (zh) 一种基于机器视觉的船舶液体管道泄漏监测方法
CN113470015B (zh) 一种基于图像处理的水体晃动检测分析方法及系统
CN116542963A (zh) 一种基于机器学习的浮法玻璃缺陷检测系统及检测方法
CN115546141A (zh) 一种基于多维度度量的小样本Mini LED缺陷检测方法及系统
CN114898148B (zh) 一种基于深度学习的鸡蛋污损检测方法及系统

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant