本发明属于先进测量
技术领域：
，涉及了一种基于双目视觉的轴类零件半径尺寸测量方法，尤其适用于非接触式工业检测和基于视觉的机器人导航系统。
背景技术：
：轴类零件在工业应用中随处可见，对轴类零件的半径尺寸精确测量在工业检测中具有重要的意义。尤其在现代制造技术中，自动化装配技术已经逐步取代人工技能和判断力进行复杂的装配操作来提高效率，保证产品的质量及稳定性。在零件装配过程中，轴孔装配是经常会遇到的工作情况，在机械手进行抓取轴类零件以及进行装配时，必须准确测量出轴类零件的半径尺寸，才能合理控制机器手的张开角度以及避免零件间发生碰撞而产生缺陷，进而顺利完成装配任务。随着机器视觉技术的发展，目前国内外已经开展了相关的研究及应用试验，将机器视觉系统引入到轴孔尺寸测量中来，但是目前大多数研究集中在空间圆孔及其椭圆孔的识别与参数测量，对于轴类零件尺寸测量的研究相对较少。崔彦平等人研究了一种基于双目视觉的回转体目标空间三维姿态的测量方法，该方法能够在不进行特征点匹配的情况下实现回转体三维姿态的测量，但由于轴类零件表面的高反光性以及实际应用中光照不均的问题，使用该种方法在提取像面母线的亚像素直线方程时将导致精度下降，并且该种方法只能测量回转体的位姿而无法测量半径尺寸。Sun等人在“shaftdiametermeasurementusingadigitalimage”一文中提出一种基于图像处理的轴类测量方法，这种基于单目视觉原理而脱离立体几何原理的测量方法必然过分依赖图像质量及其处理效果，缺乏良好的稳定性，并且该方法要想达到较高的精度需要利用已知直径的轴类尺寸作为先验知识进行重新标定，应用场合具有一定限制。此外，西安理工大学张凯提出的目标工件轴线计算方法由于在图像处理过程中需要采用灰度值阈值和RGB彩色阈值，使得在采用普通黑白摄像机的情况下可能无法识别，进行计算失败。综上所述，提出一种无需人工介入而且能够精确测量出轴类零件半径尺寸的方法具有较大的应用价值。技术实现要素：本发明的目的在于针对现有技术的不足，提出一种基于双目视觉的轴类零件半径尺寸测量方法，能够准确计算出空间轴类工件的中轴线位姿和轴半径，精度高、速度快，设计的检验模型与实验方法简单实用，节约成本，提高效率，为双目视觉算法验证工作提供了一种新思路。本发明的上述目的是通过下述技术方案的步骤实现的：1)根据张正友标定法分别对双目视觉系统的左右摄像机进行标定，由左右摄像机构成双目视觉系统，获得左右摄像机的内参、外参和畸变参数；根据标定得到的内参和外参计算出左右摄像机的投影矩阵；2)将轴类零件放在左右摄像机的公共视野内，通过左右摄像机对轴类零件进行拍照分别获得左右图像的两幅图像(左图像和右图像)，对两幅图像依次进行高斯滤波、灰度转换和阈值处理，然后利用步骤1)获得的左右摄像机的径向畸变参数矩阵对左右图像进行矫正处理；3)对矫正后的左右图像进行处理获得各自的轴轮廓；4)利用左右图像的轴轮廓处理获得轴类零件的半径。所述的步骤1)获得的左右摄像机的内部参数、畸变参数和投影矩阵包括左摄像机的内参数矩阵Al、左摄像机的径向畸变参数矩阵Kl、左摄像机的投影矩阵Ml、右摄像机的内参数矩阵Ar、右摄像机的径向畸变参数矩阵Kr和右摄像机的投影矩阵Mr。所述步骤3)具体为：对矫正后的左右图像和模板图像进行轮廓检测，再将检测到的左右图像中的所有轮廓分别与模板图像检测获得的模板轮廓进行匹配，获得左右图像中各自与模板轮廓匹配的轮廓，并作为轴轮廓，即与轴类零件对应的轮廓；所述的模板图像中具有且仅有一个与轴类零件匹配的矩形框。矩形框与轴类零件在尺寸和形状上相匹配，具体实施矩形框置于模板图像的中间。所述步骤4)具体为：4.1)用最小旋转外接矩形逼近左右图像中的轴轮廓，提取获得每个轴轮廓的长边，左右图像中两个轴轮廓的长边共计有四条长边ll1、ll2、lr1、lr2，ll1、ll2分别表示左图像中轴轮廓的两条长边，lr1、lr2分别表示右图像中轴轮廓的两条长边；4.2)采用以下公式计算通过左右摄像机中心并与轴类零件的轴轮廓相切的四个空间平面SOAB、SOCD、S0′GH、S0′IJ，进而获得四个空间平面的法向量NOAB、NOCD、N0′GH、N0′IJ：ll1MlSOAB=0,ll2MlSOcD=0,lr1MrS0′GH=0,lr2MrS0′IJ=0式中，ll1、ll2分别表示左图像中轴轮廓的两条长边，lr1、lr2分别表示右图像中轴轮廓的两条长边；Ml、Mr分别表示左右摄像机的投影矩阵；SOAB、SOCD分别表示通过左摄像机中心并分别与轴类零件的轴轮廓两条长边相切的两个空间平面，S0′GH、S0′IJ分别表示通过右摄像机中心并分别与轴类零件的轴轮廓两条长边相切的两个空间平面；4.3)根据空间平面SOAB、SOCD及其各自的法向量NOAB、NOCD计算空间平面SOAB和空间平面SOCD之间的空间角平分平面Sl，根据空间平面S0′GH、S0′IJ及其各自的法向量N0′GH、N0′IJ计算空间平面S0′GH和空间平面S0′IJ之间的空间角平分平面Sr，取两个空间角平分平面Sl与Sr的交线作为轴类零件的中轴线；4.4)利用下式计算中轴线上任意一点p(x0，y0，z0)分别到四个空间平面SOAB、SOCD、S0′GH、S0′IJ的垂直距离li，取四个垂直距离的平均值作为所求轴类零件的半径：式中，i＝1,2,3,4时的li分别表示p点到四个空间平面SOAB、SOCD、S0＇GH、S0′IJ的距离，ai,bi,ci,di分别表示四个空间平面SOAB、SOCD、S0′GH、S0′IJ对应方程的第一、第二、第三、第四系数，x0、y0、z0为中轴线上任意一点的三维坐标，通常的选取x0＝0。本发明实施中还提出了一种新的检验模型，该模型半径已知，并对中轴线进行可视化，将该条空间直线拟合、重建结果与本发明方法所得结果进行对比，验证本方法半径尺寸测量结果的有效性与精确度。与现有方法相比，本发明的有益效果是：本发明能高精度地测量获得空间轴类工件的半径，可以有效求取较大范围的空间轴半径，测量精度高，结果可靠，方法简单实用，节约成本。并通过自行设计的检验模型对该方法进行了验证，可同时对空间轴的中轴线和半径尺寸重建结果进行检验，避免了高价的检测台，降低成本。同时，本方法具有非接触式的优点，在用传统方法无法测量的场合具有很高的应用价值，尤其适用于非接触式工业检测和基于视觉的机器人导航系统。附图说明图1为本发明方法原理示意图；图2为本发明方法流程示意图；图3为实施例的新模型验证方式示意图。具体实施方式为了更好的理解本发明，下面结合附图和实施例对本发明的技术方案作详细的描述。如图1所示为双目立体视觉系统。OlXlYlZl和OrXrYrZr分别为左右摄像机坐标系，olulvl和orurvr分别为以像素为单位的左右图像坐标系，OwXwYwZw为世界坐标系，其中Z轴指向左摄像机坐标系原点Ol。图1中EF表示轴类零件的中轴线，abcd表示轴类零件在左摄像机成像平面上的轮廓，ghij表示轴类零件在右摄像机成像平面上的轮廓。下面详细叙述本发明方法的实施步骤：1.对左右摄像机，利用张正友标定法(AFlexibleNewTechniqueforCameraCalibration.ZhengyouZhang,December,2,1998.)，确定出左摄像机的内参数矩阵Al、左摄像机的外参数矩阵Rl和Tl、左摄像机的径向畸变参数矩阵Kl、右摄像机的内参数矩阵Ar、右摄像机的外参数矩阵Rr和Tr、右摄像机的径向畸变参数矩阵Kr；并计算左摄像机对应的投影矩阵Ml、右摄像机对应的投影矩阵Mr。其中，内参数矩阵的形式为：其中，αl和βl分别表示左摄像机x轴方向的有效焦距和y轴方向的焦距，uol和vol表示左摄像机成像平面的主点坐标，γl表示左摄像机坐标轴倾斜参数，理想情况下为0；αr和βr分别表示右摄像机x轴方向的有效焦距和y轴方向的焦距，uo和vor表示右摄像机成像平面的主点坐标，γr表示右摄像机坐标轴倾斜参数，理想情况下为0。本实例中摄像机标定与计算结果为：Kl＝[-0·0114-0·0676-1·4004]Kr＝[-0·0253-0.1129-0·3821]2.将双目视觉系统放置在待检测的轴类零件附近，确保该轴类零件在左右摄像机的公共视场范围内，同时使背景尽量简单、使轴类零件的轴线方向尽可能平行于左右摄像机坐标系原点的连线，即双目视觉系统的基线。利用左右摄像机同时拍摄轴类零件，从而使在左摄像机得到一幅包含轴类零件影像的图像，相应地，在右摄像机也得到一幅包含轴类零件影像的图像。利用左摄像机的径向畸变参数矩阵Kl对左图像进行畸变矫正，得到不含有畸变信息的左图像，记为planel。同时，利用右摄像机的径向畸变参数矩阵Kr对右图像进行畸变矫正，得到不含有畸变信息的右图像，记为planer。具体的矫正过程是：对左边的图像，设某包含畸变信息的图像点在以像素为单位的图像坐标系下的坐标为其归一化的图像坐标为它们对应的不含畸变信息的图像点坐标分别为(u,v)和(x,y)。根据文献(D.C.Brown，Close-rangecameracalibration,Photogram-metricEngineering,37(8):855-866,1971)，有：利用坐标变化公式：其中，Kl为左摄像机的径向畸变参数矩阵，Al为摄像机的内参数矩阵，都可以通过摄像机的单目标定确定。可以得到：由于以上方程是非线性方程组，为了简化求解过程，上述方程组可近似为：利用上面两式可以对左右图像上的每一个图像点进行畸变矫正，从而得到不含有畸变信息的图像planel。对于右边的图像，矫正方法与左图像的矫正方法完全相同，不再赘述。3.对矫正后的左图像planel和右图像planer进行自适应阈值将图像二值化，具体的自适应阈值采用以下方法计算：设变量t在灰度值范围(0～255)内依此取整数值(共256个灰度值)，每次取值将左图像分为背景和前景两部分，同时计算下列两式：u＝w0*u0+w1*u1g＝w0*(u-u0)2+w1*(u-u1)2其中，w0为图像背景像素点占整幅图像的比例，u0为图像背景像素点的平均灰度，w1为前景像素点占整幅图像的比例，u1为前景像素点的平均灰度，u表示整幅图像的平均灰度，计算结果g表示前景和背景图像灰度值的方差。比较所得的256个g值，g值最大时变量t的值为最佳阈值，根据最佳阈值进行图像二值化分割。然后进行轮廓检测，并将检测到的所有轮廓与模板轮廓进行匹配，计算匹配相似度得分，根据得分高低获得左右图像中目标轴的轮廓；由于轴类零件的轴线方向基本平行于左右摄像机坐标系原点的连线，可知轴在左右成像平面上的影像近似于矩形，故本文采用最小旋转外接矩形逼近左右图像中匹配到的轴轮廓，进行计算得到左图像planel和右图像planer上轴轮廓的4条轴线方向的边ll1、ll2、lr1、lr2，采用最小旋转外接矩形的方法还能弥补由于轴类零件表面高光以及光照不均匀造成的图像轮廓分割误差。本实例计算的ll1、ll2、lr1、lr2值为：ll1：[463.23004,-1358.9934,847808]ll2：[463.23004,-1358.9935,292788.81]lr1：[459.90479,-1310.1235,987763.63]lr2：[459.90479,-1310.1235,432691.72]4.根据左图像planel和右图像planer上4条轴线方向的轮廓边ll1、ll2、lr1、lr2以及左摄像机投影矩阵Ml和右摄像机投影矩阵Mr计算通过左右摄像机中心并与空间轴类零件侧面相切的4个平面SOAB、SOCD、S0′GH、S0′IJ。已知图1中空间轴类零件的母线AB和CD在左摄像机成像平面上的像为ab和cd，母线GH和IJ在右摄像机成像平面上的像为gh和ij，则平面SOAB、SOCD、S0′GH、S0′IJ在世界坐标系下的方程为：ll1/MlSOAB＝0，ll2MlSOCD＝0，lr1/MrS0′GH=0,lr2MrS0′IJ=0同时获得4个空间平面的法向量NOAB、NOCD、N0′GH、N0′IJ。本实例中计算得到：SOAB：[1744445.6,-5117738.5,235642.83,0]SOCD：[1744445.6,-5117738.5,-319376.47,0]S0′GH：[1731923.3,-4933702.5,415048.59,-1.0369749e+008]S0′IJ：[1731923.3,-4933702.5,-140023.33,-1.0369749e+008]NOAB：[0.32232851,-0.94562596,0.043540712]NOCD：[0.32207313,-0.94487667,-0.058965769]S0′GH：[0.33018529,-0.94059366,0.07912761]S0′IJ：[0.33110517,-0.94321412,-0.026769344]5.通过基本几何原理，由空间平面SOAB、SOCD及其相应的法向量NOAB、NOCD可计算得到SOAB、SOCD的空间角平分平面Sl；根据空平面S0′GH、S0′IJ及其相应的法向量N0′GH、N0′IJ计算S0′GH、S0′IJ的空间角平分平面Sr；则两个平分平面Sl与Sr的交线为空间轴的中轴线。本实例中计算的得到：Sl：[0.64440167,-1.8905027,-0.015425056,0]Sr：[0.66129047,-1.8838078,0.052358266,-39.594223]6.利用下式分别计算中轴线上任意一点p(x0，y0，z0)到4个切平面SOAB、SOCD、S0′GH、S0′IJ的垂直距离，取四个垂直距离的平均值为所求轴类零件的半径：式中，i＝1,2,3,4时，li分别表示p点到四个空间平面SOAB、SOCD、S0′GH、S0′IJ的距离，ai,bi,ci,di分别表示四个空间平面SOAB、SOCD、S0′GH、S0′IJ对应方程的系数，x0、y0、z0为中轴线上任意一点的三维坐标。在本实例中，选取x0＝0，计算得到p点距离4个平面的距离分别为：29.9645mm、29.8936mm、30.0243mm、30.4308mm，平均值为30.0783mm，误差为0.26％。本实施例最后设计一种检验模型，如图3所示，该模型设计为I和II两部分。第I部分为一段半径为30mm，长度为180mm的光轴，其作用是通过本发明所提出的方法对该段光轴进行半径测量，并与实际半径比较，检验本发明所提出测量方法的精确度。第II部分为一段光轴的四分之一，与第I部分光轴同轴且半径相同，长度为100mm，模型的这一部分对中轴线进行了可视化，通过对可视化中轴线的拟合结果与本发明所提出方法对第I部分光轴的中轴线测量结果进行对比，可以验证本发明方法过程中求取中轴线的精确度。本实施例中对检验模型可见直线的拟合结果以及本发明测量结果如下表1(空间直线表示为两个空间平面的交线)：表1平面1平面2轴线方向矢量拟合结果[0.6511,-1.8910,-0.0155,0][0.6673,-1.8333,0.0576,-39.9524][-0.1303,-0.0478,0.0345]算法重建[0.6544,-1.8905,-0.0154,0][0.6613,-1.8838,0.0524,-39.5942][-0.1280,-0.0439,0.0362]本实施例对半径的重建结果与实际半径对比结果为：表2重建值实际值误差误差比率轴类零件30.0783mm30.0000mm0.0783mm0.26％由此可以看出，本文所提出的基于双目视觉的轴类零件半径测量方法可以到达较高的精度，为实现轴类零件的半径尺寸测量实现自动化提供保证。当前第1页1 2 3