CN102005048A - 基于平面上固定特征点的大视场摄像机现场标定方法 - Google Patents
基于平面上固定特征点的大视场摄像机现场标定方法 Download PDFInfo
- Publication number
- CN102005048A CN102005048A CN2010105460696A CN201010546069A CN102005048A CN 102005048 A CN102005048 A CN 102005048A CN 2010105460696 A CN2010105460696 A CN 2010105460696A CN 201010546069 A CN201010546069 A CN 201010546069A CN 102005048 A CN102005048 A CN 102005048A
- Authority
- CN
- China
- Prior art keywords
- point
- calibration
- points
- unique points
- plane
- 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
- Image Processing (AREA)
Abstract
基于平面上固定特征点的大视场摄像机现场标定方法,本发明涉及一种摄像机标定方法。用于数据标定。解决了传统的标定加工难度大、加工精度不易保证、模板制作成本高和使用受限制的问题。它包括如下步骤:利用标定平面上位于两两相交的三条直线交点处的特征点,分别通过任意两个特征点再做射线,从而把每两条射线的交点处作为扩展特征点;求取从标定平面上摄录到摄像机中的图像因镜头所造成的实物坐标系与图像坐标系之间的一阶径向畸变系数;利用获得的一阶径向畸变系数修正在图像坐标系中的坐标步骤一构造的九个特征点的图像坐标及扩展标定点的图像坐标;把得到的九个特征点及扩展标定点的图像坐标作为标定点进行摄像机标定。
Description
技术领域
本发明涉及一种摄像机标定方法。
背景技术
准确的摄像机标定是精确测量目标物体位姿的必要前提。进行摄像机标定时,通常希望标定点能布满整个标定平面,传统的标定模板一般都要设置几十个甚至上百个的高精度标定点,这在工程实际中将遇到以下问题:(1)加工难度大、加工精度不易保证:由于加工的不确定性,对大量特征点的加工一致性和稳定性不易保证,从而难以保证整个标定模板的精度,同时加工大量的标定点也增加了加工难度;(2)模板制作成本高:为了在模板中放置大量的特征点,需要扩大模板的尺寸或标定点的密度,增加了模板的制作复杂程度,进而增加了制作成本;(3)使用受限制:扩大模板尺寸,不便于放置标定模板,而增加标定点的密度,在摄像机分辨率较低或角度较大时,又容易出现分辨不清的情况,使得特征点不易提取,不利于标定过程的自动化。对于大视场测量系统,视场宽度通常能达到十几米甚至几十米,在这么大的视场中安装几十个甚至上百个的标定点,不但难度大,而且精度不易保证,整个标定过程费时费力。在大视场标定中,通常使用激光跟踪仪在幕墙上设置标定点,从而大大增加了摄像机标定的成本。
发明内容
本发明的目的是提供一种基于平面上固定特征点的大视场摄像机现场标定方法,以解决传统的标定模板需要设置几十个甚至上百个的高精度标定点,在工程实际中将遇到加工难度大、加工精度不易保证、模板制作成本高和使用受限制的问题。它包括如下步骤:
一、利用标定平面上位于两两相交的三条直线交点处的一号特征点1、二号特征点2、三号特征点3、四号特征点4、五号特征点5、六号特征点6、七号特征点7、八号特征点8和九号特征点9,分别通过任意两个特征点再做射线,从而把每两条射线的交点处作为扩展特征点;上述九个特征点的坐标为已知;
二、求取从标定平面上摄录到摄像机中的图像因镜头所造成的实物坐标系与图像坐标系之间的一阶径向畸变系数k1;
三、利用步骤二获得的一阶径向畸变系数k1修正在图像坐标系中的坐标步骤一构造的九个特征点的图像坐标及扩展标定点的图像坐标;
四、把步骤三得到的九个特征点的图像坐标及扩展标定点的图像坐标作为标定点进行摄像机标定。
本发明的标定方法考虑了摄像机畸变对特征点构造的影响,使用九个满足特定几何约束(即两两相交的三条直线)的原始特征点,标定结果接近经典的Tsai两步法,在物距为13m的条件下,重构精度在3mm以内,可满足较高的精度要求。本发明可以大大降低大视场标定的复杂度,降低标定成本,提高标定效率和可靠性,降低了模板制作成本,使用上不受限制,具有较高的应用价值。
附图说明
图1是本发明方法第一步骤的示意图,虚线的交点为扩展特征点;图2为实验中使用87个特征点的传统Tsai方法标定的结果示意图;图3为实验中仅使用九个标定点的传统Tsai方法标定的结果示意图;图4为实验中使用本发明方法得到的标定结果示意图,图5为实施方式一中计算一阶径向畸变系数
具体实施方式
具体实施方式一:下面结合图1至4具体说明本实施方式。本实施方式包括如下步骤:
一、利用标定平面上位于两两相交的三条直线交点处的一号特征点1、二号特征点2、三号特征点3、四号特征点4、五号特征点5、六号特征点6、七号特征点7、八号特征点8和九号特征点9,分别通过任意两个特征点再做射线,从而把每两条射线的交点处作为扩展特征点;上述九个特征点的坐标为已知;
二、求取从标定平面上摄录到摄像机中的图像因镜头所造成的实物坐标系与图像坐标系之间的一阶径向畸变系数k1;
三、利用步骤二获得的一阶径向畸变系数k1修正在图像坐标系中的坐标步骤一构造的九个特征点的图像坐标及扩展标定点的图像坐标;
四、把步骤三得到的九个特征点的图像坐标及扩展标定点的图像坐标作为标定点进行摄像机标定。
性能实验结果如下:图2为实际实验中使用87个特征点的传统Tsai方法标定的结果;图3为实际实验中使用9个标定点的传统Tsai方法标定的结果;图4为实际实验中使用本实施方式方法得到的标定结果。
从上述实验结果可以得到如下结论:使用发明所述的方法标定,在物距为13m的条件下,重构误差在3mm以内,与使用245个点的Tsai方法标定结果接近,并且与只使用9个点的Tsai方法标定的结果相比,重构精度显著提高。本发明使用较少的特征点,取得了接近经典两步法的标定精度,降低了标定模板的加工难度,可以应用于对精度有较高要求的场合。特别是在大视场这种物理空间较大,不适宜或不方便放置较多特征点的场合,使用本方法可以大大降低靶标加工安置的难度,降低标定成本。
具体实施方式二:下面结合图5具体说明本实施方式。本实施方式与实施方式一的不同在于具体限定了步骤二中:计算步骤一的特征点中所有共线四点的交比,应用交比不变性原理求取一阶径向畸变系数k1;具体过程如下:
设空间直线L上的四个点Pw1、Pw2、Pw3、Pw4的交比为(Pw1,Pw2;Pw3,Pw4),其在L1和L2上的对应点分别为P1、P2、P3、P4和P`1、P`2、P`3、P`4,由交比不变性原理得
(Pw1,Pw2;Pw3,Pw4)=(P1,P2;P3,P4) (1)
其中,P1、P2、P3、P4的交比定义为
又,L1上的点Pi与L2上点P′i的对应关系为
将式(1)和(2)联立方程组即可求出一阶径向畸变系数k1。如果已知多组共线点坐标及其交比,利用上面的方程组可以求得一阶径向畸变系数k1的最小二乘解。
Claims (3)
1.基于平面上固定特征点的大视场摄像机现场标定方法,其特征在于它包括如下步骤:
一、利用标定平面上位于两两相交的三条直线交点处的一号特征点(1)、二号特征点(2)、三号特征点(3)、四号特征点(4)、五号特征点(5)、六号特征点(6)、七号特征点(7)、八号特征点(8)和九号特征点(9),分别通过任意两个特征点再做射线,从而把每两条射线的交点处作为扩展特征点;上述九个特征点的坐标为已知;
二、求取从标定平面上摄录到摄像机中的图像因镜头所造成的实物坐标系与图像坐标系之间的一阶径向畸变系数(k1);
三、利用步骤二获得的一阶径向畸变系数(k1)修正在图像坐标系中的坐标步骤一构造的九个特征点的图像坐标及扩展标定点的图像坐标;
四、把步骤三得到的九个特征点的图像坐标及扩展标定点的图像坐标作为标定点进行摄像机标定。
2.根据权利要求1所述的基于平面上固定特征点的大视场摄像机现场标定方法,其特征在于计算步骤一的特征点中所有共线四点的交比,应用交比不变性原理求取一阶径向畸变系数(k1)。
3.根据权利要求2所述的基于平面上固定特征点的大视场摄像机现场标定方法,其特征在于求取一阶径向畸变系数(k1)的具体过程如下:
设空间直线L上的四个点Pw1、Pw2、Pw3、Pw4的交比为(Pw1,Pw2;Pw3,Pw4),其在L1和L2上的对应点分别为P1、P2、P3、P4和P`1、P`2、P`3、P`4,由交比不变性原理得
(Pw1,Pw2;Pw3,Pw4)=(P1,P2;P3,P4) (1)
其中,P1、P2、P3、P4的交比定义为
又,L1上的点Pi与L2上点P′i的对应关系为
将式(1)和(2)联立方程组即可求出一阶径向畸变系数(k1);如果已知多组共线点坐标及其交比,利用上面的方程组可以求得一阶径向畸变系数(k1)的最小二乘解。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010105460696A CN102005048A (zh) | 2010-11-16 | 2010-11-16 | 基于平面上固定特征点的大视场摄像机现场标定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010105460696A CN102005048A (zh) | 2010-11-16 | 2010-11-16 | 基于平面上固定特征点的大视场摄像机现场标定方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN102005048A true CN102005048A (zh) | 2011-04-06 |
Family
ID=43812387
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2010105460696A Pending CN102005048A (zh) | 2010-11-16 | 2010-11-16 | 基于平面上固定特征点的大视场摄像机现场标定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102005048A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108921899A (zh) * | 2018-07-13 | 2018-11-30 | 哈尔滨工业大学 | 一种基于像素阈值求解基本矩阵的室内视觉定位方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1801953A (zh) * | 2005-12-15 | 2006-07-12 | 上海交通大学 | 仅用一幅平面标定物图像的摄像机标定方法 |
CN101419709A (zh) * | 2008-12-08 | 2009-04-29 | 北京航空航天大学 | 一种用于摄像机标定的平面靶标特征点的自动匹配方法 |
-
2010
- 2010-11-16 CN CN2010105460696A patent/CN102005048A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1801953A (zh) * | 2005-12-15 | 2006-07-12 | 上海交通大学 | 仅用一幅平面标定物图像的摄像机标定方法 |
CN101419709A (zh) * | 2008-12-08 | 2009-04-29 | 北京航空航天大学 | 一种用于摄像机标定的平面靶标特征点的自动匹配方法 |
Non-Patent Citations (5)
Title |
---|
《IEEE JOURNAL OF ROBOTICS AND AUTOMATION》 19870831 ROGER Y. TSAI A Versatile Camera Calibration Techniaue for High-Accuracy 3D Machine Vision Metrology Using Off-the-shelf TV Cameras and Lenses 第323~344页 1-3 第RA-3卷, 第4期 2 * |
《IEEE TRANSACTIONS ON PATTERN ANALYSIS AND MACHINE INTELLIGENCE》 20001130 Zhengyou Zhang A Flexible New Technique for Camera Calibration 第1330~1334页 1-3 第22卷, 第11期 2 * |
《仪器仪表学报》 20041031 贺俊吉,等 基于交比不变性的镜头畸变参数标定方法 第597~599页 1-2 第25卷, 第5期 2 * |
《北京航空航天大学学报》 20090331 余明浪,等 基于柔性平面靶标的摄像机现场标定方法 第347~350页 1-2 第35卷, 第3期 2 * |
《工业仪表与自动化装置》 20071231 张磊,等 一种考虑畸变的摄像机线性标定简化方法 第3~10页 1-2 , 第1期 2 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108921899A (zh) * | 2018-07-13 | 2018-11-30 | 哈尔滨工业大学 | 一种基于像素阈值求解基本矩阵的室内视觉定位方法 |
CN108921899B (zh) * | 2018-07-13 | 2021-05-28 | 哈尔滨工业大学 | 一种基于像素阈值求解基本矩阵的室内视觉定位方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110969663B (zh) | 一种相机外部参数的静态标定方法 | |
CN103353388B (zh) | 一种具摄像功能的双目体式显微成像系统标定方法及装置 | |
CN100432628C (zh) | 太阳敏感器测量基准转换方法与装置 | |
CN105160702A (zh) | 基于LiDAR点云辅助的立体影像密集匹配方法及系统 | |
CN112819903A (zh) | 基于l型标定板的相机和激光雷达联合标定的方法 | |
CN106354094A (zh) | 基于空间标准球的机床随动激光扫描坐标标定方法 | |
CN107044852B (zh) | 不平整状态下全站仪测量方法 | |
CN107610178A (zh) | 一种工业摄影测量系统相机参数移动式标定方法 | |
CN105021124A (zh) | 一种基于深度图的平面零件三维位置和法向量计算方法 | |
CN104406607A (zh) | 一种多视场复合光学敏感器的标定装置及方法 | |
CN103530880A (zh) | 基于投影高斯网格图案的摄像机标定方法 | |
CN102693543B (zh) | 室外环境下变焦云台摄像机全自动标定方法 | |
CN104050650B (zh) | 基于坐标变换的完整成像的图像拼接方法 | |
CN106248014A (zh) | 一种基于单相片的三维坐标测量方法及装置 | |
CN101118158A (zh) | 离散点云整合系统及方法 | |
CN103900504A (zh) | 纳米尺度下的实时三维视觉信息反馈方法 | |
CN113379668B (zh) | 光伏板拼接方法及装置、电子设备及存储介质 | |
CN104615880A (zh) | 一种三维激光雷达点云匹配的快速icp方法 | |
CN108036791A (zh) | 一种部件间高精度自动对接的位姿检测方法 | |
CN108845663A (zh) | 设计施工系统 | |
CN106405581A (zh) | 卫星结构变形引起的多类型载荷间耦合指向精度评估方法 | |
CN103424087A (zh) | 一种大尺度钢板三维测量拼接系统及方法 | |
CN103413319A (zh) | 一种工业摄像机参数现场标定方法 | |
Chen et al. | An innovative detection method of high-speed railway track slab supporting block plane based on point cloud data from 3D scanning technology | |
CN102034234B (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 | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20110406 |