基于灰度的亚像素插值视觉测量方法

戴宗贤 莫洪波 周勇 张震 尹爱军 梁子晓

摘要：针对现有视觉测量中的检测代价高，精度低和速度慢问题，该文提出一种基于计算机机器视觉的紧密内插值亚像素测量方法。该方法基于线性插值算法的原理，结合常规边缘检测方法和图像的灰度曲线图，利用阈值分割和标准长度进行亚像素自适应阈值选择。为验证该方法的有效性，对标准量块长度进行测量实验，并分析测量系统的误差影响因素。相比于传统的Canny算子方法检测结果，该方法的平均测量准确度提升46.2%。实验结果表明该算法的测量精度较高，可以快速、精确地测量出物体的几何尺寸。

关键词：视觉测量;尺寸测量;边缘检测;亚像素插值算法

中图分类号：TP391.4 文献标志码：A 文章编号：1674-5124（2019）05-0033-05

收稿日期：2018-03-30;收到修改稿日期：2018-04-24

基金项目：国防基础科研计划（JCKY2016209B008）

作者簡介：戴宗贤（1986-），男，安徽六安市人，工程师，硕士，研究方向为测量与测试技术，图像处理。

0 引言

几何尺寸测量在现代工业中有着重要作用，是保证产品质量的关键。随着科技与工业生产的发展，对零件几何尺寸的测量精度及效率提出了更高要求。传统的接触式测量主要利用探头接触待测物体表面，从而获得所需参数信息，如卡尺、千分尺等。但该方法易受设备和人为因素影响，只适合小批量和一般精度的测量。在高精度测量中，如三坐标测量机、光学影像测量机等检测设备复杂昂贵难以广泛应用。近年来兴起的非接触测量技术因其高效、无损、高分辨率的检测优点而逐渐受到重视[1-3]，随着图像处理等技术的发展，视觉测量作为一种非接触测量方法已得到广泛研究与应用。

国外对几何尺寸视觉检测方面的研究起步较早，目前已经达到较高水平，Jan Fischer等[4]研究了线性插值的亚像元位置测量系统，针对传感器不同质像素敏感性分布的缺点，提出一种线性插值方法。Eric等[5]研究了高精度机械位移与应变的视觉测量，提出视觉敏感性概念。国内杨剑等[6]进行了大尺寸视觉测量精度的理论和实验研究，提出了一种基于距离约束的大尺寸三维视觉测量算法，并应用在大型设备曲面测量上。申晓彦等[7]研究了轴对称构件几何尺寸的图像测量技术，构建一组适合工程应用的图像测量系统。

本文在现有的几何尺寸测量算法基础上，首先介绍边缘检测的常规方法，并探讨了基于边缘点插值测量尺寸的数学原理;在此基础上，融合常规边缘检测方法，提出一种针对标准薄片零件几何尺寸的亚像素插值视觉测量方法。

1 边缘提取与亚像素测量方法

1.1 常规边缘提取方法

视觉测量大多是通过处理图像目标区域边缘而获得物体的几何参数。边缘提取方法主要分为灰度值提取和灰度梯度值提取[8-9]。基于灰度值的边缘提取一般是根据阈值对图像进行二值化处理[10-11]，这种方法简单，计算量小，但是提取精度不高，边缘定位的准确性依赖于阈值的计算方法。

基于灰度值梯度的边缘检测根据图像边缘的变化情况，检测出图像中的边缘点[12]再将边缘点连接成轮廓，从而实现边缘提取。灰度的变化率和变化方向通常作为图像边缘像素点的评价指标，分别以梯度向量的幅值和方向来表示。平面区域D内，坐标为（x，y）的像素点f（x，y）在θ方向沿x的梯度定义为

对于连续图像f（x，y）来说，其方向导数将在边缘（法线）方向上有局部最大值。因此，边缘检测就是求f（x，y）梯度的局部最大值。根据式（1），当■达到局部最大值时，有，于是有

f_xsinθ_g+f_ycosθ_g=0（2）其中，梯度最大，一般称之为梯度模;灰度的变化方向为，实际计算中一般以微分算子的形式表示。

常用的微分算子有Roberts算子、Sobel算子和Canny算子等[13]，其中Roberts算子和Sobel算子是根据图像边缘处梯度最大（正的或负的）的特性来进行边缘检测。这两种方法，Roberts边缘定位准，但对噪声敏感;而Soble加权平均，适用边宽≥2像素的检测。Canny算子是最优的阶梯型边缘（step edge）检测算子，对受到白噪声影响的图像检测准确。

1.2 亚像素测量方法

传统边缘检测算法的定位精度只能达到像素的整数级，为了进一步提高检测精度，基于亚像素的边缘检测与提取方法得到了深入研究。常用的亚像素方法有拟合法[14]、插值法[15]、灰度矩法[16]和一些组合算法[17-18]等。德国MVTec的Halcon，加拿大Adept的Hexsight，美国Cognex的VisionPro等图像处理软件，其测量精度都已达到亚像素级。赵敏等[19]应用概率统计的方法提高了图像测量分辨率，娄训志等[20]研究了基于小波的亚像素检测方法。一般常用的亚像素算法的精度为0.1～0.5个像素，一些算法在理想的情况下可以达到0.01个像素。

2 一种新的亚像素高精度测量方法

由于视觉测量中一般需要复杂的光学系统，且要求在特殊场合、特殊环境下使用，精度较低、效率和实时性较差。针对该问题，本文提出灰度插值的亚像素精确测量方法，基本过程包括：利用标准长度进行亚像素自适应阈值选择、基于阈值分割的亚像素测量。

2.1 亚像素阈值的标定

系统首先利用标准值进行标定，即确定测量阈值。本文采用自适应阀值的线性插值算法，如图1所示，f为焦距;D为物距;h为CCD靶面尺寸（高）;ν为CCD靶面尺寸（宽）。

由图可知，镜头离被测物的高度为D时，由线阵CCD传感器的靶面尺寸可知：系统所能测量的最大宽度V，最大高度H范围（量程）为

V=νD/f H=hD/f（3）

则单位像长对应的实际长度μ为

μ=V/ν=H/h=D/f（4）

而每个像素对应的实际宽度β_ν、实际高度β_h为

R_ν=k_ν×μ=k_νD/f β_h=k_h×μ=k_hD/f（5）其中k_ν×k_h为像元尺寸。

其阈值的确定过程如图2所示，其中X表示边缘像素点的位置，Y表示对应像素点的灰度值，设标准件长为l，则可知其理论像素个数N=[l/β_ν]，其中[]表示取整数。

首先利用Canny算子对图像进行边缘提取（粗定位），得到边缘点的像素坐标（x，j），设第j处两邊缘点为C₁（x_t1，y_t1），C₂（x_t2，y_t2），边缘间的像素数为N_j。

若N>N，则C₁点向左移动，C₂点向右移动;若Nj

由式（5）、（11）、（12）、（14）可得：

目标物体的尺寸为

3 实验结果及误差分析

亚像素精确尺寸测量方法设计和实验如图4所示，镜头离被测物的距离为2.93cm，系统亚像素定位精度理论为0.0019像素，传感器像元个数M=1024，镜头角度理论精度约为0.012μm。相机采用Baumer公司的TXG03c，工业镜头H0514-MP，镜头的焦距为5mm，CCD分辨率为656×490像素，为了减少外界的干扰，提高测量准确性，在测量时采用平行光光源，对三等标准量块进行灰度阈值的标定。现场实验如图4所示。

3.1 标准量块测量实验

表1是在上述实验条件下所测的不同规格的量块尺寸，并与用Canny算子进行边缘提取所得数据进行了对比。从表中可以看出，本方法测量误差均比Canny方法测量误差值小。8次测量平均绝对误差值，本方法为0.000 438mm，Canny方法为0.000813mm，本方法误差为Canny方法誤差的53.8%。从表中可以看出，本系统的测量结果更接近标准量块实际值，实验表明本算法测量结果可靠且能够实现快速无接触的测量。

3.2 误差分析

从上述实验数据可以看出，本检测系统存在一定的测量误差，主要由以下4方面引起：1）镜头到被测物的高度即物距由于存在测量误差，引起β_ν的变化，从而导致检测系统存在系统误差。2）由于本实验所选用镜头的焦距短，镜头畸变相对较大，造成图像的非线性失真，从而影响目标的边缘定位，导致测量结果的偏差。3）系统光源的影响，线阵光源的准直性和稳定性将直接影响到测量结果。理想模型中，线阵光线应该是平行投射的，但在实际中，由于光源发散性的影响，将不能得到稳定的平行光源，从而产生了误差。4）被测物体放置的影响，只有当被测物体和传感器完全垂直的时候，成像才能反映出物体的真实轮廓。上述结果表明，经过本文提出的插值方法后，可以得到准确的结果，对于批量检验薄片零件是否合格有应用价值。

4 结束语

基于视觉系统的几何尺寸测量技术已得到了广泛的研究与应用。本文介绍了图像测量的基本原理，包括测量光学成像系统和测量结构;深入研究了基于光强线性插值的亚像素自适应阈值的测量原理和方法，分析了测量系统的误差影响因素。实验结果表明该算法的测量精度较高，可以快速、精确地测量出物体的几何尺寸。

参考文献

[1]王鹏.一种局部二值化方法及其应用[D].长春：吉林大学，2017.

[2]徐晶晶.低质量文档图像的二值化算法研究[D].武汉：湖北工业大学，2017.

[3]张梅，王涛.一种工业恶劣环境下的图像二值化方法及FPGA实现[J].微型电脑应用，2017，33（7）：69-72.

[4]FISCHER J，PRIBULA O.Precise subpixel positionmeasurement with linear interpolation of CMOS sensor imagedata[C]// The 6IEEE International Conference on IntelligentData Acquisition and Advanced Computing system：Data，Technology and Applications，2011.

[5]ERIC P.Visual hyper acuity representation and computationof high precisio position information[J].Computer Vision，Graphics and Image Processing，1986，33：99-115.

[6]杨剑，大尺寸视觉测量精度的理论和实验研究[D].北京：北京邮电大学，2010

[7]申晓彦.轴对称构件几何尺寸的图形测量[D].太原：中北大学，2010

[8]董鸿燕.边缘检测的若干技术研究[D].长沙：国防科学技术大学，2008

[9]季虎，孙即祥，邵晓芳，等.图像边缘提取方法及展望[J].计算机工程与应用，2004，14：70-73.

[10]李了了.工业现场字符识别方法的研究[D].合肥：合肥工业大学，2003.

[11]周衍超，罗浩彬.机器视觉技术在工业机器人的应用研究[J].装备制造技术，2015（6）：193-195.

[12]刘力双，孙双花，吕奶光.视觉系统中调焦函数的选择方法研究[J].北京机械工业学院学报， 2008，23（1）：1-4.

[13]张广军.机器视觉[M].北京：科学出版社，2004：33-42

[14]HUECKEL M H.An operator which locate edges in digitizedpictures[J].Journal of the Association for ComputingMachinery，1971，18（1）：113-125.

[15]李艳，彭嘉雄.基于D2样条插值和LOG算子的亚像素边缘检测[J].华中理工大学学报，2000，28（3）：77-79.

[16]邹宁，柳剑，周曼丽.基于Zernike矩的边缘检测方法[J].华中理工大学学报，1999，27（5）：1-2.

[17]梁燕，刘文耀.一种新型非线性边缘检测方法[J].光电子.激光，2003，14（12）：1332-1335.

[18]ZHANG H T，DUAN F J，YE S H.Study on a fast method forsub-pixel edge detection[J].Acta Metrologica Sinica，2002，23（4）：263-265，270.

[19]赵敏，邱宗明，黄秋红，等.提高图像测量分辨率的一种方法[J].西安理工大学学报，2002，18（1）：80-83.

[20]娄训志.基于小波变换的图像亚像素检测应用研究[D].武汉：湖北工业大学，2011：22-38

（编辑：莫婕）