CN1912532A - 一种三维视觉测量量值传递方法和平面靶标 - Google Patents
一种三维视觉测量量值传递方法和平面靶标 Download PDFInfo
- Publication number
- CN1912532A CN1912532A CNA2005100901386A CN200510090138A CN1912532A CN 1912532 A CN1912532 A CN 1912532A CN A2005100901386 A CNA2005100901386 A CN A2005100901386A CN 200510090138 A CN200510090138 A CN 200510090138A CN 1912532 A CN1912532 A CN 1912532A
- Authority
- CN
- China
- Prior art keywords
- target
- transit
- coordinate system
- plane
- coordinate
- 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.)
- Pending
Links
Images
Landscapes
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明属于三维视觉测量技术,涉及对三维视觉现场全局测量量值传递方法的改进。本发明方法的操作步骤如下:将视觉传感器先移动到某一位置k,进行三维测量;放置平面靶标,拍摄靶标图像,提取至少4个非共线图像特征点;计算出从靶标坐标系ObXbYbZb到摄像机坐标系OcXcYcZc的量值传递变换矩阵;将视觉传感器移到位置k+1处,用经纬仪瞄准靶标上至少4个且非共线特征点,计算出从靶标坐标系ObXbYbZb到经纬仪坐标系OtXtYtZt的量值传递变换矩阵;计算视觉传感器在第k个位置处时从OcXcYcZc到OtXtYtZt的量值传递变换矩阵;重复操作步骤,将视觉传感器在位置k(k=1,...,N)处测得的任意一点pc ki=(xc,yc,zc)T统一到OtXtYtZt下。本发明方法只需一台经纬仪就实现了三维视觉现场全局测量量值传递,减少了工作量和处理环节,降低了成本。
Description
技术领域
本发明属于三维视觉测量技术,涉及对三维视觉现场全局测量量值传递方法的改进。
背景技术
三维视觉测量技术广泛应用在汽车、飞机、船舶、航天器的外形测量,大型模具测量,无缝钢管等超大构件测量,大型天线安装及变形监测等。单个视觉传感器无法通过一次测量就完成整个大型物体的测量,实际测量中往往采用由多个视觉传感器构成的柔性多视觉检测***,或者通过单个视觉传感器在大型物体周围的流动测量实现。因而,三维视觉测量的现场全局测量量值传递是关键所在,它决定了***的最终测量精度。
目前,三维视觉测量的现场全局量值传递一般依赖于外部三坐标测量设备,如电子经纬仪对***,三坐标测量机或者数码相机等。
电子经纬仪对由两台测角精度很高的电子经纬仪构成,可以实现大范围(几米甚至几十米)的三坐标测量,三坐标测量机的测量范围一般在几米。实际测量中先利用外部三坐标测量***建立好总体测量坐标系,之后把带有一定特征点的三维立体靶放置在测量***测量空间的不同位置,外部三坐标测量***和视觉传感器同时测量三维立体靶上的特征点,通过这些特征点实现整个***的全局量值传递。(K.Feng,M.Luo,Calibration of structured-light sensor for multi-visualinspection system.SPIE Vol.2350,1994;卢荣胜,视觉准直在线测量技术研究.天津大学博士后研究工作报告,2000;J.West,Noncontactmachine vision metrology within a CAD coordinate system.Proceedings of Autofact’2002 Conference,2002)这种方式的成本较高,高精度的三维立体靶加工较困难。
有的时候,在物体表面粘贴具有一定特征的标记点。这种方式除了可用电子经纬仪对和三坐标测量机外,还可以利用数码相机。数码相机单次可视范围较大,通过有限几次移动,就可以测出所有标记点的三维坐标,用于***中视觉传感器的全局量值传递计算。如德国GOM公司研制的ATOS激光三维扫描测量***。这种贴标的方式显然容易出现变形且标有一定的厚度而难以保证精度。另外,利用数码相机的方式数码相机本身的流动统一也损失了精度。
还有一种方式是只利用三维立体块或者粘贴的标记。这种方式适应于单传感器的流动式三维视觉测量。传感器每移动的两个相邻的位置至少有三个公共点,这样,后续测量位置的传感器坐标系通过公共点进行量值传递,都统一到起始位置的传感器坐标系,从而实现全局测量的量值传递。这种方式存在累积误差,靠后的测量位置的误差较大。如德国GOM公司研制的ATOS激光三维扫描测量***。
另外一种方式是单个视觉传感器移动到被测物周围,逐次完成被测物的局部测量,然后通过高精度的运动平台获得传感器移动的空间方位变换关系,完成量值传递。这种方式的成本较高,且受运动平台移动范围的限制,现场使用的灵活性较差。
发明内容
本发明的目的是:针对现有方法存在的不足,提出一种操作简便、测量范围大、不存在累积误差、成本相对较低的三维视觉现场全局测量量值传递方法与装置,以提高三维视觉测量的量值传递效率和精度。同时,提供一种用于三维视觉现场全局测量量值传递方法的平面靶标。
本发明的技术方案是:一种三维视觉现场全局测量量值传递方法,其特征在于,使用一台经纬仪和一个平面靶标实现三维视觉现场全局测量的量值传递,摄像机模型采用针孔模型,经纬仪的光学***视为理想的透视投影模型,设OtXtYtZt为经纬仪坐标系,Ot为经纬仪的光学中心,XtOtZt平面与水平面度盘平行,Xt轴方向由用户定义,Yt轴为通过原点的竖轴,方向竖直向下,Zt轴方向根据右手法则确定,设在Z=f处有一成像平面,并设f=1,则经纬仪坐标系与像平面坐标系之间有如下关系:
其中(u,v)为像平面上点的坐标,(xt,yt,zt)为空间点在经纬仪坐标系下的坐标。将式(1)用齐次坐标表示:
由于经纬仪瞄准目标时其测量数据是目标的水平角α和垂直角β,下式将角度转化为像平面的坐标值:
设全局坐标系与经纬仪坐标系重合,即OtXtYtZt,在整个测量过程中经纬仪不移动,视觉传感器坐标系建立在摄像机坐标系下,即OcXcYcZc,靶标坐标系ObXbYbZb的XbObYb平面与靶标平面重合,Zb轴方向符合右手法则,量值传递操作步骤如下:
1、求得摄像机内部参数,固定被测物,安装经纬仪,使经纬仪能看到被测物全貌;
2、将视觉传感器连接到计算机,打开传感器电源,启动计算机,运行测量***软件;
3、视觉传感器先移动到某一位置k,对被测物局部区域进行三维测量;
4、在被测区域附件放置平面靶标,拍摄靶标图像,保存在计算机中,采用Harris角点探测方法提取至少4个非共线图像特征点;
5、根据透视几何,利用靶标上特征点坐标及对应的像点坐标,计算出从靶标平面到摄像机像平面的3×3单射矩阵H,由于摄像机内部参数已知,则内部参数矩阵A已知,由H和A可以计算出从靶标坐标系ObXbYbZb到摄像机坐标系OcXcYcZc的量值传递变换矩阵:
其中,Rk cb和tk cb分别为传感器和靶标在第k个位置处从ObXbYbZb到OcXcYcZc旋转矩阵和平移矢量;
6、将视觉传感器移到位置k+1处,靶标不动,用经纬仪瞄准靶标上至少4个且非共线特征点,记下水平角α和垂直角β,根据式(3)计算出靶标特征点对应的经纬仪像点坐标,用与5同样的方法,计算出从靶标坐标系ObXbYbZb到经纬仪坐标系OtXtYtZt的量值传递变换矩阵:
其中,Rk tb和tk tb分别为靶标在第k个位置处从ObXbYbZb到OtXtYtZt旋转矩阵和平移矢量;
7、计算视觉传感器在第k个位置处时从OcXcYcZc到OtXtYtZt的量值传递变换矩阵为
其中
8、重复步骤3到7,视觉传感器在位置k(k=1,...,N)处测得的任意一点
由式(6)统一到OtXtYtZt下,即:
一种用于三维视觉现场全局测量量值传递方法的平面靶标,其特征在于,在靶标的一个平面上有正方形的黑白相间的方格,每个方格的边长为8~10mm,精度为0.003~0.006mm,方格的行数为8~10,方格的列数为8~10。
本发明的优点是:
1、与传统方法采用两台电子经纬仪实现现场全局三维测量量值传递的方法相比,该方法只需一台经纬仪就实现了三维视觉现场全局测量量值传递,减少了工作量和处理环节,降低了成本;
2、与依靠公共区域特征点逐步进行量值传递方法相比,该方法消除了累积误差;
3、不需运动平台,测量范围大,传感器移动自由、方便,特别适合大型物体的三维视觉测量的现场全局量值传递。
附图说明
图1是本发明平面靶标结构示意图。
图2是经纬仪模型的示意图。
图3是本发明的量值传递示意图。
图4是实施例中传感器在位置1和位置2处的测量结果显示示意图。坐标系为传感器坐标系。
图5是实施例中量值传递结果显示示意图。坐标系为经纬仪坐标系。
具体实施方式
下面对本发明做进一步详细说明。本发明首次利用一台经纬仪和一个平面靶标,实现了三维视觉现场全局测量的量值传递。本发明设计如图1所示的平面靶标。黑白方格角点为特征点,黑白方块是边长为9mm的正方形,加工精度为0.005mm。每行均有黑白方块共11个,周围靠近靶标边缘的角点在实际测量中并不使用,这样,可以使用中间100个角点作为特征点来优化求解量值传递矩阵。
本发明三维视觉现场全局测量量值传递方法,其特征在于,使用一台经纬仪和一个平面靶标实现三维视觉现场全局测量的量值传递,摄像机模型采用针孔模型,参见马颂德,张正友。计算机视觉计算理论与算法基础,科学出版社。经纬仪的光学***视为理想的透视投影模型,如图2所示,设OtXtYtZt为经纬仪坐标系,Ot为经纬仪的光学中心,XtOtZt平面与水平面度盘平行,Xt轴方向由用户定义,Yt轴为通过原点的竖轴,方向竖直向下,Zt轴方向根据右手法则确定,设在Z=f处有一成像平面,并设f=1,则经纬仪坐标系与像平面坐标系之间有如下关系:
其中(u,v)为像平面上点的坐标,(xt,yt,zt)为空间点在经纬仪坐标系下的坐标。将式(1)用齐次坐标表示:
由于经纬仪瞄准目标时其测量数据是目标的水平角α和垂直角β,下式将角度转化为像平面的坐标值:
可以看出,经纬仪模型与摄像机模型完全一致,内部参数矩阵为一单位阵。单个经纬仪实际上只是一个二维测量设备。
参见图3,这是基于单台经纬仪和平面靶标的量值传递过程示意图。设全局坐标系与经纬仪坐标系重合,即OtXtYtZt,在整个测量过程中经纬仪不移动,视觉传感器坐标系建立在摄像机坐标系下,即OcXcYcZc,靶标坐标系ObXbYbZb的XbObYb平面与靶标平面重合,Zb轴方向符合右手法则,量值传递操作步骤如下:
1、求得摄像机内部参数,固定被测物,安装经纬仪,使经纬仪能看到被测物全貌;
2、将视觉传感器连接到计算机,打开传感器电源,启动计算机,运行测量***软件;
3、视觉传感器先移动到某一位置k,对被测物局部区域进行三维测量;
4、在被测区域附件放置平面靶标,拍摄靶标图像,保存在计算机中,采用Harris角点探测方法(参见C.Harris,M.Stephens.A combined andedge detector.In Proceedings Alvey Conference,PP189-192,1988)提取至少4个非共线图像特征点;
5、根据透视几何,利用靶标上特征点坐标(在靶标坐标系下)及对应的像点坐标(在图像坐标系下),计算出从靶标平面到摄像机像平面的3×3单射矩阵H,由于摄像机内部参数已知,则内部参数矩阵A已知,由H和A可以计算出从靶标坐标系ObXbYbZb到摄像机坐标系OcXcYcZc的量值传递变换矩阵(参见Z.Zhang,A Flexible New Technique for CameraCalibration,Technical Report,MSR-TR-98-71 Microsoft Research,December 1998):
其中,Rk cb和tk cb分别为传感器和靶标在第k个位置处从ObXbYbZb到OcXcYcZc旋转矩阵和平移矢量;
6、将视觉传感器移到位置k+1处,靶标不动,用经纬仪瞄准靶标上至少4个且非共线特征点,记下水平角α和垂直角β,根据式(3)计算出靶标特征点对应的经纬仪像点坐标,用与5同样的方法,计算出从靶标坐标系ObXbYbZb到经纬仪坐标系OtXtYtZt的量值传递变换矩阵:
其中,Rk tb和tk tb分别为靶标在第k个位置处从ObXbYbZb到OtXtYtZt旋转矩阵和平移矢量;
7、计算视觉传感器在第k个位置处时从OcXcYcZc到OtXtYtZt的量值传递变换矩阵为
其中
8、重复步骤3到7,视觉传感器在位置k(k=1,...,N)处测得的任意一点
可由式(6)统一到OtXtYtZt下,即:
实施例
用本发明对一平面物体在两个位置进行了测量。求得摄像机内部参数为
归一化焦距: fx=2008.169767
fy=2005.994692
偏斜因子: Sx=-0.162852
主点坐标: u0=400.782364
v0=289.420154
径向畸变系数:k1=0.178490
k2=8.055620
传感器在位置1处靶标到摄像机的坐标变换矩阵为
靶标到经纬仪的坐标变换矩阵为
则可计算出传感器在位置1处从摄像机到经纬仪的坐标变换矩阵为
传感器在位置2处靶标到摄像机的坐标变换矩阵为
靶标到经纬仪的坐标变换矩阵为
同样也可计算出传感器在位置2时从摄像机到经纬仪的坐标变换矩阵为
传感器在位置1和位置2处的测量结果显示在传感器坐标系下,如图4所示,量值传递结果显示在经纬仪坐标系下,如图5所示,图中圆点表示传感器在位置1处的测量结果,十字表示在位置2处的测量结果。从图5可以看出,传感器在两个位置的测量结果经过量值传递,实现了三维视觉测量的全局测量。
Claims (2)
1、一种三维视觉现场全局测量量值传递方法,其特征在于,使用一台经纬仪和一个平面靶标实现三维视觉现场全局测量的量值传递,摄像机模型采用针孔模型,经纬仪的光学***视为理想的透视投影模型,设OtXtYtZt为经纬仪坐标系,Ot为经纬仪的光学中心,XtOtZt平面与水平面度盘平行,Xt轴方向由用户定义,Yt轴为通过原点的竖轴,方向竖直向下,Zt轴方向根据右手法则确定,设在Z=f处有一成像平面,并设f=1,则经纬仪坐标系与像平面坐标系之间有如下关系:
其中(u,v)为像平面上点的坐标,(xt,yt,zt)为空间点在经纬仪坐标系下的坐标。将式(1)用齐次坐标表示:
由于经纬仪瞄准目标时其测量数据是目标的水平角α和垂直角β,下式将角度转化为像平面的坐标值:
设全局坐标系与经纬仪坐标系重合,即OtXtYtZt,在整个测量过程中经纬仪不移动,视觉传感器坐标系建立在摄像机坐标系下,即OcXcYcZc,靶标坐标系ObXbYbZb的XbObYb平面与靶标平面重合,Zb轴方向符合右手法则,量值传递操作步骤如下:
1.1、求得摄像机内部参数,固定被测物,安装经纬仪,使经纬仪能看到被测物全貌;
1.2、将视觉传感器连接到计算机,打开传感器电源,启动计算机,运行测量***软件;
1.3、视觉传感器先移动到某一位置k,对被测物局部区域进行三维测量;
1.4、在被测区域附件放置平面靶标,拍摄靶标图像,保存在计算机中,采用Harris角点探测方法提取至少4个非共线图像特征点;
1.5、根据透视几何,利用靶标上特征点坐标及对应的像点坐标,计算出从靶标平面到摄像机像平面的3×3单射矩阵H,由于摄像机内部参数已知,则内部参数矩阵A已知,由H和A可以计算出从靶标坐标系ObXbYbZb到摄像机坐标系OcXcYcZc的量值传递变换矩阵:
其中,Rk cb和tk cb分别为传感器和靶标在第k个位置处从ObXbYbZb到OcXcYcZc旋转矩阵和平移矢量;
1.6、将视觉传感器移到位置k+1处,靶标不动,用经纬仪瞄准靶标上至少4个且非共线特征点,记下水平角α和垂直角β,根据式(3)计算出靶标特征点对应的经纬仪像点坐标,用与1.5同样的方法,计算出从靶标坐标系ObXbYbZb到经纬仪坐标系OtXtYtZt的量值传递变换矩阵:
其中,Rk tb和tk tb分别为靶标在第k个位置处从ObXbYbZb到OtXtYtZt旋转矩阵和平移矢量;
1.7、计算视觉传感器在第k个位置处时从OcXcYcZc到OtXtYtZt的量值传递变换矩阵为
其中
1.8、重复步骤1.3到1.7,视觉传感器在位置k(k=1,…,N)处测得的任意一点
由式(6)统一到OtXtYtZt下,即:
2、一种用于三维视觉现场全局测量量值传递方法的平面靶标,其特征在于,在靶标的一个平面上有正方形的黑白相间的方格,每个方格的边长为8~10mm,精度为0.003~0.006mm,方格的行数为8~10,方格的列数为8~10。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNA2005100901386A CN1912532A (zh) | 2005-08-11 | 2005-08-11 | 一种三维视觉测量量值传递方法和平面靶标 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNA2005100901386A CN1912532A (zh) | 2005-08-11 | 2005-08-11 | 一种三维视觉测量量值传递方法和平面靶标 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN1912532A true CN1912532A (zh) | 2007-02-14 |
Family
ID=37721533
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNA2005100901386A Pending CN1912532A (zh) | 2005-08-11 | 2005-08-11 | 一种三维视觉测量量值传递方法和平面靶标 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN1912532A (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101968347A (zh) * | 2010-08-31 | 2011-02-09 | 苏州大学 | 一种扁平粒状物表面平整度的提取方法 |
CN101614539B (zh) * | 2008-06-27 | 2011-05-04 | 北京航天计量测试技术研究所 | 分离式方位角垂直传递装置 |
CN102269587A (zh) * | 2010-11-02 | 2011-12-07 | 中国海洋大学 | 基于可控光平面的水下三维重绘装置及重绘方法 |
CN103759669A (zh) * | 2014-01-03 | 2014-04-30 | 河南科技大学 | 一种大型零件的单目视觉测量方法 |
CN105631853A (zh) * | 2015-11-06 | 2016-06-01 | 湖北工业大学 | 车载双目相机标定及参数验证方法 |
CN109945824A (zh) * | 2017-12-21 | 2019-06-28 | 中国科学院长春光学精密机械与物理研究所 | 一种载机航向基准测量与传递的方法和*** |
CN111637851A (zh) * | 2020-05-15 | 2020-09-08 | 哈尔滨工程大学 | 一种基于Aruco码的平面旋转角度视觉测量方法及装置 |
-
2005
- 2005-08-11 CN CNA2005100901386A patent/CN1912532A/zh active Pending
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101614539B (zh) * | 2008-06-27 | 2011-05-04 | 北京航天计量测试技术研究所 | 分离式方位角垂直传递装置 |
CN101968347A (zh) * | 2010-08-31 | 2011-02-09 | 苏州大学 | 一种扁平粒状物表面平整度的提取方法 |
CN102269587A (zh) * | 2010-11-02 | 2011-12-07 | 中国海洋大学 | 基于可控光平面的水下三维重绘装置及重绘方法 |
CN102269587B (zh) * | 2010-11-02 | 2012-10-31 | 中国海洋大学 | 基于可控光平面的水下三维重绘装置的重绘方法 |
CN103759669A (zh) * | 2014-01-03 | 2014-04-30 | 河南科技大学 | 一种大型零件的单目视觉测量方法 |
CN103759669B (zh) * | 2014-01-03 | 2016-11-23 | 河南科技大学 | 一种大型零件的单目视觉测量方法 |
CN105631853A (zh) * | 2015-11-06 | 2016-06-01 | 湖北工业大学 | 车载双目相机标定及参数验证方法 |
CN105631853B (zh) * | 2015-11-06 | 2018-01-30 | 湖北工业大学 | 车载双目相机标定及参数验证方法 |
CN109945824A (zh) * | 2017-12-21 | 2019-06-28 | 中国科学院长春光学精密机械与物理研究所 | 一种载机航向基准测量与传递的方法和*** |
CN109945824B (zh) * | 2017-12-21 | 2022-03-22 | 中国科学院长春光学精密机械与物理研究所 | 一种载机航向基准测量与传递的方法和*** |
CN111637851A (zh) * | 2020-05-15 | 2020-09-08 | 哈尔滨工程大学 | 一种基于Aruco码的平面旋转角度视觉测量方法及装置 |
CN111637851B (zh) * | 2020-05-15 | 2021-11-05 | 哈尔滨工程大学 | 一种基于Aruco码的平面旋转角度视觉测量方法及装置 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108844459B (zh) | 一种叶片数字化样板检测***的标定方法及装置 | |
CN1259542C (zh) | 一种空间圆几何参数的视觉测量方法 | |
CN1250942C (zh) | 一种基于平面靶标的结构光视觉传感器标定方法 | |
CN102607457B (zh) | 基于惯性导航技术的大尺寸三维形貌测量装置及方法 | |
CN1912532A (zh) | 一种三维视觉测量量值传递方法和平面靶标 | |
CN1975324A (zh) | 一种双传感器激光视觉三维测量***校准方法 | |
CN105043259B (zh) | 基于双目视觉的数控机床旋转轴误差检测方法 | |
CN104315995B (zh) | 基于虚拟多立方体标准目标的tof深度相机三维坐标标定装置和方法 | |
CN107971831B (zh) | 数控机床高动态大范围任意轮廓误差单目六维测量方法 | |
CN103759670B (zh) | 一种基于数字近景摄影的物体三维信息获取方法 | |
CN201199234Y (zh) | 分层摄影检测*** | |
CN1259543C (zh) | 轮胎多几何参数的激光视觉在线自动测量方法 | |
CN1888814A (zh) | 三维主动视觉传感器的多视点姿态估计和自标定方法 | |
CN106767913B (zh) | 一种基于单led发光点和二维转台的复眼***标定装置及标定方法 | |
KR20150129701A (ko) | 비파괴 검사를 위한 자립식 홀로노믹 추적 방법 및 장치 | |
CN105252341A (zh) | 五轴数控机床动态误差视觉测量方法 | |
CN104316083B (zh) | 一种虚拟多球体球心定位的tof深度相机三维坐标标定装置和方法 | |
CN108596979A (zh) | 一种用于激光雷达和深度相机的标定装置和方法 | |
CN111811395A (zh) | 基于单目视觉的平面位姿动态测量方法 | |
CN105115560A (zh) | 一种船舱舱容的非接触测量方法 | |
CN103162623A (zh) | 一种双摄像机垂直安装的立体测量***及标定方法 | |
CN105800464A (zh) | 一种基于自动吊钩***的定位方法 | |
CN110763136A (zh) | 一种高精度三维工件尺寸检测***及方法 | |
JP2017524122A (ja) | モバイルプラットフォームの変位を計測する方法及び装置 | |
CN113702384A (zh) | 一种回转构件表面缺陷检测装置、检测方法及标定方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
AD01 | Patent right deemed abandoned | ||
C20 | Patent right or utility model deemed to be abandoned or is abandoned |