双目视觉转轴标定方法技术

本发明专利技术属于利用视觉原理进行检测的技术，特别涉及双目视觉转轴标定方法。为提供一种操作方便，标定的过程比较简单的双目视觉转轴标定方法，发明专利技术采用的技术方案是：双目视觉转轴标定方法，包括下列步骤：１）获取球形靶标边缘；２）求取空间球心坐标；３）标定计算转轴参数。本发明专利技术主要用于利用两个摄像机模拟人眼的结构对空间特征点进行测量。

【技术实现步骤摘要】

本专利技术属于利用视觉原理进行检测的技术，特别涉及。
技术介绍
双目测量技术是计算机视觉检测技术的一个重要组成部分，利用两个摄像机模拟人眼的结构对空间特征点进行测量。为了扩大双目测量系统的测量范围，需要在测量系统中加入一个旋转台，这样就可以完成对被测物体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)，得到空间直线l0 其中 分别为摄像机坐标系到测量系统世界坐标系的旋转矩阵和平移矩阵，r1，r2…r9为旋转矩阵的旋转分量，代表摄像机坐标系与测量系统世界系之间每个坐标轴的角度关系，tx，ty，tz分别为坐标两系之间的平移分量，是两个坐标系坐标原点之间的相互位置。它们通过对双目测量系统的标定获取。 分别设置左右两个摄像机的摄像机坐标系，它们与测量系统世界坐标系Ow-XwYwZw之间的相互位置关系通过双目视觉测量系统标定获取，分别对左右摄像机像面上球形靶标边缘进行筛选，它们在空间上对应标定球形靶标上的两个大圆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)可以认为是转轴与球心轨迹平面的交点，因此该点也是转轴上的一点，利用转轴的方向矢量和转轴上的一点就可以将旋转台转轴的方位在空间中表示出来，即完成了旋转台转轴的标定。 所述的通过灰度值跳变像素的寻找粗略搜寻出靶标的边缘是指，首先从搜寻出边缘的最高点开始分别沿搜寻出的边缘向两边进行边缘剔出。 本专利技术具备以下技术效果 本专利技术主要设计了在双目视觉测量系统中利用一个未知半径的球形靶标进行旋转台转轴标定的计算方法。利用在若干组旋转台旋转的不同位置上拍摄的标定球形靶标图像，通过图像处理，空间球心计算和转轴参数计算可以获取转轴在双目测量系统中的位置。这种方法与现有的转轴标定方法比较不需要在每个标定位置上完整测出标定靶标的表面轮廓，因此操作方便，标定的过程比较简单，并且很好地避免了双目视觉测量系统测量物体表面轮廓数据难得问题。而且对标定过程中使用的球形靶标没有半径的精度要求，节约了标定成本。 附图说明 图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对应的圆心坐标Oo根据摄像机小孔成像模型，摄像机像平面上的圆心坐标O在摄像机坐标系下的坐标为(uo，vo，-f)，f为摄像机镜头的焦距，通过摄像机标定获取。连接O与摄像机镜头的光心Oc(0，0，0)，得到空间直线l0 其中 分别为摄像机坐标系到测量系统世界坐标系的旋转矩阵和平移矩阵，r1，r2…r9为旋转矩阵的旋转分量，代表摄像机坐标系与测量系统世界系之间每个坐标轴的角度关系，tx，ty，tz分别为坐标两系之间的平移分量，是两个坐标系坐标原点之间的相互位置。它们通过对双目测量系统的标定获取。通过对双目测量系统的标定获取。根据物体小孔成像的摄像机投影原理该直线一定通过标定球的球心。 通过球形靶标边缘求取空间球心坐标的方法示意图如图4所示。Oc1-Xc1Yc1Zc1和Oc2-Xc2Yc2Zc2分别为左右两个摄像机的摄像机坐标系，它们与测量系统世界坐标系Ow-XwYwZw之间的相互位置关系通过双目视觉测量系统标定获取。lc1和lc2分别为左右摄像机像面上球形靶标边缘筛选结果，它们在空间上对应标定球形靶标上的两个大圆l1和l2。通过方程组1可以分别得到lc1和lc2在各自摄像机坐标系下的圆心坐标(uo1，vo1，-f1)和(uo2，vo2，-f2)，则圆心空间直线l01和l02的交点即为标定球靶标在当前位置上的球心坐标。 3)转轴参数标定计算 在旋转台旋转的若干个位置上分别求取标定球靶标的球心坐标，得到一组围本文档来自技高网...

【技术保护点】
一种双目视觉转轴标定方法，其特征在于，包括下列步骤：１）球形靶标边缘获取：通过灰度值跳变像素的寻找粗略搜寻出靶标的边缘，然后根据边缘上相邻象素点之间的曲率变化对边缘进行筛选，得到准确的球形靶标边缘点；２）空间球心坐标求取：　 　设球形靶标边缘ｌ↓［ｃ］上点的图像坐标为（ｕ↓［ｉ］，ｖ↓［ｉ］），则其对应的圆心坐标（ｕ↓［ｏ］，ｖ↓［ｏ］）应满足：（ｕ↓［ｉ］－ｕ↓［ｏ］）↑［２］＋（ｖ↓［ｉ］－ｖ↓［ｏ］）↑［２］＝ｒ↑［２］ｉ＝１，２，…，ｎ　 　（１）使用最小二乘法求解方程组（１）即可得到ｌ↓［ｃ］对应的圆心坐标Ｏ，根据摄像机小孔成像模型，摄像机像平面上的圆心坐标Ｏ在摄像机坐标系下的坐标为（ｕ↓［ｏ］，ｖ↓［ｏ］，－ｆ），ｆ为摄像机镜头的焦距，通过摄像机标定获取，连接Ｏ与 摄像机镜头的光心Ｏ↓［ｃ］（０，０，０），得到空间直线ｌ↓［０］：＊＊＊（２）其中：＊＊＊分别为摄像机坐标系到测量系统世界坐标系的旋转矩阵和平移矩阵，ｒ↓［１］，ｒ↓［２］…ｒ↓［９］为旋转矩阵的旋转分量 ，代表摄像机坐标系与测量系统世界系之间每个坐标轴的角度关系，ｔ↓［ｘ］，ｔ↓［ｙ］，ｔ↓［ｚ］分别为坐标两系之间的平移分量，是两个坐标系坐标原点之间的相互位置，它们通过对双目测量系统的标定获取；分别设置左右两个摄像机的摄像机坐标系， 它们与测量系统世界坐标系Ｏ↓［ｗ］－Ｘ↓［ｗ］Ｙ↓［ｗ］Ｚ↓［ｗ］之间的相互位置关系通过双目视觉测量系统标定获取，分别对左右摄像机像面上球形靶标边缘进行筛选，它们在空间上对应标定球形靶标上的两个大圆ｌ↓［１］和ｌ↓［２］，通过方程组１可以分别得到ｌ↓［ｃ１］和ｌ↓［ｃ２］在各自摄像机坐标系下的圆心坐标（ｕ↓［ｏ１］，ｖ↓［ｏ１］，－ｆ↓［１］）和（ｕ↓［ｏ２］，ｖ↓［ｏ２］，－ｆ↓［２］），则根据公式（２）得到的圆心空间直线ｌ↓［０１］和ｌ↓［０２］的交点即为标定球靶标在当前位置上的球心坐标；３）转轴参数标定计算：在若干个位置上分别求取标定球靶标的球心坐标，得到一组围绕转轴中心轴线旋转的点Ｏ↓［ｎ］（ｘ↓［ｎ］，ｙ↓［ｎ］，ｚ↓［ｎ］）ｎ＝１，２，…，Ｎ，将这些球心坐标代入空间平面方程，构造线性 方程组：Ａｘ↓［ｎ］＋Ｂｙ↓［ｎ］＋Ｃｚ↓［ｎ］＋Ｄ＝０ｎ＝１，２，…，Ｎ（３）其中Ａ、Ｂ、Ｃ、Ｄ为平面方程参数，利用最小二乘...

【技术特征摘要】

【专利技术属性】
技术研发人员：王鹏，孙长库，
申请(专利权)人：天津大学，
类型：发明
国别省市：12[中国|天津]