CN114155197A - 一种元器件缺陷检测的方法及装置 - Google Patents
一种元器件缺陷检测的方法及装置 Download PDFInfo
- Publication number
- CN114155197A CN114155197A CN202111314595.4A CN202111314595A CN114155197A CN 114155197 A CN114155197 A CN 114155197A CN 202111314595 A CN202111314595 A CN 202111314595A CN 114155197 A CN114155197 A CN 114155197A
- Authority
- CN
- China
- Prior art keywords
- image
- target component
- algorithm
- ray
- module
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 68
- 230000007547 defect Effects 0.000 title claims abstract description 36
- 238000005286 illumination Methods 0.000 claims abstract description 36
- 238000001514 detection method Methods 0.000 claims abstract description 33
- 238000012545 processing Methods 0.000 claims abstract description 28
- 230000009467 reduction Effects 0.000 claims abstract description 25
- 238000003709 image segmentation Methods 0.000 claims abstract description 15
- 238000012937 correction Methods 0.000 claims description 12
- 230000011218 segmentation Effects 0.000 claims description 12
- 230000008030 elimination Effects 0.000 claims description 5
- 238000003379 elimination reaction Methods 0.000 claims description 5
- 230000008569 process Effects 0.000 description 10
- 238000010586 diagram Methods 0.000 description 7
- 230000000694 effects Effects 0.000 description 6
- 230000003044 adaptive effect Effects 0.000 description 5
- 230000008859 change Effects 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 239000000203 mixture Substances 0.000 description 3
- 238000000354 decomposition reaction Methods 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 238000007781 pre-processing Methods 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 208000003464 asthenopia Diseases 0.000 description 1
- 230000000903 blocking effect Effects 0.000 description 1
- 230000007797 corrosion Effects 0.000 description 1
- 238000005260 corrosion Methods 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000009413 insulation Methods 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 238000002372 labelling Methods 0.000 description 1
- 230000004807 localization Effects 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 238000001465 metallisation Methods 0.000 description 1
- 229910052755 nonmetal Inorganic materials 0.000 description 1
- 239000005022 packaging material Substances 0.000 description 1
- 238000004806 packaging method and process Methods 0.000 description 1
- 238000003908 quality control method Methods 0.000 description 1
- 238000007789 sealing Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0004—Industrial image inspection
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
- G01N23/02—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
- G01N23/04—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/22—Matching criteria, e.g. proximity measures
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/10—Image enhancement or restoration using non-spatial domain filtering
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/194—Segmentation; Edge detection involving foreground-background segmentation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/70—Determining position or orientation of objects or cameras
- G06T7/73—Determining position or orientation of objects or cameras using feature-based methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10116—X-ray image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20048—Transform domain processing
- G06T2207/20064—Wavelet transform [DWT]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30108—Industrial image inspection
- G06T2207/30148—Semiconductor; IC; Wafer
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Data Mining & Analysis (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Artificial Intelligence (AREA)
- Quality & Reliability (AREA)
- Evolutionary Computation (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Image Processing (AREA)
Abstract
本发明公开了一种元器件缺陷检测方法及装置,所述方法包括:获取目标元器件的X射线图像,并对所述X射线图像进行降噪处理,得到第一图像;消除所述第一图像的非均匀光照背景,得到第二图像;基于模板匹配定位所述第二图像中待检测的目标元器件图像;通过预设图像分割算法,将目标元器件图像从所述第二图像中分割出来;基于局部信息熵结合区域生长算法,对所述X射线图像中的多余物进行检测,以检测到所述目标元器件的缺陷。本申请实施例提供的元器件缺陷检测方法,高效且可节省大量人力资源,所得检测结果准确、可靠。
Description
技术领域
本发明图像处理技术领域,尤其涉及一种元器件缺陷检测方法及装置。
背景技术
在元器件封装过程中,受封装材料、环境条件和工艺参数等因素的影响,外来物质很容易进入封装元器件内部。这些外来物质的存在会导致元器件可靠性的降低,如金属或半导体多余物会导致电路不同部位之间的短路,尺寸和重量比较大的非金属多余物在振动情况下可能会导致键合丝移位或断裂,移位的键合丝可能和其它部位短路;纤维状的非金属多余物会吸潮,并可能导致器件内部不同部位绝缘降低;带有腐蚀性成份的多余物可能导致电路内部发生金属化腐蚀。因此,准确检测出元器件内部的多余物等缺陷对于元器件的质量控制具有重要意义。
现有的检测元器件缺陷的流程如下:人工将待检元器件放置在X射线检测设备的检测工作台上,然后操作设备对准待检器件,再对器件进行拍照,拍照后将图片传输到计算机上。再然后,人工浏览X射线所成像的一张张图片,对元器件内部的多余物进行检测。全部工作过程中自动化程度低,基本依赖人工操作和人为检测判断。在大批量的元器件质量检测过程中,不仅会耗费检测人员大量的时间,还存在因视觉疲劳或经验判断不足造成的人为判断失误的问题,元器件质量难以有效保证,给元器件在现场的安全运行带来风险与隐患。
发明内容
本发明实施例提供了一种元器件缺陷检测方法,以解决现有技术中存在的消耗人力资源多、检测耗时长、且检测结果准确度低的问题。
为了解决上述技术问题,本发明公开了一种元器件缺陷检测方法,所述方法包括:获取目标元器件的X射线图像,并对所述X射线图像进行降噪处理,得到第一图像;消除所述第一图像的非均匀光照背景,得到第二图像;基于模板匹配定位所述第二图像中待检测的目标元器件图像;通过预设图像分割算法,将目标元器件图像从所述第二图像中分割出来;基于局部信息熵结合区域生长算法,对所述X射线图像中的多余物进行检测,以检测到所述目标元器件的缺陷。
可选的,所述获取目标元器件的X射线图像,并对所述X射线图像进行降噪处理,得到第一图像的步骤,包括:获取目标元器件的X射线图像;基于小波阈值算法,对所述X射线图像进行降噪处理,得到第一图像。
可选的,消除所述第一图像的非均匀光照背景,得到第二图像的步骤,包括:基于二维伽马函数的自适应亮度校正方法,对所述第一图像进行处理,得到消除非均匀光照背景的第二图像。
可选的,基于模板匹配定位所述第二图像中待检测的目标元器件图像的步骤,包括:采用归一化相关系数匹配法,将模板与所述第二图像进行比对,定位到第二图像中的目标元器件图像。
可选的,所述预设图像分割算法为Grabcut算法。
为了解决上述技术问题,本发明还公开了一种元器件缺陷检测装置,其中,所述装置包括:降噪模块,用于获取目标元器件的X射线图像,并对所述X射线图像进行降噪处理,得到第一图像;消除模块,用于消除所述第一图像的非均匀光照背景,得到第二图像;定位模块,用于基于模板匹配定位所述第二图像中待检测的目标元器件图像;分割模块,用于通过预设图像分割算法,将目标元器件图像从所述第二图像中分割出来;异物检测模块,用于基于局部信息熵结合区域生长算法,对所述X射线图像中的多余物进行检测,以检测到所述目标元器件的缺陷。
可选的,所述降噪模块包括:第一子模块,用于获取目标元器件的X射线图像;第二子模块,用于基于小波阈值算法,对所述X射线图像进行降噪处理,得到第一图像。
可选的,所述消除模块具体用于:基于二维伽马函数的自适应亮度校正方法,对所述第一图像进行处理,得到消除非均匀光照背景的第二图像。
可选的,所述定位模块具体用于:采用归一化相关系数匹配法,将模板与所述第二图像进行比对,定位到第二图像中的目标元器件图像。
可选的,所述预设图像分割算法为Grabcut算法。
本发明实施例提供的元器件缺陷检测方案,获取目标元器件的X射线图像,并对X射线图像进行降噪处理,得到第一图像;消除第一图像的非均匀光照背景,得到第二图像;基于模板匹配定位第二图像中待检测的目标元器件图像;通过预设图像分割算法,将目标元器件图像从第二图像中分割出来;基于局部信息熵结合区域生长算法,对X射线图像中的多余物进行检测,以检测到目标元器件的缺陷,该种元器件缺陷检测方法,高效且可节省大量人力资源,所得检测结果准确、可靠。
附图说明
图1为本发明实施例所提供的元器件缺陷检测方法的步骤流程图;
图2为本发明实施例所提供的区域生成算法流程示意图;
图3为本发明实施例所提供的多余物检测结果示意图;
图4为本发明实施例所提供的元器件缺陷检测装置的结构框图。
具体实施方式
下面根据具体的实施例,结合附图针对本发明进行详细说明。应当理解,此处所述的具体实施例仅用于解释本发明,并不用于限定本发明。
如图1所示,本发明实施例的元器件缺陷检测方法的步骤流程图,该方法包括如下步骤:
步骤101:获取目标元器件的X射线图像,并对X射线图像进行降噪处理,得到第一图像。
一种可选的对X射线图像进行降噪处理的方式可以为:在获取目标元器件的X射线图像后,基于小波阈值算法,对X射线图像进行降噪处理,得到第一图像。
由于X射线图像混合噪声明显且缺陷对比度低,采用常规的降噪算法在滤除图像噪声的同时,也会平滑图像的边缘,影响缺陷的检测。小波阀值去噪的基本思想是将信号通过小波变换(采用Mallat算法)后,信号产生的小波系数含有信号的重要信息,将信号经小波分解后小波系数较大,噪声的小波系数较小,并且噪声的小波系数要小于信号的小波系数,通过选取一个合适的阀值,大于阀值的小波系数被认为是有信号产生的,应予以保留,小于阀值的则认为是噪声产生的,置为零从而达到去噪的目的。
步骤102:消除第一图像的非均匀光照背景,得到第二图像。
一种可行性的方式为:基于二维伽马函数的自适应亮度校正方法,对所述第一图像进行处理,得到消除非均匀光照背景的第二图像。
考虑到X光设备采集图像时,在提取出场景的光照分量后,就可以根据光照分量的分布特性构造光照不均匀校正函数,对光照不均匀图像进行校正处理,降低光照过强区域的亮度值,提高光照过低区域的亮度值。为了实现上述目标,本项目提出了一种基于二维伽马函数的自适应亮度校正方法,利用图像的光照分量的分布特性自适应地调整二维伽马函数的参数,实现提高光照不均匀图像整体质量的目的。
步骤103:基于模板匹配定位第二图像中待检测的目标元器件图像。
一种可行性的方式为:采用归一化相关系数匹配法,将模板与所述第二图像进行比对,定位到第二图像中的目标元器件图像。
在获取X射线图像时,由于实际检测设备工作机理的原因,在拍摄目标待检测元器件的时候,会拍摄到其他元器件结构的一部分,会对后续的检测造成影响。本申请实施例提出基于模板匹配定位到中心的待检测的目标元器件,得到该元器件的位置并用于后续的分割提取。对于模板匹配,采用归一化相关系数匹配法,该匹配方法精度较高,可准确定位到目标元器件。
步骤104:通过预设图像分割算法,将目标元器件图像从第二图像中分割出来。
预设图像分割算法可以为Grabcut算法,Grabcut算法通过简单的交互信息就能获取图像当中人们感兴趣的区域,它是在Graphcut的基础上主要做了一下改进:(1)在图像的前景和背景区域,用GMM(高斯混合模型)替代了灰度直方图模型;(2)用迭代的方式取代了一次最小估计求解最优值,提高了分割结果的可靠性;(3)用户只需在图像周围画一个简单的矩形框就能实现交互操作,减少了用户的工作量。综上,结合S3模板匹配得到目标元器件的矩形定位框的位置,输入到Grabcut的算法中,可以将目标元器件准确的分割出来。
步骤105:基于局部信息熵结合区域生长算法,对X射线图像中的多余物进行检测,以检测到目标元器件的缺陷。
图像的信息熵表达的是图像整体的信息量,无法详细度量图像某个区域的信息量。多余物目标占图像的整体部分较少,小目标引起的灰度值变化通常不会对图像整体部分产生较大的影响,可能会导致检测的时候造成目标的漏检。但是小目标会对局部区域的灰度值产生很大的影响,局部信息熵值会产生较大的变化。因此我们将图像信息熵转换为局部信息熵,图像局部信息熵表达的是图像局部区域灰度的起伏。
本申请实施例提供的元器件缺陷检测方法,获取目标元器件的X射线图像,并对X射线图像进行降噪处理,得到第一图像;消除第一图像的非均匀光照背景,得到第二图像;基于模板匹配定位第二图像中待检测的目标元器件图像;通过预设图像分割算法,将目标元器件图像从第二图像中分割出来;基于局部信息熵结合区域生长算法,对X射线图像中的多余物进行检测,以检测到目标元器件的缺陷,该种元器件缺陷检测方法,高效且可节省大量人力资源,所得检测结果准确、可靠。
下面结合图2~图3对本申请的元器件缺陷检测方法进行说明。该方法包括如下流程:
S1图像预处理阶段
图像预处理阶段包括:(1)基于小波阈值算法的图像降噪处理和(2)基于二维伽马函数的自适应亮度校正两个阶段,下面分别对这两个阶段进行详细说明。
(1)基于小波阈值算法的图像降噪处理
X射线图像混合噪声明显且缺陷对比度低,采用小波阈值算法,其在滤除噪声、平滑图像的基础上,能够做到很好的边缘保持效果。
小波域阈值法去噪处理的过程可简述为:首先把含有噪声的信号在己确定分解层数内对其进行各个尺度上的小波分解,即利用小波变换把自然图像变换到小波域。然后对所处大尺度中低分辨率下的所有小波系数予以保留,同时针对所有尺度中高分辨率下的小波系数,可以根据预先所设定的阈值,将其中幅值小于该阈值的小波系数统统都置为零,而将其中那些幅值大于阈值的小波系数则进行保留或者适当缩小。最终再把所处理后得到的小波系数进行相应的小波逆变换,就可以得原始信号的恢复,这个信号就是去噪后的图像。
在以上过程中,小波基和分解层数的选择,阈值的选取规则,和阈值函数的设计,均为影响最终去噪效果的关键因素。
1)小波基的选择:通常用户希望所选取的小波满足以下条件:正交性、高消失矩、紧支性、对称性或反对称性。本申请实施例中选用Daubechies小波作为小波基,一般简写成dbN,N是小波的阶数。小波函数Ψ(t)和尺度函数中的支撑区为2N-1,Ψ(t)的消失矩为N。dbN小波具有较好的正则性,即该小波作为稀疏基所引入的光滑误差不容易被察觉,使得信号重构过程比较光滑。dbN小波的特点是随着阶次(序列N)的增大消失矩阶数越大,其中消失矩越高光滑性就越好,频域的局部化能力就越强,频带的划分效果越好。
2)阀值的选择:直接影响去噪效果的一个重要因素就是阀值的选取,不同的阀值选取将有不同的去噪效果。目前主要有通用阀值(VisuShrink)、SureShrink阀值、Minimax阀值、BayesShrink阀值等。本方法选用SureShrink阈值,下面介绍该阈值方法。
SureShrink阈值是一种基于Stein的无偏似然估计(Stein’s Unbiased RiskEstimation,SURE)基本思想的自适应阈值方法,因此也称可为Stein无偏风险阈值。此基本思想建立在均方差准则的无偏估计,它是专门针对于软阈值去噪函数得出的结论,是一种软件估计器,这使得SureShrink阈值趋近于理想阈值。首先给定一个的阈值λ,计算出它的似然估计,再通过最小化非似然估计λ,就得到所选的阈值。
首先令xi~N(μi,1),j=1,2,…,k,得到Stein无偏风险估计为:
则SureShrink阈值为:
式中,∧符号表示在两个数之间取小者,xi为小波系数,并且这时噪声标准方差σn=1,否则就要对其进行标准化。
3)阈值函数的选择:在小波域阈值去噪法中,阈值函数就对小波系数进行获取的函数,它表现出来的是对高于和低于阈值的小波系数所采用不同的处理规则以及不同的估计方法。本申请实施例中采用介于常见的硬阈值函数和软阈值函数之间的折中阈值函数,这种阈值函数既使函数有连续性又能保留信号中一些较大的小波系数,使得X光影像去噪具有良好的效果。折中阈值函数其表达式如下:
式中wj.k表示的是含有噪声图像的小波变换系数,阈值用λ来表示,β是处于0和1之间得实数。可以看出它的计算量就要比半软阈值函数要简单的多,所以在小波阈值去噪中是一种比较常用的阈值函数。
基于以上的选择,我们对X光图像进行去噪处理,可以有效地去除X图像上的混合噪声。
(2)基于二维伽马函数的自适应亮度校正:
考虑到X光设备采集图像时,在提取出场景的光照分量后,就可以根据光照分量的分布特性构造光照不均匀校正函数,对光照不均匀图像进行校正处理,降低光照过强区域的亮度值,提高光照过低区域的亮度值。为了实现上述目标,本申请实施例提出了一种基于二维伽马函数的自适应亮度校正方法,利用图像的光照分量的分布特性自适应地调整二维伽马函数的参数,实现提高光照不均匀图像整体质量的目的。对于输入的图像,假设提取出的光照分量为,构造了一种新的二维伽马函数,其表达式如下:
式中:O(x,y)为校正后的输出图像的亮度值;γ为用于亮度增强的指数值,其中包含了图像的光照分量特性;m为光照分量的亮度均值。
S2基于模板匹配定位X射线图像中待检测的元器件图像
在获取X射线图像时,由于实际检测设备工作机理的原因,在拍摄目标待检测元器件的时候,会拍摄到其他元器件结构的一部分,会对后续的检测造成影响。本方法提出基于模板匹配定位到中心的待检测的目标元器件,得到该元器件的位置并用于后续的分割提取。采用模板匹配的方式,从搜索图像f(x,y)的原点开始,考虑一个二维图形F=(fi,j)M×N和一个模板图形T=(ti,j)P×Q,通过模板在图像上滑动,来寻找最佳匹配值γi,j最大的位置。
本申请实施例中结合OpenCV库,调用函数matchTemplate,选择归一化相关系数匹配的方式来实现,将模版对其均值的相对值与图像对其均值的相关值进行匹配,提高了匹配的精度,有效减少了密封圈提取不准确的情况,相关系数匹配计算公式如下:
式中,R(x,y)为输出图像结果,T′(x′,y′)为模板,I(x,y)为输入图像,表示位置信息。
S3结合Grabcut算法对目标元器件进行分割
基于S2获得了目标元器件的定位矩形位置,将得到的定位矩形位置输入到grabcut函数中作为前景区域标注矩形框,Grabcut算法原理如下:
GrabCut算法是在GraphCut算法基础上的改进和扩展。GraphCut算法把待分割图像构建为一个无向网络图G(V,E)。V是所有节点的集合,除了图中每个像素外还有两个终端节点:源点s(代表前景)和汇点t(代表背景)。E代表图的边集,由n-links和t-links组成;n-links是图中每个像素之间的连接,t-links是每个像素节点到终端节点的连接。n-links和t-links分别由不同的计算方式生成权值,从而构建能量函数;然后再用最大流/最小割算法对网络进行分割得到能量函数的最小值从而实现图像的分割。GrabCut算法首先它改变原有GraphCut算法需要手动标记前景区域和背景区域的方法为只需要在前景区域标注矩形框即可,大大降低了用户交互的复杂性;其次它采用高斯混合模型替代GraphCut中的灰度直方图模型来描述前景和背景像素的分布从而实现了对彩色图像的分割;最后,使用迭代更新的方式代替原有一次分割方法的来实现对所构建网络图能量函数的最小化获得更高的分割精度,但由于GMM(高斯混合模型)参数和迭代更新的复杂性导致算法效率降低。GrabCut用来分割的能量函数如下:
E(α,k,θ,z)=U(α,k,θ,z)+V(α,z)
式中:α为透明度系数,为0是背景,为1是前景;k∈{1,2,…,K},K为GMM分量数目,通常取5;θ={πk,uk,∑k}为每个GMM分量对应的比例,均值,协方差;E是能量函数,当E取最小值便得到最佳分割;U为区域项,用来衡量透明度系数α对于像素z的匹配程度,当所有像素被正确分类到前景或背景时取得最小值;V为边界项,反映相邻像素之间的相似性,在图像的边界处取得最小值。
S4基于局部信息熵结合区域生长算法对多余物进行检测
S4包括下述(1)、(2)以及(3)三部分,下面分别对这三部分进行详细说明。
(1)视觉显著性模型
图像的信息熵表达的是图像整体的信息量,无法详细度量图像某个区域的信息量。多余物目标占图像的整体部分较少,小目标引起的灰度值变化通常不会对图像整体部分产生较大的影响,可能会导致检测的时候造成目标的漏检。但是小目标会对局部区域的灰度值产生很大的影响,局部信息熵值会产生较大的变化。因此本申请实施例中将图像信息熵转换为局部信息熵,图像局部信息熵表达的是图像局部区域灰度的起伏。局部信息熵的表达式为:
图像局部信息熵反映图像中灰度信息分布的离散程度,图像局部信息熵值越大,局部信息量越大,图像局部灰度变化越丰富,图像局部就越有可能包含目标等有用信息。相反地,图像局部信息熵值越小,图像局部信息量越小,图像这局部灰度级变化小,包含有用信息的可能性越低。我们通过对图像的信息熵进行分块计算,得到图像的信息熵图。对于得到的信息熵图设定合适的阈值进行分割,表示为:
其中,c为常数。获得基于局部信息熵的显著性图的二值图,确定了X射线图像中显著区域性的范围。
(2)区域生长算法及实现流程
区域生长的基本思想是将具有相似性质的像素集合起来构成区域。首先对每个需要分割的区域找出一个种子像素作为生长的起点,然后将种子像素周围邻域中与种子有相同或相似性质的像素(根据事先确定的生长或相似准则来确定)合并到种子像素所在的区域中。而新的像素继续做种子向四周生长,直到再没有满足条件的像素可以包括进来,一个区域就生长而成了。
区域生长算法是从单个像素出发,按照一定的生长准则,将具有相似性质的像素合并起来构成同一区域。具体做法如下:
1)选出图像中要分割的目标物体内的一点或者说生长起始点。
2)在种子区域的基础上不断将其周围的像素点以一定的规则加入其中。
3)最终将代表该物体的所有像素点结合成一个区域。
在所有图像中,同一个子区域内部像素会有相似的特性,不同子区域之间像素特性差异一般较大。区域生长算法从种子点开始生长,直到所有符合条件的像素点加入为止。区域生长算法原理很简单,其中最为核心的是生长准则,一般可以表示为一个灰度区间:[M-F·σ,M+F·σ]其中M和σ分别表示当前己生长区域的所有像素点灰度值平均值和标准差,F是自定义的系数,如果邻域中像素点的灰度值处于这个区间中则被接纳,否则就被拒绝。
其中,区域生长算法程序流程图如图2所示。如图2所示,整幅图像中所有的像素点分为两类:己检测过的像素点和未被检测过的像素点。首先选择目标区域中一个或多个点(种子区域)加入到已生长区域,并以此作为生长起点,计算已生长区域中所有像素点灰度值平均值和标准差,接着判断目前区域的邻域中是否有符合生长准则的像素点,如果存在就将其划分入已生长区域,这样就完成一次迭代。第一次迭代完成后,重复上述步骤,直到没有满足条件的邻域像素点划分入己生长区域为止,区域生长算法结束。目前来说,由于图像的多样性和复杂性,生长准则没有一定的标准。
区域生长的关键是选择合适的生长或相似准则,大部分区域生长准则会使用图像的局部性质。生长准则可以根据不同原理制定,而使用不同的生长准则会影响区域生长的过程。常用的生长准则和方法有三种:基于区域灰度差的、基于区域内灰度分布统计性质的和基于区域形状的。
本方法采用基于区域内灰度分布特性准则,考虑以灰度分布相似性作为生长准则来决定区域的合并,具体步骤为:
步骤一,把图像分成互不重叠的小区域;
步骤二:比较邻接区域的累积灰度直方图,根据灰度分布的相似性进行区域合并;
步骤三:设定终止准则,通过反复进行步骤二中的操作将各个区域依次合并直到满足终止准则,生长过程结束。
(3)多余物检测结果与分析
其中,多余物检测结果图如图3所示,图3中的(a)图为局部熵结果图,(b)图为区域生长结果图,所检测出的多余物则为元器件存在的缺陷。
图4为本发明实施例提供的一种元器件缺陷检测装置的结构框图。
如图4所示,本发明实施例的元器件缺陷检测装置包括如下模块:
降噪模块401,用于获取目标元器件的X射线图像,并对所述X射线图像进行降噪处理,得到第一图像;
消除模块402,用于消除所述第一图像的非均匀光照背景,得到第二图像;
定位模块403,用于基于模板匹配定位所述第二图像中待检测的目标元器件图像;
分割模块404,用于通过预设图像分割算法,将目标元器件图像从所述第二图像中分割出来;
异物检测模块405,用于基于局部信息熵结合区域生长算法,对所述X射线图像中的多余物进行检测,以检测到所述目标元器件的缺陷。
可选的,所述降噪模块包括:
第一子模块,用于获取目标元器件的X射线图像;
第二子模块,用于基于小波阈值算法,对所述X射线图像进行降噪处理,得到第一图像。
可选的,所述消除模块具体用于:
基于二维伽马函数的自适应亮度校正方法,对所述第一图像进行处理,得到消除非均匀光照背景的第二图像。
可选的,所述定位模块具体用于:
采用归一化相关系数匹配法,将模板与所述第二图像进行比对,定位到第二图像中的目标元器件图像。
可选的,所述预设图像分割算法为Grabcut算法。
本申请实施例提供的元器件缺陷检测装置,可用于执行上述元器件缺陷检测方法的相关流程。二者的相关之处,参照上述元器件缺陷检测方法实施例中的相关说明即可,在此不再赘述。
需要说明的是,以上所述,仅为本发明最佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。
本发明说明书中未作详细描述的内容属于本领域技术人员的公知技术。
Claims (10)
1.一种元器件缺陷检测方法,其特征在于,所述方法包括:
获取目标元器件的X射线图像,并对所述X射线图像进行降噪处理,得到第一图像;
消除所述第一图像的非均匀光照背景,得到第二图像;
基于模板匹配定位所述第二图像中待检测的目标元器件图像;
通过预设图像分割算法,将目标元器件图像从所述第二图像中分割出来;
基于局部信息熵结合区域生长算法,对所述X射线图像中的多余物进行检测,以检测到所述目标元器件的缺陷。
2.根据权利要求1所述的方法,其特征在于,所述获取目标元器件的X射线图像,并对所述X射线图像进行降噪处理,得到第一图像的步骤,包括:
获取目标元器件的X射线图像;
基于小波阈值算法,对所述X射线图像进行降噪处理,得到第一图像。
3.根据权利要求1所述的方法,其特征在于,消除所述第一图像的非均匀光照背景,得到第二图像的步骤,包括:
基于二维伽马函数的自适应亮度校正方法,对所述第一图像进行处理,得到消除非均匀光照背景的第二图像。
4.根据权利要求1所述的方法,其特征在于,基于模板匹配定位所述第二图像中待检测的目标元器件图像的步骤,包括:
采用归一化相关系数匹配法,将模板与所述第二图像进行比对,定位到第二图像中的目标元器件图像。
5.根据权利要求1所述的方法,其特征在于,所述预设图像分割算法为Grabcut算法。
6.一种元器件缺陷检测装置,其特征在于,所述装置包括:
降噪模块,用于获取目标元器件的X射线图像,并对所述X射线图像进行降噪处理,得到第一图像;
消除模块,用于消除所述第一图像的非均匀光照背景,得到第二图像;
定位模块,用于基于模板匹配定位所述第二图像中待检测的目标元器件图像;
分割模块,用于通过预设图像分割算法,将目标元器件图像从所述第二图像中分割出来;
异物检测模块,用于基于局部信息熵结合区域生长算法,对所述X射线图像中的多余物进行检测,以检测到所述目标元器件的缺陷。
7.根据权利要求6所述的装置,其特征在于,所述降噪模块包括:
第一子模块,用于获取目标元器件的X射线图像;
第二子模块,用于基于小波阈值算法,对所述X射线图像进行降噪处理,得到第一图像。
8.根据权利要求6所述的装置,其特征在于,所述消除模块具体用于:
基于二维伽马函数的自适应亮度校正方法,对所述第一图像进行处理,得到消除非均匀光照背景的第二图像。
9.根据权利要求6所述的装置,其特征在于,所述定位模块具体用于:
采用归一化相关系数匹配法,将模板与所述第二图像进行比对,定位到第二图像中的目标元器件图像。
10.根据权利要求6所述的装置,其特征在于,所述预设图像分割算法为Grabcut算法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111314595.4A CN114155197A (zh) | 2021-11-08 | 2021-11-08 | 一种元器件缺陷检测的方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111314595.4A CN114155197A (zh) | 2021-11-08 | 2021-11-08 | 一种元器件缺陷检测的方法及装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN114155197A true CN114155197A (zh) | 2022-03-08 |
Family
ID=80459160
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111314595.4A Pending CN114155197A (zh) | 2021-11-08 | 2021-11-08 | 一种元器件缺陷检测的方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114155197A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115115616A (zh) * | 2022-07-26 | 2022-09-27 | 南通浩盛汽车科技有限公司 | 基于可变光源的汽车冲压配件缺陷检测方法 |
CN116777871A (zh) * | 2023-06-20 | 2023-09-19 | 无锡日联科技股份有限公司 | 一种基于x射线的缺陷检测方法、装置、设备及介质 |
-
2021
- 2021-11-08 CN CN202111314595.4A patent/CN114155197A/zh active Pending
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115115616A (zh) * | 2022-07-26 | 2022-09-27 | 南通浩盛汽车科技有限公司 | 基于可变光源的汽车冲压配件缺陷检测方法 |
CN115115616B (zh) * | 2022-07-26 | 2022-11-18 | 南通浩盛汽车科技有限公司 | 基于可变光源的汽车冲压配件缺陷检测方法 |
CN116777871A (zh) * | 2023-06-20 | 2023-09-19 | 无锡日联科技股份有限公司 | 一种基于x射线的缺陷检测方法、装置、设备及介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Bui et al. | Single image dehazing using color ellipsoid prior | |
Celik | Spatial entropy-based global and local image contrast enhancement | |
CN109840913B (zh) | 一种乳腺x线图像中肿块分割的方法和系统 | |
JP4590471B2 (ja) | 背景色を推定する方法およびシステム | |
CN111768411B (zh) | 冠脉中心线提取方法、装置、计算机设备和存储介质 | |
CN114155197A (zh) | 一种元器件缺陷检测的方法及装置 | |
Hosotani et al. | Image denoising with edge-preserving and segmentation based on mask NHA | |
Yuan et al. | Image haze removal via reference retrieval and scene prior | |
Zheng et al. | No-reference quality assessment for screen content images based on hybrid region features fusion | |
CN105261006B (zh) | 基于傅里叶变换的医学图像分割算法 | |
CN112581483B (zh) | 基于自学习的植物叶片叶脉分割方法和装置 | |
CN111325728B (zh) | 产品缺陷检测方法、装置、设备及存储介质 | |
CN110186929A (zh) | 一种实时的产品瑕疵定位方法 | |
CN111861918A (zh) | 一种基于sar图像的海上溢油检测方法 | |
CN115719310A (zh) | 一种眼底图像数据集的预处理方法及眼底图像训练模型 | |
Guan et al. | Quality estimation based multi-focus image fusion | |
CN110472472A (zh) | 基于sar遥感图像的机场检测方法与装置 | |
Dhar et al. | A novel method for underwater image segmentation based on M-band wavelet transform and human psychovisual phenomenon (HVS) | |
CN117274216B (zh) | 基于水平集分割的超声下颈动脉斑块检测方法及系统 | |
CN112417177B (zh) | 一种基于调研数据的知识图谱构建方法 | |
Saxena et al. | Performance Analysis of Single Image Fog Expulsion Techniques | |
Ozturk et al. | Efficient and natural image fusion method for low-light images based on active contour model and adaptive gamma correction | |
Subashini et al. | A combined preprocessing scheme for texture based ice classification from SAR images | |
Bouzaraa et al. | A Naive Bayes approach to improve histogram matching quality | |
Halomoan et al. | Combination of Adaptive Thresh old Algorithm, Watershed and Top-Hat Transform on Image Segmentation Uneven Lighting |
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 |