CN102353340A - 缸盖毛坯加工尺寸识别方法及装置 - Google Patents

缸盖毛坯加工尺寸识别方法及装置 Download PDF

Info

Publication number
CN102353340A
CN102353340A CN201110151060XA CN201110151060A CN102353340A CN 102353340 A CN102353340 A CN 102353340A CN 201110151060X A CN201110151060X A CN 201110151060XA CN 201110151060 A CN201110151060 A CN 201110151060A CN 102353340 A CN102353340 A CN 102353340A
Authority
CN
China
Prior art keywords
ccd video
beta
video camera
coordinate
point
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.)
Granted
Application number
CN201110151060XA
Other languages
English (en)
Other versions
CN102353340B (zh
Inventor
王天友
鲁祯
张瑞峰
李勤
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tianjin University
Original Assignee
Tianjin University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Tianjin University filed Critical Tianjin University
Priority to CN 201110151060 priority Critical patent/CN102353340B/zh
Publication of CN102353340A publication Critical patent/CN102353340A/zh
Application granted granted Critical
Publication of CN102353340B publication Critical patent/CN102353340B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

本发明公开了一种缸盖毛坯加工尺寸识别方法及装置。具体是由2个或2个以上CCD摄像机识别接收标靶图像,CCD摄像机镜头视野包含整个标靶平面。对缸盖毛坯加工尺寸的测取定位主要包括:CCD摄像机标定;特征信息提取;立体匹配以及三维重建。本发明以计算机和CCD摄像机作为主要硬件,由摄像机拍摄工件图像后将数据传输至计算机,通过提取、计算图像的边缘、关键线面及缸盖毛坯各气道口的空间相对位置等信息,与标准模型数据进行比对,做出该毛坯气道位置是否合格的判定。本发明在机加工之前就可判定气道位置和形状尺寸,及时发现并剔出不合格的气缸盖,防止该毛坯件进入机加工工序,避免经过一系列加工后而所造成工件不合格的经济损失。

Description

缸盖毛坯加工尺寸识别方法及装置
技术领域
[0001] 本发明涉及到一种对内燃机缸盖毛坯进行加工时,确定其尺寸的定位方法。
技术背景
[0002] 目前在铸造、机加工等生产过程中很难做到精确定位的加工,例如在发动机气缸盖制造的几个环节,如沙型,铸造,机械加工等,可能会因为定位的误差造成气道尺寸或位置的偏差,而发动机气道的尺寸直接影响着发动机缸内气体的流动和燃烧状况。当前进行气道稳流试验都是在机加工完成后进行,如果此时气道参数不能满足要求,将直接影响发动机性能,造成材料及加工成本的损失。因此,如果能在缸盖铸造成型尚未进入机加工生产线时就能够测量出气道的形位尺寸,并确定出气道参数,及时发现并剔出不合格的气缸盖,这对于加工生产过程是非常重要的。可以节约机加工成本,并且避免不合格产品流入下一道工序,同时将缺陷信息反馈给铸造工艺,使其对模具及铸造做出相应调整,实现闭环控制。
[0003] 鉴于此,本发明旨在开发一种基于机器视觉及人工神经网络技术的识别系统,在缸盖毛坯进行机加工前即可测取其气道位置和几何尺寸。
发明内容
[0004] 本发明的目的在于,提供一种适用于对内燃机缸盖毛坯件尺寸进行预先判定的方法及装置。
[0005] 以下对本发明的技术方法及计算过程予以说明。缸盖毛坯加工尺寸识别方法及装置,主要包括CXD摄像机、托板、支撑架、标靶以及计算机等。其设备组成是:由2个或2个以上CCD摄像机构成识别接收器,识别接收器接收标靶图像,CCD摄像机镜头视野包含整个标靶平面,标靶是平面阵列的圆孔图案,或者是黑白平面的方块阵列。
[0006] 对缸盖毛坯加工尺寸的测取定位包括如下步骤:
[0007] (I)CCD摄像机标定
[0008] 获取两个CXD摄像机的内部参数,
Figure CN102353340AD00061
[0009] 首先对两个CXD摄像机之间的旋转平移矩阵进行标定,2D点表示为m = [u,ν]τ, 3D点表示为M = [X,Y,Z]T,2D和3D的增广矩阵:
[0010]
Figure CN102353340AD00062
[0011]
[0012] 3D点M与它投影图像点之间的关系为:
[0013][0015]
Figure CN102353340AD00071
[0014] 其中,M是世界坐标系中的三维坐标,s是任意标准矢量,(R,t)定义为外部参数, 表示世界坐标系和CCD摄像机坐标系间的旋转与平移关系,CCD摄像机内部参数用A矩阵表不:
Figure CN102353340AD00072
[0016] 其中,(uQ,v0)是坐标系原点,a和β是图像上u和ν坐标轴的标准矢量,Y表示两坐标的垂直度,由于采用平面标靶,所以假设Z = 0,则由表达式(d)得到:
[0017]
Figure CN102353340AD00073
[0018] 其中,议MHHF是标靶平面上的齐次坐标,抓=[Hi]『为标靶平面上的点投影到图像平面上对应点的齐次坐标,[ri r2 r3]和t分别是C⑶摄像机坐标系相对于世界坐标系的旋转矩阵及平移向量,用一个3X3的映射矩阵:
[0019] H= Lh1 h2 h3] = AAtr1 r2 t]
[0020] 得到内部参数矩阵A的两个约束条件
[0021] Zz1H1Zz2=O (g)
[0022] ^A-tA'1 H1 = hT2A-TA-\ (h)
[0023]设:
[0024]
Figure CN102353340AD00074
[0025]
[0026] B是对称矩阵,表示一个六位向量b = [B11, B12,B22, B13,B23, B33]τ,基于绝对二次曲线原理求出B,并对B求逆,利用Cholesky分解,从B中倒出内部参数矩阵Α,由内部参数矩阵A和映射矩阵H计算每幅图像相对于平面标靶的CCD摄像机外部参数矩阵R和平移向量 t :
[0027] Γ! = λ A^h1
[0028] r2 = λ A_1h2
[0029] r3 = T1Xr2 (j)[0030]
t = λ
[0031] 将(j)作为初值,Γι、r2、r3> t表示世界坐标系和CXD摄像机坐标系间的旋转与平移关系,通过最大似然估计优化求解其参数,目标函数为:
[0032]
Figure CN102353340AD00081
[0033] 其中为第i幅图像第j个点图像的齐次坐标,m为第i幅图像第j个点映射出的图像的齐次坐标,A为内参,Ri和Ti是第i幅图像对于CCD摄像机坐标系的旋转和平移矩阵,Mj为第i幅图像的第j个世界坐标系的三维点,目标函数趋于0,计算出CCD摄像机的准确参数值。径向畸变参数求解,由径向畸变参数模型:
[0034]
Figure CN102353340AD00082
[0037] 其中(x,y)为理想的无畸变的图像坐标,(二^为实际的畸变图像坐标,ki、!^为径向畸变参数;
[0038]
[0039]
[0040]
[0041]
Y表示两坐标的垂直度,故设Y = 0,由:
Figure CN102353340AD00083
[0042] 得到
[0043]
[0044]
[0045]
[0046]
得如下方程
[0047]
Figure CN102353340AD00084
其中(u,ν)是(k)式求出的内参,计算出的理想投影点的像素坐标,由(U,V),获
[0048] 对于η幅图片的m个点,共有2XnXm个方程,表示为Dk = d
[0049]其中:[0050]
[0051]
[0052]
k =
D =
'K
_众2_
(u-u0)(x2+y2)(u-u0)(x2+ y2)2 (V-V0Xx2+/)(v-v0)(x2+/)2
d =
u-u
V-V
[0053] k= (DTD) _1DTd
[0054] 目标函数为:
[0055]
Figure CN102353340AD00091
[0056] 由此标定出CXD摄像机的内参和径向畸变参数;
[0057] 立体标定依赖于查找两台C⑶摄像机之间的旋转矩阵R和平移向量T,由于该装置具有两个CCD摄像机,由q式可分别获得到左CCD摄像机的旋转矩阵Ii1和平移向量T1以及右CCD摄像机的旋转矩阵艮和平移向量八
[0058]
Figure CN102353340AD00092
[0059]
[0060] 由(r)式得到两个(XD摄像机的位置关系R禾口 T ;
[0061] (2)特征信息提取
[0062] 采用canny算法进行边缘检测、获取缸盖毛坯图形的轮廓信息,然后采用均值滤波计算排除掉噪声点,识别出关键点线,对经过边缘检测后的缸盖毛坯图像再应用均值滤波计算,滤掉其中的杂点,得到有效的关键点线;
[0063] (3)立体匹配
[0064] 1)消除畸变:采用步骤(1)中的(L)式消除径向方向上的镜头畸变,输出无畸变的图像;
[0065] 2)(XD摄像机校正:在物理意义上将CXD摄像机在数学意义上对准到同一个观察平面上,方法是通过数学的算法来实现使左右两幅图像的y坐标相等,这样左右两幅图像的视差只存在χ轴上,输出行对准的校正图像;
[0066] 3)图像匹配:查找两个CXD摄像机视场中的相同特征,算出视差;
[0067] (4)三维重建
[0068] 由两个CXD摄像机拍摄分别获取两张图像,由视差计算出缸盖毛坯进气端面轮廓特征图像点的世界坐标值,根据特征信息提取所得到的有效的关键点线,用最小二乘法做曲线拟合,获取缸盖毛坯进气端面特征轮廓的函数方程,将设计的缸盖进气端面加工位置点带入该函数方程,与轮廓特征点进行对比,确定被测缸盖毛坯加工位置以及相应尺寸。
[0069] 其过程原理是,通过缸盖端面的位置形状信息能推测出气道的尺寸及位置是否出现偏差,为了获取进缸盖端面的位置形状信息,用两个CCD摄像机并排连接构建双目测量系统。根据CCD摄像机所拍摄的两幅二维图像重建待测工件的三维几何形状,并实现对关键线面的精确测量,然后将测量数据与标准模型进行比对得出合格与否的结论。系统所涉及的关键技术包括:(XD摄像机标定,特征信息提取,立体匹配,三维重建。
[0070] 如果将CXD摄像机所拍摄到的图像还原为空间中实际的物体,那么CXD摄像机拍摄的图像与空间中的物体之间,存在一种线性关系:
[0071][像]=M[物]
[0072] 这里,矩阵M可看成是CXD摄像机成像的几何模型。M中的参数就是CXD摄像机参数,也称内参数。通常,这些参数是要通过实验与计算来得到的。这个求解参数的过程就称为CCD摄像机标定。标定的目的是为了获取由三维空间点(世界坐标系)到二维平面点 (图像坐标系)的投影变换关系,为此将涉及CCD摄像机模型的内、外参数。完成标定后以矩阵形式存储在系统中用作由两维数据来恢复三维立体图像过程中的转换标度。采取的技术方案是,先分别获得到两个CXD摄像机的内、外参数及畸变参数;再通过同一世界坐标中的一组定标点来建立两个CCD摄像机之间的位置关系。具体方法是将标靶(具有标准平面圆孔图案阵列,也可以采用标准平面黑白方块整列)置于世界坐标系中的已知位置,然后根据靶标角点(标靶方块图案交接点)与其在图像投影之间的映射关系计算出描述CCD摄像机光学特性和空间位置的各种参数。
[0073] 特征信息提取是对两个CXD摄像机所拍摄二维图像的数字处理,主要基于边缘检测和形态学图像处理技术。由于立体视觉系统是基于视差原理实现三维重建的,所以为了计算视差需要将三维空间中测量关注特征点,在不同图像中的映像点对应起来实现匹配。 为了提高在线状态的匹配速度,需要在离线状态下对两个CCD摄像机的空间几何结构通过图像校正变换成标准的外极线几何结构。立体匹配包括消除畸变、图像匹配、CCD摄像机校正,其中C⑶摄像机校正通过3个基本步骤完成:(1)从立体图像对中的一幅图像选择与实际物理结构相应的图像特征;(2)在另一幅图像确定出同一物理结构的对应图像特征,这是实现匹配的关键;C3)确定这两个特征之间的相对位置得到视差。
[0074] 获取空间任一点(与两个图像中对应)的坐标和满足两个摄像机的参数矩阵,即可进行该点的三维重建。
[0075] 根据双目立体视觉原理,用两个CCD摄像机模仿双眼,与被测物构成三角形,被测物在两个像面上成立体像对,然后根据匹配的同名像点,依据立体视差原理获取被测物的三维轮廓数据信息。本发明依据两CCD摄像机的内外参数,通过两个成像平面上的极线约束关系,建立对应点之间的联系,并由此联立方程,计算出图像点在实际物理空间的世界坐标值。在此基础上,测算出描述待测工件关键点线的离散数据,例如进气道入口端面内具有 4个方孔和3个圆孔,出口底面内具有的16个圆孔(如图6),然后采用三次样条插值的曲线拟合方法计算出相应的数学方程,用来与标准模型进行比对给出合格与否的判定。
[0076] 本发明的特点以及由此产生的积极效果是,在机加工之前就可判定气道位置和形状尺寸,及时发现并剔出不合格的气缸盖,防止该毛坯件进入机加工工序,避免经过一系列加工后而所造成工件不合格的经济损失。
附图说明
[0077] 图1是CXD摄像装置组件图。
[0078] 图2是标靶示意简图。[0079] 图3获取的缸盖进气端面轮廓特征示意图。
[0080] 图4获取的缸盖进气端面的真实轮廓参数点示意图。
[0081] 图5标准无缺陷的缸盖进气端面轮廓特征示意图。
[0082] 图6缸盖外形结构示意简图。
具体实施方式
Figure CN102353340AD00111
[0083] 以下结合附图和实施例对本发明的方法过程做进一步的说明。需要说明的是该实施例是叙述性的,而不是限定性的,不能以此实施例限定本发明的保护范围。
[0084] 缸盖毛坯加工尺寸识别方法及装置具有:(XD摄像机1、托板2、支撑架3、靶标4以及计算机等(如图1)。由2个或2个以上CCD摄像机构成识别接收器,识别接收器接收标靶图像,CCD摄像机镜头视野包含整个标靶平面,标靶是平面阵列的圆孔图案,或者是黑白平面的方块阵列(如图幻。两个CCD摄像机并列固定在托板上,托板焊接在支撑架的顶部, 支撑架的高度可以调节。
[0085] 缸盖毛坯加工尺寸的测取定位包括如下步骤:
[0086] (I)CCD摄像机标定
[0087] 获取两个CXD摄像机的内部参数,A =
[0088] 首先对两个CXD摄像机之间的旋转平移矩阵进行标定,2D点表示为m = [u,ν]τ, 3D点表示为M = [X,Y,Z]T,2D和3D的增广矩阵:
Figure CN102353340AD00112
^[XJ,Ζ,If 3D点M与它投影图像点之间的关系为
Figure CN102353340AD00113
[0093] 其中,M是世界坐标系中的三维坐标,s是任意标准矢量,(R,t)定义为外部参数, 表示世界坐标系和CCD摄像机坐标系间的旋转与平移关系,CCD摄像机内部参数用A矩阵表不:
Figure CN102353340AD00114
[0095] 其中,(uQ,v0)是坐标系原点,a和β是图像上u和ν坐标轴的标准矢量,Y表示两坐标的垂直度,由于采用平面标靶,所以假设Z = 0,则由表达式(d)得到:
Figure CN102353340AD00115
[0097] 其中,议MHHF是标靶平面上的齐次坐标,抓=[Hi]『为标靶平面上的点投影到图像平面上对应点的齐次坐标,[巧r2 r3]和t分别是C⑶摄像机坐标系相对于世界坐标系的旋转矩阵及平移向量,用-
-个3X3的映射矩阵
[0098] H= Lh1 h2 h3] = AAtr1 r2 t]
[0099] 得到内部参数矩阵A的两个约束条件
[0100] Zz1H1Zz2=O (g)
[0101] ^A-tA'1 H1 = hT2A-TA-\ (h)
[0102]设:
[0103]
Figure CN102353340AD00121
[0105] B是对称矩阵,表示一个六位向量b = [B11, B12,B22, B13,B23, Β33]τ,基于绝对二次曲线原理求出B,并对B求逆,利用Cholesky分解,从B中倒出内部参数矩阵A,由内部参数矩阵A和映射矩阵H计算每幅图像相对于平面标靶的CCD摄像机外部参数矩阵R和平移向量 t :
[0106]
[0107]
[0108] [0109]
Figure CN102353340AD00122
[0110] 将(j)作为初值,ri、r2、r3、t表示世界坐标系和CXD摄像机坐标系间的旋转与平移关系,通过最大似然估计优化求解其参数,目标函数为:
[0111]
Figure CN102353340AD00123
[0112] 其中Hiij为第i幅图像第j个点图像的齐次坐标,m为第i幅图像第j个点映射出图像的齐次坐标,A为内参,Ri和Ti是第i幅图像对于CCD摄像机坐标系的旋转和平移矩阵,Mj为第i幅图像的第j个世界坐标系的三维点,目标函数趋于0,计算出CCD摄像机的准确参数值。径向畸变参数求解,由径向畸变参数模型:
[0113]
Figure CN102353340AD00124
[0114] (L)
[0115]
Figure CN102353340AD00131
[0116] 其中(x,y)为理想的无畸变的图像坐标,(ij)为实际的畸变图像坐标,1^、1ί2为径向畸变参数;
0125]
得如下方程
0126]
0117]
0118]
0119]
0120]
Y表示两坐标的垂直度,故设Y = 0,由:
Figure CN102353340AD00132
0121] 得到
0122]
0123]
0124]
Figure CN102353340AD00133
其中(u,ν)是(k)式求出的内参,计算出的理想投影点的像素坐标,由(U,V),获
Figure CN102353340AD00134
0127] 对于η幅图片的m个点,共有2XnXm个方程,表示为Dk = d
0128] 其中:
0129]
Figure CN102353340AD00135
0132]
0133]
0134]
0135]
k = (DtD) _1DTd 目标函数为:
由此标定出CCD摄像机的内参和径向畸变参数;[0136] 由q式可分别获得到左CXD摄像机的旋转矩阵Ii1和平移向量T1以及右CXD摄像机的旋转矩阵艮和平移向量Ί;
[0137] R = Rr(R1)1 (r)
[0138] T = Tr-RT1
[0139] 由(r)式得到两个(XD摄像机的位置关系R和T。
[0140] (2)特征信息提取
[0141] 采用canny算法进行边缘检测、获取缸盖毛坯图形的轮廓信息,然后采用均值滤波计算排除掉噪声点,识别出关键点线,对经过边缘检测后的缸盖毛坯图像再应用均值滤波计算,滤掉其中的杂点,得到有效的关键点线;
[0142] (3)立体匹配
[0143] 1)消除畸变:采用步骤(1)中的(L)式消除径向方向上的镜头畸变,输出无畸变的图像;
[0144] 2)(XD摄像机校正:在物理意义上将CXD摄像机在数学意义上对准到同一个观察平面上,方法是通过数学的算法来实现使左右两幅图像的y坐标相等,这样左右两幅图像的视差只存在χ轴上,输出行对准的校正图像;
[0145] 3)图像匹配:查找两个C⑶摄像机视场中的相同特征,算出视差。
[0146] (4)三维重建
[0147] 由两个CXD摄像机拍摄分别获取两张图像,由视差计算出缸盖毛坯进气端面轮廓特征图像点的世界坐标值,根据特征信息提取所得到的有效的关键点线,用最小二乘法作曲线拟合,获取缸盖毛坯进气端面特征轮廓的函数方程,将设计的缸盖进气端面加工位置点带入该函数方程,与轮廓特征点进行对比,确定被测缸盖毛坯加工位置以及相应尺寸。
[0148] 本实施例以缸盖毛坯进气端面加工尺寸的测取定位为例。
[0149] 首先对左右两个CXD摄像机分别进行标定,将黑白平面的方块阵列标靶(如图2) 置于CCD摄像机前,摄像机的视野包含整个标靶平面。然后根据靶标角点(标靶方块图案交接点)与其在图像投影之间的映射关系,分别计算出描述左右两个CCD摄像机的内参数、 外参数和径向畸变参数:
r M0
;外参数Lr1 r2 r3 t];径向畸变参数Lk1 k2]。
[0151] 在分别得到左侧摄像机的旋转矩阵Ii1和平移向量T1以及右侧摄像机的旋转矩阵 Rr和平移向量I后,利用(r)式,得到两个C⑶摄像机的位置关系R和Τ:
[0152] R = Rr(R1)1 (r)
[0153] T = Tr-RT1
[0154] 即可完成标定。此时两个C⑶摄像机之间的相对位置不能发生改变。
[0155] 然后采用canny算法进行边缘检测、获取缸盖毛坯图形的轮廓信息;采用均值滤波计算排除掉噪声点,识别出关键点线,对经过边缘检测后的缸盖毛坯图像再应用均值滤波计算,滤掉其中的杂点,得到有效的关键点线(图3所示)。
[0156] 采用步骤⑴中的(L)式消除径向方向上的镜头畸变,输出无畸变的图像;通过数学的算法使两个CCD摄像机将视场中的景物置于同一个观察平面,即:使(摄像机获取的)
[0150]
Figure CN102353340AD00141
左右两幅图像的1坐标相等,这样两幅图像仅在X轴上存在视差,获取输出行对准的校正图像。图3是获取的缸盖进气端面的轮廓特征,由于CCD摄像机沿镜头方向与被测缸盖的进气端面并不垂直,所以图3并不能反映真实的缸盖进气端面的轮廓特征点参数。所以需进行三维重建。根据已获取的视差,可计算出缸盖进气端面的轮廓特征图像点在实际物理空间的世界坐标值。本实施例获取的缸盖进气端面的真实轮廓参数点如图4所示,将图4与标准的无缺陷的缸盖进气端面轮廓特征点参数(如图幻进行对比,由此来判断缸盖毛坯件是否合格。

Claims (2)

1.缸盖毛坯加工尺寸识别方法及装置,具有CCD摄像机(1)、托板O)、支撑架(3)、标靶以及计算机,其特征是由2个或2个以上CCD摄像机构成识别接收器,识别接收器接收标靶图像,CCD摄像机镜头视野包含整个标靶平面,标靶是平面阵列的圆孔图案,或者是黑白平面的方块阵列。
2.缸盖毛坯加工尺寸识别方法,其特征是对缸盖毛坯加工尺寸的测取定位包括如下步骤:(I)CCD摄像机标定获取两个CCD摄像机的内部参数,A =首先对两个CCD摄像机之间的旋转平移矩阵进行标定,2D点表示为m = [u, ν]τ, 3D点表示为M = [X,Y,Z]T,2D和3D的增广矩阵:
Figure CN102353340AC00021
3D点M与它投影图像点之间的关系为:
Figure CN102353340AC00022
其中,M是世界坐标系中的三维坐标,s是任意标准矢量,(R,t)定义为外部参数,表示世界坐标系和CCD摄像机坐标系间的旋转与平移关系,CCD摄像机内部参数用A矩阵表示:
Figure CN102353340AC00023
其中,(u。,v。)是坐标系原点,a和β是图像上u和ν坐标轴的标准矢量,Y表示两坐标的垂直度,由于采用平面标靶,所以假设Z = 0,则由表达式(d)得到:
Figure CN102353340AC00024
其中,是标靶平面上的齐次坐标,/Τ。= [nif为标靶平面上的点投影到图像平面上对应点的齐次坐标,[ri r2 r3]和t分别是C⑶摄像机坐标系相对于世界坐标系的旋转矩阵及平移向量,用一个3X3的映射矩阵: H = Di1 h2 h3] = AAtr1 r2 t] (f)得到内部参数矩阵A的两个约束条件:
Figure CN102353340AC00025
Figure CN102353340AC00031
B是对称矩阵,表示一个六位向量b = [B11, B12,B22, B13,B23,も3]\基于绝对二次曲线原 理求出B,并对B求逆,利用Cholesky分解,从B中倒出内部參数矩阵A,由内部參数矩阵A 和映射矩阵H计算每幅图像相对于平面标靶的CCD摄像机外部參数矩阵R和平移向量t :
Figure CN102353340AC00032
将(j)作为初值,巧、r2、r3、t表示世界坐标系和CXD摄像机坐标系间的旋转与平移关系,通过最大似然估计优化求解其參数,目标函数为:
Figure CN102353340AC00033
其中为第i幅图像第j个点图像的齐次坐标,m为第i幅图像第j个点映射出图像 的齐次坐标,A为内參,Ri和Ti是第i幅图像对于CXD摄像机坐标系的旋转和平移矩阵,Mj 为第i幅图像的第j个世界坐标系的三维点,目标函数趋于0,计算出CXD摄像机的准确參 数值。径向畸变參数求解,由径向畸变參数模型:
Figure CN102353340AC00034
其中(x,y)为理想的无畸变的图像坐标,(X力为实际的畸变图像坐标,I^lc2为径向畸 变參数;Y表示两坐标的垂直度,故设Y = 0,由:
Figure CN102353340AC00035
得到:
Figure CN102353340AC00036
Figure CN102353340AC00041
其中(u,ν)是(k)式求出的内参,计算出的理想投影点的像素坐标,由(u,V),获得如下方程:(O)(U -W0X^ -U0Xx 2+/)厂 'K V u—u(V- -V0Xx2 + /Xv- -V0Xx2 + /)2 _ Λ. V V-V对于η幅图片的m个点,共有2 X η Xm个方程,表示为Dk = d ; 其中:
Figure CN102353340AC00042
目标函数为:由此标定出CCD摄像机的内参和径向畸变参数;由q式可分别获得到左CCD摄像机的旋转矩阵Ii1和平移向量T1以及右CCD摄像机的旋转矩阵艮和平移向量Ί;R = Rr(R1)1 (1)T = Tr-RT1由(r)式得到两个CXD摄像机的位置关系R和T ;(2)特征信息提取采用carmy算法进行边缘检测、获取缸盖毛坯图形的轮廓信息,然后采用均值滤波计算排除掉噪声点,识别出关键点线,对经过边缘检测后的缸盖毛坯图像再应用均值滤波计算,滤掉其中的杂点,得到有效的关键点线;(3)立体匹配1)消除畸变:采用步骤(1)中的(L)式消除径向方向上的镜头畸变,输出无畸变的图像;2) CXD摄像机校正:在物理意义上将CXD摄像机在数学意义上对准到同一个观察平面上,方法是通过数学的算法来实现使左右两幅图像的y坐标相等,这样左右两幅图像的视差只存在χ轴上,输出行对准的校正图像;3)图像匹配:查找两个C⑶摄像机视场中的相同特征,算出视差;(4)三维重建由两个CXD摄像机拍摄分别获取两张图像,由视差计算出缸盖毛坯进气端面轮廓特征图像点的世界坐标值,根据特征信息提取所得到的有效的关键点线,用最小二乘法作曲线拟合,获取缸盖毛坯进气端面特征轮廓的函数方程,将设计的缸盖进气端面加工位置点带入该函数方程,与轮廓特征点进行对比,确定被测缸盖毛坯加工位置以及相应尺寸。
CN 201110151060 2011-06-08 2011-06-08 缸盖毛坯加工尺寸识别方法及装置 Active CN102353340B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201110151060 CN102353340B (zh) 2011-06-08 2011-06-08 缸盖毛坯加工尺寸识别方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201110151060 CN102353340B (zh) 2011-06-08 2011-06-08 缸盖毛坯加工尺寸识别方法及装置

Publications (2)

Publication Number Publication Date
CN102353340A true CN102353340A (zh) 2012-02-15
CN102353340B CN102353340B (zh) 2013-08-28

Family

ID=45576961

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201110151060 Active CN102353340B (zh) 2011-06-08 2011-06-08 缸盖毛坯加工尺寸识别方法及装置

Country Status (1)

Country Link
CN (1) CN102353340B (zh)

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103759638A (zh) * 2014-01-10 2014-04-30 北京力信联合科技有限公司 一种零件检测方法
CN103954213A (zh) * 2014-04-17 2014-07-30 北京力信联合科技有限公司 一种分析零件的实测图的方法
CN104655018A (zh) * 2015-03-16 2015-05-27 武汉大学 一种花键套端面尺寸检测系统及方法
CN105261010A (zh) * 2015-09-18 2016-01-20 北京林业大学 一种不需控制点坐标测量的相机定标方法
CN106918306A (zh) * 2017-04-22 2017-07-04 许晟明 基于光场单相机的工业产品三维形貌实时检测系统
CN107204017A (zh) * 2017-06-08 2017-09-26 爱佩仪中测(成都)精密仪器有限公司 一种三维测量中的单个摄像机标定方法
CN107256569A (zh) * 2017-06-08 2017-10-17 爱佩仪中测(成都)精密仪器有限公司 基于双目视角的三维测量双摄像机标定方法
CN107588721A (zh) * 2017-08-28 2018-01-16 武汉科技大学 一种基于双目视觉的零件多尺寸的测量方法及系统
CN108109179A (zh) * 2017-12-29 2018-06-01 天津科技大学 基于针孔摄像机模型的摄像机姿态校正方法
CN108536151A (zh) * 2018-05-06 2018-09-14 长春北方化工灌装设备股份有限公司 一种视觉导引的闭环执行系统及视觉导引方法
CN109335649A (zh) * 2018-09-26 2019-02-15 重庆巴奥科技有限公司 一种用于智能货物摆放的无规则分布抓取提炼方法
CN111156872A (zh) * 2019-12-31 2020-05-15 东风汽车零部件(集团)有限公司活塞轴瓦分公司 汽车双离合换挡毂三坐标测量夹具及测量方法
CN111862048A (zh) * 2020-07-22 2020-10-30 浙大城市学院 基于关键点检测和深度卷积神经网络的鱼体姿态与长度自动分析方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1566900A (zh) * 2003-06-11 2005-01-19 北京航空航天大学 一种空间圆几何参数的视觉测量方法
WO2006120759A1 (ja) * 2005-05-12 2006-11-16 Techno Dream 21 Co., Ltd. 3次元形状計測方法およびその装置
WO2009056050A1 (fr) * 2007-10-24 2009-05-07 Shenzhen Huawei Communication Technologies Co. , Ltd. Procédé et dispositif d'étalonnage de caméra vidéo
CN101876532A (zh) * 2010-05-25 2010-11-03 大连理工大学 测量系统中的摄像机现场标定方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1566900A (zh) * 2003-06-11 2005-01-19 北京航空航天大学 一种空间圆几何参数的视觉测量方法
WO2006120759A1 (ja) * 2005-05-12 2006-11-16 Techno Dream 21 Co., Ltd. 3次元形状計測方法およびその装置
WO2009056050A1 (fr) * 2007-10-24 2009-05-07 Shenzhen Huawei Communication Technologies Co. , Ltd. Procédé et dispositif d'étalonnage de caméra vidéo
CN101876532A (zh) * 2010-05-25 2010-11-03 大连理工大学 测量系统中的摄像机现场标定方法

Cited By (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103759638B (zh) * 2014-01-10 2019-04-02 北京力信联合科技有限公司 一种零件检测方法
CN103759638A (zh) * 2014-01-10 2014-04-30 北京力信联合科技有限公司 一种零件检测方法
CN103954213A (zh) * 2014-04-17 2014-07-30 北京力信联合科技有限公司 一种分析零件的实测图的方法
CN104655018A (zh) * 2015-03-16 2015-05-27 武汉大学 一种花键套端面尺寸检测系统及方法
CN104655018B (zh) * 2015-03-16 2017-07-25 武汉大学 一种花键套端面尺寸检测系统及方法
CN105261010B (zh) * 2015-09-18 2017-12-15 北京林业大学 一种不需控制点坐标测量的相机定标方法
CN105261010A (zh) * 2015-09-18 2016-01-20 北京林业大学 一种不需控制点坐标测量的相机定标方法
CN106918306A (zh) * 2017-04-22 2017-07-04 许晟明 基于光场单相机的工业产品三维形貌实时检测系统
CN107204017A (zh) * 2017-06-08 2017-09-26 爱佩仪中测(成都)精密仪器有限公司 一种三维测量中的单个摄像机标定方法
CN107256569A (zh) * 2017-06-08 2017-10-17 爱佩仪中测(成都)精密仪器有限公司 基于双目视角的三维测量双摄像机标定方法
CN107588721A (zh) * 2017-08-28 2018-01-16 武汉科技大学 一种基于双目视觉的零件多尺寸的测量方法及系统
CN108109179B (zh) * 2017-12-29 2021-05-18 天津科技大学 基于针孔摄像机模型的摄像机姿态校正方法
CN108109179A (zh) * 2017-12-29 2018-06-01 天津科技大学 基于针孔摄像机模型的摄像机姿态校正方法
CN108536151A (zh) * 2018-05-06 2018-09-14 长春北方化工灌装设备股份有限公司 一种视觉导引的闭环执行系统及视觉导引方法
CN109335649A (zh) * 2018-09-26 2019-02-15 重庆巴奥科技有限公司 一种用于智能货物摆放的无规则分布抓取提炼方法
CN111156872A (zh) * 2019-12-31 2020-05-15 东风汽车零部件(集团)有限公司活塞轴瓦分公司 汽车双离合换挡毂三坐标测量夹具及测量方法
CN111862048A (zh) * 2020-07-22 2020-10-30 浙大城市学院 基于关键点检测和深度卷积神经网络的鱼体姿态与长度自动分析方法
CN111862048B (zh) * 2020-07-22 2021-01-29 浙大城市学院 基于关键点检测和深度卷积神经网络的鱼体姿态与长度自动分析方法

Also Published As

Publication number Publication date
CN102353340B (zh) 2013-08-28

Similar Documents

Publication Publication Date Title
CN102353340B (zh) 缸盖毛坯加工尺寸识别方法及装置
WO2014024579A1 (ja) 光学データ処理装置、光学データ処理システム、光学データ処理方法、および光学データ処理用プログラム
CN104537707B (zh) 像方型立体视觉在线移动实时测量系统
Ortin et al. Indoor robot motion based on monocular images
CN103278138B (zh) 一种复杂结构薄部件三维位置及姿态的测量方法
Chatterjee et al. Algorithms for coplanar camera calibration
CN101354796B (zh) 基于泰勒级数模型的全向立体视觉三维重建方法
CN102005039A (zh) 基于泰勒级数模型的鱼眼相机立体视觉深度测量方法
CN107816942A (zh) 一种基于十字结构光视觉系统的平面二维尺寸测量方法
CN110288656A (zh) 一种基于单目摄像头的目标定位方法
Ye et al. An accurate 3D point cloud registration approach for the turntable-based 3D scanning system
Orteu et al. Camera calibration for 3D reconstruction: application to the measurement of 3D deformations on sheet metal parts
CN106447709A (zh) 一种快速高精度双目视差匹配方法
Su et al. A simple rectification method of stereo image pairs with calibrated cameras
Ann et al. Study on 3D scene reconstruction in robot navigation using stereo vision
Pachidis et al. Pseudo-stereo vision system: a detailed study
CN109579695B (zh) 一种基于异构立体视觉的零件测量方法
Lee et al. CCD camera calibrations and projection error analysis
Cui et al. Epipolar geometry for prism-based single-lens stereovision
Zhang et al. Measuring propeller blade width using binocular stereo vision
Pedersini et al. Calibration and self-calibration of multi-ocular camera systems
Cheng et al. Light pen calibration for a monocular-vision-based coordinate measuring system
Xu et al. A real-time ranging method based on parallel binocular vision
Liu et al. A calibration method of binocular vision system for large forging dimension measurement
Yamazoe et al. Multiple camera calibration with bundled optimization using silhouette geometry constraints

Legal Events

Date Code Title Description
PB01 Publication
C06 Publication
SE01 Entry into force of request for substantive examination
C10 Entry into substantive examination
GR01 Patent grant
C14 Grant of patent or utility model