一种用于pcb视觉定位的检测方法
【专利摘要】一种用于PCB视觉定位的检测方法。本发明涉及PCB的视觉定位与检测领域。本发明是为了解决传统PCB检测方法定位精度低和通用性差的问题。本发明主要技术为采用贴片机视觉系统获取PCB定位标志图像，进行阀值分割得到二值化预处理后的图像，对PCB上的定位标志进行形状判别，然后采用先粗后精的边缘检测方法得到定位标志边缘点的亚像素坐标，计算定位标志中心点的坐标。本发明主要用于PCB上定位标志的质心检测。
【专利说明】-种用于PCB视觉定位的检测方法
【技术领域】
[0001] 本发明涉及PCB的视觉定位和检测领域，具体是一种用于PCB视觉定位的检测方 法。
【背景技术】
[0002] PCB印刷电路板定位是通过对母板上的定位标志进行检测，从而得到PCB坐标系 与贴片机机器坐标系之间的转换关系。快速、准确地完成PCB上定位标志的对准对提高贴 片机的工作效率和精度起着至关重要的作用。已有的对PCB的检测，基于Hough变换的方 法来对十字形定位标志的四条长边缘直线进行抽取，根据长边缘直线的4个交点的坐标平 均值可以得到十字形的中心点。但该方法采用的Hough变换计算速度较慢，采用的Sobel 边缘检测算子检测精度只能达到像素级。也有算法针对定位标志的特点，提出了结合图像 分割、几何特征识别与模式匹配的PCB定位方法。首先采用联通域分割算法对定位标志进 行图像分割以减少匹配模板的面积，然后根据几何特征识别方法进行相应的粗定位，以确 定定位标志中心点的大致坐标，最后将模板在近似值附近邻域内进行匹配以寻找到目标图 像与模板相似度最大的点的位置。当定位标志图像较大时，模板匹配时耗时较长，该方法的 实时性将变得较差。
【发明内容】
[0003] 本发明是为了解决传统的PCB检测方法定位精度低，研究对象单一，通用性差，而 提出一种用于PCB视觉定位的检测方法。
[0004] 本发明的技术方案如下：
[0005] -种用于PCB视觉定位的检测方法按以下步骤实现：
[0006] 步骤一：采用贴片机视觉系统获取PCB的图像。
[0007] 步骤二：将步骤一得到的图像采用最大类间方差法进行阈值分割，得到二值化预 处理后的图像。
[0008] 步骤三：如果定位标志图像是中空的，采用形态学膨胀运算对中空区域进行填充。 然后对图像中的黑点进行最小二乘圆拟合，从而计算得到面积率A和圆形度S。
[0009] 区域进行填充的步骤为：
[0010] 初始化：随机选取原始图像中某个连通区域内的一个亮点B。。
[0011] 循环：采用十字形结构元S对仏进行膨胀运算，膨胀后的结果与原始图像的补图 像相交得到&，重复进行膨胀运算与相交运算，直到Bi+1 ==Bp最终的膨胀结果与原始图 像取并集得到实心图像。
[0012] 面积率入和圆形度S的定义为：
【权利要求】
1. 一种用于PCB视觉定位的检测方法，其特征在于它包括下述步骤： 步骤一：采用贴片机视觉系统获取PCB的图像； 步骤二：将步骤一得到的图像采用最大类间方差法进行阈值分割，得到二值化预处理 后的图像； 步骤三：如果定位标志图像是中空的，采用形态学膨胀运算对中空区域进行填充。然后 对图像中的黑点进行最小二乘圆拟合，从而计算得到面积率A和圆形度S; 区域进行填充的步骤为： 初始化：随机选取原始图像中某个连通区域内的一个亮点B。； 循环：采用十字形结构元S对^进行膨胀运算，膨胀后的结果与原始图像的补图像相 交得到&，重复进行膨胀运算与相交运算，直到Bi+1 ==Bp最终的膨胀结果与原始图像取 并集得到实心图像； 面积率A和圆形度S的定义为：
其中，n为拟合圆内的黑点数，N为目标图像的总黑点数，S为二值图像中的黑点总数，L为二值图像中黑色区域的边界长度； 步骤四：判断定位标志的形状，判断顺序依次为：圆形、正方形、三角形、格子形以及 十字形。若入>0.95且6>1，则为圆形；若入>0.8且S>〇. 7,则为正方形；若入>0.6且 S>0. 7,则为三角形；若入>0. 4且S>〇. 3,则为格子；若入>0. 3且S>〇. 1，则为十字形； 步骤五：采用Sobel边缘检测算子对定位标志图像进行边缘检测，得到定位标志边缘 点的像素级坐标； 步骤六：采用基于空间矩的亚像素边缘检测法对步骤五得到的像素级坐标进行精确定 位，得到定位标志边缘点的亚像素级坐标； 二维函数f(x，y)的p+q空间矩定义为： Mpq = / /xpy9f(x,y)dydx其中中，P、q为非负整数。 步骤七：对步骤六得到的亚像素坐标点进行最小二乘圆拟合，得到定位标志的中心点 坐标； 步骤八：利用步骤七得到的中心点坐标定位PCB，从而计算出PCB坐标系和机器坐标系 之间的转换关系。
2. 根据权利要求1所述的一种用于PCB视觉定位的检测方法，其特征在于步骤六所述 的基于空间矩的亚像素边缘检测法，具体按照以下步骤实现： (1)根据二级理想边缘模型，利用5X5的模板计算每个像素级边缘点的6个空间矩M〇〇>M01,M10,Mg2,M2〇,Mn ； ⑵计算旋转后的6个空间矩M'QQ，M'Q1，M' 1Q，M' Q2,M' 2Q，M'n ; (3)计算二级理想边缘模型的边缘参数；
(4) 采用最大类间方差法求得最佳阀值it、kt。若k>ktn1彡it，则该像素点为亚像 素边缘点； (5) 计算边缘点的亚像素坐标；
其中，（xs，ys)为亚像素边缘点的坐标，（x，y)为像素级边缘点的坐标； (6) 计算边缘交点的判别值e;
选取适当的阀值T,若|e| >T,则边缘点为交点；若|e|〈T,则边缘点为普通边缘点； (7) 若为边缘交点，对像素间的亚像素进行双线性插值，得到亚像素级的灰度值；采用 共轭梯度法对边缘交点的亚像素坐标进行寻优；经过多次迭代即可准确地得到交点的亚像 素坐标； 目标函数的表达式为：
其中，（1，7)的初始值为交点，(1。!￡(1，7)、(1。7(1，7)是函数(：(1，7)的一阶偏导数 ； (8) 若不为交点附近的边缘点，根据三级理想边缘模型，利用5X5的模板计算每个像 素级边缘点的6个空间矩M'QQ，M'Q1，M' 1Q，M'Q2，M' 2Q，M'n ; 计算旋转后的6个空间矩M"QQ，M"Q1，M%，M" Q2,M" 2Q，M"n ; 计算角度为P的边缘的亚像素位置le;
式中，X=Ak/k,f(l) = (l-l2)1.5; 计算原理误差E; E=l-le (10)； 利用式（10)对式（2)计算出的1进行补偿； (9) 计算边缘点的亚像素坐标；
(10) 若为交点附近的边缘点，根据交点处的边缘模型，利用5X5的模板计算每个像素 级边缘点旋转后的6个空间矩M'QQ，M'Q1，M' 1Q，M' Q2,M' 2Q，M'n; 在交点处，6个空间矩满足下述关系： C(x,y) = 2(M,20+M,02)-M,00 = 0 (式 12);
在普通边缘点出，6个空间矩满足下述关系：
得到边缘交点的判断公式为：；
选取适当的阀值T,若| 5 |彡T,则边缘点为交点；若| 5 |〈T,则边缘点为普通边缘点； 对像素间的亚像素进行双线性插值，得到亚像素级的灰度值。采用共轭梯度法对交点 的亚像素坐标进行寻优，目标函数的表达式为：
式中，（x，y)的初始值为交点，y)是函数C(x，y)的一阶偏导数；经过 多次迭代即可准确地得到交点的亚像素坐标。
【文档编号】G06K9/00GK104408722SQ201410693065
【公开日】2015年3月11日 申请日期:2014年11月26日 优先权日:2014年11月26日
【发明者】黄筱调, 王祖进, 丁爽 申请人:南京工业大学