技术领域
[0001] 本
发明涉及PCB板的自动在线缺陷检测，具体涉及一种基于投影法和连通域分析的PCB板缺陷检测方法。
背景技术
[0002] 印制
电路板(Printed Circuit Board，简称PCB),是各种
电子元器件的载体，在电子通信等许多领域中都有着广泛的应用，需求量巨大。随着工艺的提高，电子技术的发展，
电路板向小元件、高
密度、细间距方向发展。由于其结构复杂、生产工序繁多，出现断路、
短路、线细等缺陷不可避免。传统的人工目测、电子检测等方法因其检测速度慢，易漏检，且
接触法检测易损伤PCB等缺点，已经不能满足生产的需要。
自动光学检测技术(AOI)在PCB的缺陷检测上应用越来越多，也出现了很多商用的检测设备。但大都价格昂贵。如何经济有效的实现PCB的自动检测，成为中小企业迫切需要解决的问题。
[0003]
机器视觉检测
算法是自动光学检测的关键技术，国内外学者也提出了不少算法，主要包括三类：参考比较法、设计准则校验法和混合法。
[0004] ⑴参考比较法是将待检测PCB图像和标准PCB图像逐点进行比较，或者将待检测PCB图像上提取的特征与标准PCB图像所提取的特征进行比较，有差异即认为有缺陷，常见的算法是通过差影、异或运算识别缺陷，通过统计待
检测区域的连通域个数及面积识别缺陷等。这种方法速度快，
硬件实现容易，但光照和
定位的要求比较高，若待检测图像和参考图像没有精确配准，将产生虚假警报。
[0005] ⑵设计准则校验法是通过检测待检测PCB板子是否满足预先定义的设计准则来判断其是否有缺陷。这种方法的优点是不需要标准PCB图像的信息，可以降低对存储空间的需求，不需要配准就能进行常见缺陷的检测，但是运算复杂，计算量大，还必须根据设计准则设计出相应的图案数据结构，只能检测出不符合设计准则的缺陷。
[0006] ⑶混合法是将以上两种方法结合起来对PCB进行检测，希望综合两者的优点来实现检测。
发明内容
[0007] 本发明的目的是对尚未安装电子元器件的具有重复单元的PCB进行缺陷检测。对于大规格的PCB板，为了满足低成本高
精度的检测要求，需要多个面阵相机同时进行
图像采集。若先完成图像拼接再整板进行缺陷检测，数据量大，运算时间长，无法满足在线检测的要求。因此，考虑将PCB上的重复单元分割出来，逐个进行缺陷检测，以减少图像拼接。由于各单元
块的间隙很小，实现单元块的分割是本发明要解决的一个重要问题。完成单元分割后，对单元块进行缺陷检测。PCB板上单元块的电路类型多样，要实现准确快速的在线检测，就要求检测算法准确且稳定可靠，但又不能太过复杂，具有较好的通用性。
[0008] 本发明对采集到的PCB图像，先进行预处理，再通过投影法实现单元块的分割，最后与标准单元块进行比较，实现缺陷检测，主要检测开路和短路等缺陷。
[0009] 本发明提供一种基于投影法和连通域分析的PCB板缺陷检测方法，主要包括以下步骤：
[0010] 步骤1，采集待检测PCB板灰度图像，并对其进行预处理；
[0011] 步骤2，利用投影法对预处理后的图像进行单元分割，即将PCB板上重复的单元块独立分割开；
[0012] 步骤3，缺陷检测，统计无缺陷模板中单元块的连通域数量Num及面积，并将各连通域面积按由大到小排列，记作序列Area，将连通域数量Num和面积序列Area作为特征，再统计待检测图像中各个单元块的相应参数，将其逐个与模板的标准参数比对，判断是否有缺陷及缺陷类型。
[0013] 进一步的，所述步骤3缺陷检测的具体实现方式如下，
[0014] 步骤3.1，统计无缺陷模板中单元块的连通域数量Num及面积，并将各连通域面积按由大到小排列，记作序列Area；
[0015] 步骤3.2，统计待检测图像中各个单元块的相应参数NumD和AreaD，将其逐个与模板的标准参数Num和Area比对，判断是否有缺陷及缺陷类型；判断方式如下，[0016] ①若待检测单元块的连通域数量大于标准参数，即NumD＞Num，则缺陷为断路；
[0017] ②若待检测单元块的连通域数量小于标准参数，即NumD＜Num，则缺陷为短路；
[0018] ③若待检测单元块的连通域数量等于标准参数，即NumD＝Num，继续对比面积，若各连通域面积的相对误差不超过±δ，则判定无缺陷，否则，也标记为有缺陷。
[0019] 进一步的，所述步骤1中的预处理包括如下子步骤，
[0020] 步骤1.1，对采集到的灰度图像f进行伽
马变换得到fg，
[0021] fg(x,y)＝c·f(x,y)γ
[0022] 其中，x,y表示
像素点坐标，c和γ为常数参数；
[0023] 步骤1.2，对伽马变换后的图像用Otsu的方法计算得到全局
阈值Th，得到二值图像fb，
[0024]
[0025] 进一步的，所述步骤2中利用投影法对预处理后的图像进行单元分割的具体实现方式如下，
[0026] 步骤2.1，采用张正友相机标定法对系统进行标定，获得相机的内外参数，利用这些参数对二值化后的图像fb进行畸变校正得到fc；
[0027] 步骤2.2，对校正后的图像fc减半
采样，得到fch，
[0028] fch(m,n)＝fc(2m-1,2n-1)
[0029] 其中，m＝1,2,...,x/2；n＝1,2,...,y/2，x,y表示像素点坐标；
[0030] 步骤2.3，对图像fch做Radon变换求出待检测PCB板的倾斜
角度θ，然后根据倾斜角度θ将二值图像fc旋正，得到fbr；
[0031] 步骤2.4，对旋正后的二值图像fbr求其行投影V(i)和列投影H(j)；
[0032]
[0033] 步骤2.5，设置合适的阈值ThV和ThH，然后通过截底将位于底部较低的投影截去，[0034]
[0035] 步骤2.6，搜索截底后的投影结果VN(i)和HN(j)，找到的波谷
位置即为图像fbr中单元块间隙的位置，实现单元块的分割。
[0036] 进一步的，所述步骤2.3中获得倾斜角度θ的具体实现方式为，在-1°～1°的范围内，每隔0.1°对图像fch进行Radon变换，获得从-1°～1°，每隔0.1°，共21个角度对应的投影图，统计各投影图中像素值为0的点的数量，数量最多的投影图对应的角度即为PCB板的倾斜角度θ。
[0037] 进一步的，所述的步骤3.2中δ＝0.8％。
[0038] 进一步的，所述的步骤1.1中的参数c＝1，γ＝0.5。
[0039] 与
现有技术相比，本发明的优点和有益效果如下：
[0040] 1.在进行缺陷检测之前，先对采集到的图像进行单元分割，解决了大规格PCB板检测时先拼接再进行缺陷检测，数据量大，运算时间长的问题；
[0041] 2.本发明的单元分割采用投影法，算法简单，速度快，可靠性好；
[0042] 3.投影之前，先对图像减半采样，再通过Randon变换实现PCB板子的倾斜矫正，提高检测速度；
[0043] 4.缺陷检测时采用先对比待检测
图像分割单元块与无缺陷模板单元块中连通域的数量，再对比连通域的面积序列，检测PCB板上的短路和断路缺陷，速度快，鲁棒性好。
[0044] 实验表明，该发明对于由多个重复单元组成的PCB裸板上的短路及开路等缺陷具有很好的检测效果，能检测到的最小缺陷为0.15mm，检测速度最快可达10m/min，检测精度高、速度快，成本低，具有很好的实际应用价值。
附图说明
[0046] 图2为本发明实施例中采集到的PCB板图像；
[0047] 图3为本发明实施例中预处理后图像；
[0048] 图4为本发明实施例中倾斜校正后图像的行投影和列投影，其中，(a)行投影，(b)列投影；
[0049] 图5为本发明实施例中截底后的行投影和列投影，其中，(a)行投影，(b)列投影；
[0050] 图6为本发明实施例中PCB板单元分割结果；
[0051] 图7为本发明实施例中标准单元块；
[0052] 图8为本发明实施例中待检测单元块，其中，(a)断路，图中椭圆线框内为缺陷，图下方为缺陷的放大图；(b)短路，图中椭圆线框内为缺陷，图下方为缺陷的放大图；(c)无缺陷，(d)有缺陷，图中椭圆线框内为缺陷，图下方为缺陷的放大图。
具体实施方式
[0053] 下面结合附图和实施例对本发明的技术方案作进一步说明。
[0054] 参见图1，本发明所描述的一种基于投影法和连通域分析的PCB板缺陷检测方法，包括以下步骤：
[0055] 步骤1，采集待检测PCB板灰度图像f，并对其进行预处理。其过程如下：
[0056] 步骤1.1，首先进行图像采集，所用的是华睿A3A20CG8CMOS面阵相机，该相机
分辨率为1200万像素，能够满足检测精度的要求，采集到的原始PCB板图像如图2所示。对采集到的灰度图像f进行伽马变换得到fg：
[0057] fg(x,y)＝c·f(x,y)γ
[0058] 其中，x,y表示像素点坐标，参数c＝1，γ＝0.5。
[0059] 步骤1.2，对伽马变换后的图像用Otsu的方法计算得到全局阈值Th，得到二值图像fb：
[0060]
[0061] 伽马变换矫正了幂律响应现象，增强了图像
对比度和细节特征，之后再使用最大类间方差法(Otsu)实现了图像的二值化，将电路从PCB板上分割出来，效果如图3所示。
[0062] 步骤2，对预处理后的图像fb进行单元分割，即将PCB板上的各个重复单元块独立的分割开，本实施例中采用投影法来实现单元分割。
[0063] 图像投影是将二值图像在某一方向上做线性积分，对于数字图像即为累加求和。通过
水平和垂直方向的投影，可将二维图像的像素分布特征简化为x轴和y轴上的两个一维函数，可以获得图像在这两个方向上的特性。
[0064] 投影法就是根据这两个一维函数反映的特征来判定图像的切割位置坐标，再用这个坐标来切割原图像，最终实现单元分割的。
[0065] 本发明中，PCB板上的重复单元呈整齐的行列排列，其二值化后的图像背景为黑色，目标物为白色。若将二值化后的图像旋正，再分别求取行投影和列投影，那么投影的波谷处应恰好为单元块间的间隙位置。
[0066] 具体实施步骤如下：
[0067] 步骤2.1，由于采集到的图像有畸变，在分割前先采用张正友相机标定法对系统进行标定，获得相机的内外参数，利用这些参数对二值化后的图像fb进行畸变校正得到fc，这部分为现有技术，本发明不予撰述，具体可参见文献Zhang,Z.(2000)A Flexible New Technique for Camera Calibration.IEEE Transactions on PatternAnalysis and Machine Intelligence,22,1330-1334.http://dx.doi.org/10.1109/34.888718；
[0068] 步骤2.2，为了提高效率，在求板子倾斜角度前先对校正后的图像fc减半采样,得到fch：
[0069] fch(m,n)＝fc(2m-1,2n-1)
[0070] 其中，m＝1,2,...,x/2；n＝1,2,...,y/2。
[0071] 步骤2.3，倾斜校正。对图像fch做Radon变换求出待检测板子的倾斜角度θ，然后将二值图像fc旋正，得到fbr。
[0072] 经过机械矫正装置的
姿态调整，传送带上的PCB板在到达相机视场下时的倾斜角度在-1°～1°之间，那么在这个范围内，每隔0.1°对图像fch进行Radon变换，即可获得从-1°～1°，每隔0.1°，共21个角度对应的投影图，统计各投影图中像素值为0的点的数量，数量最多的投影图对应的角度即为板子的倾斜角度θ。之后，将图像fc转过-θ得到fbr，板子图像即被旋正。本例中，求出θ＝-1°。
[0073] 步骤2.4，对旋正后的二值图像fbr求其行投影V(i)和列投影H(j)；
[0074]
[0075] 得到的投影图如图4所示，其中(a)为行投影，横坐标为图像的行，纵坐标为图像每行各点的像素值之和；(b)为列投影，横坐标为图像的列，纵坐标为图像每列各点的像素值之和。
[0076] 步骤2.5，由于单元块的缝隙处可能出现噪声、余
铜和标记等，从图4中可以看到缝隙处灰度值为零的点可能不存在或是数量较少，这将给投影法分割带来困难。因此，在分割前还需对投影图进行处理，即设置合适的阈值ThV和ThH(本实施例中ThV和ThH由实验得到的，ThV＝ThH＝100)，然后通过截底将位于底部较低的投影截去：
[0077]
[0078] 截底后投影图如图5所示，其中(a)为行投影，横坐标为图像的行，纵坐标为图像每行各点的像素值之和；(b)为列投影，横坐标为图像的列，纵坐标为图像每列各点的像素值之和。截底之后，若某行或某列的像素值之和小于阈值ThV或ThH，即认为该行或该列出现的是噪声、余铜或标记，而不是单元块的内容，将投影值置为零，则消除了这些因素的影响，便于单元块的分割。
[0079] 步骤2.6，搜索截底后的投影结果VN(i)和HN(j)，找到的波谷位置即为图像fbr中单元块间隙的位置，实现单元块的分割，分割结果如图6所示。
[0080] 步骤3，缺陷检测。将分割出的单元块与标准单元块进行比较，实现缺陷检测。
[0081] 步骤3.1，统计无缺陷模板中单元块(如图7)的连通域数量Num和面积，并将各连通域面积按由大到小排列，记作Area。本例中，连通域数量Num＝10，面积序列Area＝[172699,129637，53646，53039，38744，29135，22785，21611，13814，13498]；
[0082] 步骤3.2，统计待检测图像中各个单元块的相应参数NumD和AreaD，将其逐个与模板的标准参数Num和Area比对，判断是否有缺陷及缺陷类型。
[0083] ①若待检测单元块的连通域数量大于标准参数(如图8(a)，NumD1＝11)，即NumD1＞Num，则缺陷为断路；
[0084] ②若待检测单元块的连通域数量小于标准参数(如图8(b)，NumD2＝9)，即NumD2＜Num，则缺陷为短路；
[0085] ③若待检测单元块的连通域数量等于标准参数(如图8(c)和(d)，Num D3＝Num D4＝10)，即NumD3＝NumD4＝Num，继续对比其面积；若其面积的相对误差不超过±0.8％，则判定无缺陷。否则，也标记为有缺陷。本例图8(c)中Area3＝[171491,129118，53377，52720，38589，29018，22693，21546，13747，13450]，计算得各连通域面积的相对误差δ3max＝0.7，即δ3max＜δ，因此单元块3无缺陷。
[0086] 图8(d)中，Area4＝[155843,129563，53392，53015，38637，29125，27384，22727，21609，16778]，计算得各连通域面积的相对误差δ4max＞0.8，即δ4max＞δ，因此单元块4有缺陷。
[0087] 本文中所描述的具体实施例仅仅是对本发明精神作举例说明。本发明所属技术领域的技术人员可以对所描述的具体实施例做各种各样的
修改或补充或采用类似的方式替代，但并不会偏离本发明的精神或者超越所附
权利要求书所定义的范围。
