CN113028990A - 一种基于加权最小二乘的激光跟踪姿态测量系统及方法 - Google Patents
一种基于加权最小二乘的激光跟踪姿态测量系统及方法 Download PDFInfo
- Publication number
- CN113028990A CN113028990A CN202110250220.XA CN202110250220A CN113028990A CN 113028990 A CN113028990 A CN 113028990A CN 202110250220 A CN202110250220 A CN 202110250220A CN 113028990 A CN113028990 A CN 113028990A
- Authority
- CN
- China
- Prior art keywords
- coordinate system
- target
- attitude angle
- measurement
- acquiring
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/002—Measuring arrangements characterised by the use of optical techniques for measuring two or more coordinates
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
- G01S5/16—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using electromagnetic waves other than radio waves
- G01S5/163—Determination of attitude
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明提供了一种基于加权最小二乘的激光跟踪姿态测量系统及方法,包括:激光跟踪测量单元、单目视觉测量单元以及靶标,靶标包括角锥棱镜及光电位置传感器;激光跟踪测量单元和单目视觉测量单元分别固定在预设位置,靶标固定放置在被测物上;综合利用角锥棱镜与光电位置传感器测量方位角与俯仰角时的高分辨率、单目视觉测量横滚角时的高灵敏度,通过加权最小二乘,基于激光光束向量在不同坐标系之间的转换关系完成姿态角的测量。该测量系统及方法提高了大尺寸范围内激光跟踪姿态测量的整体精度,克服了单目姿态测量精度随距离增加线性递减的问题,同时改善了利用纵向投影比进行姿态测量时近距离精度较低的问题。
Description
技术领域
本发明实施例涉及精密测量技术领域,更具体地,涉及一种基于 加权最小二乘融合的激光跟踪姿态测量系统及方法。
背景技术
随着我国国民经济和制造业的不断发展,诸如大型装备建造、航 空航天、汽车等领域,对于飞行器对接、汽车装配、机器人跟踪测量 等方面的姿态测量要求越来越高,大尺寸范围内高精度姿态测量技术 已经成为这些领域重要研究方向之一。
国内外都有对于姿态测量方面的研究,目前国外对于姿态测量方 法的研究较为成熟,已有产品进行商用的公司主要有三家:Leica、 FARO、API。Leica公司开发了T-Probe、T-MAC、T-Scan等一系列姿 态测量装置,采用驻机定位技术,配合激光跟踪仪、高精度相机实现 姿态角的测量;FARO公司提出了使用6Probe探头的方案,该探头上 布置有11个LED发光点,通过激光跟踪仪转台的双目系统进行跟踪 测量,可以覆盖隐藏区域和小型特征,即使在障碍物后边也能跟随移 动的目标,并可以在目标稳定时锁定它;API公司设计了STS特殊反 射器,配有2个编码器和1个水平传感器,根据编码器与水平传感器 确定位置和方向信息;国内华中科技大学、上海交通大学、天津大学 等基于盾构机施工中的导向系统,对激光跟踪姿态测量系统展开了研 究,通过全站仪、CCD相机、角锥棱镜、倾角仪等实现了姿态测量; 湖北工业大学的张慧娟提出了一种单目视觉自动姿态测量方法,将 EPNP与SoftPOSIT相结合,充分利用两种算法的特性实现姿态测量; 湖北工业大学的闫坤将光电位置传感器与单目视觉进行组合,基于特 定靶点的纵向投影比获取合作目标的滚动角度,通过坐标系之间的转 换关系获取姿态角。
现有技术中以激光跟踪设备和单目视觉进行姿态测量时,测量精 度随着距离的增加呈现线性递减的趋势;利用纵向投影比进行姿态测 量时,测量精度不随距离的变化而变化,但近距离测量精度较低。
发明内容
本发明实施例提供了一种克服上述问题或者至少部分地解决上述 问题的基于加权最小二乘融合的激光跟踪姿态测量系统及方法。
一方面,本发明实施例提供了一种基于加权最小二乘融合的激光 跟踪姿态测量系统,其特征在于,包括:激光跟踪测量单元、单目视 觉测量单元以及靶标;其中,
所述激光跟踪测量单元和所述单目视觉测量单元分别固定在预设 位置,所述靶标固定放置在被测物上;所述靶标包括角锥棱镜、光电 位置传感器以及多个靶点;所述角锥棱镜在顶点处开孔,所述光电位 置传感器的感光面与角锥棱镜开孔面平行,所述角锥棱镜的顶点与光 电位置传感器之间预留一定间距,所述激光跟踪测量单元所发出激光 通过所述角锥棱镜开孔打在所述光电位置传感器的感光面上形成光 斑,所述多个靶点的数量不少于4个,且所述不少于4个的靶点非共 面;
所述激光跟踪测量单元用于获取所述角锥棱镜顶点在第一坐标系 中的坐标;其中,所述第一坐标系为所述激光跟踪测量单元的自身坐 标系;
所述单目视觉测量单元用于通过所述多个靶点获取所述靶标在第 二坐标系中的第一姿态角;其中,所述第二坐标系为所述单目视觉测 量单元的自身坐标系;
所述光电位置传感器用于获取所述光斑在第三坐标系中的坐标; 其中,所述第三坐标系为所述光电位置传感器的自身坐标系;
所述第三坐标系下的坐标可结合所述角锥棱镜的顶点与光电位置 传感器之间间距转换为第四坐标系下的坐标;其中,所述第四坐标系 的X、Y轴与所述第二坐标系X、Y轴平行,Z轴方向指向所述光电位 置传感器感光面中心,原点位于角锥棱镜顶点;
根据所述角锥棱镜顶点在所述第一坐标系中的坐标,获取所述激 光束在所述第一坐标系中对应的第一空间向量;根据所述光斑在所述 第四坐标系中的坐标,获取所述激光束在所述第四坐标系中对应的第 二空间向量;根据所述第一空间向量、所述第二空间向量以及所述靶 标在第二坐标系中的第一姿态角,获取所述靶标在所述第一坐标系中 的第二姿态角,进而获得所述被测物的姿态角。
进一步地,所述激光跟踪测量单元为全站仪或激光跟踪仪,所述 单目视觉测量单元为工业相机。
进一步地,所述不少于4个的非共面靶点由反光标志或红外LED 制成。
另一方面本发明实例提供了一种基于加权最小二乘的激光跟踪姿 态测量的方法,其特征在于,包括:
将所述靶标固定放置在被测物上,将所述激光跟踪测量单元与所 述单目视觉测量单元分别固定在预设位置上,通过标定获取第一坐标 系与第二坐标系之间的转换矩阵;
利用所述光电位置传感器获取所属光斑在所述第三坐标系下的坐 标,利用所述单目视觉测量单元通过所述不少于4个的非共面靶点获 取所述靶标在所述第二坐标系下的第一姿态角,利用所述激光跟踪测 量单元获取所述角锥棱镜的顶点在第一坐标系下的坐标;
根据所述光斑在所述第三坐标系下的坐标和所述角锥棱镜的顶点 与光电位置传感器之间间距,获取所述光斑在所述第四坐标系下的坐 标;
根据所述角锥棱镜的顶点在所述第一坐标系下的坐标,获取所述 激光在所述第一坐标系中对应的第一空间向量,根据所述光斑在所述 第四坐标系下的坐标,获取所述激光在所述第四坐标系中对应的第二 空间向量;
根据所述靶标在所述第二坐标系下的第一姿态角、所述第一空间 向量、所述第二空间向量以及所述第一坐标系与第二坐标系之间的转 换矩阵,获取所述靶标在所述第一坐标系下的第二姿态角;
将所述第一姿态角、所述第二姿态角转换到同一坐标系下,通过 加权最小二乘融合获得最终姿态角。
进一步地,利用所述单目视觉测量单元通过所述不少于4个的非 共面靶点获取所述靶标在所述第二坐标系下的第一姿态角,具体包括:
在进行姿态角测量时,首先利用所述单目视觉测量单元获取所述 靶标的像素坐标,通过EPNP算法进行二维点到三维点的匹配,获取 所述第二坐标系与所述第四坐标系之间的转换矩阵与位移向量;
将所述转换矩阵与位移向量作为初始值,通过SoftPOSIT算法完 成后续实时状态下的姿态解算,获取单目视觉测量时姿态角,即所述 靶标在所述第二坐标系下的第一姿态角。
进一步地,根据所述光斑在所述第三坐标系下的坐标和所述角锥 棱镜的顶点与光电位置传感器之间间距,获取所述光斑在所述第四坐 标系下的坐标,具体包括:
将所述光斑在所述第三坐标系下的横坐标与纵坐标,分别作为所 述光斑在第四坐标系下的横坐标与纵坐标,将所述角锥棱镜的顶点与 光电位置传感器之间间距作为所述光斑的Z轴坐标,即所述光斑在所 述第四坐标系下的坐标。
进一步地,根据所述角锥棱镜的顶点在所述第一坐标系下的坐标, 获取所述激光在所述第一坐标系中对应的第一空间向量,根据所述光 斑在所述第四坐标系下的坐标,获取所述激光在所述第四坐标系中对 应的第二空间向量,具体包括:
将所述第一坐标系的原点作为起点,所述角锥棱镜的顶点作为终 点,获取所述第一空间向量,将所述第四坐标系的原点作为起点,所 述光斑作为终点,获取所述第二空间向量。
进一步地,根据所述靶标在所述第二坐标系下的第一姿态角、所 述第一空间向量、所述第二空间向量以及所述第一坐标系与第二坐标 系之间的转换矩阵,获取所述靶标在所述第一坐标系下的第二姿态角, 具体包括:
对所述第一空间向量和所述第二空间向量分别进行单位化处理, 得到第一单位向量和第二单位向量;
根据所述靶标在所述第二坐标系下的第一姿态角、所述第一单位 向量、所述第二单位向量以及所述第一坐标系与第二坐标系之间的转 换矩阵,基于光束向量在不同坐标系下的转换关系,获取靶标在所述 第一坐标系下的第二姿态角;
进一步地,将所述第一姿态角、所述第二姿态角转换到同一坐标 系下,通过加权最小二乘融合获得最终姿态角,具体包括:
将所述第一姿态角通过所述第一坐标系与第二坐标系之间的转换 矩阵转换到第一坐标系下,获取第三姿态角;
假设所述第二姿态角、所述第三姿态角为传感器所测结果,将两个 传感器对于姿态角的观测方程设置为:
Y=Hx+e (7)
式中:x为一维状态量;Y为二维测量向量,设Y=[y1,y2]T;e为二 维测量噪声向量,包含传感器的内部噪声和环境干扰噪声,设 e=[e1,e2]T,若传感器的测量方差为Ri,i=1,2,则E[ei 2]=E[(x-yi)2]=Ri;H 为已知的二维常向量,设H=[1,1]T。
x=x-x (9)
由式(8)和式(9),可得估计方差为:
对式(10)求极小值,取wi的偏导数,并令其等于零得
将(11)代入式(10),可得估计误差为
由式(12)不难看出,采用加权融合的估计方差比任何一个传感 器的测量方差都小。
通过上述过程进行加权最小二乘融合,获取最终姿态角。
本发明实施例提供了一种基于加权最小二乘融合的激光跟踪姿态 测量系统及方法,通过使用加权最小二乘融合的方法,综合利用角锥 棱镜与光电位置传感器测量方位角与俯仰角时的高分辨率,单目视觉 测量横滚角时的高灵敏度,基于激光光束向量在不同坐标系之间的转 换关系完成姿态角的测量。该测量系统及方法提高了大尺寸范围内激 光跟踪姿态测量的整体精度,克服了单目姿态测量时精度随距离增加 线性递减的问题,同时改善了利用纵向投影比进行姿态测量时近距离 精度较低的问题。
附图说明
为了更清楚的说明本发明实施例或现有技术中的技术方案,下面 将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而 易见地,下面描述中的附图是本发明的一些实施例,对于本领域普通 技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图 获得其他的附图。
图1为本发明实施例提供的一种基于加权最小二乘融合的激光跟 踪姿态测量系统示意图;
图2为本发明实施例提供的一种靶标的剖面图;
图3为本发明实施例提供的一种利用上述基于加权最小二乘融合 的激光跟踪姿态测量系统进行姿态角测量的方法的流程图。
图4为本发明实施例提供的以单目视觉自动测量、基于纵向投影 比结合光束向量在不同坐标系下的转换关系解算、EPNP与SoftPOSIT 联合解算并融合光束向量在不同坐标系下的转换关系解算三种方法的 仿真结果对比图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结 合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚地 描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的 实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创 造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
图1为本发明实例提供了一种基于加权最小二乘融合的激光跟踪 姿态测量系统,包括:激光跟踪测量单元1、单目视觉测量单元2以及 靶标3;其中,
所述激光跟踪测量单元1和所述单目视觉测量单元2分别固定在 预设位置,所述靶标3固定放置在被测物上;如图2所示,所述靶标 包括角锥棱镜4、光电位置传感器5以及多个靶点6;所述角锥棱镜4 在顶点处开孔,所述光电位置传感器5的感光面与所述角锥棱镜4开 孔面平行,所述角锥棱镜4的顶点与所述光电位置传感器5之间预留 一定间距,所述激光跟踪测量单元1所发出激光通过所述角锥棱镜4 开孔打在所述光电位置传感器5的感光面上形成光斑,所述多个靶点6 的数量不少于4个,且所述不少于4个的靶点6非共面;
所述激光跟踪测量单元用于获取所述角锥棱镜顶点在第一坐标系 中的坐标;其中,所述第一坐标系为所述激光跟踪测量单元的自身坐 标系;
所述单目视觉测量单元用于通过所述多个靶点获取所述靶标在第 二坐标系中的第一姿态角;其中,所述第二坐标系为所述单目视觉测 量单元的自身坐标系;
所述光电位置传感器用于获取所述光斑在第三坐标系中的坐标; 其中,所述第三坐标系为所述光电位置传感器的自身坐标系;
所述第二坐标系下的坐标可结合所述角锥棱镜的顶点与光电位置 传感器之间间距转换为第四坐标系下的坐标;其中,所述第四坐标系 的X、Y轴与所述第二坐标系X、Y轴平行,Z轴方向指向所述光电位 置传感器感光面中心,原点位于角锥棱镜顶点。
根据所述角锥棱镜顶点在所述第一坐标系中的坐标,获取所述激 光束在所述第一坐标系中对应的第一空间向量;根据所述光斑在所述 第四坐标系中的坐标,获取所述激光束在所述第四坐标系中对应的第 二空间向量;根据所述第一空间向量、所述第二空间向量以及所述靶 标在第二坐标系中的第一姿态角,获取所述靶标在所述第一坐标系中 的第二姿态角,进而获得所述被测物的姿态角。
其中,第四坐标系可以理解为对应于靶标的坐标系,第四坐标系 在靶标运动(即被测物运动)时也跟随靶标做相应的运动,故第四坐 标系相对于第二坐标系的姿态角即为靶标在第二坐标系中的第一姿态 角。激光跟踪测量单元所发出激光通过角锥棱镜开孔打在光电位置传 感器上,该激光光束可以在第一坐标系与第四坐标系下分别表示,即 第一空间向量和第二空间向量,根据向量在不同坐标系之间的转换关 系即可以得到第四坐标系相对于第一坐标系的第二姿态角,进而得到 靶标在激光跟踪测量单元中的姿态角。
具体地,根据光斑在第三坐标系中的坐标结合预设距离,获取光 斑在第四坐标系中的坐标;根据角锥棱镜的顶点在第一坐标系中的坐 标,获取激光束在第三坐标系中对应的第一向量,根据光斑在第四坐 标系中的坐标,获取激光束在第四坐标系中对应的第二向量;根据第 一坐标向量,第二坐标向量以及靶标在第二坐标系中的第一姿态角, 获取靶标在第一坐标系的第二姿态角,将第一姿态角、第二姿态角转 换到同一坐标系下,通过加权最小二乘融合获得最终姿态角。
本发明实施例提供了一种基于加权最小二乘融合的激光跟踪姿态 测量系统,通过使用加权最小二乘融合的方法,综合利用角锥棱镜与 光电位置传感器测量方位角与俯仰角时的高分辨率,单目视觉测量横 滚角时的高灵敏度,基于激光光束向量在不同坐标系之间的转换关系 完成姿态角的测量。该测量系统及方法提高了大尺寸范围内激光跟踪 姿态测量的整体精度,克服了单目姿态测量时精度随距离增加线性递 减的问题,同时改善了利用纵向投影比进行姿态测量时近距离精度较 低的问题。
在上述实施例中,所述激光跟踪测量单元为全站仪或激光跟踪仪, 所述单目视觉测量单元为工业相机。
在上述实施例中,所述多个靶点的数量不少于4个且非共面,且 每个靶点由反光标志或红外LED制成。
具体地,每个靶点由反光标志或红外LED制成,在实际测量中, 可以根据具体需求进行选用。
图3为本发明实施例提供的一种利用上述基于加权最小二乘融合 的激光跟踪姿态测量系统进行姿态角测量的方法的流程图,如图3所 示,包括:
S301,将所述靶标固定放置在被测物上,将所述激光跟踪测量单 元与所述单目视觉测量单元分别固定在预设位置上,通过标定获取第 一坐标系与第二坐标系之间的转换矩阵;
S302,利用所述光电位置传感器获取所属光斑在所述第三坐标系 下的坐标,利用所述单目视觉测量单元通过所述不少于4个的非共面 靶点获取所述靶标在所述第二坐标系下的第一姿态角,利用所述激光 跟踪测量单元获取所述角锥棱镜的顶点在第一坐标系下的坐标;
S303,根据所述光斑在所述第三坐标系下的坐标和所述角锥棱镜 的顶点与光电位置传感器之间间距,获取所述光斑在所述第四坐标系 下的坐标;
S304,根据所述角锥棱镜的顶点在所述第一坐标系下的坐标,获 取所述激光在所述第一坐标系中对应的第一空间向量,根据所述光斑 在所述第四坐标系下的坐标,获取所述激光在所述第四坐标系中对应 的第二空间向量;
S305,根据所述靶标在所述第二坐标系下的第一姿态角、所述第 一空间向量、所述第二空间向量以及所述第一坐标系与第二坐标系之 间的转换矩阵,获取所述靶标在所述第一坐标系下的第二姿态角;
S306,将所述第一姿态角、所述第二姿态角转换到同一坐标系下, 通过加权最小二乘融合获得最终姿态角。
在步骤S301中,将所述靶标固定放置在被测物上,将所述激光跟 踪测量单元与所述单目视觉测量单元分别固定在预设位置上,通过标 定获取第一坐标系与第二坐标系之间的转换矩阵,其中标定方法可以 理解为现有技术,在此不再赘述。
在步骤S302-S305中,第四坐标系可以理解为对应于靶标的坐标 系,第四坐标系在靶标运动(即被测物运动)时也跟随靶标做相应的 运动,故第四坐标系相对于第二坐标系的姿态角即为靶标在第二坐标 系中的第一姿态角。激光跟踪测量单元所发出激光通过角锥棱镜开孔 打在光电位置传感器上,该激光光束可以在第一坐标系与第四坐标系 下分别表示,即第一空间向量和第二空间向量,根据向量在不同坐标 系之间的转换关系即可以得到第四坐标系相对于第一坐标系的第二姿 态角,进而得到靶标在激光跟踪测量单元中的姿态角。
本发明实施例提供了一种姿态角测量方法,通过使用加权最小二 乘融合的方法,综合利用角锥棱镜与光电位置传感器测量方位角与俯 仰角时的高分辨率,单目视觉测量横滚角时的高灵敏度,基于激光光 束向量在不同坐标系之间的转换关系完成姿态角的测量。该测量系统 及方法提高了大尺寸范围内激光跟踪姿态测量的整体精度,克服了单 目姿态测量时精度随距离增加线性递减的问题,同时改善了利用纵向 投影比进行姿态测量时近距离精度较低的问题。
在上述实施例中,所述利用所述单目视觉测量单元通过所述不少 于4个的非共面靶点获取所述靶标在所述第二坐标系下的第一姿态角, 具体包括:
在进行姿态角测量时,首先利用所述单目视觉测量单元获取所述 靶标的像素坐标,通过EPNP算法进行二维点到三维点的匹配,获取 所述第二坐标系与所述第四坐标系之间的转换矩阵与位移向量;
具体地,利用所述单目视觉测量单元采集所述靶标在第一状态下 的所述多个靶点的图像,利用预设图像处理算法对第一状态下的图像 进行处理,得到所述多个靶点的像素坐标(ui,vi),i=1,2,3,…,n。假设靶标 中靶点的坐标为右上标t用来标注在第四坐标系下的坐 标,虚拟控制点在第四坐标系下的非齐次坐标为靶标靶点的坐标与虚拟控制点的坐标关系如 式(13)。由线性关系在欧式变换下的不变性可知,存在式(14),右 上标c用来标注在相机坐标系下的坐标。
式中,fx,fy,ux,uy为相机的内部参数,当有n个空间和图像 点对应时,可得到含2n个方程的线性方程组,记为矩阵形式Mx=0, 其中M是2n×12的矩阵,向量为12×1的向量, 包含了n个虚拟控制点在第二坐标系下的非齐次坐标。根据欧式变换 的保距性,由n个靶标靶点之间的距离已知,即可求出n个靶标靶点 在第二坐标系下的坐标,进而求得第二坐标系与所述第四坐标系之间 的转换矩阵与位移向量。
将所述转换矩阵与位移向量作为初始值,通过SoftPOSIT算法完 成后续实时状态下的姿态解算,获取单目视觉测量时姿态角,即所述 靶标在所述第二坐标系下的第一姿态角。
具体地,由于现实中镜头畸变的存在,空间点P的像点会偏离实 际位置,假设偏离前后像点坐标分别为Q(xi,yi)与QL(xi',yi'),位姿求解 问题可转变为求解两者距离最小值的问题,可表示为:
设旋转矩阵R=[R1 R2 R3]T,其中为矩阵R的行向量,平 移矩阵T=[Tx TyTz]T,f为相机的焦距,w为比例系数,P0为靶标坐 标系中心,空间点P的齐次坐标为P0P=(Xi,Yi,Zi,1)T,图像坐标为(xi,yi), 则有
设w=f/Tz,代入式(17)并除以Tz可得
将式(18)带入式(16),设M1=s(R1,Tx),M2=s(R2,Ty),Si=(Xi,Yi,Zi,1)T, 即可得
位姿估计的过程转变为求解D的最小值,满足D最小值的位姿为 SoftPOSIT最后一次求解得到的位姿信息,位姿信息包括位置信息与姿 态信息,求解所得姿态信息即靶标在所述第二坐标系下的第一姿态角, 具体包括方位角α1、俯仰角β1、横滚角γ1。
在上述实施例中,所述根据所述光斑在所述第三坐标系下的坐标 和所述角锥棱镜的顶点与光电位置传感器之间间距,获取所述光斑在 所述第四坐标系下的坐标,具体包括:
将所述光斑在所述第三坐标系下的横坐标与纵坐标,分别作为所 述光斑在第四坐标系下的横坐标与纵坐标,将所述角锥棱镜的顶点与 光电位置传感器之间间距作为所述光斑的Z轴坐标,即所述光斑在所 述第四坐标系下的坐标。
具体地,由于第三坐标系与第四坐标系X、Y轴方向平行,第四 坐标系的Z轴方向指向第三坐标系原点,故光斑在第三坐标系下的横 坐标与纵坐标及其在第四坐标系下的横坐标与纵坐标,光斑在第四坐 标系下的Z轴方向坐标即为预设间距。
预设间距的具体值可结合所需被测物相对于激光跟踪测量单元转 动角度确定,假设间距为h,光电位置传感器感光面边长为L,转动角 度为λ,则
在上述实施例中,所述根据所述角锥棱镜的顶点在所述第一坐标 系下的坐标,获取所述激光在所述第一坐标系中对应的第一空间向量, 根据所述光斑在所述第四坐标系下的坐标,获取所述激光在所述第四 坐标系中对应的第二空间向量,具体包括:
将所述第一坐标系的原点作为起点,所述角锥棱镜的顶点作为终 点,获取所述第一空间向量,将所述第四坐标系的原点作为起点,所 述光斑作为终点,获取所述第二空间向量。
具体地,根据第一坐标系原点的坐标值与角锥棱镜的顶点的坐标 值,可以计算得到第一空间向量的坐标值;根据第四坐标系的原点的 坐标值与光斑的坐标值,可以计算得到第二空间向量的坐标值。
在上述实施例中,所述根据所述靶标在所述第二坐标系下的第一 姿态角、所述第一空间向量、所述第二空间向量以及所述第一坐标系 与第二坐标系之间的转换矩阵,获取所述靶标在所述第一坐标系下的 第二姿态角,具体包括:
对所述第一空间向量和所述第二空间向量分别进行单位化处理, 得到第一单位向量和第二单位向量;
根据所述靶标在所述第二坐标系下的第一姿态角、所述第一单位 向量、所述第二单位向量以及所述第一坐标系与第二坐标系之间的转 换矩阵,基于光束向量在不同坐标系下的转换关系,获取靶标在所述 第一坐标系下的第二姿态角。
将单目视觉解算横滚角γ1与式(22)代入式(21)即可求解其余两 个角度α2、β2,将这三个角通过第一坐标系与第二坐标系之间的关系转 化到第一坐标系下,获取第二姿态角α'2、β′2、γ'2。
在上述实施例中,所述将所述第一姿态角、所述第二姿态角转换 到同一坐标系下,通过加权最小二乘融合获得最终姿态角,具体包括:
将所述第一姿态角通过所述第一坐标系与第二坐标系之间的转换 矩阵转换到第一坐标系下,获取第三姿态角α′1、β′1、γ′1;
将第二姿态角与第三姿态角的测量方差代入式(5),即可得到各 个姿态角归一化之前的权值;
本发明的优点可以通过以下仿真实验进行进一步说明:
假定激光跟踪设备使用激光跟踪仪,测点误差为0.01mm/m,单目 视觉测量单元分辨率为2590×2048pixel,像元尺寸大小为4.8um,工 业相机位于激光跟踪仪Z方向上0.3m处,光电位置传感器测点误差为 0.005mm,图像处理误差为0.1pixel。预设姿态角为[0°,0°,0°],激光跟 踪仪与靶标之间间距为3~15m。
本专利将单目视觉自动测量(简称为单目视觉法)、基于纵向投影 比结合光束向量在不同坐标系下的转换关系解算(简称为纵向投影比 法)、EPNP与SoftPOSIT联合解算并融合光束向量在不同坐标系下的 转换关系解算(简称为加权融合法)三种方法进行对比,以激光跟踪 仪与靶标之间距离为约束,姿态评价的方式采用均方差的方式,假设 姿态均方差为ζi,i=1,2,3,则:
仿真结果如图4所示,从图中可以看出,在距离为3~8m时,使用 单目视觉法精度较高,在距离为9~15m时纵向投影比法表现较好,而 加权融合法在各个位置都能取得较好的效果。
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而 非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领 域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技 术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修 改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方 案的精神和范围。
Claims (9)
1.一种基于加权最小二乘的激光跟踪姿态测量系统,其特征在于,包括:激光跟踪测量单元、单目视觉测量单元以及靶标;其中,
所述激光跟踪测量单元和所述单目视觉测量单元分别固定在预设位置,所述靶标固定放置在被测物上;所述靶标包括角锥棱镜、光电位置传感器以及多个靶点;所述角锥棱镜在顶点处开孔,所述光电位置传感器的感光面与角锥棱镜开孔面平行,所述角锥棱镜的顶点与光电位置传感器之间预留一定间距,所述激光跟踪测量单元所发出激光通过所述角锥棱镜开孔打在所述光电位置传感器的感光面上形成光斑,所述多个靶点的数量不少于4个,且所述不少于4个的靶点非共面;
所述激光跟踪测量单元用于获取所述角锥棱镜顶点在第一坐标系中的坐标;其中,所述第一坐标系为所述激光跟踪测量单元的自身坐标系;
所述单目视觉测量单元用于通过所述多个靶点获取所述靶标在第二坐标系中的第一姿态角;其中,所述第二坐标系为所述单目视觉测量单元的自身坐标系;
所述光电位置传感器用于获取所述光斑在第三坐标系中的坐标;其中,所述第三坐标系为所述光电位置传感器的自身坐标系;
所述第三坐标系下的坐标可结合所述角锥棱镜的顶点与光电位置传感器之间间距转换为第四坐标系下的坐标;其中,所述第四坐标系的X、Y轴与所述第二坐标系X、Y轴平行,Z轴方向指向所述光电位置传感器感光面中心,原点位于角锥棱镜顶点;
根据所述角锥棱镜顶点在所述第一坐标系中的坐标,获取所述激光束在所述第一坐标系中对应的第一空间向量;根据所述光斑在所述第四坐标系中的坐标,获取所述激光束在所述第四坐标系中对应的第二空间向量;根据所述第一空间向量、所述第二空间向量以及所述靶标在第二坐标系中的第一姿态角,获取所述靶标在所述第一坐标系中的第二姿态角,进而获得所述被测物的姿态角。
2.根据权利要求1所述基于加权最小二乘融合的激光跟踪姿态测量系统,其特征在于,所述激光跟踪测量单元为全站仪或激光跟踪仪,所述单目视觉测量单元为工业相机。
3.根据权利要求1所述基于加权最小二乘融合的激光跟踪姿态测量系统,其特征在于,所述不少于4个的非共面靶点由反光标志或红外LED制成。
4.一种利用权利要求1-3所述基于加权最小二乘融合的激光跟踪姿态测量的方法,其特征在于,包括:
将所述靶标固定放置在被测物上,将所述激光跟踪测量单元与所述单目视觉测量单元分别固定在预设位置上,通过标定获取第一坐标系与第二坐标系之间的转换矩阵;
利用所述光电位置传感器获取所属光斑在所述第三坐标系下的坐标,利用所述单目视觉测量单元通过所述不少于4个的非共面靶点获取所述靶标在所述第二坐标系下的第一姿态角,利用所述激光跟踪测量单元获取所述角锥棱镜的顶点在第一坐标系下的坐标;
根据所述光斑在所述第三坐标系下的坐标和所述角锥棱镜的顶点与光电位置传感器之间间距,获取所述光斑在所述第四坐标系下的坐标;
根据所述角锥棱镜的顶点在所述第一坐标系下的坐标,获取所述激光在所述第一坐标系中对应的第一空间向量,根据所述光斑在所述第四坐标系下的坐标,获取所述激光在所述第四坐标系中对应的第二空间向量;
根据所述靶标在所述第二坐标系下的第一姿态角、所述第一空间向量、所述第二空间向量以及所述第一坐标系与第二坐标系之间的转换矩阵,获取所述靶标在所述第一坐标系下的第二姿态角;
将所述第一姿态角、所述第二姿态角转换到同一坐标系下,通过加权最小二乘融合获得最终姿态角。
5.根据权利要求4所述方法,其特征在于,所述利用所述单目视觉测量单元通过所述不少于4个的非共面靶点获取所述靶标在所述第二坐标系下的第一姿态角,具体包括:
在进行姿态角测量时,首先利用所述单目视觉测量单元获取所述靶标的像素坐标,通过EPNP算法进行二维点到三维点的匹配,获取所述第二坐标系与所述第四坐标系之间的转换矩阵与位移向量;
将所述转换矩阵与位移向量作为初始值,通过SoftPOSIT算法完成后续实时状态下的姿态解算,获取单目视觉测量时姿态角,即所述靶标在所述第二坐标系下的第一姿态角。
6.根据权利要求4所述方法,其特征在于,所述根据所述光斑在所述第三坐标系下的坐标和所述角锥棱镜的顶点与光电位置传感器之间间距,获取所述光斑在所述第四坐标系下的坐标,具体包括:
将所述光斑在所述第三坐标系下的横坐标与纵坐标,分别作为所述光斑在第四坐标系下的横坐标与纵坐标,将所述角锥棱镜的顶点与光电位置传感器之间间距作为所述光斑的Z轴坐标,即所述光斑在所述第四坐标系下的坐标。
7.根据权利要求4所述方法,其特征在于,所述根据所述角锥棱镜的顶点在所述第一坐标系下的坐标,获取所述激光在所述第一坐标系中对应的第一空间向量,根据所述光斑在所述第四坐标系下的坐标,获取所述激光在所述第四坐标系中对应的第二空间向量,具体包括:
将所述第一坐标系的原点作为起点,所述角锥棱镜的顶点作为终点,获取所述第一空间向量,将所述第四坐标系的原点作为起点,所述光斑作为终点,获取所述第二空间向量。
8.根据权利要求4所述方法,其特征在于,所述根据所述靶标在所述第二坐标系下的第一姿态角、所述第一空间向量、所述第二空间向量以及所述第一坐标系与第二坐标系之间的转换矩阵,获取所述靶标在所述第一坐标系下的第二姿态角,具体包括:
对所述第一空间向量和所述第二空间向量分别进行单位化处理,得到第一单位向量和第二单位向量;
根据所述靶标在所述第二坐标系下的第一姿态角、所述第一单位向量、所述第二单位向量以及所述第一坐标系与第二坐标系之间的转换矩阵,基于光束向量在不同坐标系下的转换关系,获取靶标在所述第一坐标系下的第二姿态角。
9.根据权利要求4所述方法,其特征在于,所述将所述第一姿态角、所述第二姿态角转换到同一坐标系下,通过加权最小二乘融合获得最终姿态角,具体包括:
将所述第一姿态角通过所述第一坐标系与第二坐标系之间的转换矩阵转换到第一坐标系下,获取第三姿态角;
假设所述第二姿态角、所述第三姿态角为传感器所测结果,将两个传感器对于姿态角的观测方程设置为:
Y=Hx+e (1)
式中:x为一维状态量;Y为二维测量向量,设Y=[y1,y2]T;e为二维测量噪声向量,包含传感器的内部噪声和环境干扰噪声,设e=[e1,e2]T,若传感器的测量方差为Ri,i=1,2,则E[ei 2]=E[(x-yi)2]=Ri;H为已知的二维常向量,设H=[1,1]T。
x=x-x (3)
由式(2)和式(3),可得估计方差为:
对式(4)求极小值,取wi的偏导数,并令其等于零得
将(5)代入式(4),可得估计误差为
由式(6)不难看出,采用加权融合的估计方差比任何一个传感器的测量方差都小。
通过上述过程进行加权最小二乘融合,获取最终姿态角。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110250220.XA CN113028990B (zh) | 2021-03-08 | 2021-03-08 | 一种基于加权最小二乘的激光跟踪姿态测量系统及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110250220.XA CN113028990B (zh) | 2021-03-08 | 2021-03-08 | 一种基于加权最小二乘的激光跟踪姿态测量系统及方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113028990A true CN113028990A (zh) | 2021-06-25 |
CN113028990B CN113028990B (zh) | 2022-11-18 |
Family
ID=76466698
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110250220.XA Active CN113028990B (zh) | 2021-03-08 | 2021-03-08 | 一种基于加权最小二乘的激光跟踪姿态测量系统及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113028990B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114322886A (zh) * | 2022-01-10 | 2022-04-12 | 深圳市中图仪器股份有限公司 | 具有多传感器的姿态探头 |
CN115493616A (zh) * | 2022-08-03 | 2022-12-20 | 湖北工业大学 | 一种激光跟踪姿态角现场精度的评定方法 |
CN116165673A (zh) * | 2022-12-12 | 2023-05-26 | 深圳市中图仪器股份有限公司 | 基于反向跟踪的六维探头 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101363714A (zh) * | 2008-09-26 | 2009-02-11 | 浙江大学 | 一种机翼水平位姿测量与评估方法 |
CN101363715A (zh) * | 2008-09-26 | 2009-02-11 | 浙江大学 | 基于激光跟踪仪的飞机机身姿态计算方法 |
CN106312397A (zh) * | 2016-10-12 | 2017-01-11 | 华南理工大学 | 一种激光视觉引导的焊接轨迹自动跟踪系统及方法 |
US20180216941A1 (en) * | 2015-07-31 | 2018-08-02 | Tianjin University | Indoor mobile robot position and posture measurement system based on photoelectric scanning and measurement method |
CN109448055A (zh) * | 2018-09-20 | 2019-03-08 | 中国科学院光电研究院 | 单目视觉姿态确定方法及系统 |
CN109540058A (zh) * | 2018-11-21 | 2019-03-29 | 北京航天新风机械设备有限责任公司 | 一种基于点集测量的刚体位置与姿态测量不确定度评估方法 |
CN109737913A (zh) * | 2018-11-23 | 2019-05-10 | 湖北工业大学 | 一种激光跟踪姿态角测量系统及方法 |
CN110017810A (zh) * | 2019-05-16 | 2019-07-16 | 湖北工业大学 | 一种光电位置传感器与单目视觉组合姿态测量系统及方法 |
CN111993159A (zh) * | 2020-08-27 | 2020-11-27 | 江苏科技大学 | 一种轴类工件在位非接触检测方法 |
-
2021
- 2021-03-08 CN CN202110250220.XA patent/CN113028990B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101363714A (zh) * | 2008-09-26 | 2009-02-11 | 浙江大学 | 一种机翼水平位姿测量与评估方法 |
CN101363715A (zh) * | 2008-09-26 | 2009-02-11 | 浙江大学 | 基于激光跟踪仪的飞机机身姿态计算方法 |
US20180216941A1 (en) * | 2015-07-31 | 2018-08-02 | Tianjin University | Indoor mobile robot position and posture measurement system based on photoelectric scanning and measurement method |
CN106312397A (zh) * | 2016-10-12 | 2017-01-11 | 华南理工大学 | 一种激光视觉引导的焊接轨迹自动跟踪系统及方法 |
CN109448055A (zh) * | 2018-09-20 | 2019-03-08 | 中国科学院光电研究院 | 单目视觉姿态确定方法及系统 |
CN109540058A (zh) * | 2018-11-21 | 2019-03-29 | 北京航天新风机械设备有限责任公司 | 一种基于点集测量的刚体位置与姿态测量不确定度评估方法 |
CN109737913A (zh) * | 2018-11-23 | 2019-05-10 | 湖北工业大学 | 一种激光跟踪姿态角测量系统及方法 |
CN110017810A (zh) * | 2019-05-16 | 2019-07-16 | 湖北工业大学 | 一种光电位置传感器与单目视觉组合姿态测量系统及方法 |
CN111993159A (zh) * | 2020-08-27 | 2020-11-27 | 江苏科技大学 | 一种轴类工件在位非接触检测方法 |
Non-Patent Citations (2)
Title |
---|
李丽娟等: "加权整体最小二乘在激光跟踪仪转站中的应用", 《光学精密工程》 * |
杨剑等: "加权最小二乘算法在机器视觉系统中的应用", 《光学精密工程》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114322886A (zh) * | 2022-01-10 | 2022-04-12 | 深圳市中图仪器股份有限公司 | 具有多传感器的姿态探头 |
CN114322886B (zh) * | 2022-01-10 | 2024-03-22 | 深圳市中图仪器股份有限公司 | 具有多传感器的姿态探头 |
CN115493616A (zh) * | 2022-08-03 | 2022-12-20 | 湖北工业大学 | 一种激光跟踪姿态角现场精度的评定方法 |
CN116165673A (zh) * | 2022-12-12 | 2023-05-26 | 深圳市中图仪器股份有限公司 | 基于反向跟踪的六维探头 |
CN116165673B (zh) * | 2022-12-12 | 2023-09-22 | 深圳市中图仪器股份有限公司 | 基于反向跟踪的六维探头 |
Also Published As
Publication number | Publication date |
---|---|
CN113028990B (zh) | 2022-11-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113028990B (zh) | 一种基于加权最小二乘的激光跟踪姿态测量系统及方法 | |
CN108571971B (zh) | 一种agv视觉定位系统及方法 | |
CN110116407B (zh) | 柔性机器人位姿测量方法及装置 | |
CN108012325B (zh) | 一种基于uwb和双目视觉的导航定位方法 | |
CN109737913B (zh) | 一种激光跟踪姿态角测量系统及方法 | |
CN111964694B (zh) | 一种用于三维测量的激光测距仪标定方法 | |
CN108594245A (zh) | 一种目标物运动监测系统及方法 | |
CN104315995B (zh) | 基于虚拟多立方体标准目标的tof深度相机三维坐标标定装置和方法 | |
CN106871787B (zh) | 大空间线扫描成像三维测量方法 | |
Shang et al. | Measurement methods of 3D shape of large-scale complex surfaces based on computer vision: A review | |
CN104316083B (zh) | 一种虚拟多球体球心定位的tof深度相机三维坐标标定装置和方法 | |
CN109855822B (zh) | 一种基于无人机的高铁桥梁竖向动扰度测量方法 | |
CN111220126A (zh) | 一种基于点特征和单目相机的空间物体位姿测量方法 | |
CN109212497A (zh) | 一种空间六自由度车载雷达天线位姿偏差测量及对接方法 | |
CN113048938B (zh) | 一种合作目标设计及姿态角测量系统及方法 | |
CN112045655B (zh) | 用于大尺度多站点场景的移动机器人位姿测量方法及系统 | |
US20190339071A1 (en) | Marker, and Posture Estimation Method and Position and Posture Estimation Method Using Marker | |
CN114413958A (zh) | 无人物流车的单目视觉测距测速方法 | |
Liu et al. | A high-accuracy pose measurement system for robotic automated assembly in large-scale space | |
CN109813509A (zh) | 基于无人机实现高铁桥梁竖向动扰度测量的方法 | |
Peng et al. | A measuring method for large antenna assembly using laser and vision guiding technology | |
CN208350997U (zh) | 一种目标物运动监测系统 | |
CN113340272B (zh) | 一种基于无人机微群的地面目标实时定位方法 | |
Li et al. | Extrinsic calibration of non-overlapping multi-camera system with high precision using circular encoded point ruler | |
CN111207688B (zh) | 在载运工具中测量目标对象距离的方法、装置和载运工具 |
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 |