CN108225371A - 一种惯导/相机安装误差标定方法 - Google Patents

一种惯导/相机安装误差标定方法 Download PDF

Info

Publication number
CN108225371A
CN108225371A CN201611150018.5A CN201611150018A CN108225371A CN 108225371 A CN108225371 A CN 108225371A CN 201611150018 A CN201611150018 A CN 201611150018A CN 108225371 A CN108225371 A CN 108225371A
Authority
CN
China
Prior art keywords
coordinate system
camera
inertial navigation
turntable
calibration
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
CN201611150018.5A
Other languages
English (en)
Other versions
CN108225371B (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 Automation Control Equipment Institute BACEI
Original Assignee
Beijing Automation Control Equipment Institute BACEI
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 Automation Control Equipment Institute BACEI filed Critical Beijing Automation Control Equipment Institute BACEI
Priority to CN201611150018.5A priority Critical patent/CN108225371B/zh
Publication of CN108225371A publication Critical patent/CN108225371A/zh
Application granted granted Critical
Publication of CN108225371B publication Critical patent/CN108225371B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C25/00Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass

Landscapes

  • Engineering & Computer Science (AREA)
  • Manufacturing & Machinery (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Gyroscopes (AREA)
  • Navigation (AREA)

Abstract

本发明属于视觉导航与组合导航技术领域,具体涉及一种惯导/相机安装误差标定方法。本方法包括如下步骤:步骤一、系统安装及其连接方式;步骤二、转台旋转并记录数据;步骤三、构建12个未知量的方程组并计算;步骤四、输出结果。本发明提出的标定方法无需激光扫描仪与差分卫星导航系统等繁重的设备,而是利用转台与标定板,在室内即可完成惯导与相机之间安装误差的标定,提高了工作效率,降低了标定成本。

Description

一种惯导/相机安装误差标定方法
技术领域
本发明属于视觉导航与组合导航技术领域,具体涉及一种惯导/相机安装误差标定方法。
背景技术
视觉导航以传感器获取的客观事物的视觉图像为基础,通过计算机来模拟人的视觉功能,对图像进行识别和理解,进而获取载体的导航信息。
视觉导航的传感器即光学相机。将惯导与相机捷联安装在一起,形成惯性/视觉组合导航系统,可以将惯性导航的绝对导航信息与视觉导航的相对导航信息进行融合,特别适用于无人机自主着陆、智能机器人导航、工业测量等领域。而惯性/视觉组合导航是建立在惯导/相机安装误差的精密标定基础上的。
惯导与相机捷联安装在一起,惯导坐标系与相机坐标系之间存在一个姿态转转与平移的关系,这个关系称之为相对姿态与相对位置,也称之为姿态安装误差与位置安装误差,统称为安装误差。
发明内容
根据上述问题,本发明提供一种惯导/相机安装误差标定方法,完成惯性/视觉组合导航系统中,惯导与相机姿态安装误差与位置安装误差的精密标定问题。
为了实现这一目的,本发明采取的技术方案是:
一种惯导/相机安装误差标定方法,该方法包括如下步骤:
步骤一、系统安装及其连接方式;步骤二、转台旋转并记录数据;步骤三、构建12个未知量的方程组并计算;步骤四、输出结果;
所述系统安装及其连接方式,将惯导与相机安装在转台上,将标定板放置在转台附近,且处于相机的视场范围内。
一种惯导/相机安装误差标定方法,步骤二、旋转转台并记录数据:
系统安装完成后,旋转转台4次,旋转时要求保证:
1、惯导与相机之间的相对位置与相对姿态不变;
2、标定板位置保持不变,从而保证转台与标定板之间的相对位置与相对姿态不变;
3、转台的XYZ三个轴同时旋转一定的角度(三个轴旋角的大小与方向不要求一致);
4、旋转时要求保证标定板不超出相机的视场范围;
记录4个状态下惯导与转台之间的标定矩阵相机与标定板之间的标定矩阵其中通过光瞄与转台旋转参数获得,通过高精度光学测量方法获得。
一种惯导/相机安装误差标定方法,步骤三、构建12个未知量的方程组并计算:
假设某空间点P由坐标系a中的坐标Pa转换为坐标系b中的坐标Pb,其过程可表示为,
公式(1),其中Pa=[xa,ya,za]T为空间点P在坐标系a中的坐标;是一个平移向量,其中tx为由坐标系a向坐标系b平移过程中x方向的平移量,ty为由坐标系a向坐标系b平移过程中y方向的平移量,tz为由坐标系a向坐标系b平移过程中x方向的平移量;是一个旋转矩阵,描述了从坐标系a向坐标系b的旋转过程,即坐标系a按的旋转顺序,分别旋转后,其三轴将与坐标系b平行,
T中的6个参数(α,β,γ,tx,ty,tz)描述了坐标系a与坐标系b之间的相对位置与相对姿态关系,安装误差标定就是要把这6个参数求解出来;
公式(9)可以整理成以下形式
其中
分别称为空间点P在坐标系a向坐标系b下的齐次坐标,中的元素由标定参数计算而来,称之为标定矩阵;
将惯导与相机安装在转台上,将标定板放置在转台附近,且处于相机的视场范围内,转台、惯导、相机、标定板四者之间的关系可写成归一化矩阵连乘的形式:
式中:
——4×4矩阵,惯导与转台之间的标定矩阵,通过光瞄与转台旋转参数方法求得,已知量;
——4×4矩阵,惯导与相机之间的标定矩阵,未知量,待求;
——4×4矩阵,相机与标定板之间的标定矩阵,通过光学测量方法求得,已知量;
——4×4矩阵,标定板与转台之间的标定矩阵,未知量;
在上面的连乘公式中,这两个未知的4×4矩阵中可以由共12个未知量表示;
将式(14)写成以下形式:
可以整理出3个以为变量的方程,分别表示为:
因此,得到该状态下12个未知量的3个方程;
转台旋转4次,每次都获得该状态下的3个方程,一共可以获得4个状态下12个未知量的12个方程,表示如下:
采用通用的矩阵分解法、最小二乘法求解以上线性方程组,得到线性方程组的解,获得共12个未知量的值。
一种惯导/相机安装误差标定方法,步骤四、输出结果,将步骤三中的作为惯导与相机之间的安装误差输出。
本发明的有益效果为:
本发明提出的标定方法无需激光扫描仪与差分卫星导航系统等繁重的设备,而是利用转台与标定板,在室内即可完成惯导与相机之间安装误差的标定,提高了工作效率,降低了标定成本。
具体实施方式
本发明具体实施例如下:
一种惯导/相机安装误差标定方法,包括下述步骤:
(1)系统安装及其连接方式;(2)转台旋转并记录数据;(3)构建12个未知量的方程组并计算;(4)输出结果,完成惯导/相机安装误差标定。
步骤一:系统安装及其连接方式
将惯导与相机安装在转台上,将标定板放置在转台附近,且处于相机的视场范围内。
步骤二:旋转转台并记录数据
系统安装完成后,旋转转台4次,旋转时要求保证:
(1)惯导与相机之间的相对位置与相对姿态不变;
(2)标定板位置保持不变,从而保证转台与标定板之间的相对位置与相对姿态不变;
(3)转台的XYZ三个轴同时旋转一定的角度(三个轴旋角的大小与方向不要求一致);
(4)旋转时要求保证标定板不超出相机的视场范围。
记录4个状态下惯导与转台之间的标定矩阵相机与标定板之间的标定矩阵其中通过光瞄与转台旋转参数获得,通过高精度光学测量方法获得。
步骤三:构建12个未知量的方程组并计算
假设某空间点P由坐标系a中的坐标Pa转换为坐标系b中的坐标Pb,其过程可表示为
公式(1),其中Pa=[xa,ya,za]T为空间点P在坐标系a中的坐标;是一个平移向量,其中tx为由坐标系a向坐标系b平移过程中x方向的平移量,ty为由坐标系a向坐标系b平移过程中y方向的平移量,tz为由坐标系a向坐标系b平移过程中x方向的平移量;是一个旋转矩阵,描述了从坐标系a向坐标系b的旋转过程,即坐标系a按的旋转顺序,分别旋转后,其三轴将与坐标系b平行,
T中的6个参数(α,β,γ,tx,ty,tz)描述了坐标系a与坐标系b之间的相对位置与相对姿态关系,安装误差标定就是要把这6个参数求解出来。
公式(9)可以整理成以下形式
其中
分别称为空间点P在坐标系a向坐标系b下的齐次坐标,中的元素由标定参数计算而来,称之为标定矩阵。
在本发明中,需要用到转台、惯导、相机、标定板四个装置,假设这4个装置分别定义了四个坐标系,分别称之为坐标系base、坐标系imu、坐标系cam、坐标系cal。
将惯导与相机安装在转台上,将标定板放置在转台附近,且处于相机的视场范围内,转台、惯导、相机、标定板四者之间的关系可写成归一化矩阵连乘的形式:
式中:
——4×4矩阵,惯导与转台之间的标定矩阵,通过光瞄与转台旋转参数方法求得,已知量;
——4×4矩阵,惯导与相机之间的标定矩阵,未知量,待求;
——4×4矩阵,相机与标定板之间的标定矩阵,通过光学测量方法求得,已知量;
——4×4矩阵,标定板与转台之间的标定矩阵,未知量。
在上面的连乘公式中,这两个未知的4×4矩阵中可以由共12个未知量表示。
将式(14)写成以下形式:
可以整理出3个以为变量的方程,分别表示为:
因此,得到该状态下12个未知量的3个方程。
转台旋转4次,每次都获得该状态下的3个方程,一共可以获得4个状态下12个未知量的12个方程,表示如下:
采用通用的矩阵分解法、最小二乘法求解以上线性方程组,得到线性方程组的解,获得共12个未知量的值。
步骤四:输出结果
将步骤三中的作为惯导与相机之间的安装误差输出。

Claims (4)

1.一种惯导/相机安装误差标定方法,该方法包括如下步骤:
步骤一、系统安装及其连接方式;步骤二、转台旋转并记录数据;步骤三、构建12个未知量的方程组并计算;步骤四、输出结果;
其特征在于:所述系统安装及其连接方式,将惯导与相机安装在转台上,将标定板放置在转台附近,且处于相机的视场范围内。
2.如权利要求1所述的一种惯导/相机安装误差标定方法,其特征在于,步骤二、旋转转台并记录数据:
系统安装完成后,旋转转台4次,旋转时要求保证:
1、惯导与相机之间的相对位置与相对姿态不变;
2、标定板位置保持不变,从而保证转台与标定板之间的相对位置与相对姿态不变;
3、转台的XYZ三个轴同时旋转一定的角度(三个轴旋角的大小与方向不要求一致);
4、旋转时要求保证标定板不超出相机的视场范围;
记录4个状态下惯导与转台之间的标定矩阵相机与标定板之间的标定矩阵其中通过光瞄与转台旋转参数获得,通过高精度光学测量方法获得。
3.如权利要求1所述的一种惯导/相机安装误差标定方法,其特征在于,步骤三、构建12个未知量的方程组并计算:
假设某空间点P由坐标系a中的坐标Pa转换为坐标系b中的坐标Pb,其过程可表示为,
公式(1),其中Pa=[xa,ya,za]T为空间点P在坐标系a中的坐标;是一个平移向量,其中tx为由坐标系a向坐标系b平移过程中x方向的平移量,ty为由坐标系a向坐标系b平移过程中y方向的平移量,tz为由坐标系a向坐标系b平移过程中x方向的平移量;是一个旋转矩阵,描述了从坐标系a向坐标系b的旋转过程,即坐标系a按的旋转顺序,分别旋转后,其三轴将与坐标系b平行,
T中的6个参数(α,β,γ,tx,ty,tz)描述了坐标系a与坐标系b之间的相对位置与相对姿态关系,安装误差标定就是要把这6个参数求解出来;
公式(1)可以整理成以下形式
其中
分别称为空间点P在坐标系a向坐标系b下的齐次坐标,中的元素由标定参数计算而来,称之为标定矩阵;
将惯导与相机安装在转台上,将标定板放置在转台附近,且处于相机的视场范围内,转台、惯导、相机、标定板四者之间的关系可写成归一化矩阵连乘的形式:
式中:
——4×4矩阵,惯导与转台之间的标定矩阵,通过光瞄与转台旋转参数方法求得,已知量;
——4×4矩阵,惯导与相机之间的标定矩阵,未知量,待求;
——4×4矩阵,相机与标定板之间的标定矩阵,通过光学测量方法求得,已知量;
——4×4矩阵,标定板与转台之间的标定矩阵,未知量;
在上面的连乘公式中,这两个未知的4×4矩阵中可以由共12个未知量表示;
将式(6)写成以下形式:
可以整理出3个以为变量的方程,分别表示为:
因此,得到该状态下12个未知量的3个方程;
转台旋转4次,每次都获得该状态下的3个方程,一共可以获得4个状态下12个未知量的12个方程,表示如下:
采用通用的矩阵分解法、最小二乘法求解以上线性方程组,得到线性方程组的解,获得共12个未知量的值。
4.如权利要求1所述的一种惯导/相机安装误差标定方法,其特征在于,步骤四、输出结果,将步骤三中的作为惯导与相机之间的安装误差输出。
CN201611150018.5A 2016-12-14 2016-12-14 一种惯导/相机安装误差标定方法 Active CN108225371B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201611150018.5A CN108225371B (zh) 2016-12-14 2016-12-14 一种惯导/相机安装误差标定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201611150018.5A CN108225371B (zh) 2016-12-14 2016-12-14 一种惯导/相机安装误差标定方法

Publications (2)

Publication Number Publication Date
CN108225371A true CN108225371A (zh) 2018-06-29
CN108225371B CN108225371B (zh) 2021-07-13

Family

ID=62638144

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201611150018.5A Active CN108225371B (zh) 2016-12-14 2016-12-14 一种惯导/相机安装误差标定方法

Country Status (1)

Country Link
CN (1) CN108225371B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109238309A (zh) * 2018-10-31 2019-01-18 中航航空服务保障(天津)有限公司 一种捷联航姿安装板校准装置及其校准方法
CN110132309A (zh) * 2019-06-05 2019-08-16 西京学院 一种采煤机摇臂惯性/视觉组合定姿装置标定方法
CN111486867A (zh) * 2020-03-19 2020-08-04 天津大学 一种视觉和惯性混合跟踪组件安装参数的标定装置及方法
CN112611361A (zh) * 2020-12-08 2021-04-06 华南理工大学 一种无人机机载测绘吊舱的摄像头安装误差测定方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102435188A (zh) * 2011-09-15 2012-05-02 南京航空航天大学 一种用于室内环境的单目视觉/惯性全自主导航方法
CN103630138A (zh) * 2013-12-09 2014-03-12 天津工业大学 基于摄像机头标定算法的无人机视觉导航方法
CN105067011A (zh) * 2015-09-15 2015-11-18 沈阳飞机工业(集团)有限公司 一种基于视觉标定及坐标转换的测量系统整体校准方法
US20160005164A1 (en) * 2013-02-21 2016-01-07 Regents Of The University Of Minnesota Extrinsic parameter calibration of a vision-aided inertial navigation system
CN105509765A (zh) * 2014-09-23 2016-04-20 北京自动化控制设备研究所 一种惯性/dvl/usbl安装误差标定方法
CN105606127A (zh) * 2016-01-11 2016-05-25 北京邮电大学 一种双目立体相机与惯性测量单元相对姿态标定方法
CN106153041A (zh) * 2015-03-30 2016-11-23 北京自动化控制设备研究所 一种基于多景深信息的视觉里程计测速方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102435188A (zh) * 2011-09-15 2012-05-02 南京航空航天大学 一种用于室内环境的单目视觉/惯性全自主导航方法
US20160005164A1 (en) * 2013-02-21 2016-01-07 Regents Of The University Of Minnesota Extrinsic parameter calibration of a vision-aided inertial navigation system
CN103630138A (zh) * 2013-12-09 2014-03-12 天津工业大学 基于摄像机头标定算法的无人机视觉导航方法
CN105509765A (zh) * 2014-09-23 2016-04-20 北京自动化控制设备研究所 一种惯性/dvl/usbl安装误差标定方法
CN106153041A (zh) * 2015-03-30 2016-11-23 北京自动化控制设备研究所 一种基于多景深信息的视觉里程计测速方法
CN105067011A (zh) * 2015-09-15 2015-11-18 沈阳飞机工业(集团)有限公司 一种基于视觉标定及坐标转换的测量系统整体校准方法
CN105606127A (zh) * 2016-01-11 2016-05-25 北京邮电大学 一种双目立体相机与惯性测量单元相对姿态标定方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
KELLY.J ET AL: "Fast Relative Pose Calibration for Visual and Inertial Sensors", 《11TH INTERNATIONAL SYMPOSIUM ON EXPERIMENTAL ROBOTICS (ISER)》 *
刘畅等: "一种基于惯性/视觉信息融合的无人机自主着陆导航算法", 《导航定时与授时》 *
陈林等: "INS/Vision组合导航中视觉系统动态定位方法研究", 《传感技术学报》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109238309A (zh) * 2018-10-31 2019-01-18 中航航空服务保障(天津)有限公司 一种捷联航姿安装板校准装置及其校准方法
CN109238309B (zh) * 2018-10-31 2023-12-26 中航航空服务保障(天津)有限公司 一种捷联航姿安装板校准装置及其校准方法
CN110132309A (zh) * 2019-06-05 2019-08-16 西京学院 一种采煤机摇臂惯性/视觉组合定姿装置标定方法
CN110132309B (zh) * 2019-06-05 2023-04-25 西京学院 一种采煤机摇臂惯性/视觉组合定姿装置标定方法
CN111486867A (zh) * 2020-03-19 2020-08-04 天津大学 一种视觉和惯性混合跟踪组件安装参数的标定装置及方法
CN111486867B (zh) * 2020-03-19 2023-08-11 天津大学 一种视觉和惯性混合跟踪组件安装参数的标定装置及方法
CN112611361A (zh) * 2020-12-08 2021-04-06 华南理工大学 一种无人机机载测绘吊舱的摄像头安装误差测定方法

Also Published As

Publication number Publication date
CN108225371B (zh) 2021-07-13

Similar Documents

Publication Publication Date Title
CN109696183B (zh) 惯性测量单元的标定方法及装置
WO2021180128A1 (zh) 一种基于rgbd传感器和imu传感器的定位方法
CN109813336B (zh) 惯性测量单元标定方法
CN105518486B (zh) 追踪可移动目标物体的方位的系统和方法
CN106052584B (zh) 一种基于视觉及惯性信息融合的轨道空间线形测量方法
CN110057352A (zh) 一种相机姿态角确定方法及装置
CN108225371A (zh) 一种惯导/相机安装误差标定方法
CN110987021B (zh) 一种基于转台基准的惯性视觉相对姿态标定方法
CN109460052A (zh) 一种可拼组飞行器的控制方法
CN103527620A (zh) 可实现回转角度测量的球铰链及测量方法
CN105241478B (zh) 一种单轴调制激光陀螺惯性导航系统多坐标系标定方法
CN108253966A (zh) 无人机飞行三维模拟显示方法
CN106643489A (zh) 一种基于激光跟踪仪的地面二维码标定方法及装置
CN110728716B (zh) 一种标定方法、装置及飞行器
CN107782309A (zh) 非惯性系视觉和双陀螺仪多速率ckf融合姿态测量方法
CN113267794B (zh) 一种基线长度约束的天线相位中心校正方法及装置
Grivon et al. Development of an innovative low-cost MARG sensors alignment and distortion compensation methodology for 3D scanning applications
CN106017452A (zh) 双陀螺抗扰动寻北方法
CN109093620B (zh) 一种双目相机辅助的空间非合作目标动力学参数辨识方法
CN115371665A (zh) 一种基于深度相机和惯性融合的移动机器人定位方法
Xian et al. Fusing stereo camera and low-cost inertial measurement unit for autonomous navigation in a tightly-coupled approach
MacArthur et al. Unmanned ground vehicle state estimation using an unmanned air vehicle
CN112489118B (zh) 一种无人机机载传感器组外参快速标定方法
CN105928513B (zh) 一种基于位置姿态测量系统的机载合成孔径雷达运动参数测量方法
CN103995239A (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