CN102646278A - 一种基于x光图像的元器件封装气泡检测方法 - Google Patents

一种基于x光图像的元器件封装气泡检测方法 Download PDF

Info

Publication number
CN102646278A
CN102646278A CN2012100712346A CN201210071234A CN102646278A CN 102646278 A CN102646278 A CN 102646278A CN 2012100712346 A CN2012100712346 A CN 2012100712346A CN 201210071234 A CN201210071234 A CN 201210071234A CN 102646278 A CN102646278 A CN 102646278A
Authority
CN
China
Prior art keywords
bubble
ray image
image
edge
pixel
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.)
Granted
Application number
CN2012100712346A
Other languages
English (en)
Other versions
CN102646278B (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.)
South China University of Technology SCUT
Original Assignee
South China University of Technology SCUT
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 South China University of Technology SCUT filed Critical South China University of Technology SCUT
Priority to CN201210071234.6A priority Critical patent/CN102646278B/zh
Publication of CN102646278A publication Critical patent/CN102646278A/zh
Application granted granted Critical
Publication of CN102646278B publication Critical patent/CN102646278B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Analysing Materials By The Use Of Radiation (AREA)

Abstract

本发明公开了一种基于X光图像的元器件封装气泡检测方法,包括以下步骤:(1)采集被封装元器件的X光图像,得到X光图像矩阵;(2)对所述X光图像矩阵进行预处理;(3)初步确定气泡的边缘:通过Canny边缘检测法对经预处理的X光图像进行检侧,得到二值图像矩阵A;(4)确定气泡的边缘轮廓、质心以及面积。本发明实现了对元器件封装中的气泡进行精确定位,并且具有检测速度快的优点。

Description

一种基于X光图像的元器件封装气泡检测方法
技术领域
本发明涉及在元器件封装中的气泡检测方法,特别涉及一种基于X光图像的元器件封装气泡检测方法。
背景技术
封装,是指把硅片上的电路管脚用导线接引到外部接头处,以便与其它器件连接。封装形式是指安装半导体集成电路芯片用的外壳,它不仅起着安装、固定、密封、保护芯片及增强电热性能等方面的作用,而且还通过芯片上的接点用导线连接到封装外壳的引脚上,这些引脚又通过印刷电路板上的导线与其他器件相连接,从而实现内部芯片与外部电路的连接。因为芯片必须与外界隔离,以防止空气中的杂质对芯片电路的腐蚀而造成电气性能下降。另一方面,封装后的芯片也更便于安装和运输。由于封装技术的好坏还直接影响到芯片自身性能的发挥和与之连接的PCB(印制电路板)的设计和制造,因此它是至关重要的。往往在电子元器件封装过程中在基板与芯片之间以及不同的基板之间常常出现气泡,气泡的出现将严重影响电子元器件的电气特性,机械特性以及使用寿命等。因此元器件检测中气泡检测成为一个非常重要的指标。通过检测单个气泡的大小位置以及所有气泡的面积来确定电子元器件的等级。由于X光图像的低对比度、低灰度、高噪声,气泡形状的任意性、位置的随机性气泡检测技术已成为封装电子元器件检测的重点和难点。对于气泡检测是影响电子元器件检测的关键因素,也是本发明解决的问题。
虽然,国内研究人员已在电子元器件无损检测方面做了大量工作,但气泡检测不仅受生产过程中现场噪声和工况的影响,而且还要根据X光和电子元器件自身形状特点来考虑电子元器件的气泡检测。目前已有的检测方法不能很好的检测到气泡的轮廓,而且确定的气泡的面积质心准确性不高,已有检测方法检测速度慢不能满足现场的实时性要求。因此,为了提高基于X光图像的电子元器件气泡检测技术,就必须发明新的方法,以提高检测速度与准确性。
发明内容
为了克服现有技术的上述不足,本发明的目的在于提供一种基于X光图像的元器件封装气泡检测方法。
本发明的目的通过以下技术方案实现:
一种基于X光图像的元器件封装气泡检测方法,包括以下步骤:
(1)采集被封装元器件的X光图像,得到X光图像矩阵;
(2)对所述X光图像矩阵进行预处理;
(3)初步确定气泡的边缘:通过Canny边缘检测法对经预处理的X光图像进行检测,得到气泡的边缘,并得到图像的二值矩阵A;
(4)确定气泡的边缘轮廓、质心以及面积,具体为:
(4-1)设定膨胀腐蚀模板矩阵B,先通过二值图像矩阵A被B形态学开运算去除孤立点,再通过A被B进行形态学闭运算连接非封闭边缘,得到二值图像的封闭区域;
(4-2)对二值图像的封闭区域进行气泡填充:
设Xk为气泡内的点,令X0=1,则气泡填充过程为:
根据
Figure BDA0000144051500000021
进行迭代,直至Xk=Xk-1;{Xk| k=0,1、2...}为组成气泡的点集;
其中 C = 0 1 0 1 1 1 0 1 0 ;
(4-3)对气泡轮廓进行标记:
设Yk为气泡轮廓上的点,Y0为气泡轮廓上的一个已知点,对二值图像中的连通分量进行提取,其过程如下:
根据
Figure BDA0000144051500000023
进行迭代,直至Yk=Yk-1;{Yk| k=0,1、2...}为标志气泡轮廓的点集;
(4-4)根据(4-3)标记的气泡轮廓确定气泡的质心及面积。
步骤(2)所述对所述X光图像矩阵进行预处理,包括以下步骤:
(2-1)通过高斯滤波器进行高斯滤波;
(2-2)对高斯滤波后的图像进行对比度拉伸。
步骤(3)所述通过Canny边缘检测法对经预处理的X光图像进行检侧,包括以下步骤:
(3-1)使用带有指定标准偏差σ的高斯滤波器对X光图像进行平滑处理;
(3-2)计算X光图像每一点的局部梯度
Figure BDA0000144051500000024
和边缘方向α(x,y)=arctan(ly/lx),并确定边缘点,得到梯度幅度图像;其中边缘点定义为梯度方向上其强度局部最大点;
(3-3)对梯度幅度图像进行非最大值抑制处理,得到脊像素;
(3-4)用阈值T1和T2对脊像素进行阈值处理,值大于T2的脊像素为强边缘像素,T1与T2之间的脊像素为弱边缘像素;阈值T1和T2根据实际情况确定;
(3-5)将8连接的弱边缘像素集成到强边缘像素,执行边缘连接。
步骤(2-1)所述通过高斯滤波器进行高斯滤波,具体为:
假设f(x,y)为输入图像,h(s,t)为高斯滤波模板,g(x,y)为滤波后的输出图像,根据以下公式进行高斯滤波过程:
g ( x , y ) = Σ s = - a a Σ s = - b b h ( s , t ) * f ( x + s , y + t ) Σ s = - a a Σ t = - b b h ( s , t ) .
步骤(2-2)所述对高斯滤波后的图像进行对比度拉伸,具体为:
对比度拉伸后图像l(x,y)=c*log(1+(g(x,y))),其中c为对比度拉伸系数,由用户根据具体需要确定。
步骤(3-3)所述对梯度幅度图像进行非最大值抑制处理,具体为:追踪梯度幅度图像中脊的顶部,将所有不在脊的顶部的像素设为零。
与现有技术相比,本发明具有以下优点和有益效果:本发明针对低对比度、低灰度、高噪声的X光图像,对图像中的任意形状、任意位置、任意大小的气泡均可实现检测,并且具有更高的准确率,检测速度更快。
附图说明
图1为本发明基于X光图像的元器件封装气泡检测方法的流程图。
具体实施方式
下面结合实施例及附图,对本发明作进一步地详细说明,但本发明的实施方式不限于此。
实施例
如图1所示,本发明所述基于X光图像的元器件封装气泡检测方法,包括以下步骤:
(1)采集被封装元器件的X光图像,得到X光图像矩阵;
(2)对所述X光图像矩阵进行预处理,包括以下步骤:
(2-1)通过高斯滤波器进行高斯滤波,具体为:
假设f(x,y)为输入图像,H(s,t)为二维高斯函数,g(x,y)为滤波后的输出图像,高斯滤波过程为:
H ( s , t ) = e - ( s 2 + t 2 ) 2 δ 2
g ( x , y ) = f ( x , y ) ⊕ H ( s , t )
可由高斯函数H(s,t)离散化确定高斯滤波模板h(s,t),再根据以下公式进行高斯滤波:
g ( x , y ) = Σ s = - a a Σ s = - b b h ( s , t ) * f ( x + s , y + t ) Σ s = - a a Σ t = - b b h ( s , t ) ;
(2-2)对高斯滤波后的图像进行对比度拉伸,以增加气泡与图像背景的对比度,具体为:
对比度拉伸后图像l(x,y)=c*log(1+(g(x,y)));其中c为对比度拉伸系数,由用户根据具体需要确定。
(3)初步确定气泡的边缘:通过Canny边缘检测法对经预处理的X光图像进行检侧,得到气泡的边缘,并得到二值图像矩阵A;
所述通过Canny边缘检测法对经预处理的X光图像进行检侧,包括以下步骤:
(3-1)使用带有指定标准偏差σ的高斯滤波器对X光图像进行平滑处理;
(3-2)计算X光图像每一点的局部梯度
Figure BDA0000144051500000044
和边缘方向α(x,y)=arctan(ly/lx),并确定边缘点,得到梯度幅度图像;其中边缘点定义为梯度方向上其强度局部最大点;
(3-3)对梯度幅度图像进行非最大值抑制处理,得到脊像素;
(3-4)用阈值T1和T2对脊像素进行阈值处理,值大于T2的脊像素为强边缘像素,T1与T2之间的脊像素为弱边缘像素;阈值T1和T2根据实际情况确定;
(3-5)将8连接的弱边缘像素集成到强边缘像素,执行边缘连接。
(4)确定气泡的边缘轮廓、质心以及面积,具体为:
(4-1)设定膨胀腐蚀模板矩阵为B,先通过二值图像矩阵A被B形态学开运算去除孤立点,再通过A被B进行形态学闭运算连接非封闭边缘,得到二值图像的封闭区域;
(4-2)对二值图像的封闭区域进行气泡填充:
设Xk为气泡内的点,令X0=1,则气泡填充过程为:
根据进行迭代,直至Xk=Xk-1;{Xk| k=0,1、2...}为组成气泡的点集;
其中 C = 0 1 0 1 1 1 0 1 0 ;
(4-3)对气泡轮廓进行标记:
设Yk为气泡轮廓上的点,Y0为气泡轮廓上的一个已知点,对二值图像中的连通分量进行提取,其过程如下:
根据
Figure BDA0000144051500000053
进行迭代,直至Yk=Yk-1;{Yk| k=0,1、2...}为标志气泡轮廓的点集;
(4-4)根据(4-3)标记的气泡轮廓确定气泡的质心及面积。
上述实施例为本发明较佳的实施方式,但本发明的实施方式并不受所述实施例的限制,其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化,均应为等效的置换方式,都包含在本发明的保护范围之内。

Claims (6)

1.一种基于X光图像的元器件封装气泡检测方法,其特征在于,包括以下步骤:
(1)采集被封装元器件的X光图像,得到X光图像矩阵;
(2)对所述X光图像矩阵进行预处理;
(3)初步确定气泡的边缘:通过Canny边缘检测法对经预处理的X光图像进行检侧,得到气泡的边缘,并得到二值图像矩阵A;
(4)确定气泡的边缘轮廓、质心以及面积,具体为:
(4-1)设定膨胀腐蚀模板矩阵为B,先通过二值图像矩阵A被B形态学开运算去除孤立点,再通过A被B进行形态学闭运算连接非封闭边缘,得到二值图像的封闭区域;
(4-2)对二值图像的封闭区域进行气泡填充:
设Xk为气泡内的点,令X0=1,则气泡填充过程为:
根据进行迭代,直至Xk=Xk-1;{Xk| k=0,1、2...}为组成气泡的点集;
其中 C = 0 1 0 1 1 1 0 1 0 ;
(4-3)对气泡轮廓进行标记:
设Yk为气泡轮廓上的点,Y0为气泡轮廓上的一个已知点,对二值图像中的连通分量进行提取,其过程如下:
根据
Figure FDA0000144051490000013
进行迭代,直至Yk=Yk-1;{Yk| k=0,1、2...}为标志气泡轮廓的点集;
(4-4)根据(4-3)标记的气泡轮廓确定气泡的质心及面积。
2.根据权利要求1所述的一种基于X光图像的元器件封装气泡检测方法,其特征在于,步骤(2)所述对所述X光图像矩阵进行预处理,包括以下步骤:
(2-1)通过高斯滤波器进行高斯滤波;
(2-2)对高斯滤波后的图像进行对比度拉伸。
3.根据权利要求1所述的一种基于X光图像的元器件封装气泡检测方法,其特征在于,步骤(3)所述通过Canny边缘检测法对经预处理的X光图像进行检侧,包括以下步骤:
(3-1)使用带有指定标准偏差σ的高斯滤波器对X光图像进行平滑处理;
(3-2)计算X光图像每一点的局部梯度
Figure FDA0000144051490000021
和边缘方向α(x,y)=arctan(ly/lx),并确定边缘点,得到梯度幅度图像;其中边缘点定义为梯度方向上其强度局部最大点;
(3-3)对梯度幅度图像进行非最大值抑制处理,得到脊像素;
(3-4)用阈值T1和T2对脊像素进行阈值处理,值大于T2的脊像素为强边缘像素,T1与T2之间的脊像素为弱边缘像素;阈值T1和T2根据实际情况确定;
(3-5)将8连接的弱边缘像素集成到强边缘像素,执行边缘连接。
4.根据权利要求2所述的一种基于X光图像的元器件封装气泡检测方法,其特征在于,步骤(2-1)所述通过高斯滤波器进行高斯滤波,具体为:
假设f(x,y)为输入图像,h(s,t)为高斯滤波模板,g(x,y)为滤波后的输出图像,根据以下公式进行高斯滤波过程:
g ( x , y ) = Σ s = - a a Σ s = - b b h ( s , t ) * f ( x + s , y + t ) Σ s = - a a Σ t = - b b h ( s , t ) .
5.根据权利要求4所述的一种基于X光图像的元器件封装气泡检测方法,其特征在于,步骤(2-2)所述对高斯滤波后的图像进行对比度拉伸,具体为:
对比度拉伸后图像l(x,y)=c*log(1+(g(x,y))),其中c为对比度拉伸系数,由用户根据具体需要确定。
6.根据权利要求4所述的一种基于X光图像的元器件封装气泡检测方法,其特征在于,步骤(3-3)所述对梯度幅度图像进行非最大值抑制处理,具体为:追踪梯度幅度图像中脊的顶部,将所有不在脊的顶部的像素设为零。
CN201210071234.6A 2012-03-16 2012-03-16 一种基于x光图像的元器件封装气泡检测方法 Active CN102646278B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210071234.6A CN102646278B (zh) 2012-03-16 2012-03-16 一种基于x光图像的元器件封装气泡检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210071234.6A CN102646278B (zh) 2012-03-16 2012-03-16 一种基于x光图像的元器件封装气泡检测方法

Publications (2)

Publication Number Publication Date
CN102646278A true CN102646278A (zh) 2012-08-22
CN102646278B CN102646278B (zh) 2014-07-16

Family

ID=46659090

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210071234.6A Active CN102646278B (zh) 2012-03-16 2012-03-16 一种基于x光图像的元器件封装气泡检测方法

Country Status (1)

Country Link
CN (1) CN102646278B (zh)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103020392A (zh) * 2012-12-29 2013-04-03 上海华碧检测技术有限公司 一种空洞率的测算方法
CN103500459A (zh) * 2013-09-24 2014-01-08 华南理工大学 一种基于规则化参数的印刷电路板缺陷检测方法
WO2017025001A1 (zh) * 2015-08-07 2017-02-16 深圳市瀚海基因生物科技有限公司 单分子定位装置、方法及系统
CN106530287A (zh) * 2016-10-24 2017-03-22 武汉新芯集成电路制造有限公司 一种基于晶圆内部缺陷检测的图像自动识别系统
CN106780493A (zh) * 2017-01-24 2017-05-31 广东省特种设备检测研究院珠海检测院 一种基于焊缝x射线数字图像的圆形缺陷识别方法
CN108090892A (zh) * 2016-11-22 2018-05-29 通用汽车环球科技运作有限责任公司 接头质量的自动确定
CN109166083A (zh) * 2018-09-03 2019-01-08 哈尔滨工业大学 一种用于去除水下图像气泡噪声的方法
CN110361400A (zh) * 2019-07-01 2019-10-22 创新奇智(合肥)科技有限公司 一种铸铁工件的气泡检测方法及电子设备
CN114331923A (zh) * 2022-03-11 2022-04-12 中国空气动力研究与发展中心低速空气动力研究所 一种基于改进Canny算法的冰结构中气泡轮廓提取方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030099330A1 (en) * 2000-02-05 2003-05-29 Domingo Mery Method for automatically detecting casting defects in a test piece
CN102175700A (zh) * 2011-01-20 2011-09-07 山东大学 数字x射线图像焊缝分割和缺陷检测方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030099330A1 (en) * 2000-02-05 2003-05-29 Domingo Mery Method for automatically detecting casting defects in a test piece
CN102175700A (zh) * 2011-01-20 2011-09-07 山东大学 数字x射线图像焊缝分割和缺陷检测方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
何新英等: "基于数学形态学和Canny算子的边缘提取方法", 《计算机应用》 *
赵熹等: "X射线检测图像中缺陷部分的自动提取", 《无损探伤》 *
高红霞等: "基于X光检测的元器件任意形状气泡提取方法", 《半导体技术》 *

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103020392A (zh) * 2012-12-29 2013-04-03 上海华碧检测技术有限公司 一种空洞率的测算方法
CN103500459A (zh) * 2013-09-24 2014-01-08 华南理工大学 一种基于规则化参数的印刷电路板缺陷检测方法
CN103500459B (zh) * 2013-09-24 2017-01-18 华南理工大学 一种基于规则化参数的印刷电路板缺陷检测方法
WO2017025001A1 (zh) * 2015-08-07 2017-02-16 深圳市瀚海基因生物科技有限公司 单分子定位装置、方法及系统
CN106530287A (zh) * 2016-10-24 2017-03-22 武汉新芯集成电路制造有限公司 一种基于晶圆内部缺陷检测的图像自动识别系统
CN108090892A (zh) * 2016-11-22 2018-05-29 通用汽车环球科技运作有限责任公司 接头质量的自动确定
CN108090892B (zh) * 2016-11-22 2022-02-15 通用汽车环球科技运作有限责任公司 自动评定接头质量的方法和组件
CN106780493A (zh) * 2017-01-24 2017-05-31 广东省特种设备检测研究院珠海检测院 一种基于焊缝x射线数字图像的圆形缺陷识别方法
CN109166083A (zh) * 2018-09-03 2019-01-08 哈尔滨工业大学 一种用于去除水下图像气泡噪声的方法
CN110361400A (zh) * 2019-07-01 2019-10-22 创新奇智(合肥)科技有限公司 一种铸铁工件的气泡检测方法及电子设备
CN114331923A (zh) * 2022-03-11 2022-04-12 中国空气动力研究与发展中心低速空气动力研究所 一种基于改进Canny算法的冰结构中气泡轮廓提取方法
CN114331923B (zh) * 2022-03-11 2022-05-13 中国空气动力研究与发展中心低速空气动力研究所 一种基于改进Canny算法的冰结构中气泡轮廓提取方法

Also Published As

Publication number Publication date
CN102646278B (zh) 2014-07-16

Similar Documents

Publication Publication Date Title
CN102646278B (zh) 一种基于x光图像的元器件封装气泡检测方法
CN106066919B (zh) 应用于近/亚阈值数字电路的统计静态时序分析方法
CN103884980A (zh) 基于电源电流的硬件木马检测方法和系统
CN102049732A (zh) 一种硅片边缘膜厚测量方法
CN104950246A (zh) 基于延时的硬件木马检测方法和系统
CN104199276A (zh) 基于fpga的信号时差测量方法及时间数字转换器
CN103702516B (zh) 印制电路板系数计算方法及计算系统
CN104933252A (zh) 一种面向多抽象层次电路的可靠性评估方法
WO2010011448A3 (en) Predictive modeling of interconnect modules for advanced on-chip interconnect technology
CN104484747A (zh) 利用截尾样本确定产品合格率的方法
CN103632342A (zh) 一种用于集成电路封装中的x射线图像模糊增强方法
CN104062008A (zh) 一种考虑整体度量的实测光谱曲线中异常光谱的剔除方法
WO2007027758A3 (en) Functional cells for automated i/o timing characterization of an integrated circuit
CN110135161B (zh) 一种硬件木马的在片检测方法
CN106777523B (zh) 一种电阻子电路噪声模型结构及其建模方法
CN112581447A (zh) 一种基于全局缺陷和局部缺陷的fpc软板线路检测方法
Lin A feasible method to predict thin package actual warpage based on an FEM model integrated with empirical data
CN115690064A (zh) 对准方法、装置、电子设备及计算机可读存储介质
CN105975646A (zh) 分析半导体装置的性能的系统、方法和计算机程序产品
Cai et al. Improving the measurement accuracy of an absolute imaging position encoder via a new edge detection method
CN101398864B (zh) 电路板制作仿真系统及方法
CN106250594A (zh) 一种PCB中自动检测Place_Bound的方法
CN102214252A (zh) 一种对半导体器件进行提参建模的方法
CN107908912A (zh) 一种自动检查线路电源与地连接的方法
CN105699734A (zh) 检测信号延迟时间的装置及方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant