CN109974608B - 基于机器视觉技术的零件孔径尺寸测量方法 - Google Patents
基于机器视觉技术的零件孔径尺寸测量方法 Download PDFInfo
- Publication number
- CN109974608B CN109974608B CN201910201755.0A CN201910201755A CN109974608B CN 109974608 B CN109974608 B CN 109974608B CN 201910201755 A CN201910201755 A CN 201910201755A CN 109974608 B CN109974608 B CN 109974608B
- Authority
- CN
- China
- Prior art keywords
- coordinate system
- world coordinate
- small circular
- tubes
- hole
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/08—Measuring arrangements characterised by the use of optical techniques for measuring diameters
- G01B11/12—Measuring arrangements characterised by the use of optical techniques for measuring diameters internal diameters
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/70—Determining position or orientation of objects or cameras
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Theoretical Computer Science (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明涉及一种基于机器视觉技术的零件孔径尺寸测量方法,属于机器视觉测量技术领域。能够实现以上零件孔径尺寸的非接触测量。在标定CCD摄像机的内部参数和畸变系数后,将三根尺寸已知的相同小圆管插入零件孔内并拍摄图像,拟合小圆管端面的内、外边缘曲线;然后,用相邻小圆管重建世界坐标系,求出小圆管端面中心点的世界坐标,并据此解出世界坐标系和摄像机坐标系的变换关系;最后,由于零件孔边缘倒角会影响检测零件孔的边缘图像点,所以不采用拟合零件孔边缘曲线的方法,而是利用零件孔和小圆管的几何关系求出零件孔径尺寸。本发明利用机器视觉测量技术,采用在零件孔内放置三根尺寸相同小圆管的方法,实现零件孔径尺寸的非接触测量。
Description
技术领域
本发明涉及机器视觉测量技术领域,特别涉及一种基于机器视觉技术的零件孔径尺寸测量方法。
背景技术
在重工、农机、汽车等机械行业中,有很多零件都带有孔,例如汽车的变速箱、机床的主轴箱、轴承座等等一般都带有不同尺寸的孔。带孔的零件在加工制造、表面热处理等过程会不可避免地出现不同程度的变形,从而使孔径尺寸出现较大的误差,这些变形或过大的误差若不能被有效精准地测量,在后续的精加工或装配过程中将会带来各种各样的问题,会直接影响产品的质量,甚至直接影响产品的废品率。所以,为了降低零件的废品率、提高成品率、节约加工和原材料的成本以及保证工件能够达到公差要求,快速、准确、有效地测量孔的直径尺寸是带孔零件加工过程中保证产品质量的重要工序,也是产品加工后进行质量检测的必要手段。
工业上现有主要的孔径测量方法按照测量装置是否与被测零件相接触进行划分,一般可分为接触测量方法和非接触测量方法。接触测量方法主要包括:卡尺测量法、螺旋测微器测量法、内径千分尺测量法和三坐标测量机测量法等方法,接触测量方法不仅工作强度大,而且效率低,不利于大范围推广。非接触测量方法一般是指基于光学原理的测量方法,主要包括机器视觉测量法、激光准直仪测量法、结构光测量法等方法。非接触测量方法与传统的接触测量方法相比,具有很多优点。非接触测量方法可以实时地将零件的信息反馈到控制系统,易于信息的集成和管理,能够实现智能测量。此外,非接触测量方法可用于长时间恶劣条件下的测量,也可用于无法进行接触情况下的测量。
机器视觉测量技术主要是利用物体的图像作为测量的载体,被测尺寸所在的物体表面就是成像平面,可以利用机器视觉中的边缘检测技术和摄像机的标定技术相结合完成高精度测量。用机器视觉测量技术测量孔径时,由于孔的特殊性使孔表面的外部参数不易标定,通常采用在孔表面粘贴面积较大标定纸的方法或用工装夹具使标定板和孔表面处于同一平面来进行标定,这些方法都会使标定的精度有所降低。在求解孔径尺寸时,一般采用拍摄孔表面图像,拟合孔边缘曲线的方法求出,但由于孔表面多带有倒角,被测孔边缘倒角会对检测孔边缘图像点产生影响,从而影响测量精度。
发明内容
本发明的目的在于提供一种基于机器视觉技术的零件孔径尺寸测量方法,解决现有零件孔表面不易标定且精度不高,而且零件孔边缘存在的倒角结构使其边缘在成像照片中边界不清,导致零件孔径照片无法直接通过机器视觉中的边缘检测技术实现准确提取,从而造成现有应用机器视觉技术测量零件孔径尺寸的过程繁琐、效率低、精度差的技术问题。本发明通过测量三根小管端面圆中心点坐标的方法标定孔表面的外部参数,解决了孔表面的外部参数不易标定且精度不高的问题,再利用孔内各小管间的几何关系解出孔的直径尺寸,消除了被测孔边缘倒角对检测孔边缘图像点产生的影响,改善了测量精度,实现孔径尺寸的非接触测量。
本发明的上述目的通过以下技术方案实现:
基于机器视觉技术的零件孔径尺寸测量方法,采用在零件孔内放置三根尺寸已知的相同小圆管的方式,实现了零件孔径尺寸的非接触测量,具体过程包括以下步骤:
步骤(1)设定世界坐标系I,建立孔径测量模型;
(1.1)设定理论图像平面和坐标系,先标定CCD摄像机的内部参数和畸变系数,然后保持CCD摄像机的位置不变,将零件孔水平放置并安装好;
(1.2)设定世界坐标系I,让世界坐标系I的原点与摄像机坐标系的原点重合,Z轴垂直于零件孔表面,X-Y平面平行于零件孔表面;
(1.3)在不标定零件孔表面外部参数的情况下,利用设定世界坐标系I的形成的条件,建立世界坐标系I和摄像机坐标系的变换关系方程,然后建立孔径测量模型。
步骤(2)将三根尺寸已知的相同小圆管插入零件孔内使四者端面共面并拍摄端面图像,在理论图像平面拟合三根小圆管端面的内、外边缘曲线,由于没有标定零件孔表面的外部参数,所以先根据相邻小圆管重建世界坐标系,然后在新世界坐标系下根据相邻小圆管的位置关系求出三根小圆管端面中心点O1、O2、O3的世界坐标,具体步骤如下:
(2.1)准备三根尺寸已知的相同小圆管A、B、C,且其尺寸能够保证三根小圆管轻松放置在零件孔内;
(2.2)将三根小圆管插入待检测的零件孔内,且三根小圆管和零件孔四者的端面共面,其中A为底部中间的小圆管,其两侧小圆管B和小圆管C由于重力作用均与小圆管A相外切;
(2.3)用CCD摄像机采集放有三根小圆管的待检测零件孔端面图像一幅,并检测小圆管端面图像的内、外边缘点的亚像素坐标,然后在理论图像平面拟合三根小圆管的端面内、外边缘曲线方程;
(2.4)取两根相邻小圆管A和B,在以其端面中心点O1和O2的连线为轴线的无穷多个平面束中,一定存在一个平面与理论图像平面垂直,将该平面与世界坐标系I的X-Y平面的交线作为世界坐标系II的X轴方向,其余条件与世界坐标系I相同,新建世界坐标系II;
(2.5)利用在世界坐标系II下小圆管A和B相外切的关系以及二次曲线不变量的特点,建立方程组求解两根小圆管A和B端面中心点O1和O2的世界坐标;
(2.6)重复步骤(2.4)和(2.5),以两根相邻的小圆管A和C同理建立世界坐标系III,并求得小圆管C的端面中心点O3在世界坐标系III下的坐标。
步骤(3)利用三根小圆管端面中心点的世界坐标和小圆管端面边缘曲线不变量的特点,求出世界坐标系和摄像机坐标系的变换关系,具体步骤如下:
(3.1)在世界坐标系II下,为求出世界坐标系II和摄像机坐标系的变换矩阵Q,利用小圆管端面边缘曲线的系数和不变量的特点,建立小圆管端面中心点坐标和变换矩阵Q中元素之间的关系;
(3.2)利用步骤(3.1)中建立的关系,求解世界坐标系II和摄像机坐标系的变换矩阵Q;
(3.3)通过步骤(3.2)得到变换矩阵Q的多组解,利用二次曲线不变量的表达式和零件孔与摄像机的相对位置对解集进行筛选,得到变换矩阵Q的唯一解,从而得到世界坐标系II和摄像机坐标系的变换关系;
(3.4)重复步骤(3.1)至步骤(3.3),同理求出世界坐标系III和摄像机坐标系的变换关系。
步骤(4)利用零件孔与三根小圆管的几何关系,求出零件孔径尺寸,具体步骤如下:
(4.1)根据世界坐标系II、III和摄像机坐标系之间的变换关系,将处于不同世界坐标系下三根小圆管端面中心点O1、O2、O3的坐标变换到同一世界坐标系II下;
(4.2)由于零件孔与三根小圆管相内切,三根小圆管之间相外切,所以,两组相邻小圆管端面中心点的连接线段O1O2和O1O3的垂直平分线相交于零件孔表面圆的中心点O0;
(4.3)利用步骤(4.2)的几何关系,建立方程组求出零件孔表面圆中心点O0在世界坐标系II下的坐标;
(4.4)在世界坐标系II下,根据零件孔表面圆中心点的坐标和三根小圆管中任意一根的端面中心点的坐标,利用零件孔与三根小圆管相内切的位置关系,求出零件的孔径尺寸。
本发明的有益效果在于:本发明能够实现φ30以上零件孔径尺寸的非接触测量。基于平面标定法标定CCD摄像机的内部参数和镜头畸变系数,通过标定优化有效地修正了因为畸变而产生的图像失真;根据孔的特征设定世界坐标系,利用设定世界坐标系形成的条件,建立世界坐标系和摄像机坐标系的变换关系方程,从而建立孔径测量模型;在理论图像平面,采集一幅放有三根小管的被测孔端面图像,通过曲线拟合,得到小管端面内、外边缘曲线方程的系数,根据相邻小管建立新的世界坐标系,解出不同世界坐标系下三根小管的端面圆中心点坐标;利用小管端面圆中心点的坐标和二次曲线不变量的关系,解出世界坐标系和摄像机坐标系之间的变换关系,完成孔表面外部参数的标定,有效地解决了现有方法不易标定孔表面外部参数且精度不高的问题;由于实际的孔边缘多带有倒角,在视觉检测中,倒角会对检测孔边缘的图像点产生影响,所以,利用孔和小管的几何关系与位置关系,解出被测孔的直径尺寸,消除了孔边缘倒角对检测孔边缘图像点产生的影响,省略了提取孔边缘图像位置和拟合孔边缘曲线方程的过程,减少了由于图像成像和曲线拟合带来的误差,简化了测量过程,改善了测量精度。
附图说明
此处所说明的附图用来提供对本发明的进一步理解,构成本申请的一部分,本发明的示意性实例及其说明用于解释本发明,并不构成对本发明的不当限定。
图1为本发明的标定CCD摄像机的内部参数和畸变系数的结构示意图;
图2为本发明的标定CCD摄像机内部参数和畸变系数的理论图像平面和相关坐标系;
图3为本发明设定的世界坐标系I;
图4为本发明的三根相同小圆管在零件孔内的安装示意图;
图5为本发明的世界坐标系II的X-Y轴;
图6为本发明零件孔内三根小圆管的几何关系。
图中:1、背景光源;2、标定板;3、镜头;4、CCD摄像机;5、计算机;6、支架;{Ow,Xw,Yw,Zw}为标定用世界坐标系;{Oc,Xc,Yc,Zc}为摄像机坐标系;{Oi,x,y}为理论图像坐标系;{Op,u,v}为图像像素坐标系。
具体实施方式
下面结合附图进一步说明本发明的详细内容及其具体实施方式。
参见图1至图6所示,本发明的基于机器视觉技术的零件孔径尺寸测量方法,利用设定世界坐标系I形成的条件,建立世界坐标系I和摄像机坐标系的变换关系方程,从而建立孔径测量模型;本方法首创性地提出,在零件孔内放置三根尺寸已知的相同小圆管并采集零件孔端面图像,通过在理论图像平面的曲线拟合,得到小圆管端面内、外边缘的曲线方程,然后根据相邻小圆管建立新的世界坐标系II、III,在不同世界坐标系下求出三根小圆管端面中心点的世界坐标,再利用三根小圆管端面中心点的世界坐标和小圆管端面边缘曲线不变量的特点,求出世界坐标系和摄像机坐标系之间的变换关系,该求解世界坐标系和摄像机坐标系变换关系的方式完全规避了旧有方法采用在孔表面粘贴面积较大标定纸或通过工装夹具使用标定板等辅助方式标定零件孔表面的外部参数,然后求解世界坐标系和摄像机坐标系的变换关系时额外引入的误差,有效地解决了现有方法不易标定零件孔表面外部参数且精度不高的问题;在测量零件孔径尺寸时,根据世界坐标系II、III和摄像机坐标系之间的变换关系,先将三根小圆管端面中心点坐标变换到同一世界坐标系,然后根据零件孔与三根小圆管相内切,三根小圆管之间相外切的关系求出零件孔表面圆中心点的世界坐标,再结合任意一根小圆管端面中心点的世界坐标求出零件的孔径尺寸,该求解零件孔径尺寸的方法相对于现有公知的拍摄零件孔表面图像,通过拟合零件孔边缘曲线求出零件孔径尺寸的方法而言,完全避免了在提取零件孔边缘图像位置时,由于零件孔边缘倒角而产生的误差,从而忽略了零件孔边缘倒角对零件孔径尺寸测量结果的不利影响,简化了测量过程,改善了测量精度。
本发明基于机器视觉的测量模型建立了零件孔径的测量模型,采用测量尺寸已知的三根相同小圆管端面中心点坐标的方法,标定零件孔表面的外部参数,并利用零件孔内各小圆管间的几何关系计算零件孔的直径尺寸,结合边缘检测技术、角点检测技术、二次曲线不变量算法、曲线拟合技术等方法实现了零件孔径尺寸的非接触测量,具体过程包括以下步骤:
步骤1.设定世界坐标系I,建立孔径测量模型。具体过程包括以下步骤:
1.1)设定理论图像平面和相关坐标系,先标定CCD摄像机4的内部参数和畸变系数,然后保持CCD摄像机4的位置不变,将零件孔水平放置并安装好。
主要是基于张正友提出的摄像机平面标定算法,利用高精度平面标定板不同姿态的图像,标定出摄像机内部参数和镜头3的畸变系数。
如图1所示,先把CCD摄像机4固定在支架6上,并与计算机5相连,然后调节背景光源1的亮度,照射标定板2,利用CCD摄像机4采集9幅不同位姿的标定板图像。在标定CCD摄像机4的内参及畸变系数的过程中,主要涉及四个坐标系,分别为世界坐标系{Ow,Xw,Yw,Zw}、摄像机坐标系{Oc,Xc,Yc,Zc}、理论图像坐标系{Oi,x,y}和图像像素坐标系{Op,u,v}。具体坐标系的设定和位置关系如图2所示;
检测出标定板中平面方块角点的亚像素坐标,根据张正友提出的摄像机平面标定算法,利用检测得到的角点坐标求解摄像机内部参数和畸变系数,最后利用不同位姿的标定板提取到的角点坐标,对摄像机的内部参数矩阵和畸变系数进行优化。
1.2)设定世界坐标系I,让世界坐标系I的原点与摄像机坐标系的原点重合,Z轴垂直于零件孔表面,X-Y平面平行于零件孔表面。
在摄像机坐标系下,设孔表面平面方程为:
A1x+B1y+C1z+1=0 (1)
则Z轴的方向向量为(A1,B1,C1)。
设与孔表面平面相交的平面H的方程表达式为:
A2x+B2y+C2z+1=0 (2)
则以平面(1)、(2)的交线为轴线的平面束为:
λ(A2x+B2y+C2z+1)+(A1x+B1y+C1z+1)=0 (3)
在(3)中,一定存在一个平面J与理论图像平面垂直,利用垂直关系可以解出则可以解出λ=-C1/C2,代入公式(3)可以解出平面J的空间表达式:
(A1C2-C1A2)x+(B1C2-C1B2)y+(C2-C1)=0 (4)
以平面(1)和(4)的交线作为世界坐标系I的X轴方向,其方向向量为Ax=C1(C1B2-B1C2),Bx=C1(A1C2-C1A2),Cx=A1(B1C2-C1B2)-B1(A1C2-C1A2)。
按空间笛卡尔直角坐标系的建立规则,X轴与Z轴叉乘的方向为Y轴方向,如图3所示。
在世界坐标系I下,设零件孔表面圆中心点O0的坐标为(X0,Y0,Z0),孔的半径为r,则零件孔的空间表达式为:
1.3)在不标定零件孔表面外部参数的情况下,利用设定世界坐标系I的形成的条件,建立世界坐标系I和摄像机坐标系的变换关系方程,然后建立孔径测量模型。
在摄像机坐标系下,根据求得的世界坐标系I各坐标轴的方向向量可以得到各坐标轴的方向余弦:
其中:
(e11,e21,e31)是世界坐标系I的X轴在摄像机坐标系的方向余弦;
(e12,e22,e32)是世界坐标系I的Y轴在摄像机坐标系的方向余弦;
(e13,e23,e33)是世界坐标系I的Z轴在摄像机坐标系的方向余弦;
式(7)中,Ax=C1(C1B2-B1C2),Bx=C1(A1C2-C1A2),Cx=A1(B1C2-C1B2)-B1(A1C2-C1A2);
式(8)中,ey1=(e23e31-e33e21),ey2=(e33e11-e13e31),ey3=(e13e21-e23e11)。
由公式(6)至公式(8)得到世界坐标系I和摄像机坐标系的变换关系:
由于世界坐标系I和摄像机坐标系的原点重合,所以[T]=0。由于旋转矩阵E是对称矩阵,所以有:
设世界坐标系I下点的坐标表示为(X,Y,Z),摄像机坐标系下点的坐标表示为(x,y,z)。
在世界坐标系I下,以零件孔表面圆作为底圆,建立以世界坐标系I的原点作为顶点的斜圆锥面方程为:
斜圆锥面方程(11)与Z=Z0的平面的交线是世界坐标系I下的一个平面圆R,其方程为:
(X-X0)2+(Y-Y0)2=r2 (12)
将公式(12)写成矩阵的形式为:
将世界坐标系I下的斜圆锥面方程(11)写成矩阵的形式为:
将世界坐标系I下的自变量转换为摄像机坐标系下的自变量:
根据式(13)、(16)可知,在摄像机坐标系下的斜圆锥面方程和理论图像平面的交线即为理论图像坐标系下的零件孔边缘图像投影曲线方程S:
式(17)中,
它们是零件孔边缘图像投影曲线方程S一般表达式的系数。
则在理论图像平面,零件孔边缘图像投影曲线方程S的二次曲线不变量为:
步骤2.将三根尺寸已知的相同小圆管插入零件孔内使四者端面共面并拍摄端面图像,在理论图像平面拟合三根小圆管端面的内、外边缘曲线,由于没有标定零件孔表面的外部参数,所以先根据相邻小圆管重建世界坐标系,然后在新世界坐标系下根据相邻小圆管的位置关系求出三根小圆管端面中心点O1、O2、O3的世界坐标。具体过程包括以下步骤:
2.1)准备三根尺寸已知的相同小圆管A、B、C,且其尺寸能够保证三根小圆管轻松放置在零件孔内。
2.2)将三根小圆管插入待检测的零件孔内,并使其四者的端面共面,其中A为底部中间的小圆管,其两侧小圆管B和小圆管C由于重力作用均与小圆管A相外切。设小圆管A、B间的切点为f,小圆管A、C间的切点为g,三根小圆管在零件孔内的安装示意图如图4所示。
2.3)用CCD摄像机采集放有三根小圆管的待检测零件孔端面图像一幅,并检测小圆管端面图像的内、外边缘点的亚像素坐标,然后在理论图像平面拟合三根小圆管的端面内、外边缘曲线方程。
根据采集的边缘点亚像素坐标,通过曲线拟合,可以得到在理论图像平面下的每根小圆管端面外边缘二次曲线方程的系数mx_1=[ax_1 bx_1 cx_1 dx_1 ex_1 fx_1]和内边缘二次曲线方程的系数:mx_2=[ax_2 bx_2 cx_2 dx_2 ex_2 fx_2]。
2.4)取两根相邻小圆管A和B,在以其端面中心点O1和O2的连线为轴线的无穷多个平面束中,一定存在一个平面与理论图像平面垂直,将该平面与世界坐标系I的X-Y平面的交线作为世界坐标系II的X轴方向,其余条件与世界坐标系I相同,新建世界坐标系II。
让世界坐标系II的原点与摄像机坐标系的原点重合,X-Y平面仍与零件孔表面平行,世界坐标系II的Z轴方向仍垂直于零件孔表面。在以两根小圆管A、B的端面中心点O1和O2的连线为轴线的无穷多个平面束中,一定存在一个平面与理论图像平面垂直,将该平面与世界坐标系I的X-Y平面的交线作为世界坐标系II的X轴方向,利用空间笛卡尔直角坐标系的建立规则,X轴与Z轴叉乘的方向为Y轴的方向,如图5所示,建立新的世界坐标系II。这样做可以为求解小圆管端面中心点的世界坐标提供约束条件。
2.5)利用在世界坐标系II下小圆管A和B相外切的关系以及二次曲线不变量的特点,建立方程组求解两根小圆管A和B端面中心点O1和O2的世界坐标;
设世界坐标系II下,小圆管A和B的端面中心点O1和O2的坐标分别为(X0_1,Y0_1,Z0_1)和(X0_2,Y0_2,Z0_2),根据图5可知,O1和O2的坐标满足:
式(20)中,r1_1和r2_1是两根小圆管A和B的端面外圆的半径。
在理论图像坐标系下,根据公式(18)与拟合曲线和实际曲线之间的二次曲线不变量的关系,可以得到关于小圆管A的等式:
方程组(22)中:r1_2是小圆管A的端面内圆半径,k1_1和k1_2是比例系数,D1_1和D1_2分别是小圆管A的端面外边缘曲线不变量和内边缘曲线不变量,即
整理方程组(21)、(22)可得:
解(23)式得:
将(24)式代入方程组(21)、(22)可得:
再将(25)式代入方程组(21)可得:
同理,对小圆管B做相同处理,可以得到:
由于Z0_1=Z0_2,所以联立公式(20)、(27)和(28),可以得到求解两根小圆管A、B端面中心点O1、O2坐标的方程组:
通过解方程组(29)就可得到两根小圆管A、B端面中心点O1、O2在世界坐标系II下的坐标。
2.6)重复步骤2.4)和2.5),以两根相邻小圆管A和C同理建立世界坐标系III,并求得小圆管C的端面中心点O3在世界坐标系III下的坐标。
步骤3.利用三根小圆管端面中心点的世界坐标和小圆管端面边缘曲线不变量的特点,求出世界坐标系和摄像机坐标系的变换关系。具体过程包括以下步骤:
3.1)在世界坐标系II下,为求出世界坐标系II和摄像机坐标系的变换矩阵Q,利用小圆管端面边缘曲线的系数和不变量的特点,建立小圆管端面中心点坐标和变换矩阵Q中元素之间的关系。
可以任取小圆管A和B中的一根,以小圆管A为例,在理论图像坐标系下,通过曲线拟合,得到小圆管A端面外边缘曲线的拟合系数为mn=[an bn cn dn en fn],在理论图像坐标系下,小圆管A的端面外边缘曲线实际系数的表达式如式(30)所示。
在理论图像坐标系下,根据拟合曲线和实际曲线的二次曲线不变量的关系可知,拟合的曲线系数mn=[an bn cn dn en fn]和实际的曲线系数mx=[ax bx cx dx ex fx]之间相差一个比例系数k1,即
mx=[ax bx cx dx ex fx]=k1·[an bn cn dn en fn] (31)
由于小圆管A的部分相关参数已经求出,所以通过式(21)可以求出该比例系数k1。
通过曲线拟合,得到小圆管A的端面外边缘曲线方程系数an-r1和内边缘曲线方程系数an-r2,再根据式(21)、(22)求得比例系数k1、k2,然后可以得到小圆管A的端面内、外边缘曲线的实际系数:
ax-r1=k1an-r1,ax-r2=k2an-r2 (32)
3.2)利用步骤3.1中建立的关系,求解世界坐标系II和摄像机坐标系的变换矩阵Q。
根据式(30)可得
联立式(33)、(34)可得:
通过联立式(35)和方程组(36)解得的解不唯一,共有四组解。
同理,将式(37)、(38)按相同的方式求解。
可以得到e21、e22、e23和e31、e32、e33的各四组解。
3.3)通过步骤3.2)可以得到变换矩阵Q的多组解,利用二次曲线不变量的表达式和零件孔与摄像机的相对位置对解集进行筛选,可以得到变换矩阵Q的唯一解,从而得到世界坐标系II和摄像机坐标系的变换关系。
由于得到的坐标系变换矩阵Q共有64组解,所以在解集中挑选同时满足式(30)和通过式(30)变换得到的式(39)、(40)、(41)的解。
通过筛选,可以剩下两组解Q1和Q2。
将采集的零件孔中小圆管A的端面图像数据分别按照两组坐标变换矩阵的解Q1和Q2进行坐标变换,可以得到两组数据,通过数据的分布可以看出,利用Q1转换的数据点分布在世界坐标系II下Z轴的正方向,而利用Q2转换的数据点分布在世界坐标系II下Z轴的负方向。根据世界坐标系II的建立方式和零件孔与摄像机的实际位置关系可知,采集的图像数据点均应在世界坐标系II的Z轴正方向,所以,排除坐标变换矩阵Q2,得到唯一的坐标系变换矩阵Q1,求得了世界坐标系II和摄像机坐标系的变换关系。
3.4)重复步骤3.1)、3.2)、3.3),同理求出世界坐标系III和摄像机坐标系的变换关系。
步骤4.利用零件孔与三根小圆管的几何关系,求出零件孔径尺寸。具体过程包括以下步骤:
4.1)根据世界坐标系II、III和摄像机坐标系之间的变换关系,将处于不同世界坐标系下三根小圆管端面中心点O1、O2、O3的坐标变换到同一世界坐标系II下。
根据已经求得的坐标系变换关系,设世界坐标系II表示为{Ow,Xw1,Yw1,Zw1},世界坐标系III表示为{Ow,Xw2,Yw2,Zw2},两个世界坐标系与摄像机坐标系的变换关系为:
将式(43)、(44)联立求解两个世界坐标系的变换关系如式(45)所示,根据式(45)可以将不同世界坐标系下三根小圆管的端面中心点O1、O2、O3的坐标变换到世界坐标系II。
4.2)由于零件孔与三根小圆管相内切,三根小圆管之间相外切,所以,两组相邻小圆管端面中心点的连接线段O1O2和O1O3的垂直平分线相交于零件孔表面圆的中心点O0。
如图6所示,根据三根小圆管和零件孔的几何关系可知,相邻小圆管端面中心点连接线段O1O2和O1O3的垂直平分线的交点是零件孔表面圆的中心点O0。
4.3)利用步骤4.2)的几何关系,建立方程组求出零件孔表面圆中心点O0在世界坐标系II下的坐标。
在世界坐标系II下,设零件孔表面中心点的坐标为O0(xow,yow,zow),三根小圆管A、B、C的端面中心点坐标分别为O1(xo1,yo1,zo1),O2(xo2,yo2,zo2),O3(xo3,yo3,zo3),根据直线的点斜式方程,求出O1O2和O1O3的垂直平分线的方程组:
根据式(46)、(47)可以解出零件孔表面中心点O0的世界坐标:
4.4)在世界坐标系II下,根据零件孔表面圆中心点的坐标和任意一根小圆管端面中心点的坐标,利用零件孔与三根小圆管相内切的位置关系,求出零件的孔径尺寸。
根据图6可以看出,在世界坐标系II下,零件孔表面的中心点与任意一根小圆管端面中心点的距离加上小圆管的半径即为零件孔的半径,其二倍为零件孔的内径。以小圆管A为例:
式(50)中,d为零件孔的内径,r1-1为小圆管A的外径,孔表面中心点O0的坐标为(xow,yow,zow),小圆管A端面中心点O1的坐标为(xo1,yo1,zo1)。
根据式(50)即可求出零件孔径尺寸。
以上所述仅为本发明的优选实例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡对本发明所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (3)
1.一种基于机器视觉技术的零件孔径尺寸测量方法,其特征在于:采用在零件孔内放置三根尺寸已知的相同小圆管的方式,实现了零件孔径尺寸的非接触测量,具体过程包括以下步骤:
步骤(1)设定世界坐标系I,建立孔径测量模型;
(1.1)设定理论图像平面和坐标系,先标定CCD摄像机的内部参数和畸变系数,然后保持CCD摄像机的位置不变,将零件孔水平放置并安装好;
(1.2)设定世界坐标系I,让世界坐标系I的原点与摄像机坐标系的原点重合,Z轴垂直于零件孔表面,X-Y平面平行于零件孔表面;
(1.3)在不标定零件孔表面外部参数的情况下,利用设定世界坐标系I的形成的条件,建立世界坐标系I和摄像机坐标系的变换关系方程,然后建立孔径测量模型;
步骤(2)将三根尺寸已知的相同小圆管插入零件孔内使四者端面共面并拍摄端面图像,在理论图像平面拟合三根小圆管端面的内、外边缘曲线,由于没有标定零件孔表面的外部参数,所以先根据相邻小圆管重建世界坐标系,然后在新世界坐标系下根据相邻小圆管的位置关系求出三根小圆管端面中心点O1、O2、O3的世界坐标,具体步骤如下:
(2.1)准备三根尺寸已知的相同小圆管A、B、C,且其尺寸能够保证三根小圆管轻松放置在零件孔内;
(2.2)将三根小圆管插入待检测的零件孔内,且三根小圆管和零件孔四者的端面共面,其中A为底部中间的小圆管,其两侧小圆管B和小圆管C由于重力作用均与小圆管A相外切;
(2.3)用CCD摄像机采集放有三根小圆管的待检测零件孔端面图像一幅,并检测小圆管端面图像的内、外边缘点的亚像素坐标,然后在理论图像平面拟合三根小圆管的端面内、外边缘曲线方程;
(2.4)取两根相邻小圆管A和B,在以其端面中心点O1和O2的连线为轴线的无穷多个平面束中,一定存在一个平面与理论图像平面垂直,将该平面与世界坐标系I的X-Y平面的交线作为世界坐标系II的X轴方向,其余条件与世界坐标系I相同,新建世界坐标系II;
(2.5)利用在世界坐标系II下小圆管A和B相外切的关系以及二次曲线不变量的特点,建立方程组求解两根小圆管A和B端面中心点O1和O2在世界坐标系II下的坐标;
(2.6)重复步骤(2.4)和(2.5),以两根相邻的小圆管A和C同理建立世界坐标系III,并求得小圆管C的端面中心点O3在世界坐标系III下的坐标;
步骤(3)利用三根小圆管端面中心点的世界坐标和小圆管端面边缘曲线不变量的特点,求出世界坐标系II、III和摄像机坐标系的变换关系;
步骤(4)利用零件孔与三根小圆管的几何关系,求出零件孔径尺寸。
2.根据权利要求1所述的基于机器视觉技术的零件孔径尺寸测量方法,其特征在于:步骤(3)所述的利用三根小圆管端面中心点的世界坐标和小圆管端面边缘曲线不变量的特点,求出世界坐标系和摄像机坐标系的变换关系,具体步骤如下:
(3.1)在世界坐标系II下,为求出世界坐标系II和摄像机坐标系的变换矩阵Q,利用小圆管端面边缘曲线的系数和不变量的特点,建立小圆管端面中心点坐标和变换矩阵Q中元素之间的关系;
(3.2)利用步骤(3.1)中建立的关系,求解世界坐标系II和摄像机坐标系的变换矩阵Q;
(3.3)通过步骤(3.2)得到变换矩阵Q的多组解,利用二次曲线不变量的表达式和零件孔与摄像机的相对位置对解集进行筛选,得到变换矩阵Q的唯一解,从而得到世界坐标系II和摄像机坐标系的变换关系;
(3.4)重复步骤(3.1)至步骤(3.3),同理求出世界坐标系III和摄像机坐标系的变换关系。
3.根据权利要求1所述的基于机器视觉技术的零件孔径尺寸测量方法,其特征在于:步骤(4)所述的利用零件孔与三根小圆管的几何关系,求出零件孔径尺寸,具体步骤如下:
(4.1)根据世界坐标系II、III和摄像机坐标系之间的变换关系,将处于不同世界坐标系下三根小圆管端面中心点O1、O2、O3的坐标变换到同一世界坐标系II下;
(4.2)由于零件孔与三根小圆管相内切,三根小圆管之间相外切,所以,两组相邻小圆管端面中心点的连接线段O1O2和O1O3的垂直平分线相交于零件孔表面圆的中心点O0;
(4.3)利用步骤(4.2)的几何关系,建立方程组求出零件孔表面圆中心点O0在世界坐标系II下的坐标;
(4.4)在世界坐标系II下,根据零件孔表面圆中心点的坐标和三根小圆管中任意一根的端面中心点的坐标,利用零件孔与三根小圆管相内切的位置关系,求出零件的孔径尺寸。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910201755.0A CN109974608B (zh) | 2019-03-18 | 2019-03-18 | 基于机器视觉技术的零件孔径尺寸测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910201755.0A CN109974608B (zh) | 2019-03-18 | 2019-03-18 | 基于机器视觉技术的零件孔径尺寸测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109974608A CN109974608A (zh) | 2019-07-05 |
CN109974608B true CN109974608B (zh) | 2020-11-03 |
Family
ID=67079191
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910201755.0A Active CN109974608B (zh) | 2019-03-18 | 2019-03-18 | 基于机器视觉技术的零件孔径尺寸测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109974608B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110553586B (zh) * | 2019-09-11 | 2020-12-25 | 吉林大学 | 基于机器视觉检测物体表面数据点集的方法 |
CN113063351B (zh) * | 2021-03-23 | 2023-03-14 | 江南造船(集团)有限责任公司 | 管尺寸偏差自动测量系统及测量方法 |
CN113639635B (zh) * | 2021-08-09 | 2023-02-10 | 温州市工业科学研究院 | 一种基于图像识别技术的零件测量方法 |
CN114279325B (zh) * | 2021-12-16 | 2024-02-09 | 天津大学 | 视觉测量模块测量坐标系空间位置关系标定系统及方法 |
CN116188472B (zh) * | 2023-05-04 | 2023-07-07 | 无锡康贝电子设备有限公司 | 一种数控机床零件的在线视觉检测方法 |
CN118328917B (zh) * | 2024-06-12 | 2024-10-08 | 中国航空工业集团公司金城南京机电液压工程研究中心 | 一种带镀层的波纹管类零件几何尺寸自动化测量方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH11176891A (ja) * | 1997-12-11 | 1999-07-02 | Canon Inc | ボンディングワイヤ検査方法およびその装置 |
CN1566900A (zh) * | 2003-06-11 | 2005-01-19 | 北京航空航天大学 | 一种空间圆几何参数的视觉测量方法 |
KR101295330B1 (ko) * | 2011-12-20 | 2013-08-12 | 충북대학교 산학협력단 | 레이저 조사장치를 이용한 직경 측정장치 및 측정방법 |
CN106152956A (zh) * | 2016-09-23 | 2016-11-23 | 苏州蓝王机床工具科技有限公司 | 一种激光内径测量仪 |
CN206847565U (zh) * | 2016-11-29 | 2018-01-05 | 四川凌峰航空液压机械有限公司 | 非接触测量深小孔径尺寸的测量装置 |
CN108180851A (zh) * | 2017-12-22 | 2018-06-19 | 中国航空工业集团公司北京航空精密机械研究所 | 一种用于测量气膜孔形位参数的五轴影像测量装置 |
CN108562233A (zh) * | 2018-01-15 | 2018-09-21 | 吉林大学 | 利用二次曲线不变量的轴类零件直径尺寸在线测量方法 |
CN108871216A (zh) * | 2018-07-12 | 2018-11-23 | 湘潭大学 | 一种基于视觉引导的机器人多孔接触式自动测量方法 |
-
2019
- 2019-03-18 CN CN201910201755.0A patent/CN109974608B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH11176891A (ja) * | 1997-12-11 | 1999-07-02 | Canon Inc | ボンディングワイヤ検査方法およびその装置 |
CN1566900A (zh) * | 2003-06-11 | 2005-01-19 | 北京航空航天大学 | 一种空间圆几何参数的视觉测量方法 |
KR101295330B1 (ko) * | 2011-12-20 | 2013-08-12 | 충북대학교 산학협력단 | 레이저 조사장치를 이용한 직경 측정장치 및 측정방법 |
CN106152956A (zh) * | 2016-09-23 | 2016-11-23 | 苏州蓝王机床工具科技有限公司 | 一种激光内径测量仪 |
CN206847565U (zh) * | 2016-11-29 | 2018-01-05 | 四川凌峰航空液压机械有限公司 | 非接触测量深小孔径尺寸的测量装置 |
CN108180851A (zh) * | 2017-12-22 | 2018-06-19 | 中国航空工业集团公司北京航空精密机械研究所 | 一种用于测量气膜孔形位参数的五轴影像测量装置 |
CN108562233A (zh) * | 2018-01-15 | 2018-09-21 | 吉林大学 | 利用二次曲线不变量的轴类零件直径尺寸在线测量方法 |
CN108871216A (zh) * | 2018-07-12 | 2018-11-23 | 湘潭大学 | 一种基于视觉引导的机器人多孔接触式自动测量方法 |
Non-Patent Citations (3)
Title |
---|
Superresolution Restoration of an Image Sequence: Adaptive Filtering Approach;Michael Elad等;《IEEE TRANSACTIONS ON IMAGE PROCESSING》;19991231;第8卷(第3期);全文 * |
基于机器视觉的零件尺寸测量;田原嫄等;《激光与光电子学进展》;20101231;全文 * |
基于线结构光视觉技术的微小直径高精度测量系统;刘斌等;《仪器仪表学报》;20141231;第35卷(第12期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN109974608A (zh) | 2019-07-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109974608B (zh) | 基于机器视觉技术的零件孔径尺寸测量方法 | |
CN110057295B (zh) | 一种免像控的单目视觉平面距离测量方法 | |
CN103499302B (zh) | 基于结构光视觉成像系统的凸轮轴直径尺寸在线测量方法 | |
CN111369630A (zh) | 一种多线激光雷达与相机标定的方法 | |
CN103471531B (zh) | 轴类零件直线度在线非接触测量方法 | |
CN110455198B (zh) | 基于线结构光视觉的矩形花键轴键宽及直径测量方法 | |
CN105716542B (zh) | 一种基于柔性特征点的三维数据拼接方法 | |
CN109242908B (zh) | 用于水下双目视觉测量系统的标定方法 | |
CN108562233B (zh) | 利用二次曲线不变量的轴类零件直径尺寸在线测量方法 | |
CN111536902A (zh) | 一种基于双棋盘格的振镜扫描系统标定方法 | |
CN109870111B (zh) | 基于机器视觉的渐开线圆柱齿轮齿距累积误差测量方法 | |
CN103615980B (zh) | 一种板件上圆孔参数的测量方法及系统 | |
CN110455225B (zh) | 基于结构光视觉的矩形花键轴同轴度及键位置度测量方法 | |
CN109974618B (zh) | 多传感器视觉测量系统的全局标定方法 | |
CN103278138A (zh) | 一种复杂结构薄部件三维位置及姿态的测量方法 | |
CN107339935B (zh) | 用于全视角扫描测量系统的靶标空间交会测量方法 | |
WO2024138916A1 (zh) | 一种2d面阵相机与线激光3d传感器联合标定方法及装置 | |
CN112258583B (zh) | 基于等畸变量分区的近景图像畸变标定方法 | |
CN114331924B (zh) | 大工件多相机视觉测量方法 | |
CN111024047A (zh) | 基于正交双目视觉的六自由度位姿测量装置及方法 | |
CN115235379A (zh) | 一种单目线激光三维视觉传感器参数在位标定装置及方法 | |
CN110428471B (zh) | 一种针对光学自由曲面子孔径偏折测量的精确自定位方法 | |
CN114001651B (zh) | 一种基于双目视觉测量和先验检测数据的大型细长筒类构件位姿原位测量方法 | |
CN111968182B (zh) | 一种双目相机非线性模型参数的标定方法 | |
CN114279325A (zh) | 视觉测量模块测量坐标系空间位置关系标定系统及方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |