基于亚像素边缘算法的平面零件尺寸测量方法
【专利摘要】本发明涉及基于亚像素边缘算法的平面零件尺寸测量方法，利用CCD设备作为测量平台，包括下列步骤：1、钣金零件成像特性分析；2、Canny边缘粗定位；3、基于多项式拟合的亚像素检测；4、区分被测物上下边缘；5三角法求取被测物准确边缘；6、被测物尺寸恢复，突破了目前基于CCD的尺寸测量技术中多用于小尺寸物体测量的局限，提出了亚像素精度边缘检测算法，首先，结合光学成像理论分析了背光光源下不同厚度钣金零件的边缘特性，提出面向相机一侧有上下边缘，然后根据钣金零件图像的边缘分布特征，采用Canny边缘检测和多项式拟合算法实现亚像素边缘提取，提出平均距离法区分零件上下边缘，进而利用三角法解算获取零件真实边缘，最终利用成像原理和直线与零件表面交点确定零件尺寸。
【专利说明】基于亚像轰边缘算法的平面零件尺寸测量方法
【技术领域】
[0001] 本发明属于图像处理技术（机器视觉）领域，涉及一种基于亚像素边缘算法的平 面零件尺寸测量方法，可W实现对加工后的机械零件进行高精度、高效率的检测。
【背景技术】
[0002] 检测技术是现代制造业的基础技术之一，是保证产品质量的关键。在饭金加工领 域，零件尺寸的在线高精度检测已经成为一个重要研究课题。机器视觉测量技术具有非接 触、速度快、高精度和环境适应性好等优点，能够很好地满足现代工业在线测量的要求，已 经成为先进制造领域内最先进的在线测量技术之一，而边缘检测是机器视觉检测技术的基 础和关键，边缘的定位精度直接影响到尺寸检测的精度。
[0003] 传统的边缘检测算法仅能达到像素级精度，不能很好地满足测量精度要求，因此 对像素级的精度再细分即亚像素级的边缘定位技术研究十分必要。在视觉检测系统中， 亚像素边缘的检测和定位精度是其最终测量精度的关键。亚像素边缘检测技术最早由 化eckel提出，现已发展为插值法、矩法和拟合法等多种方法。它们的边缘检测精度、抗噪声 能力W及算法复杂度各不相同。插值法定位依赖于相邻像素间的灰度、坐标关系和插值函 数的选取，计算时间短，但定位精度低。基于矩的亚像素技术利用的是积分算子，对噪声不 敏感，但是通常矩方法缺少明确的准则来区分像素是否属于边缘，另外对包含梯度变化的 像素集都会产生响应并且只在边缘像素的近邻域有效，如果边缘位置比积分区域远，基于 矩的方法就会失效。王希军，陈阔等人分析了不同亚像素定位算法受各方面因素的影响情 况，分别对几类亚像素边缘定位算法的稳定性W及定位精度的评价方法进行了研究。在被 测物有一定的厚度且尺寸较大时，在距离相机一侧时会有上下边缘，而在相机另一侧时，只 有上边缘，因此边缘信息需要判别上边缘与下边缘，Antonis.Jan提出了利用法线角和参考 角的差值关系判别上下边缘，但是在差值临界区会出现判别错误，并且应用在圆形轮廓中 会出现边缘不连续现象。
[0004] 针对饭金零件亚像素边缘准确提取困难的问题，本文通过分析不同厚度饭金零件 的边缘特性，采用Canny边缘检测和多项式拟合算法实现亚像素边缘提取，提出平均距离 法区分零件上下边缘，进而利用H角法解算获取零件真实边缘，从而测量零件尺寸。并通过 实验验证了本文的亚像素准确边缘提取算法的可行性。
【发明内容】
[0005] 本发明的目的就是为克服现有技术的不足，提供一种亚像素边缘算法的平面零件 尺寸测量方法，测量精度在1X0. 75m2的视场内，系统测量精度达到了 0. 05mm，关键技术可 W提高检测机械零件的精度和效率，避免了传统人工检测时所带来的主观性、易疲劳、速度 慢、成本高、强度大等缺点。特别是对金属板材精密裁切提供理论和技术上的支持。
[0006] 本发明是通过该样的技术方案实现的；基于亚像素边缘算法的平面零件尺寸测量 方法，利用CCD设备作为测量平台，其特征在于，包括下列步骤：
[0007] (I)结合光学成像理论分析了背光光源下不同厚度饭金零件的边缘特性；
[0008] (2)根据饭金零件图像的边缘分布特征，采用Canny边缘粗定位；
[0009] (3)基于多项式拟合的亚像素检测；
[0010] (4)通过平均距离法区分被测物上下边缘；
[0011] (5)采用H角法求取饭金零件准确边缘；
[0012] (6)利用成像原理和直线与零件表面交点确定被测物尺寸。
[0013] 本发明的实质性特点是，基于亚像素边缘算法的平面零件尺寸测量方法，该模型 综合运用了基于多项式拟合的亚像素检测，平均距离法区分零件上下边缘等技术。首先， 结合光学成像理论分析了背光光源下不同厚度饭金零件的边缘特性，提出面向相机一侧有 上下边缘。然后根据饭金零件图像的边缘分布特征，采用Canny边缘检测和多项式拟合算 法实现亚像素边缘提取，提出平均距离法区分零件上下边缘，进而利用H角法解算获取零 件真实边缘，最终利用成像原理和直线与零件表面交点确定零件尺寸。实验结果表明，在 1X0. 75m2的视场内，系统测量精度达到了 0. 05mm。本发明与现有技术相比较具有如下优 占. ；、、、?
[0014] 1.结合光学成像理论分析了背光光源下不同厚度饭金零件的边缘特性，提出面向 相机一侧有上下边缘。
[0015] 2.提出平均距离法区分零件上下边缘。进而利用H角法解算获取零件真实边缘， 最终利用成像原理和直线与零件表面交点确定零件尺寸。
[0016] 3.本发明在计算边缘时采用亚像素的边缘提取，精度较高。
【专利附图】
【附图说明】
[0017] 图1 ;本发明的方法流程图；
[001引图2;边缘灰度提取图；
[001引 图3 ;H角法图。
【具体实施方式】
[0020] 本发明的流程图如图1所示，首先结合光学成像理论分析了背光光源下不同厚度 饭金零件的边缘特性，提出面向相机一侧有上下边缘。然后根据饭金零件图像的边缘分布 特征，采用Canny边缘检测和多项式拟合算法实现亚像素边缘提取，提出平均距离法区分 零件上下边缘。之后进而利用H角法解算获取零件真实边缘，最终利用成像原理和直线与 零件表面交点确定零件尺寸。下面结合附图，对本发明技术方案的具体实施过程加W说明。
[0021] 1.饭金零件的边缘特性分析
[0022] 实际的CCD成像系统中，CCD感光元不但接收照射到自身感光面的光，还接收照射 到相邻感光元的光，尤其是对边缘点，物体和背景的不同反射特性W及CCD器件的积分效 应，造成CCD器件对阶跃边缘的响应产生由明到暗或由暗到明的渐变过程，所W边缘在图 像中表征为一种灰度分布，如图2所示；另外，当被测物有一定的厚度且尺寸较大时，在距 离相机一侧存在两个梯度值较大的过渡带，而在相机另一侧时，只有上边缘；因此边缘信息 需要区分上边缘与下边缘，并只提取对相机都看见的上边缘；然后通过Canny算子进行整 像素级的边缘提取，得到整像素级的边缘后用H次多项式拟合的方法提取亚像素边缘精确 定位，再通过平均距离法确定饭金零件的上边缘，通过H角法最终确定零件的实际边缘，最 后利用成像原理和直线与零件表面交点确定被测物尺寸。
[0023] 2.Canny算子粗定位边缘
[0024] 亚像素边缘定位前必须先采用像素级边缘检测方法确定边缘点的位置，然后根据 边缘点附近的灰度分布来进行亚像素定位；选用Canny算子进行整像素级边缘提取，首先 对图像进行高斯平滑，通过高斯卷积实现；然后对平滑后的图像进行简单的2维一阶微分 操作，得到梯度大小和方向，采用W下2X2大小的模板作为对X方向和y方向偏微分的一 阶近似
【权利要求】
1. 基于亚像素边缘算法的平面零件尺寸测量方法，利用CCD设备作为测量平台，其特 征在于，包括下列步骤： (1) 结合光学成像理论分析了背光光源下不同厚度钣金零件的边缘特性； (2) 根据钣金零件图像的边缘分布特征，采用Canny边缘粗定位； (3) 基于多项式拟合的亚像素检测； (4) 通过平均距离法区分被测物上下边缘； (5) 采用三角法求取钣金零件准确边缘； (6) 利用成像原理和直线与零件表面交点确定被测物尺寸。
2. 根据权利要求1所述的基于亚像素边缘算法的平面零件尺寸测量方法，其特征在 于，步骤⑴中 实际的CCD成像系统中，CCD感光元不但接收照射到自身感光面的光，还接收照射到相 邻感光元的光，尤其是对边缘点，物体和背景的不同反射特性以及CCD器件的积分效应，造 成CCD器件对阶跃边缘的响应产生由明到暗或由暗到明的渐变过程，所以边缘在图像中表 征为一种灰度分布，如图2所示；另外，当被测物有一定的厚度且尺寸较大时，在距离相机 一侧存在两个梯度值较大的过渡带，而在相机另一侧时，只有上边缘；因此边缘信息需要区 分上边缘与下边缘，并只提取对相机都看见的上边缘；然后通过Canny算子进行整像素级 的边缘提取，得到整像素级的边缘后用三次多项式拟合的方法提取亚像素边缘精确定位， 再通过平均距离法确定钣金零件的上边缘，通过三角法最终确定零件的实际边缘，最后利 用成像原理和直线与零件表面交点确定被测物尺寸。
3. 根据权利要求1所述的基于亚像素边缘算法的平面零件尺寸测量方法，其特征在 于，步骤⑵中， 亚像素边缘定位前必须先采用像素级边缘检测方法确定边缘点的位置，然后根据边缘 点附近的灰度分布来进行亚像素定位；选用Canny算子进行整像素级边缘提取，首先对图 像进行高斯平滑，通过高斯卷积实现；然后对平滑后的图像进行简单的2维一阶微分操作， 得到梯度大小和方向，采用以下2X2大小的模板作为对X方向和y方向偏微分的一阶近 似；
0. arctan [Gy (x, y) /Gx (x, y)] 式中：Gx，Gy分别为图像像素点x方向和y方向偏微分的一阶近似，M(x，y)为此点的梯 度大小，9为该点梯度方向； 然后采用"非最大抑制"算法寻找图像中的可能边缘点，最后通过双门限值递归寻找图 像边缘点得到单像素宽度边缘图像。
4. 根据权利要求1所述的基于亚像素边缘算法的平面零件尺寸测量方法，其特征在 于，步骤（3)中， 提取亚像素边缘使用三次多项式拟合法，根据灰度分布函数的特点，利用多项式函数 对灰度分布作最小二乘拟合，并根据拟合出的函数曲线来确定灰度分布边缘点即亚像素边 缘点的位置； 采用的三次多项式的形式为f (X，y) = ax3+bx2+cx+d并且假设:P = /(X)将全部测量值 与回归直线的偏离平方和记为S，则
式中=Xi, Yi为边缘像素坐标，a, b, c, d为多项式拟合的系数。
5. 根据权利要求1所述的基于亚像素边缘算法的平面零件尺寸测量方法，其特征在 于，步骤⑷中， 使用平均距离法来确定被测物的上下边缘，首先对钣金零件图像进行运用亚像素精度 阈值分割的算法得到轮廓C，设轮廓C的亚像素坐标为集合Q (?, q2，….qi)，待检测图像边 缘点坐标为集合Pb1, P2,….Pt)，计算两边缘间的偏差量H(P, Q),计算方法为： 首先，对P中的每一个点，计算其到Q上所有点的欧式距离，距离用符号I |p_q| I表示， 然后，将得到的距离进行排序，取其中距离最小值为此点的偏差量，用h(p，Q)表示，最后， 计算h(p，Q)的集合即为两边缘间的偏差量H(P，Q);
H(P, Q) = {h(p, Q) |p G P} 式中：(xl,yl)为待检测图像边缘点坐标，（x2,y2)为亚像素轮廓坐标，q」为亚像素轮 廓的一点，B为亚像素轮廓像素个数； 计算两边缘间的偏差量H(P，Q)，集合的中间值mid，计算H(P，Q)中小于mid的所有值 的均值D，将轮廓Q向内收缩D长度，得到即为图像的边缘轮廓；
式中：f是集合H(P，Q)中小于mid的值的个数，Di为集合H(P，Q)中小于mid的所有 值。
6. 根据权利要求1所述的基于亚像素边缘算法的平面零件尺寸测量方法，其特征在 于，步骤（5)中， 根据摄像机成像原理，对于有厚度的被测物，底部轮廓会被遮挡，导致成像成比例的增 大，所以上述算法求得的边缘需要进一步求解；创建边缘点的三维模型如图3所示，图中 虚线表示从摄像机焦点到边缘点的视线，由上述的边缘提取算法得到点P的图像坐标（px， Py)，若已知物体厚度T和摄像机标定可以求得摄像机焦点的z坐标fpz，则实际边缘点p ' (px '，py ')坐标可以通过三角形相似法可以，如下式所示；
式中：(px, Py)为边缘点图像坐标，p ' (px ' , py ')为边缘点实际坐标。
7.根据权利要求1所述的基于亚像素边缘算法的平面零件尺寸测量方法，其特征在 于，步骤（6)中， 获取边缘信息后，结合标定参数，即可解算测量数据；通过图像点坐标求取对应的物方 坐标是利用成像的基本原理，根据图像上目标成像点所对应的光线直线与边缘点所在平面 的交点来确定边缘点的空间坐标，进而得到目标的位置信息；定义该光线L。的两个点分别 为摄像机的投影中心和图像坐标系中的（r，c)T点，通过摄像机标定模型，得到光线L。在摄 像机坐标系中的方程： Lc = (0, 0, 0)T+A (u, V, f)T 根据摄像机标定得到的外参得到光线L。在世界坐标系下的方程： Lc= 〇w+ X (Iw-Ow) = Ow+ 入 Dw 式中：0W = -RtT是转换后的光心点，Iw = Rt((u, V, f)t-T)是转换到成像平面上的点， Dw= (dx，dy，dz)T是光线的方向向量，IT 1是摄像机标定中旋转矩阵R的逆矩阵；入为比例系 数，u, V为图像坐标，f为相机焦距； 计算光线与测量平面Z = 0的交点坐标，即目标点的世界坐标：
式中：(ox, oy, oz)为光心点世界坐标，（dx，dy，dz)为光心的方向向量坐标。
【文档编号】G01B11/00GK104359403SQ201410675297
【公开日】2015年2月18日 申请日期:2014年11月21日 优先权日:2014年11月21日
【发明者】耿磊, 李文杰, 肖志涛, 张芳, 吴骏, 李月龙, 袁菲, 杨振杰, 苏静静 申请人:天津工业大学