CN100562707C - 双目视觉转轴标定方法 - Google Patents
双目视觉转轴标定方法 Download PDFInfo
- Publication number
- CN100562707C CN100562707C CNB2008100520547A CN200810052054A CN100562707C CN 100562707 C CN100562707 C CN 100562707C CN B2008100520547 A CNB2008100520547 A CN B2008100520547A CN 200810052054 A CN200810052054 A CN 200810052054A CN 100562707 C CN100562707 C CN 100562707C
- Authority
- CN
- China
- Prior art keywords
- coordinate
- edge
- rotating shaft
- camera
- target
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Image Analysis (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明属于利用视觉原理进行检测的技术,特别涉及双目视觉转轴标定方法。为提供一种操作方便,标定的过程比较简单的双目视觉转轴标定方法,发明采用的技术方案是:双目视觉转轴标定方法,包括下列步骤:1)获取球形靶标边缘;2)求取空间球心坐标;3)标定计算转轴参数。本发明主要用于利用两个摄像机模拟人眼的结构对空间特征点进行测量。
Description
技术领域
本发明属于利用视觉原理进行检测的技术,特别涉及双目视觉转轴标定方法。
技术背景
双目测量技术是计算机视觉检测技术的一个重要组成部分,利用两个摄像机模拟人眼的结构对空间特征点进行测量。为了扩大双目测量系统的测量范围,需要在测量系统中加入一个旋转台,这样就可以完成对被测物体360°的测量。为了保证被测物随转台旋转之后测量的准确性,需要对转台进行标定,找出旋转台的旋转轴线。
现在对于计算机视觉测量中转轴标定的方法主要标准平面或高精度半径已知的标准球来实现。通过被标定的视觉测量设备首先多个旋转位置获取标准平面或标准球的表面三维数据,然后根据三维点数据拟合出平面或球的表面方程,最后根据平面或标准球的几何特征求取出旋转台的旋转轴线的位置。这些方法需要标准平面或高精度半径已知的标准球,标定成本比较高;需要在多个位置获取平面或标准球的表面三维数据,标定效率比较低。针对双目视觉测量系统的特点,其测量被测物体的完整表面三维数据过程复杂,需要设计出一套适用于双目测量系统的转台标定方法。
发明内容
为克服现有技术的不足,本发明的目的在于:提供一种操作方便,标定的过程比较简单的双目视觉转轴标定方法及实施装置,本发明采用的技术方案是:双目视觉转轴标定方法,包括下列步骤:
1)球形靶标边缘获取:通过灰度值跳变像素的寻找粗略搜寻出靶标的边缘,然后根据边缘上相邻象素点之间的曲率变化对边缘进行筛选,得到准确的球形靶标边缘点;
2)空间球心坐标求取:
设球形靶标边缘lc上点的图像坐标为(ui,vi),则其对应的圆心坐标(uo,vo)应满足:
(ui-uo)2+(vi-vo)2=r2 i=1,2,…,n (1)
使用最小二乘法求解方程组(1)即可得到lc对应的圆心坐标O,根据摄像机小孔成像模型,摄像机像平面上的圆心坐标O在摄像机坐标系下的坐标为(uo,vo,-f),f为摄像机镜头的焦距,通过摄像机标定获取,连接O与摄像机镜头的光心Oc(0,0,0),得到空间直线lo:
其中:
分别为摄像机坐标系到测量系统世界坐标系的旋转矩阵和平移矩阵,r2,r2…r9为旋转矩阵的旋转分量,代表摄像机坐标系与测量系统世界系之间每个坐标轴的角度关系,tx,ty,tz分别为坐标两系之间的平移分量,是两个坐标系坐标原点之间的相互位置。它们通过对双目测量系统的标定获取。
分别设置左右两个摄像机的摄像机坐标系,它们与测量系统世界坐标系Ow-XwYwZw之间的相互位置关系通过双目视觉测量系统标定获取,分别对左右摄像机像面上球形靶标边缘进行筛选,它们在空间上对应标定球形靶标上的两个大圆l1和l2,通过方程组1可以分别得到lc1和lc2在各自摄像机坐标系下的圆心坐标(uo1,vo1,-f1)和(uo2,vo2,-f2),则根据公式(2)得到的圆心空间直线lO1和lO2的交点即为标定球靶标在当前位置上的球心坐标;
3)转轴参数标定计算:
在若干个位置上分别求取标定球靶标的球心坐标,得到一组围绕转轴中心轴线旋转的点On(xn,yn,zn)n=1,2,…,N,将这些球心坐标代入空间平面方程,构造线性方程组:
Axn+Byn+Czn+D=0 n=1,2,…,N (3)
其中A、B、C、D为平面方程参数。利用最小二乘法求解这个方程组,可以得到球心旋转轨迹所在的平面方程PB,该平面的法线矢量即可认为是旋转台转轴的方向矢量:
在轨迹平面内PB内利用搜寻法求取与球心坐标距离最小的点,即在约束条件:Ax+By+Cz+D=0下求取目标函数:
的最小值,点OA(xA,yA,zA)可以认为是转轴与球心轨迹平面的交点,因此该点也是转轴上的一点,利用转轴的方向矢量和转轴上的一点就可以将旋转台转轴的方位在空间中表示出来,即完成了旋转台转轴的标定。
所述的通过灰度值跳变像素的寻找粗略搜寻出靶标的边缘是指,首先从搜寻出边缘的最高点开始分别沿搜寻出的边缘向两边进行边缘剔出。
本发明具备以下技术效果:
本发明主要设计了在双目视觉测量系统中利用一个未知半径的球形靶标进行旋转台转轴标定的计算方法。利用在若干组旋转台旋转的不同位置上拍摄的标定球形靶标图像,通过图像处理,空间球心计算和转轴参数计算可以获取转轴在双目测量系统中的位置。这种方法与现有的转轴标定方法比较不需要在每个标定位置上完整测出标定靶标的表面轮廓,因此操作方便,标定的过程比较简单,并且很好地避免了双目视觉测量系统测量物体表面轮廓数据难得问题。而且对标定过程中使用的球形靶标没有半径的精度要求,节约了标定成本。
附图说明
图1转轴标定球靶标。
图2球靶标边缘提取结果。
图3球靶标边缘筛选结果。
图4空间球心坐标求取。
图5转轴参数标定。
图6转轴标定图像。
具体实施方式
下面结合附图对本发明做进一步的说明。
如图1所示是标定转轴所需的球形靶标,靶标使用白色的立柱支起一个黑色的球形靶标,并且为了靶标轮廓提取方便保证背景颜色为白色。
1)球形靶标边缘获取
对于球形靶标图像由于背景灰度值较高,而球形靶标灰度值较低,因此通过灰度值跳变像素的寻找可以粗略搜寻出靶标的边缘,如图2所示。但是由于底座的影响这个边缘并不是靶标的真实边缘,需要对搜寻出的边缘进行筛选。首先从搜寻出边缘的最高点开始分别沿搜寻出的边缘向两边进行边缘剔出,因为该点没有受到底座的影响,可以认为就是球形靶标边缘上的点,然后根据边缘中相邻点的曲率变化剔出两边曲率变化跳变点之间的搜寻边缘,得到可信的球靶标边缘像素,如图3所示。
2)空间球心坐标求取
设球形靶标边缘lc上点的图像坐标为(ui,vi),则其对应的圆心坐标(uo,vo)应满足:
(ui-uo)2+(vi-vo)2=R2 i=1,2,…,n (1)
使用最小二乘法求解方程组1即可得到lc对应的圆心坐标O。根据摄像机小孔成像模型,摄像机像平面上的圆心坐标O在摄像机坐标系下的坐标为(uo,vo,-f),f为摄像机镜头的焦距,通过摄像机标定获取。连接O与摄像机镜头的光心Oc(0,0,0),得到空间直线lo:
其中:
分别为摄像机坐标系到测量系统世界坐标系的旋转矩阵和平移矩阵,r1,r2…r9为旋转矩阵的旋转分量,代表摄像机坐标系与测量系统世界系之间每个坐标轴的角度关系,tx,ty,yz分别为坐标两系之间的平移分量,是两个坐标系坐标原点之间的相互位置。它们通过对双目测量系统的标定获取。通过对双目测量系统的标定获取。根据物体小孔成像的摄像机投影原理该直线一定通过标定球的球心。
通过球形靶标边缘求取空间球心坐标的方法示意图如图4所示。Oc1-Xc1Yc1Zc1和Oc2-Xc2Yc2Zc2分别为左右两个摄像机的摄像机坐标系,它们与测量系统世界坐标系Ow-XwYwZw之间的相互位置关系通过双目视觉测量系统标定获取。lc1和lc2分别为左右摄像机像面上球形靶标边缘筛选结果,它们在空间上对应标定球形靶标上的两个大圆l1和l2。通过方程组1可以分别得到lc1和lc2在各自摄像机坐标系下的圆心坐标(uo1,vo1,-f1)和(uo2,vo2,-f2),则圆心空间直线lO1和lO2的交点即为标定球靶标在当前位置上的球心坐标。
3)转轴参数标定计算
在旋转台旋转的若干个位置上分别求取标定球靶标的球心坐标,得到一组围绕转轴中心轴线旋转的点On(xn,yn,zn)n=1,2,…,N,如图5所示。将这些球心坐标代入空间平面方程,构造线性方程组:
Axn+Byn+Czn+D=0 n=1,2,…,N (3)
其中A、B、C、D为平面方程参数。利用最小二乘法求解这个方程组,可以得到球心旋转轨迹所在的平面方程PB,该平面的法线矢量即可认为是旋转台转轴的方向矢量:
在轨迹平面内PB内利用搜寻法求取与球心坐标距离最小的点,即在约束条件:Ax+By+Cz+D=0下求取目标函数:
的最小值,点OA(xA,yA,zA)可以认为是转轴与球心轨迹平面的交点,因此该点也是转轴上的一点。
利用转轴的方向矢量和转轴上的一点就可以将旋转台转轴的方位在空间中表示出来,即完成了旋转台转轴的标定。
在标定过程中首先在旋转台不同的旋转位置上分别利用左右摄像机拍摄标定球靶标的图像,如图6所示。在每幅图像中根据球靶标与背景之间灰度的变化提取球靶标的轮廓边缘,并利用轮廓边缘相邻像素点之间的曲率变化完成轮廓边缘的筛选,得到可信的球靶标边缘轮廓像素点并拟合出轮廓圆心所在的位置。根据双目视觉测量系统的空间几何关系,求解左右摄像机轮廓圆心投影成像直线的交点,即可得到该位置处标定球靶标球心的空间坐标,结果如表1所示。
表1球心空间坐标
0° | 60° | 120° | 180° | 240° | 300° | |
X<sub>w</sub> | -51.604 | -48.399 | -47.271 | -49.391 | -52.623 | -53.723 |
Y<sub>w</sub> | -11.006 | -11.007 | -11.017 | -11.011 | -11.007 | -10.996 |
Z<sub>w</sub> | 3.814 | 3.218 | 0.132 | -2.373 | -1.788 | 1.316 |
对这6个球心坐标使用平面拟合方法,利用最小二乘求解平面方程组,得到球心轨迹所在的平面方程,则该平面的法线矢量u=[-0.002251 -0.999997 0.001034]Г即为旋转台转轴的方向矢量。通过求解5式的目标函数,得到转轴上的点为OA(-50.502749,-11.007636,0.713734)。
本发明可以替代普通的转轴标定方法,应用到基于双目立体视觉测量原理,需要进行旋转台转轴标定的测量系统中去。
Claims (2)
1、一种双目视觉转轴标定方法,其特征在于,包括下列步骤:
1)球形靶标边缘获取:通过灰度值跳变像素的寻找粗略搜寻出靶标的边缘,然后根据边缘上相邻象素点之间的曲率变化对边缘进行筛选,得到准确的球形靶标边缘点;
2)空间球心坐标求取:
设球形靶标边缘lc上点的图像坐标为(ui,vi),则其对应的圆心坐标(uo,vo)应满足:
(ui-uo)2+(vi-vo)2=r2 i=1,2,…,n (1)
使用最小二乘法求解方程组(1)即可得到lc对应的圆心坐标O,根据摄像机小孔成像模型,摄像机像平面上的圆心坐标O在摄像机坐标系下的坐标为(uo,vo,-f),f为摄像机镜头的焦距,通过摄像机标定获取,连接O与摄像机镜头的光心Oc(0,0,0),得到空间直线l0:
其中:
分别为摄像机坐标系到测量系统世界坐标系的旋转矩阵和平移矩阵,r1,r2…r9为旋转矩阵的旋转分量,代表摄像机坐标系与测量系统世界坐标系之间每个坐标轴的角度关系,tx,ty,tz分别为坐标两系之间的平移分量,是两个坐标系坐标原点之间的相互位置,旋转分量、平移分量通过对双目测量系统的标定获取;
分别设置左右两个摄像机的摄像机坐标系,它们与测量系统世界坐标系Ow-XwYwZw之间的相互位置关系通过双目视觉测量系统标定获取,分别对左右摄像机像面上球形靶标边缘进行筛选,lc1和lc2分别为左右摄像机像面上球形靶标边缘筛选结果,lc1和lc2在空间上对应标定球形靶标上的两个大圆l1和l2,通过方程组(1)可以分别得到lc1和lc2在各自摄像机坐标系下的圆心坐标(uo1,vo1,-f1)和(uo2,vo2,-f2),则根据公式(2)得到的圆心空间直线l01和l02的交点即为标定球形靶标在当前位置上的球心坐标;
3)转轴参数标定计算:
在若干个位置上分别求取标定球形靶标的球心坐标,得到一组围绕转轴中心轴线旋转的点On(xn,yn,zn)n=1,2,…,N,将这些球心坐标代入空间平面方程,构造线性方程组:
Axn+Byn+Czn+D=0 n=1,2,…,N (3)
其中A、B、C、D为平面方程参数,利用最小二乘法求解这个方程组,可以得到球心旋转轨迹所在的平面方程PB,该平面的法线矢量即是旋转台转轴的方向矢量:
在轨迹平面内PB内利用搜寻法求取与球心坐标距离最小的点,即在约束条件:Ax+By+Cz+D=0下求取目标函数:
的最小值,点OA(xA,yA,zA)是转轴与球心轨迹平面的交点,因此该点也是转轴上的一点,利用转轴的方向矢量和转轴上的一点就可以将旋转台转轴的方位在空间中表示出来,即完成了旋转台转轴的标定。
2、根据权利要求1所述的一种双目视觉转轴标定方法,其特征在于,所述的通过灰度值跳变像素的寻找粗略搜寻出靶标的边缘,需要对搜寻出的边缘进行筛选,首先从搜寻出边缘的最高点开始分别沿搜寻出的边缘向两边进行边缘剔出。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2008100520547A CN100562707C (zh) | 2008-01-11 | 2008-01-11 | 双目视觉转轴标定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2008100520547A CN100562707C (zh) | 2008-01-11 | 2008-01-11 | 双目视觉转轴标定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101216296A CN101216296A (zh) | 2008-07-09 |
CN100562707C true CN100562707C (zh) | 2009-11-25 |
Family
ID=39622779
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2008100520547A Expired - Fee Related CN100562707C (zh) | 2008-01-11 | 2008-01-11 | 双目视觉转轴标定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100562707C (zh) |
Families Citing this family (23)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101696876B (zh) * | 2009-10-26 | 2011-05-18 | 宁波大红鹰学院 | 一种vcm磁钢的视觉检测方法 |
CN101750021B (zh) * | 2009-12-04 | 2011-05-11 | 深圳先进技术研究院 | Ct系统中几何参数的标定方法、装置 |
CN101825439B (zh) * | 2010-03-25 | 2011-06-08 | 天津大学 | 基于多相机组合的发动机缸体结合面孔组在线测量方法 |
CN102981011B (zh) * | 2012-12-10 | 2014-10-29 | 西北工业大学 | 一种基于双目视觉的空间目标旋转状态辨识方法 |
CN103162622B (zh) * | 2013-02-28 | 2016-06-29 | 西安交通大学 | 单目视觉系统及其使用的便携式球靶标及其测量方法 |
CN103745458B (zh) * | 2013-12-26 | 2015-07-29 | 华中科技大学 | 一种鲁棒的基于双目光流的空间目标旋转轴及质心估计方法 |
CN104101299A (zh) * | 2014-08-05 | 2014-10-15 | 吉林大学 | 汽车视觉检测系统的摄像机三维桁架标定靶标 |
CN104515478B (zh) * | 2014-12-11 | 2017-08-01 | 华中科技大学 | 一种高精度的航空发动机叶片自动三维测量方法和系统 |
CN105783710B (zh) * | 2014-12-24 | 2018-09-11 | 北京中电科电子装备有限公司 | 一种位置标定的方法及装置 |
CN107810382B (zh) * | 2015-06-26 | 2020-12-22 | 卡尔蔡司工业测量技术有限公司 | 用于在坐标测量机的测量空间中组装旋转设备的适配器元件 |
CN105606026B (zh) * | 2016-02-16 | 2018-12-07 | 广东工业大学 | 一种球心坐标测量装置及其测量方法 |
CN105787493B (zh) * | 2016-02-23 | 2019-10-18 | 北京九碧木信息技术有限公司 | 一种基于bim智能提取放样特征点的方法 |
CN105758343B (zh) * | 2016-03-15 | 2018-06-26 | 哈尔滨工业大学 | 基于双标准球的c轴回转中心标定的装置及方法 |
CN106500666A (zh) * | 2016-09-12 | 2017-03-15 | 东莞中子科学中心 | 一种摄影测量标志 |
CN106595517B (zh) * | 2016-11-29 | 2019-01-29 | 天津大学 | 投射条纹几何分布特征结构光测量系统标定方法 |
CN107367229B (zh) * | 2017-04-24 | 2020-05-05 | 天津大学 | 自由双目立体视觉转轴参数标定方法 |
CN108507462B (zh) * | 2018-02-05 | 2019-10-18 | 黑龙江科技大学 | 一种全息干涉四轴测量设备旋转轴的标定方法 |
CN108592787B (zh) * | 2018-03-19 | 2019-11-19 | 上海交通大学 | 3d旋转扫描系统的转轴标定方法与系统 |
CN109163657B (zh) * | 2018-06-26 | 2020-04-10 | 浙江大学 | 一种基于双目视觉三维重建的圆形目标位姿检测方法 |
CN109059761B (zh) * | 2018-08-01 | 2020-06-19 | 西安交通大学 | 一种基于eiv模型的手持靶标测头标定方法 |
CN109514554B (zh) * | 2018-11-30 | 2021-06-22 | 天津大学 | 利用机器人末端视觉系统的工具坐标系快速标定方法 |
CN110044349B (zh) * | 2019-05-13 | 2021-03-16 | 郑州轻工业学院 | 一种基于经纬仪的球体定位方法 |
CN112116667B (zh) * | 2020-09-22 | 2023-11-24 | 扬州大学 | 一种发动机表面加工孔直径测量方法 |
-
2008
- 2008-01-11 CN CNB2008100520547A patent/CN100562707C/zh not_active Expired - Fee Related
Non-Patent Citations (4)
Title |
---|
一种基于立体模版的双目视觉传感器现场标定方法. 陈刚,车仁生,叶东,黄庆成.光学 精密工程,第12卷第6期. 2004 |
一种基于立体模版的双目视觉传感器现场标定方法. 陈刚,车仁生,叶东,黄庆成.光学 精密工程,第12卷第6期. 2004 * |
计算机视觉中双目视觉综述. 沈洪宇,柴毅.科技资讯,第34期. 2007 |
计算机视觉中双目视觉综述. 沈洪宇,柴毅.科技资讯,第34期. 2007 * |
Also Published As
Publication number | Publication date |
---|---|
CN101216296A (zh) | 2008-07-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN100562707C (zh) | 双目视觉转轴标定方法 | |
CN103759669B (zh) | 一种大型零件的单目视觉测量方法 | |
CN102927908B (zh) | 机器人手眼系统结构光平面参数标定装置及方法 | |
CN102692214B (zh) | 一种狭窄空间双目视觉测量定位装置及方法 | |
CN103903260B (zh) | 一种摄像机内参数快速标定的靶标方法 | |
CN109238235B (zh) | 单目序列图像实现刚体位姿参数连续性测量方法 | |
CN105698699A (zh) | 一种基于时间转轴约束的双目视觉测量方法 | |
CN104634248B (zh) | 一种双目视觉下的转轴标定方法 | |
CN109242915A (zh) | 基于多面立体靶标的多相机系统标定方法 | |
CN105486289B (zh) | 一种激光摄影测量系统及相机标定方法 | |
CN107886546B (zh) | 利用球像及公共自极三角形标定抛物折反射摄像机的方法 | |
CN103106661B (zh) | 空间二条相交直线线性求解抛物折反射摄像机内参数 | |
CN106197265B (zh) | 一种空间自由飞行模拟器视觉精密定位方法 | |
CN103426168B (zh) | 基于一维标定杆的普通、广角、鱼眼立体摄像机的通用标定方法 | |
CN107680137B (zh) | 一种基于球体投影公切线的多相机标定及参数优化方法 | |
CN103983186B (zh) | 双目视觉系统校正方法及校正设备 | |
CN110503687B (zh) | 一种空中光电测量平台目标定位方法 | |
CN102982551B (zh) | 空间三条不平行直线线性求解抛物折反射摄像机内参数 | |
CN104835144A (zh) | 利用一个球的球心的像及正交性求解摄像机内参数 | |
CN103606147B (zh) | 多台不共视场量测相机的坐标系转换标定方法 | |
CN103106662B (zh) | 空间一条直线线性求解抛物折反射摄像机内参数 | |
CN109141226A (zh) | 单相机多角度的空间点坐标测量方法 | |
CN102136140B (zh) | 一种基于矩形图样的视频图像距离检测方法 | |
CN110470226A (zh) | 一种基于无人机系统的桥梁结构位移测量方法 | |
Gerke | Using horizontal and vertical building structure to constrain indirect sensor orientation |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20091125 Termination date: 20110111 |