CN109883443B - 一种线结构光传感器空间姿态标定方法 - Google Patents

一种线结构光传感器空间姿态标定方法 Download PDF

Info

Publication number
CN109883443B
CN109883443B CN201910122872.8A CN201910122872A CN109883443B CN 109883443 B CN109883443 B CN 109883443B CN 201910122872 A CN201910122872 A CN 201910122872A CN 109883443 B CN109883443 B CN 109883443B
Authority
CN
China
Prior art keywords
ellipse
coordinate system
light sensor
linear structure
geometric
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
Application number
CN201910122872.8A
Other languages
English (en)
Other versions
CN109883443A (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.)
Beijing University of Technology
Original Assignee
Beijing University of Technology
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 Beijing University of Technology filed Critical Beijing University of Technology
Priority to CN201910122872.8A priority Critical patent/CN109883443B/zh
Publication of CN109883443A publication Critical patent/CN109883443A/zh
Application granted granted Critical
Publication of CN109883443B publication Critical patent/CN109883443B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Length Measuring Devices By Optical Means (AREA)

Abstract

本发明公开了一种线结构光传感器空间姿态标定方法,属于精密测量领域。本发明基于线结构光传感器和标准圆柱,提出了一种线结构光传感器空间姿态标定方法。这种方法适用于线结构光测量回转轴类零件,包括圆柱,圆锥,齿轮,轴承等。该方法同时利用了代数距离法最小二乘拟合和几何距离法最小二乘拟合,同时,根据几何关系,在拟合过程中加入了约束条件,最终得到准确的空间姿态标定结果,实现了线结构光对回转类零件的测量。

Description

一种线结构光传感器空间姿态标定方法
技术领域
本发明涉及一种线结构光传感器空间自由姿态标定方法,属于测量测试领域。
背景技术
线结构光传感器是基于激光三角测量原理的一类传感器,可以在测量范围内获取传感器距离激光反射表面的距离。线结构光具有非接触、快速、精度高、结构简单、使用简便的特点,且使用人员无需掌握专业的激光知识。线结构光传感器可以一次读取被测物表面的轮廓。一般情况下,线结构光传感器输出数据为二维信息元素,结合其他定位信息,可得出被测物三维几何信息。
在测量物体表面几何尺寸前,必须标定线结构光传感器与被测物测量所在坐标系的空间位置关系,标定的方法与被测物的几何形状关系密切。而被测物所在的几何坐标系一般为直角坐标系和极坐标系,所以线结构光传感器的姿态标定也基于这两种主要的坐标系进行,其最终目的为将线结构光传感器测得的距离值转换为被测物坐标系的三维值:笛卡尔坐标值(x,y,z)或者圆柱坐标值(ρ,θ,h)。
在实际测量过程中,需要按照被测物特点选择合适标定方法,本方法是针对测量柱状零件(圆柱、齿轮等)等旋转件且该旋转件放置可旋转的轴系上的传感器姿态标定。
发明内容
为了解决线结构光传感器空间坐标标定问题,本发明专利提出了一种自由空间位置的线结构光传感器标定方法。该方法采用一个固定到回转轴系的高精度环规作为校准件,该环规的圆柱直径与被测物直径相当,该圆柱表面在线结构光的测量测量范围之内。线结构光传感器以一定角度倾斜对准该环规,手动微调线结构光传感器的位置,使线结构光传感器发射光线能够进入接受区,使其能够准确的读取测量值。线结构光的线性激光束在圆柱表面成一段圆弧,该圆弧实际为线结构光传感器输出光平面与圆柱体相交所成椭圆的一段椭圆弧,线结构光传感器读取的数据为传感器距离椭圆弧各点的距离。该椭圆有如下特点:
1.该椭圆所在的平面为激光平面
2.该椭圆中心在圆柱轴线上
3.该椭圆短轴与圆柱半径相等
将以上特点作为约束条件,通过最小二乘法拟合出该椭圆方程,进而计算出椭圆的长轴,长短轴比值为该平面的倾斜率,然后将测点转换到以椭圆中心为中心,长短轴为X,Y轴的坐标系中,最后将椭圆坐标系中的点转换到圆柱坐标系中,由此可得三维坐标,如图1所示。
一种线结构光传感器空间姿态标定方法,该方法的具体实施步骤如下:
S1.建立线结构光传感器坐标系
建立如图2所示的传感器坐标系S0(X0,Y0),根据椭圆定义,建立椭圆表达式,
F(a,x)=a·x=ax2+bxy+cy2+dx+ey+f=0
其中,a=[a b c d e f]T,a即a、b、c、d、e、f为椭圆表达式的参数向量,x=[x2 x yy2 x y 1]T,x、y为方程各项;
根据最小代数距离,利用最小二乘法拟合该椭圆,定义
Figure BDA0001972567370000021
其中D=[x1,x2,x3···xN],D中个元素为各点到椭圆中心的距离向量;
该椭圆短轴为圆柱半径R,根据椭圆公式得该椭圆短半轴为:
Figure BDA0001972567370000022
整理得
Figure BDA0001972567370000023
根据椭圆定义,其a≠c,可设
Figure BDA0001972567370000024
整理上式可得
Figure BDA0001972567370000025
其中
Figure BDA0001972567370000026
C1为系数矩阵;
又根据椭圆定义可得4ac-b2<0,设b2-4ac=1,即
Figure BDA0001972567370000027
其中
Figure BDA0001972567370000031
C2是椭圆参数关系的系数矩阵;
联立关于C1,C2等式可得
Figure BDA0001972567370000032
最后得方程组
Figure BDA0001972567370000033
其中,C=C1+C2,λ为广义特征矩阵的特征值
解此广义特征值问题,得椭圆一般表达式的各参数
a0=[a0 b0 c0 d0 e0 f0]
另,椭圆各几何参数中心坐标(Xc,Yc),长轴A,短轴B,倾斜角θ值为,
Figure BDA0001972567370000034
Figure BDA0001972567370000035
Figure BDA0001972567370000036
Figure BDA0001972567370000037
Figure BDA0001972567370000038
S2.几何距离法迭代法拟合椭圆,提高椭圆拟合精度
根据几何距离法定义,首先求椭圆上与测点距离最小的点的坐标,其距离可用如下等式表示
xi=R(Xi-XC)或者Xi=R-1xi-Xc
其中Xi=[xi,yi]T
Figure BDA0001972567370000041
为坐标旋转算子,且R-1=RT,i为点序列
由椭圆的标准表达式
Figure BDA0001972567370000042
以及给定的测量点(xi,yi)并求导可得
Figure BDA0001972567370000043
联立等式可得方程组:
Figure BDA0001972567370000044
对x,y求偏导可得:
Figure BDA0001972567370000045
则迭代公式为
Figure BDA0001972567370000046
其中Qk为第k个点的偏导值。
设置迭代初始值
Figure BDA0001972567370000047
其中
Figure BDA0001972567370000048
迭代3~5次即可得到一个比较准确的距离最近的点。
设a=(Xc Yc a b α)T将根据S1的代数距离法拟合结果a0=(Xc0 Yc0 A0 B0 α0)T作为几何距离法拟合迭代的初始值。对x求偏导可得,
Figure BDA0001972567370000049
且其雅克比矩阵为
Figure BDA00019725673700000410
对f1,f2求偏导可得
Figure BDA0001972567370000051
根据以上方程可得
Jx′i,a=(B-1Q-1B)|x=x′i
其中
B=(B1 B2 B3 B4 B5)
Figure BDA0001972567370000052
Figure BDA0001972567370000053
Figure BDA0001972567370000054
Figure BDA0001972567370000055
Figure BDA0001972567370000056
建立雅克比迭代矩阵
Figure BDA0001972567370000057
将S1中的计算结果作为初始值带入进行迭代,即可得到椭圆参数的准确拟合结果。
S3.将传感器坐标系转换到椭圆坐标系下
如图3所示,设椭圆坐标系为O1,则两坐标系的关系为
Figure BDA0001972567370000058
即可将传感器坐标系转换到椭圆坐标系。
Figure BDA0001972567370000059
为坐标转换前的值,
Figure BDA00019725673700000510
为坐标转换后的坐标值。
S4.将椭圆坐标系转换到圆柱三维坐标系
设椭圆的倾斜角为β,则
Figure BDA00019725673700000511
其中,a,b为椭圆的长轴和短轴长度;
如图4所示,以椭圆中心为中心,建立垂直于椭圆轴线的圆平面,以此平面为Z轴方向原点所在的平面,以圆柱坐标系的X,Y轴为X,Y轴建立圆柱空间直角坐标系Q2,则椭圆平面坐标系可转换圆柱坐标系的三维数据,转换矩阵如下
Figure BDA0001972567370000061
至此,完成整个标定过程。
Figure BDA0001972567370000062
为三维坐标值。
附图说明
图1标定原理整体图示。
图2传感器测量椭圆图示。
图3传感器坐标系转换为椭圆坐标系。
图4椭圆坐标系转换为圆柱三维坐标系。
具体实施方式
以下结合具体加工实例对本发明进行说明:
适当方式固定线结构光传感器,使其平面激光束从倾斜方向与圆柱相交,调整位置使激光接触区域在线结构光传感器测量有效区域内,测量数据即为传感器坐标系下的二维数据,多次测量取均值,得到二维点集合X0=(xi,yi)。
线结构光传感器位置标定步骤如下:
1.根据S1的方法,根据代数距离法最小拟合方法,得到初步的椭圆参数
a0=[Xc Yc a b α]
2.按照S2中关于几何距离拟合方法,以a0为初始值,迭代计算得出较精确的椭圆参数
a=[Xc Yc a b α]
3.在椭圆参数的基础上,结合该椭圆的几何位置关系,建立椭圆坐标,将测量点转换到椭圆坐标系下。
Figure BDA0001972567370000063
4.将椭圆坐标系中的测点转换到圆柱三维直角坐标系下。
Figure BDA0001972567370000064

Claims (1)

1.一种线结构光传感器空间姿态标定方法,其特征在于:该方法利用线结构光传感器的测量特点,其激光平面与标准圆柱环规曲面相交为一椭圆,通过精确拟合该椭圆得到线结构光传感器与圆柱环规的空间几何关系,将线结构光传感器坐标系的测点转换到椭圆坐标系下,根据椭圆与圆柱的几何关系,将测点转换到圆柱三维直角坐标系;带约束条件的代数距离最小二乘法,用于计算椭圆参数的初始值;
几何距离最小二乘法用于迭代计算精确的椭圆参数值;该线结构光传感器空间姿态标定方法包括如下步骤:
步骤1:固定线结构光传感器,使其平面激光束从倾斜方向与圆柱相交,调整位置使激光接触区域在线结构光传感器测量有效区域内,测量数据即为线结构光传感器坐标系下的二维数据,多次测量取均值,得到二维点集合X0=(xi,yi);
步骤2:根据带约束条件的代数距离最小二乘法,得到初步的椭圆参数;
步骤3:按照几何距离最小二乘法,以椭圆参数向量a0为初始值,迭代计算得出较精确的椭圆参数;
步骤4:在椭圆参数的基础上,结合该椭圆的几何位置关系,建立椭圆坐标系,将测量点转换到椭圆坐标系下;
步骤5:将椭圆坐标系中的测点转换到圆柱三维直角坐标系下;
带约束条件的代数距离最小二乘法利用椭圆短轴固定为圆柱半径的特征,建立约束条件;
Figure FDA0003050994630000011
其中,a=[a0 b0 c0 d0 e0 f0],为平面自由椭圆的参数向量;C1为系数矩阵;
利用椭圆的特殊形式为:
Figure FDA0003050994630000012
其必须满足参数关系4ac-b2<0,设b2-4ac=1,即
Figure FDA0003050994630000013
C2是椭圆参数关系的系数矩阵;
几何距离最小二乘法是利用代数距离法得到的椭圆几何参数为初始值,使用雅克比矩阵迭代出一个精确的椭圆参数结果;
Figure FDA0003050994630000021
其中
Figure FDA0003050994630000022
为迭代的雅克比矩阵算子,
Figure FDA0003050994630000023
为迭代增量,
Figure FDA0003050994630000024
为各点的二维计算值。
CN201910122872.8A 2019-02-19 2019-02-19 一种线结构光传感器空间姿态标定方法 Active CN109883443B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910122872.8A CN109883443B (zh) 2019-02-19 2019-02-19 一种线结构光传感器空间姿态标定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910122872.8A CN109883443B (zh) 2019-02-19 2019-02-19 一种线结构光传感器空间姿态标定方法

Publications (2)

Publication Number Publication Date
CN109883443A CN109883443A (zh) 2019-06-14
CN109883443B true CN109883443B (zh) 2021-06-18

Family

ID=66928586

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910122872.8A Active CN109883443B (zh) 2019-02-19 2019-02-19 一种线结构光传感器空间姿态标定方法

Country Status (1)

Country Link
CN (1) CN109883443B (zh)

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110260822B (zh) * 2019-06-18 2021-01-19 西安交通大学 一种多目结构光系统高精度标定方法
CN111046535B (zh) * 2019-11-25 2022-12-09 暨南大学 一种激光加工热分布计算方法
CN111156948A (zh) * 2019-12-29 2020-05-15 苏州赛腾精密电子股份有限公司 用于3c玻璃面板检测的三维数据坐标校正方法及装置
CN111968187B (zh) * 2020-10-21 2021-01-12 三代光学科技(天津)有限公司 一种环形结构光参数标定方法
CN112465979A (zh) * 2020-11-10 2021-03-09 国网内蒙古东部电力有限公司检修分公司 一种基于剖切面几何特征的圆柱面点云拟合方法
CN112880592B (zh) * 2021-01-20 2022-03-04 湘潭大学 一种基于芯轴的数控转台顶尖的倾斜标定方法
CN113358025B (zh) * 2021-05-21 2022-08-02 北京工业大学 一种线激光传感器空间位姿标定件及标定方法
CN113607053B (zh) * 2021-08-13 2022-06-14 西北工业大学 基于筒体内表面特征点的位姿调整装置、方法及系统
CN113847884A (zh) * 2021-09-18 2021-12-28 武汉光谷卓越科技股份有限公司 基于线扫描的精细化三维测量与建模方法
CN114719784B (zh) * 2022-04-11 2023-08-18 沈阳理工大学 一种圆柱形称重传感器倾斜角度检测装置和方法

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6178649B1 (en) * 1998-01-16 2001-01-30 Chyi-Yiing Wu Optical calibrating apparatus for emitting datum laser line
CN1384351A (zh) * 2001-04-30 2002-12-11 北京航空航天大学 结构光三维双视觉标定点发生方法及装置
WO2013102572A1 (en) * 2012-01-05 2013-07-11 Teknologian Tutkimuskeskus Vtt Arrangement for optical measurements and related method
CN103411553A (zh) * 2013-08-13 2013-11-27 天津大学 多线结构光视觉传感器的快速标定方法
CN104848801A (zh) * 2015-06-05 2015-08-19 北京航空航天大学 一种基于平行双圆柱靶标的线结构光视觉传感器标定方法
TW201621811A (zh) * 2014-12-09 2016-06-16 財團法人工業技術研究院 三維掃描器校正系統及其校正方法
CN205650975U (zh) * 2016-02-23 2016-10-19 南京中建化工设备制造有限公司 基于结构光视觉的非标件自动化焊接加工系统
CN109000572A (zh) * 2018-09-17 2018-12-14 四川大学 一种圆柱直径检测方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5316118B2 (ja) * 2009-03-12 2013-10-16 オムロン株式会社 3次元視覚センサ

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6178649B1 (en) * 1998-01-16 2001-01-30 Chyi-Yiing Wu Optical calibrating apparatus for emitting datum laser line
CN1384351A (zh) * 2001-04-30 2002-12-11 北京航空航天大学 结构光三维双视觉标定点发生方法及装置
WO2013102572A1 (en) * 2012-01-05 2013-07-11 Teknologian Tutkimuskeskus Vtt Arrangement for optical measurements and related method
CN103411553A (zh) * 2013-08-13 2013-11-27 天津大学 多线结构光视觉传感器的快速标定方法
TW201621811A (zh) * 2014-12-09 2016-06-16 財團法人工業技術研究院 三維掃描器校正系統及其校正方法
CN104848801A (zh) * 2015-06-05 2015-08-19 北京航空航天大学 一种基于平行双圆柱靶标的线结构光视觉传感器标定方法
CN205650975U (zh) * 2016-02-23 2016-10-19 南京中建化工设备制造有限公司 基于结构光视觉的非标件自动化焊接加工系统
CN109000572A (zh) * 2018-09-17 2018-12-14 四川大学 一种圆柱直径检测方法

Non-Patent Citations (7)

* Cited by examiner, † Cited by third party
Title
Calibration method for line-structured light vision sensor based on a single ball target;Zhen Liu 等;《Optics and Lasers in Engineering》;20150220;第20-28页 *
Direct least square fitting of ellipse;Andrew 等;《IEEE Transactions on Pattern Analysis and Machine Intelligence》;19991231;第476-480页 *
三维空间中线结构光与相机快速标定方法;陈文亮;《红外技术》;20170131;第39卷(第1期);第81-85页 *
刘思远.基于线结构光视觉的轴径测量算法研究.《中国博士学位论文全文数据库信息科技辑》.2016,(第08期), *
基于线结构光视觉的轴径测量算法研究;刘思远;《中国博士学位论文全文数据库信息科技辑》;20160815(第08期);第1-12、55-68页 *
线结构光标定方法综述;张曦 等;《激光与光电子学进展》;20180228(第02期);第1-11页 *
结构光三维视觉检测中光条椭圆拟合方法与仿真研究;魏振忠 等;《仪器仪表学报》;20030228;第24卷(第1期);第69-74页 *

Also Published As

Publication number Publication date
CN109883443A (zh) 2019-06-14

Similar Documents

Publication Publication Date Title
CN109883443B (zh) 一种线结构光传感器空间姿态标定方法
Chao et al. Calibration of laser beam direction for optical coordinate measuring system
CN107702651B (zh) 基于测距传感器的四轴测量系统坐标系与被测叶片模型坐标系的配准方法
CN109163675B (zh) 一种基于激光跟踪仪检测角摆轴位置精度的方法
CN103697824A (zh) 用于坐标测量机的测头的系统标定方法
Chao et al. Extrinsic calibration of a laser displacement sensor in a non-contact coordinate measuring machine
CN112325773B (zh) 一种激光位移传感器的光束方向矢量和原点位置标定方法
CN113160334B (zh) 一种基于手眼相机的双机器人系统标定方法
CN110806571B (zh) 一种多结构光传感器空间姿态标定件及其标定方法
CN106875439B (zh) 基于三维点云模型的单晶硅棒外形尺寸测量方法
Chen et al. Measuring principle and uncertainty analysis of a large volume measurement network based on the combination of iGPS and portable scanner
CN109884659B (zh) 基于激光追踪仪多站位测量系统的大型精密转台标定方法
CN109323665B (zh) 一种线结构光驱动全息干涉的精密三维测量方法
Ravishankar et al. Automated inspection of aircraft parts using a modified ICP algorithm
Zhou et al. A combined calibration method of a mobile robotic measurement system for large-sized components
US11073383B2 (en) Geometric error measurement method for feed drive system and computer-readable recording medium having program recorded for executing same
CN111562563A (zh) 激光雷达转台标定方法、装置和计算机可读存储介质
CN110211175B (zh) 准直激光器光束空间位姿标定方法
CN114485392B (zh) 基于激光跟踪仪的航天器大尺寸机械基准建立方法和系统
CN112365602B (zh) 一种基于激光点云的路面三维重构误差计算方法
CN111046584A (zh) 星载仪器设备的精密调整方法
CN113916128A (zh) 一种基于光笔式视觉测量系统的提高精度的方法
Pechenin et al. Development of a method of ICP algorithm accuracy improvement during shaped profiles and surfaces control
CN115164793A (zh) 微纳米级复合测量的坐标统一标定器及标定方法
CN114034244B (zh) 一种基于Makima算法获得CMM测量范围内任意点几何误差的方法

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