一种基于鲁棒主成分分析的ic元件焊点检测方法
【专利摘要】本发明公开了一种基于鲁棒主成分分析的IC元件焊点检测方法，包括：采集多个待检测IC元件焊点的测试图片后，形成测试集合；分别将测试集合的每张测试图片向量化成列向量后，将获得的所有列向量合并成IC焊点测试矩阵；进行优化分解处理后获得其低秩矩阵和稀疏矩阵；将稀疏矩阵的每个列向量进行逆向量化后获得对应的图片，进而对每个图片进行二值化分割得到对应的二值化图像；获取根据训练样本训练获得的频率分布图后，结合二值化图像计算每个测试图片的缺陷度；将该图片的缺陷度与训练样本的缺陷度阈值进行比对后获得IC元件焊点的检测结果。本发明可以准确、有效、快速地进行IC元件焊点缺陷检测，可广泛应用于IC元件的焊点检测行业中。
【专利说明】
一种基于鲁棒主成分分析的1C元件焊点检测方法
技术领域
[0001] 本发明涉及数字图像处理应用领域，特别是涉及一种基于鲁棒主成分分析的1C元 件焊点检测方法。
【背景技术】
[0002] 名词解释：
[0003] RPCA:全称Robust principle component analysis，鲁棒主成分分析；
[0004] PCB:全称Printed Circuit Board，印制电路板；
[0005] IC: Integrated Circuit，集成电路；
[0006] 印制电路板缺陷检测是自动光学检测（automatic optical inspection，Α0Ι)应 用的一个热点方向，近年来得到越来越多的关注。目前主要采用的检测方式是通过CCD检测 PCB元件的图像后进行处理从而实现缺陷检测。实际使用中，印制电路板遇到的情况非常复 杂，CCD采集到的PCB元件图像往往出现不同程度的变化和不规则现象，例如:光照强度不均 匀，光照角度改变，CCD摄像头采集的图像会有一定角度的偏转，元件尺寸越来越小，PCB板 内元件密度越来越大等等，这些问题使得PCB焊点缺陷检测变得相当困难。而1C元件焊点的 尺寸比一般的CHIP元件焊点的尺寸要小的多，虚焊和正常焊点在图像上非常相似，这使得 1C元件焊点的虚焊检测一直是缺陷检测中难以攻克的难题。
[0007] 现有的较为成熟的对1C元件焊点缺陷的检测方法大多数为基于特征的方法。这种 方法将缺陷检测分为两个步骤:提取特征和分类。在提取特征阶段，选择颜色梯度，区域面 积，周长，水力半径等等典型特征；在分类阶段，选择较为成熟的分类器，例如神经网络， AdaB 〇〇st，SVM等等，对提取的特征进行分类。这些方法在CHIP元件焊点取得了较好的应用 效果。但是，由于1C元件焊点尺寸小，焊点密集度大，且虚焊焊点样本难以收集，使得目前的 基于分类器的方法在1C元件焊点虚焊检测中难以得到良好的分类效果。另外，在线检测方 法对时间要求严格，神经网络这些较为成熟的分类器计算量大，难以符合在线监测在时间 上的要求。同时虽然也有人提出基于单高斯模型的像素点建模的1C元件焊点检测方法，这 种方法的检测速度快，但是准确率低，无法在实际生产中应用。
[0008] 另外，目前也出现了基于统计建模的方法，此类方法将焊点的缺陷检测问题作为 一般的目标检测问题，也就是说，将合格的样本图片当作背景，把不合格的图片当作前景。 然后，进行背景建模，最后使用该模型进行1C焊点缺陷检测。然而，目前的这类统计建模的 方法都是基于像素点的，而在目标检测当中，基于像素点的方法往往会对背景部分的考虑 过于简单，因为这种方式主要是认为背景是静止的或者是变化缓慢的，而在合格的1C引脚 图片样本中，每一张图片相比都有很大程度的不同，即背景是剧烈变化的。这会严重地抑制 训练和检测的效果，导致模型训练结果不理想，从而导致检测结果不准确。
[0009] 总的来说，目前的检测方法无法有效、准确、快速地检测出1C元件的焊点缺陷。
【发明内容】
[0010] 为了解决上述的技术问题，本发明的目的是提供一种基于鲁棒主成分分析的ic元 件焊点检测方法。
[0011] 本发明解决其技术问题所采用的技术方案是：
[0012] -种基于鲁棒主成分分析的1C元件焊点检测方法，包括步骤：
[0013] S1、采集多个待检测1C元件焊点的测试图片后，形成测试集合；
[0014] S2、分别将测试集合的每张测试图片向量化成列向量后，将获得的所有列向量合 并成1C焊点测试矩阵；
[0015] S3、对1C焊点测试矩阵进行优化分解处理后获得其低秩矩阵和稀疏矩阵；
[0016] S4、将稀疏矩阵的每个列向量进行逆向量化后获得对应的图片，进而对每个图片 进行二值化分割从而得到对应的二值化图像；
[0017] S5、获取根据训练样本训练获得的频率分布图后，结合二值化图像计算每个测试 图片的缺陷度；
[0018] S6、将该图片的缺陷度与训练样本的缺陷度阈值进行比对后获得1C元件焊点的检 测结果。
[0019] 进一步，所述步骤S1之前还包括对训练样本进行训练的步骤，具体包括：
[0020] A1、分别将训练样本的每张1C焊点训练图片向量化成列向量后，将获得的所有列 向量合并成1C焊点样本矩阵；
[0021] A2、对1C焊点样本矩阵进行优化分解处理后获得其低秩矩阵和稀疏矩阵；
[0022] A3、将稀疏矩阵的每个列向量进行逆向量化后获得对应的样本图片，进而对每个 样本图片进行二值化分割从而得到对应的二值化图像；
[0023] A4、建立频率分布图后，依次获取每个二值化图像，对频率分布图进行更新，进而 在更新完毕后计算训练样本的缺陷度阈值。
[0024]进一步，所述步骤S2,其具体为：
[0025]分别将测试集合的每张测试图片向量化成NX 1维列向量后，将获得的所有列向量 合并成NXM维的1C焊点测试矩阵;其中，Μ表示测试集合的测试图片的数量，N表示列向量的 维度，且N=WXH，W表示测试图片的宽度，Η表示测试图片的高度。
[0026] 进一步，所述步骤S3,包括：
[0027] S31、建立如下的增广拉格朗日函数作为1C焊点测试矩阵的优化函数：
[0029]其中，ML，S，X)表示优化函数，D表示1C焊点测试矩阵，L表示1C焊点测试矩阵的 低秩矩阵，S表示1C焊点测试矩阵的稀疏矩阵，X表示与1C焊点测试矩阵相同维度的拉格朗 日乘子矩阵，μ表示补偿正标量，I |L| I*表示低秩矩阵L的核范数，I |S| U表示稀疏矩阵3的1^ 范数，λ表示调节参数，I I · I If表示Frobenius范数，〈.，.>表示矩阵内积运算子；
[0030] S32、求解获得低秩矩阵L、稀疏矩阵S和拉格朗日乘子矩阵X的迭代更新公式如下：
[0032] 其中，上标k表示迭代次数，υ Σ 示矩阵（
的奇异值分解，? (.)表 示压缩算子
表示矩阵Σ中元素如果大于
，则取原值，否则，取0;
[0033] S33、根据迭代更新公式对低秩矩阵L、稀疏矩阵S和拉格朗日乘子矩阵X进行更新， 同时更新补偿正标量μ和迭代次数k: yk+1=ykP，k=k+Ι，其中，p表示迭代参数；
[0034] S34、迭代执行步骤S33直到满足收敛条件后，获取当前的低秩矩阵和稀疏矩阵。 [00 35]进一步，所述步骤S33,其具体为：
[0036] 根据迭代更新公式对低秩矩阵L、稀疏矩阵S和拉格朗日乘子矩阵X进行更新，同时 更新补偿正标量以和迭代次数1^1^ 1 =以1^，1^ = 1^+1，其中，初始化的迭代条件如下:1^ = 3()= Χ°=0,μ°>0,Ρ>1Λ=0〇
[0037] 进一步，所述步骤S34中所述收敛条件为：| |D-Lk-Sk| |F/| |D| |Κ10-7或迭代次数k 达到预设最大次数。
[0038]进一步，所述步骤S4,其具体为：
[0039]将稀疏矩阵的每个列向量进行逆向量化后获得对应的WXH的图片，进而对每个图 片进行二值化分割从而得到对应的二值化图像，W表示测试图片的宽度，Η表示测试图片的 高度。
[0040] 进一步，所述步骤S5中所述频率分布图是通过以下方式训练获得的：
[0041] 建立一个WXH的矩阵并将该矩阵的所有元素值初始化为1后作为频率分布图，进 而依次获取训练样本的每张1C焊点训练图片的稀疏矩阵对应的二值化图像后，对频率分布 图进行更新，针对二值化图像中像素值为255的像素点，将频率分布图的相应位置的元素值 加1，W表示测试图片的宽度，Η表示测试图片的高度。
[0042] 进一步，所述步骤S5中所述计算每个测试图片的缺陷度的步骤，是通过下式计算 缺陷度的：
[0044] 上式中，Em表示测试图片的缺陷度，Μ表示测试集合的测试图片个数，f(x，y)表示 频率分布图，b(x，y)表示测试图片的二值化图像，W表示测试图片的宽度，Η表示测试图片的 高度。
[0045] 进一步，所述步骤S6,其具体为：
[0046] 判断该缺陷度是否大于训练样本的缺陷度阈值，若是，则判断该1C元件焊点为虚 焊焊点，反之，判断该1C元件焊点为正常焊点。
[0047]进一步，所述步骤SI，其具体为：
[0048] 采集多个待检测1C元件焊点的图片，并将每张图片由RGB颜色空间转换为HSV空间 后，将其Η通道图像作为1C元件焊点的测试图片，进而形成测试集合。
[0049] 进一步，所述训练样本包括多个1C焊点训练图片，每个1C焊点训练图片是通过将 合格的1C引脚图片由RGB颜色空间转换为HSV空间后获得的Η通道图像。
[0050] 进一步，所述测试集合中测试图片的数量与训练样本的1C焊点训练图片的数量一 致，均为Μ。
[0051] 本发明的有益效果是：本发明的一种基于鲁棒主成分分析的1C元件焊点检测方 法，包括步骤:S1、采集多个待检测1C元件焊点的测试图片后，形成测试集合;S2、分别将测 试集合的每张测试图片向量化成列向量后，将获得的所有列向量合并成1C焊点测试矩阵； S3、对1C焊点测试矩阵进行优化分解处理后获得其低秩矩阵和稀疏矩阵;S4、将稀疏矩阵的 每个列向量进行逆向量化后获得对应的图片，进而对每个图片进行二值化分割从而得到对 应的二值化图像;S5、获取根据训练样本训练获得的频率分布图后，结合二值化图像计算每 个测试图片的缺陷度;S6、将该图片的缺陷度与训练样本的缺陷度阈值进行比对后获得1C 元件焊点的检测结果。本方法通过求解稀疏矩阵将1C元件焊点检测扩展到最优化的层面， 因此，任何属于背景部分的变化，即任何1C元件焊点图片之间的不同部分，都能够通过低秩 矩阵捕捉到以及进行约束，最终提高检测效果，因此，本方法可以准确、有效、快速地进行1C 元件焊点缺陷检测。
【附图说明】
[0052]下面结合附图和实施例对本发明作进一步说明。
[0053]图1是本发明的一种基于鲁棒主成分分析的1C元件焊点检测方法的流程图。
【具体实施方式】
[0054]参照图1，本发明提供了一种基于鲁棒主成分分析的1C元件焊点检测方法，包括步 骤：
[0055] S1、采集多个待检测1C元件焊点的测试图片后，形成测试集合；
[0056] S2、分别将测试集合的每张测试图片向量化成列向量后，将获得的所有列向量合 并成1C焊点测试矩阵；
[0057] S3、对1C焊点测试矩阵进行优化分解处理后获得其低秩矩阵和稀疏矩阵；
[0058] S4、将稀疏矩阵的每个列向量进行逆向量化后获得对应的图片，进而对每个图片 进行二值化分割从而得到对应的二值化图像；
[0059] S5、获取根据训练样本训练获得的频率分布图后，结合二值化图像计算每个测试 图片的缺陷度；
[0060] S6、将该图片的缺陷度与训练样本的缺陷度阈值进行比对后获得1C元件焊点的检 测结果。
[0061] 进一步作为优选的实施方式，所述步骤S1之前还包括对训练样本进行训练的步 骤，具体包括：
[0062] Α1、分别将训练样本的每张1C焊点训练图片向量化成列向量后，将获得的所有列 向量合并成1C焊点样本矩阵；
[0063] A2、对1C焊点样本矩阵进行优化分解处理后获得其低秩矩阵和稀疏矩阵；
[0064] A3、将稀疏矩阵的每个列向量进行逆向量化后获得对应的样本图片，进而对每个 样本图片进行二值化分割从而得到对应的二值化图像；
[0065] A4、建立频率分布图后，依次获取每个二值化图像，对频率分布图进行更新，进而 在更新完毕后计算训练样本的缺陷度阈值。
[0066] 进一步作为优选的实施方式，所述步骤S2，其具体为：
[0067] 分别将测试集合的每张测试图片向量化成NX 1维列向量后，将获得的所有列向量 合并成NXM维的1C焊点测试矩阵;其中，Μ表示测试集合的测试图片的数量，N表示列向量的 维度，且N=WXH，W表示测试图片的宽度，Η表示测试图片的高度。
[0068] 进一步作为优选的实施方式，所述步骤S3,包括：
[0069] S31、建立如下的增广拉格朗日函数作为1C焊点测试矩阵的优化函数：
[0071] 其中，lu(L，S，X)表示优化函数，D表示1C焊点测试矩阵，L表示1C焊点测试矩阵的 低秩矩阵，S表示1C焊点测试矩阵的稀疏矩阵，X表示与1C焊点测试矩阵相同维度的拉格朗 日乘子矩阵，μ表示补偿正标量，I |L| I*表示低秩矩阵L的核范数，| |S| U表示稀疏矩阵3的1^ 范数，λ表示调节参数，I I · I If表示Frobenius范数，〈.，.>表示矩阵内积运算子；
[0072] S32、求解获得低秩矩阵L、稀疏矩阵S和拉格朗日乘子矩阵X的迭代更新公式如下：
[0074] 其中，上标k表示迭代次数，示矩阵
1的奇异值分解，?(.)表 示压缩算子且
表示矩阵Σ中元素如果大于
，则取原值，否则，取0;
[0075] S33、根据迭代更新公式对低秩矩阵L、稀疏矩阵S和拉格朗日乘子矩阵X进行更新， 同时更新补偿正标量μ和迭代次数k: yk+1=ykP，k=k+Ι，其中，ρ表示迭代参数；
[0076] S34、迭代执行步骤S33直到满足收敛条件后，获取当前的低秩矩阵和稀疏矩阵。 [0077]进一步作为优选的实施方式，所述步骤S33，其具体为：
[0078] 根据迭代更新公式对低秩矩阵L、稀疏矩阵S和拉格朗日乘子矩阵X进行更新，同时 更新补偿正标量以和迭代次数1^1^ 1 =以1^，1^ = 1^+1，其中，初始化的迭代条件如下:1^ = 3()= Χ°=0,μ°>0,Ρ>1Λ=0〇
[0079] 进一步作为优选的实施方式，所述步骤S34中所述收敛条件为：| |D-Lk-Sk| |f/||D I f〈1(T7或迭代次数k达到预设最大次数。
[0080] 进一步作为优选的实施方式，所述步骤S4，其具体为：
[0081] 将稀疏矩阵的每个列向量进行逆向量化后获得对应的WXH的图片，进而对每个图 片进行二值化分割从而得到对应的二值化图像，W表示测试图片的宽度，Η表示测试图片的 高度。
[0082] 进一步作为优选的实施方式，所述步骤S5中所述频率分布图是通过以下方式训练 获得的：
[0083] 建立一个WXH的矩阵并将该矩阵的所有元素值初始化为1后作为频率分布图，进 而依次获取训练样本的每张1C焊点训练图片的稀疏矩阵对应的二值化图像后，对频率分布 图进行更新，针对二值化图像中像素值为255的像素点，将频率分布图的相应位置的元素值 加1，W表示测试图片的宽度，Η表示测试图片的高度。
[0084] 进一步作为优选的实施方式，所述步骤S5中所述计算每个测试图片的缺陷度的步 骤，是通过下式计算缺陷度的：
[0086] 上式中，Em表示测试图片的缺陷度，Μ表示测试集合的测试图片个数，f(x，y)表示 频率分布图，b(x，y)表示测试图片的二值化图像，W表示测试图片的宽度，Η表示测试图片的 高度。
[0087] 进一步作为优选的实施方式，所述步骤S6，其具体为：
[0088] 判断该缺陷度是否大于训练样本的缺陷度阈值，若是，则判断该1C元件焊点为虚 焊焊点，反之，判断该1C元件焊点为正常焊点。
[0089] 进一步作为优选的实施方式，所述步骤S1，其具体为：
[0090] 采集多个待检测1C元件焊点的图片，并将每张图片由RGB颜色空间转换为HSV空间 后，将其Η通道图像作为1C元件焊点的测试图片，进而形成测试集合。
[0091] 进一步作为优选的实施方式，所述训练样本包括多个1C焊点训练图片，每个1C焊 点训练图片是通过将合格的1C引脚图片由RGB颜色空间转换为HSV空间后获得的Η通道图 像。
[0092]进一步作为优选的实施方式，所述测试集合中测试图片的数量与训练样本的1C焊 点训练图片的数量一致，均为Μ。
[0093]以下结合具体实施例对本发明做详细说明。
[0094] 参照图1，一种基于鲁棒主成分分析的1C元件焊点检测方法，包括步骤：
[0095] S1、采集多个待检测1C元件焊点的测试图片后，形成测试集合;具体如下：采集多 个待检测1C元件焊点的图片，并将每张图片由RGB颜色空间转换为HSV空间后，将其Η通道图 像作为1C元件焊点的测试图片，进而形成测试集合。测试集合中测试图片的数量与训练样 本的1C焊点训练图片的数量一致，均为Μ。
[0096] S2、分别将测试集合的每张测试图片向量化成NX 1维列向量后，将获得的所有列 向量合并成NXM维的1C焊点测试矩阵;其中，Μ表示测试集合的测试图片的数量，N表示列向 量的维度，且N=WXH，W表示测试图片的宽度，Η表示测试图片的高度。
[0097] S3、对1C焊点测试矩阵进行优化分解处理后获得其低秩矩阵和稀疏矩阵，包括S31 ~S34:
[0098] S31、建立如下的增广拉格朗日函数作为1C焊点测试矩阵的优化函数：
[0100] 其中，ML，S，X)表示优化函数，D表示1C焊点测试矩阵，L表示1C焊点测试矩阵的 低秩矩阵，S表示1C焊点测试矩阵的稀疏矩阵，X表示与1C焊点测试矩阵相同维度的拉格朗 日乘子矩阵，μ表示补偿正标量，I |L| I*表示低秩矩阵L的核范数，I |S| U表示稀疏矩阵3的1^ 范数，λ表示调节参数，I I · I If表示Frobenius范数，〈.，.>表示矩阵内积运算子；
[0101] S32、求解获得低秩矩阵L、稀疏矩阵S和拉格朗日乘子矩阵X的迭代更新公式如下：
[0103] 其中，上标k表示迭代次数，U Σ ￥7表示矩阵
的奇异值分解，? (.)表 示压缩算子且
表示矩阵Σ中元素如果大于
则取原值，否则，取0;这里，UXVT中的 矩阵Σ是表示运用压缩算子?(.)进行运算<
后获得的矩阵；
[0104] 某个ΝΧΜ维的矩阵C的奇异值分解1]Σν%^具体含义如下：
[0105] U为ΝΧΝ矩阵，且U的列为CCT的正交特征向量，V为ΜΧΜ的矩阵，且V的列为CTC的正 交特征向量，VT表示V的转置矩阵。CCWPCTC的特征值相同，均为ai，a 2，一ar，其中r表示矩阵C 的秩。Σ为NXM的矩阵，且该矩阵中元素Σ ^的值为& : ，其余元素值为〇，元素Σ ^的值 按照大小顺序降序排列。
[0106]
>的奇异值分解υΣντ按照以上定义进行求解即可获得。S33、 根据迭代更新公式对低秩矩阵L、稀疏矩阵S和拉格朗日乘子矩阵X进行更新，同时更新补偿 正标量μ和迭代次数1^^1{+1 =以1^汰=1^1，其中，初始化的迭代条件如下:1^ = 5() = # = 〇4()> 0，P>1，k=0〇
[0107] S34、迭代执行步骤S33直到满足收敛条件后，获取当前的低秩矩阵和稀疏矩阵，其 中，I |D-Lk-sk| |F/| |D| IK10-7或迭代次数k达到预设最大次数，优选的，本实施例中，预设最 大次数为500次。
[0108] S4、将稀疏矩阵的每个列向量进行逆向量化后获得对应的WXH的图片，进而对每 个图片进行二值化分割从而得到对应的二值化图像。
[0109] S5、获取根据训练样本训练获得的频率分布图后，结合二值化图像计算每个测试 图片的缺陷度；
[0110] 本步骤中，频率分布图是通过以下方式训练获得的：
[0111] 建立一个WXH的矩阵并将该矩阵的所有元素值初始化为1后作为频率分布图，进 而依次获取训练样本的每张1C焊点训练图片的稀疏矩阵对应的二值化图像后，对频率分布 图进行更新，针对二值化图像中像素值为255的像素点，将频率分布图的相应位置的元素值 加1，W表示测试图片的宽度，Η表示测试图片的高度。
[0112]缺陷度的计算公式如下：
[0114] 上式中，Em表示测试图片的缺陷度，Μ表示测试集合的测试图片个数，f(x，y)表示 频率分布图，b(x，y)表示测试图片的二值化图像，W表示测试图片的宽度，Η表示测试图片的 高度。
[0115] S6、将该图片的缺陷度与训练样本的缺陷度阈值进行比对后获得1C元件焊点的检 测结果，具体如下：
[0116] 判断该缺陷度是否大于训练样本的缺陷度阈值，若是，则判断该1C元件焊点为虚 焊焊点，反之，判断该1C元件焊点为正常焊点。
[0117] 优选的，步骤S1之前还包括对训练样本进行训练的步骤，具体包括：
[0118] Α1、分别将训练样本的每张1C焊点训练图片向量化成列向量后，将获得的所有列 向量合并成1C焊点样本矩阵；训练样本包括多个1C焊点训练图片，每个1C焊点训练图片是 通过将合格的1C引脚图片由RGB颜色空间转换为HSV空间后获得的Η通道图像。
[0119] 本步骤与步骤S2相似，具体为:分别将训练样本的每张1C焊点训练图片向量化成Ν XI维列向量后，将获得的所有列向量合并成ΝΧΜ维的1C焊点测试矩阵;其中，Μ表示训练样 本的1C焊点训练图片的数量，Ν表示列向量的维度，且N = WXH，W表示测试图片的宽度，Η表 示测试图片的高度。
[0120] Α2、对1C焊点样本矩阵进行优化分解处理后获得其低秩矩阵和稀疏矩阵;获得低 秩矩阵和稀疏矩阵的详细过程与步骤S3中相同，两者只是处理对象略有区别，具体包括Α21 ~Α24:
[0121] Α21、建立如下的增广拉格朗日函数作为1C焊点样本矩阵的优化函数：
[0123] 其中，ML，S，X)表示优化函数，D表示1C焊点样本矩阵，L表示1C焊点样本矩阵的 低秩矩阵，S表示1C焊点样本矩阵的稀疏矩阵，X表示与1C焊点样本矩阵相同维度的拉格朗 日乘子矩阵，μ表示补偿正标量，I |L| I*表示低秩矩阵L的核范数，I |S| U表示稀疏矩阵3的1^ 范数，λ表示调节参数，I I · I If表示Frobenius范数，〈.，.>表示矩阵内积运算子；
[0124] A22、求解获得低秩矩阵L、稀疏矩阵S和拉格朗日乘子矩阵X的迭代更新公式如下：
[0126]其中，上标k表示迭代次数，U Σ ￥7表示矩阵
I的奇异值分解，? (.)表 示压缩算子，各符号的取值及含义如前所示。
[0127] A23、根据迭代更新公式对低秩矩阵L、稀疏矩阵S和拉格朗日乘子矩阵X进行更新， 同时更新补偿正标量μ和迭代次数k:y k+1 = ykP，k=k+Ι，其中，初始化的迭代条件如下:LQ = S°=X0 = 0,y°>0,P>l,k=0〇
[0128] A24、迭代执行步骤S33直到满足收敛条件后，获取当前的低秩矩阵和稀疏矩阵，其 中，收敛条件为：I |D-Lk-Sk| |F/| |D| |F〈1(T7或迭代次数k达到预设最大次数。
[0129] A3、将稀疏矩阵的每个列向量进行逆向量化后获得对应的样本图片，进而对每个 样本图片进行二值化分割从而得到对应的二值化图像；
[0130] A4、建立频率分布图后，依次获取每个二值化图像，对频率分布图进行更新，进而 在更新完毕后根据下式分别计算训练样本的每个1C焊点训练图片的缺陷度并获取最大值 作为训练样本的缺陷度阈值：
[0132] 上式中，Em表示测试图片的缺陷度，Μ表示测试集合的测试图片个数，f(x，y)表示 频率分布图，b(x，y)表示测试图片的二值化图像，W表示测试图片的宽度，Η表示测试图片的 高度。
[0133] 本发明中，优化分解获得低秩矩阵和稀疏矩阵的原理如下：
[0134] 具体优化问题为：
[0136] 其中，| |L| I*表示低秩矩阵L的核范数，即低秩矩阵L的特征值之和，| |S| U表示稀 疏矩阵S的Li范数，λ表示调节参数。
[0137] 上述问题是一个带等式约束条件的优化问题，运用基于拉格朗日乘子的迭代更新 法可以得到该问题的最优解。首先，该方法将上述优化问题转化为增广拉格朗日函数：
[0139] X表示与1C焊点测试矩阵相同维度(ΝΧΜ维）的拉格朗日乘子矩阵，μ表示补偿正标 量，I I · I |f表示Frobenius范数，〈·，·>表示矩阵内积运算子；
[0140]基于拉格朗日乘子的迭代更新法的基本思想是将其余变量当作是固定不变的情 况下（当作是常量），对另一个变量进行迭代更新，其更新方向是在缓慢减小。详细更新方法 如下：
[0145] 因此，获得低秩矩阵L的迭代更新公式。
[0146] 类似的，稀疏矩阵S的更新如下：
[0150] 拉格朗日乘子矩阵X的更新公式如下：
[0151] xk+1=xk+yk(D-Lk+1-S k+1)
[0152] 最后，当满足收敛条件：I |D-Lk-Sk| |F/| |D| |Κ1(ΤΜ、于10-7或迭代次数达到预设最 大次数时，停止迭代更新，得到最终收敛值，即1C焊点样本矩阵的低秩矩阵L和稀疏矩阵S。 根据上述原理即可获得步骤S3或Α2中的迭代更新步骤。
[0153] 本发明通过求解稀疏矩阵S将1C元件焊点检测扩展到最优化的层面，因此，任何属 于背景部分的变化，即任何1C元件焊点图片之间的不同部分，都能够通过低秩矩阵L捕捉到 以及进行约束(最终提高检测效果），最后通过对应的稀疏矩阵S计算缺陷度与缺陷度阈值 进行比对后，即可获得1C元件焊点的检测结果，本方法可以准确、有效、快速地进行焊点缺 陷检测。
[0154]以上是对本发明的较佳实施进行了具体说明，但本发明创造并不限于实施例，熟 悉本领域的技术人员在不违背本发明精神的前提下还可做出种种的等同变形或替换，这些 等同的变型或替换均包含在本申请权利要求所限定的范围内。
【主权项】
1. 一种基于鲁棒主成分分析的IC元件焊点检测方法，其特征在于，包括步骤： 51、 采集多个待检测1C元件焊点的测试图片后，形成测试集合； 52、 分别将测试集合的每张测试图片向量化成列向量后，将获得的所有列向量合并成 1C焊点测试矩阵； 53、 对1C焊点测试矩阵进行优化分解处理后获得其低秩矩阵和稀疏矩阵； 54、 将稀疏矩阵的每个列向量进行逆向量化后获得对应的图片，进而对每个图片进行 二值化分割从而得到对应的二值化图像； 55、 获取根据训练样本训练获得的频率分布图后，结合二值化图像计算每个测试图片 的缺陷度； 56、 将该图片的缺陷度与训练样本的缺陷度阔值进行比对后获得1C元件焊点的检测结 果。2. 根据权利要求1所述的一种基于鲁棒主成分分析的1C元件焊点检测方法， 其特征在于，所述步骤S1之前还包括对训练样本进行训练的步骤，具体包括： A1、分别将训练样本的每张1C焊点训练图片向量化成列向量后，将获得的所有列向量 合并成1C焊点样本矩阵； A2、对1C焊点样本矩阵进行优化分解处理后获得其低秩矩阵和稀疏矩阵； A3、将稀疏矩阵的每个列向量进行逆向量化后获得对应的样本图片，进而对每个样本 图片进行二值化分割从而得到对应的二值化图像； A4、建立频率分布图后，依次获取每个二值化图像，对频率分布图进行更新，进而在更 新完毕后计算训练样本的缺陷度阔值。3. 根据权利要求1所述的一种基于鲁棒主成分分析的1C元件焊点检测方法， 其特征在于，所述步骤S2，其具体为： 分别将测试集合的每张测试图片向量化成NX 1维列向量后，将获得的所有列向量合并 成NXM维的1C焊点测试矩阵；其中，Μ表示测试集合的测试图片的数量，N表示列向量的维 度，且N=WXH，W表示测试图片的宽度，Η表示测试图片的高度。4. 根据权利要求1所述的一种基于鲁棒主成分分析的1C元件焊点检测方法， 其特征在于，所述步骤S3，包括： 531、 建立如下的增广拉格朗日函数作为1C焊点测试矩阵的优化函数： .心(￡.，&.巧=11 LII* +λ. I间I, +昔1_正 _S Ik. + < 义，'0_王-S > 其中，l^LiSiX)表示优化函数，D表示1C焊点测试矩阵，L表示1C焊点测试矩阵的低秩矩 阵，S表示1C焊点测试矩阵的稀疏矩阵，X表示与1C焊点测试矩阵相同维度的拉格朗日乘子 矩阵，μ表示补偿正标量，II L II *表示低秩矩阵L的核范数，II S II1表示稀疏矩阵5的^范数，λ 表示调节参数，I I · IIf表示Frobenius范数，<.，.〉表示矩阵内积运算子； 532、 求解获得低秩矩阵L、稀疏矩阵S和拉格朗日乘子矩阵X的迭代更新公式如下：其中，上标k表示迭代次数，u Σ八表示矩阵的奇异值分解，? Ο表示压 缩算子且表示矩阵Σ中元素如果大于^，则取原值，否则，取0; μ* μ 533、 根据迭代更新公式对低秩矩阵L、稀疏矩阵S和拉格朗日乘子矩阵X进行更新，同时 更新补偿正标量μ和迭代次数k: yk+i=ykp，k=k+1，其中，P表示迭代参数； 534、 迭代执行步骤S33直到满足收敛条件后，获取当前的低秩矩阵和稀疏矩阵。5. 根据权利要求4所述的一种基于鲁棒主成分分析的1C元件焊点检测方法， 其特征在于，所述步骤S33，其具体为： 根据迭代更新公式对低秩矩阵L、稀疏矩阵S和拉格朗日乘子矩阵X进行更新，同时更新 补偿正标量4和迭代次数4:4^ = 4*^0，4 = 4+1，其中，初始化的迭代条件如下:1^*^ = 5*^ = 乂*^ = 〇，μ°〉〇，ρ〉1，k = 0。6. 根据权利要求4所述的一种基于鲁棒主成分分析的IC元件焊点检测方法，其特征在 于，所述步骤S34中所述收敛条件为：II D-Lk-Sk II f/Μ D II F<1(T7或迭代次数k达到预设最大 次数。7. 根据权利要求1所述的一种基于鲁棒主成分分析的1C元件焊点检测方法， 其特征在于，所述步骤S4，其具体为： 将稀疏矩阵的每个列向量进行逆向量化后获得对应的WXH的图片，进而对每个图片进 行二值化分割从而得到对应的二值化图像，W表示测试图片的宽度，Η表示测试图片的高度。8. 根据权利要求1所述的一种基于鲁棒主成分分析的1C元件焊点检测方法， 其特征在于，所述步骤S5中所述频率分布图是通过W下方式训练获得的： 建立一个WXH的矩阵并将该矩阵的所有元素值初始化为1后作为频率分布图，进而依 次获取训练样本的每张1C焊点训练图片的稀疏矩阵对应的二值化图像后，对频率分布图进 行更新，针对二值化图像中像素值为255的像素点，将频率分布图的相应位置的元素值加1， W表示测试图片的宽度，Η表示测试图片的高度。9. 根据权利要求1所述的一种基于鲁棒主成分分析的1C元件焊点检测方法，其特征在 于，所述步骤S5中所述计算每个测试图片的缺陷度的步骤，是通过下式计算缺陷度的：上式中，Em表示测试图片的缺陷度，Μ表示测试集合的测试图片个数，f(x，y)表示频率分 布图，b(x，y)表示测试图片的二值化图像，W表示测试图片的宽度，Η表示测试图片的高度。10. 根据权利要求1所述的一种基于鲁棒主成分分析的1C元件焊点检测方法， 其特征在于，所述步骤S6，其具体为： 判断该缺陷度是否大于训练样本的缺陷度阔值，若是，则判断该1C元件焊点为虚焊焊 点，反之，判断该1C元件焊点为正常焊点。
【文档编号】G01N21/956GK106093074SQ201610428079
【公开日】2016年11月9日
【申请日】2016年6月16日
【发明人】蔡念, 周杨, 叶倩, 刘根, 王晗, 翁韶伟
【申请人】广东工业大学