CN112697073A - 一种三维姿态测量方法 - Google Patents

一种三维姿态测量方法 Download PDF

Info

Publication number
CN112697073A
CN112697073A CN202011248011.3A CN202011248011A CN112697073A CN 112697073 A CN112697073 A CN 112697073A CN 202011248011 A CN202011248011 A CN 202011248011A CN 112697073 A CN112697073 A CN 112697073A
Authority
CN
China
Prior art keywords
theta
cmos imaging
imaging elements
substrate
dimensional
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
Application number
CN202011248011.3A
Other languages
English (en)
Other versions
CN112697073B (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.)
Wuhan No 2 Ship Design Institute No 719 Research Institute of China Shipbuilding Industry Corp
Original Assignee
Wuhan No 2 Ship Design Institute No 719 Research Institute of China Shipbuilding Industry Corp
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 Wuhan No 2 Ship Design Institute No 719 Research Institute of China Shipbuilding Industry Corp filed Critical Wuhan No 2 Ship Design Institute No 719 Research Institute of China Shipbuilding Industry Corp
Priority to CN202011248011.3A priority Critical patent/CN112697073B/zh
Publication of CN112697073A publication Critical patent/CN112697073A/zh
Application granted granted Critical
Publication of CN112697073B publication Critical patent/CN112697073B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00Measuring arrangements characterised by the use of optical techniques
    • G01B11/26Measuring arrangements characterised by the use of optical techniques for measuring angles or tapers; for testing the alignment of axes
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00Measuring arrangements characterised by the use of optical techniques
    • G01B11/002Measuring arrangements characterised by the use of optical techniques for measuring two or more coordinates

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Length Measuring Devices By Optical Means (AREA)

Abstract

本发明公开了一种三维姿态测量方法,通过构建虚拟单维高分辨成像元件,显著提高了光斑斜率测量的精度,进而实现高精度的滚转角测量,使一次测量过程能够同时实现秒级三维高精度姿态测量。

Description

一种三维姿态测量方法
技术领域
本发明属于光学技术领域,具体涉及一种三维姿态测量方法。
背景技术
光电自准直仪是一种利用光的自准直原理测量小角度的仪器。现有技术中,在对待测量目标进行姿态测量时,利用光电自准直仪可有效实现方位、俯仰二维角度的秒级精度测量。关于滚转角测量,公开号CN1335483A等专利提出了偏振相位法测量方案,但是上述方法都无法同时实现高精度三维姿态测量。在现有技术中,也有技术提出了基于自准直仪的三维姿态测量方法,该方法中方位、俯仰角度测量具有常规自准直仪的高精度,但滚转角测量精度受限于目前的CMOS成像元件的分辨率,最高只能达到约10″量级。具体分析如下:
当前CMOS成像元件分辨率为千万像素级别,例如4800万像素,单维像素约为几千,例如4800万像素分别率即为8000×6000(4:3)。根据现有技术中滚转角的测量原理,需对成像进行斜率计算。如图2所示,十字光斑斜率a=y/L,因此滚转角测量精度取决于y的测量精度与L的最大值。不考虑算法优化y测量精度为1个像素,L最大值为成像元件的单维像素。按照前述4800万像素成像元件的分辨率计算,其斜率测量精度为1/8000=25.8〞,由于从原理上讲滚转角为斜率变化的1/2,因此滚转角测量精度约13〞。
总之,现有技术中采用光电自准直仪测量目标三维姿态的方法,仍然存在滚转角测量精度不高的问题。
发明内容
有鉴于此,本发明提供了一种三维姿态测量方法,能够解决当前成像元件分辨率不足的问题,实现对测量目标的秒级三维姿态测量。
本发明提供的一种三维姿态测量方法,包括以下步骤:
步骤1、将待测量目标与反射面组合体固定;
步骤2、将光电自准直仪中测量滚转角的光电传感器替换为虚拟单维高分辨成像元件;所述虚拟单维高分辨成像元件包括基底和两个CMOS成像元件,所述两个CMOS成像元件分别固定于所述基底的同侧两端,所述基底的长度大于10倍的所述CMOS成像元件的尺寸;
步骤3、采用公式(1)统一所述两个CMOS成像元件的坐标系:
Figure BDA0002770694710000021
其中,(x,y)、(x′,y′)分别为两个CMOS成像元件的坐标系中的坐标,θ、x0、y0均为坐标系变换参数;
θ的计算过程为:将待测量目标任意旋转一个滚转角度,在两个CMOS成像元件的坐标系中分别采用公式(2)和公式(3)拟合求得滚转角,计算两个坐标系下滚转角的差值θ1-θ′1即为θ的取值,即θ=θ1-θ′1,其中,b1及b′1均为拟合参数;
y=tanθ1*x+b1 (2)
y′=tanθ′1*x′+b′1 (3)
采用公式(4)计算x0、y0的取值:
Figure BDA0002770694710000022
步骤4、采用步骤3得到的统一坐标系,分别计算CMOS1和CMOS2上的光强质心,由质心连线即可计算得到光斑斜率。
进一步地,所述基底的材料为熔融石英。
有益效果:
本发明通过构建虚拟单维高分辨成像元件,显著提高了光斑斜率测量的精度,进而实现高精度的滚转角测量,使一次测量过程能够同时实现秒级三维高精度姿态测量。
附图说明
图1为本发明提供的一种三维姿态测量方法所采用的结构示意图。
图2为本发明提供的一种三维姿态测量方法所采用的传感器的结构示意图。
其中,1-分光镜A,2-分光镜B,3-滤光片A,4-滤光片B,5-物镜,6-分划板,7-光源,8-光电传感器A,9-光电传感器B及10-待测目标。
具体实施方式
下面结合附图并举实施例,对本发明进行详细描述。
本发明提出的技术方案是在申请号为201910994110.7的专利中所公开的技术方案的基础上的改进。申请号为201910994110.7的专利所公开的技术方案中采用反射面组合体代替平面反射镜,采用识别方法区分出平面反射面反射光和直角反射面反射光,通过分析平面反射面反射光计算待测量目标的俯仰角、方位角,通过分析直角反射面反射光计算待测量目标的滚转角,从而实现待测量目标三维姿态的测量。
本发明提供了一种三维姿态测量方法,其核心思想是:在采用反射面组合体代替平面反射镜的基础上,将现有光电自准直仪中测量滚转角的光电传感器替换为由两个CMOS成像元件组成的虚拟单维高分辨成像元件,通过分析计算两个CMOS成像元件的光斑来测量目标的滚转角。
本发明提供了一种三维姿态测量方法,如图1所示,具体包括以下步骤:
步骤1、将待测量目标与反射面组合体固定。
步骤2、将光电自准直仪中测量滚转角的光电传感器B9替换为虚拟单维高分辨成像元件,该虚拟单维高分辨成像元件的结构如图2所示。
虚拟单维高分辨成像元件包括基底和两个CMOS成像元件,其中,基底长度为CMOS成像元件尺寸10倍以上,宽度大于CMOS成像元件尺寸即可,在基底两端分别各固定安装一个CMOS成像元件,两个CMOS成像元件相对于基底的高度保持一致最佳。由此可将光电传感器B9看作一个中间成像部分缺失的单维高分辨率(长度方向)成像元件。
步骤3、统一两个CMOS成像元件的坐标系。设COMS1面阵坐标系为XY,COMS2面阵坐标系为X′Y′。以坐标统一到COMS1为例,COMS2上某点坐标为(x′,y′),在COMS1坐标系统中的坐标为(x,y),则有两个坐标系下坐标之间的关系满足公式(1):
Figure BDA0002770694710000041
其中,θ为X′Y′相对于XY坐标系的旋转角(以逆时针为正),(x0,y0)为XY坐标系的原点在X′Y′坐标系中的坐标,由此可见,只要求出变换参数θ、x0、y0即可得到两坐标系之间的变换关系。
将装置实施测量操作,将待测量目标任意旋转一个滚转角度,使反射十字光斑的一条光斑的像处于CMOS1和CMOS2视野中。改变测量目标滚转角度,使反射十字光斑的一条光斑的像仍处于CMOS1和CMOS2视野中,重复上述操作,得到旋转后的光斑在CMOS1坐标系中的拟合直线表达式如公式(2)所示,
y=tanθ1*x+b1 (2)
在CMOS2坐标系中的拟合直线表达式如公式(3)所示。
y′=tanθ′1*x′+b′1 (3)
然后,改变测量目标滚转角度,使反射十字光斑的一条光斑的像仍处于CMOS1和CMOS2视野中,重复上述操作,得到旋转后斑在CMOS1坐标系中的拟合直线表达式如公式(4)所示:
y=tanθ2*x+b2 (4)
在CMOS2坐标系中的拟合直线表达式如公式(5)所示:
y′=tanθ′2*x′+b′2 (5)
其中,b1及b′1均为拟合参数,则,通过计算θ1与θ′1的差值或计算θ2与θ′2的差值即可得到θ的值,即如公式(6)所示:
θ=θ1-θ′1=θ2-θ′2 (6)
公式(1)中的x0、y0采用公式(7)计算可得:
Figure BDA0002770694710000051
由此,公式(1)中所需要的参数θ、x0、y0均已计算得到,将上述参数代入公式(1)即可完成两个CMOS成像元件坐标系的统一。
步骤4、基于统一坐标系计算光斑斜率。采用步骤3得到的统一坐标系,分别计算CMOS1和CMOS2上的光强质心,由质心连线即可计算得到光斑斜率。
例如,采用热膨胀系数为5×10-7/℃的熔融石英作为基底,基底尺寸250mm×50mm。选取两个分辨率4000×3000,像元尺寸1.1m的CMOS分别固定于基底两端,两传感器距离约220mm。以此构成虚拟单维高分辨成像元件,长度方向相当于分辨率约220mm/1.1m=2×105,按照前述分析,对测量目标滚转角的分辨率可达1/(2×2×105)rad=0.5〞。
由于两CMOS相对距离远大于像元尺寸,需保证相对位置的稳定,尤其是在与两传感器连线垂直方向上,1个像元的相对位移,就可能带来与分辨率相当的测量误差。两CMOS固定在一块基底上,最大的位移不稳定性来源于温度引起的基底材料的热胀冷缩,因此采用低热膨胀系数材料作为基底,按照热膨胀系数为5×10-7/℃估计,环境温度变化1摄氏度,两CMOS距离仅变化约0.1m,远小于像元尺寸,因此温度对传感器精度影响可忽略不计。
综上所述,以上仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (2)

1.一种三维姿态测量方法,其特征在于,包括以下步骤:
步骤1、将待测量目标与反射面组合体固定;
步骤2、将光电自准直仪中测量滚转角的光电传感器替换为虚拟单维高分辨成像元件;所述虚拟单维高分辨成像元件包括基底和两个CMOS成像元件,所述两个CMOS成像元件分别固定于所述基底的同侧两端,所述基底的长度大于10倍的所述CMOS成像元件的尺寸;
步骤3、采用公式(1)统一所述两个CMOS成像元件的坐标系:
Figure FDA0002770694700000011
其中,(x,y)、(x′,y′)分别为两个CMOS成像元件的坐标系中的坐标,θ、x0、y0均为坐标系变换参数;
θ的计算过程为:将待测量目标任意旋转一个滚转角度,在两个CMOS成像元件的坐标系中分别采用公式(2)和公式(3)拟合求得滚转角,计算两个坐标系下滚转角的差值θ1-θ′1即为θ的取值,即θ=θ1-θ′1,其中,b1及b′1均为拟合参数;
y=tanθ1*x+b1 (2)
y′=tanθ′1*x′+b′1 (3)
采用公式(4)计算x0、y0的取值:
Figure FDA0002770694700000012
步骤4、采用步骤3得到的统一坐标系,分别计算CMOS1和CMOS2上的光强质心,由质心连线即可计算得到光斑斜率。
2.根据权利要求1所述的方法,其特征在于,所述基底的材料为熔融石英。
CN202011248011.3A 2020-11-10 2020-11-10 一种三维姿态测量方法 Active CN112697073B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011248011.3A CN112697073B (zh) 2020-11-10 2020-11-10 一种三维姿态测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011248011.3A CN112697073B (zh) 2020-11-10 2020-11-10 一种三维姿态测量方法

Publications (2)

Publication Number Publication Date
CN112697073A true CN112697073A (zh) 2021-04-23
CN112697073B CN112697073B (zh) 2022-07-15

Family

ID=75507003

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011248011.3A Active CN112697073B (zh) 2020-11-10 2020-11-10 一种三维姿态测量方法

Country Status (1)

Country Link
CN (1) CN112697073B (zh)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100157058A1 (en) * 2008-11-20 2010-06-24 Hella Kgaa Hueck & Co. Method and Device for Compensating a Roll Angle
CN103697864A (zh) * 2013-12-27 2014-04-02 武汉大学 一种基于大虚拟相机的窄视场双相机影像拼接方法
CN105021211A (zh) * 2015-06-05 2015-11-04 中国船舶重工集团公司第七0七研究所 一种基于自准直仪的姿态测试装置及方法
CN206989960U (zh) * 2017-07-31 2018-02-09 天津大学 一种动态位姿实时测量装置
EP3435333A1 (en) * 2017-07-26 2019-01-30 Aptiv Technologies Limited Method of determining the roll angle of a vehicle mounted camera
CN109343578A (zh) * 2018-11-16 2019-02-15 北京理工大学 一种基于视觉反馈的tdi-ccd相机双目视觉系统视场对齐方法
CN111220094A (zh) * 2019-10-18 2020-06-02 中国船舶重工集团公司第七一九研究所 一种基于光电自准直仪的三维姿态测量方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100157058A1 (en) * 2008-11-20 2010-06-24 Hella Kgaa Hueck & Co. Method and Device for Compensating a Roll Angle
CN103697864A (zh) * 2013-12-27 2014-04-02 武汉大学 一种基于大虚拟相机的窄视场双相机影像拼接方法
CN105021211A (zh) * 2015-06-05 2015-11-04 中国船舶重工集团公司第七0七研究所 一种基于自准直仪的姿态测试装置及方法
EP3435333A1 (en) * 2017-07-26 2019-01-30 Aptiv Technologies Limited Method of determining the roll angle of a vehicle mounted camera
CN206989960U (zh) * 2017-07-31 2018-02-09 天津大学 一种动态位姿实时测量装置
CN109343578A (zh) * 2018-11-16 2019-02-15 北京理工大学 一种基于视觉反馈的tdi-ccd相机双目视觉系统视场对齐方法
CN111220094A (zh) * 2019-10-18 2020-06-02 中国船舶重工集团公司第七一九研究所 一种基于光电自准直仪的三维姿态测量方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
GUAN WEIPENG 等: "High Precision Indoor Visible Light Positioning Algorithm Based on Double LEDs Using CMOS Image Sensor", 《APPLIED SCIENCES-BASEL》 *
吕勇 等: "基于双CCD的滚转角测量方法研究", 《光学技术》 *
吕勇: "基于多准直光的位姿测量方法及应用研究", 《中国优秀博硕士学位论文全文数据库(博士)工程科技Ⅱ辑》 *
孙国燕 等: "三维姿态角高精度测量装置", 《光学精密工程》 *
张士涛: "音圈式大行程快速反射镜及其视轴稳定技术研究", 《中国优秀博硕士学位论文全文数据库(博士)工程科技Ⅱ辑》 *

Also Published As

Publication number Publication date
CN112697073B (zh) 2022-07-15

Similar Documents

Publication Publication Date Title
US8675077B2 (en) Alignment metrology and resolution measurement system for imaging arrays
EP2435783B1 (en) Hybrid sensor
CN110140023B (zh) 标记器和使用了标记器的姿势推定方法及位置姿势推定方法
CN107339935B (zh) 用于全视角扫描测量系统的靶标空间交会测量方法
CN113048938B (zh) 一种合作目标设计及姿态角测量系统及方法
WO2019080888A1 (zh) 多轴激光位移测量系统中干涉仪的安装偏差标定方法
CN110501026B (zh) 基于阵列星点的相机内方位元素标定装置及方法
CN109827521A (zh) 一种快速多线结构光视觉测量系统标定方法
CN105157700A (zh) 用于空间光学相机的无地面控制点定位方法和装置
CN106525054A (zh) 一种采用星上推扫遥感图像信息的单星自主测定轨方法
CN115854921A (zh) 基于结构光测量物体表面形状的方法及其测量装置
CN112729109B (zh) 一种点云数据的校正方法及装置
CN109990801B (zh) 基于铅垂线的水平仪装配误差标定方法
CN112697073B (zh) 一种三维姿态测量方法
Kim et al. Measurement of six-degree-of-freedom absolute postures using a phase-encoded pattern target and a monocular vision system
CN113916128A (zh) 一种基于光笔式视觉测量系统的提高精度的方法
Korolev et al. Use of information redundancy in optical digital measurement systems with 2D sensor
CN116105638A (zh) 一种u型转台同轴度检测系统及检测方法
CN114526693B (zh) 一种基于非标准柱面角锥镜的滚转角测量方法
US9052159B2 (en) System for determining the spatial orientation of a movable apparatus
CN105758339A (zh) 基于几何误差修正技术的光轴与物面垂直度检测方法
CN113390394B (zh) 一种带有摄影比例尺的光束法平差算法
CN114088060B (zh) 一种用于光学遥感卫星指向测量的星地相机成像系统
CN110260817B (zh) 基于虚拟标志点的复杂曲面偏折测量自定位方法
CN114413789B (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