CN109242910B - 一种基于任意已知平面形状的单目相机自标定方法 - Google Patents
一种基于任意已知平面形状的单目相机自标定方法 Download PDFInfo
- Publication number
- CN109242910B CN109242910B CN201810952873.0A CN201810952873A CN109242910B CN 109242910 B CN109242910 B CN 109242910B CN 201810952873 A CN201810952873 A CN 201810952873A CN 109242910 B CN109242910 B CN 109242910B
- Authority
- CN
- China
- Prior art keywords
- matrix
- sub
- image
- pixel
- plane shape
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 37
- 230000009466 transformation Effects 0.000 claims abstract description 20
- 239000011159 matrix material Substances 0.000 claims description 62
- 239000013598 vector Substances 0.000 claims description 28
- 125000004432 carbon atom Chemical group C* 0.000 claims description 9
- 238000007476 Maximum Likelihood Methods 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 6
- 238000000605 extraction Methods 0.000 claims description 6
- 238000006243 chemical reaction Methods 0.000 claims description 5
- 238000001514 detection method Methods 0.000 claims description 3
- 238000003708 edge detection Methods 0.000 claims description 3
- 230000033001 locomotion Effects 0.000 abstract description 7
- 230000001186 cumulative effect Effects 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 238000003754 machining Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/80—Analysis of captured images to determine intrinsic or extrinsic camera parameters, i.e. camera calibration
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/13—Edge detection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20016—Hierarchical, coarse-to-fine, multiscale or multiresolution image processing; Pyramid transform
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20092—Interactive image processing based on input by user
- G06T2207/20104—Interactive definition of region of interest [ROI]
Landscapes
- Engineering & Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Image Analysis (AREA)
- Studio Devices (AREA)
Abstract
本发明公开了一种基于任意已知平面形状的单目相机自标定方法。本发明充分利用目标中已知平面形状信息,对目标与相机间位姿关系无任何要求,适应尺度、旋转及射影变换,且对遮挡、杂波具有一定的鲁棒性;无需额外标定物及相机特殊运动,大大降低相机标定过程的复杂度;解决了自标定结果实时性、准确性差的技术问题。
Description
技术领域
本发明涉及相机标定技术领域,具体涉及一种基于任意已知平面形状的单目相机自标定方法。
背景技术
在机器视觉测量领域,为确定空间物体的三维几何位置与其在图像中对应点之间的转换关系,必须建立相机成像的几何模型,而这一转换关系由相机成像的几何模型决定。相机标定就是求解这些几何模型参数的过程,其标定结果的精度及算法的稳定性直接影响相机工作产生结果的准确性。因此,相机标定在机器视觉应用中十分关键,是后续测量工作的前提,改善标定方法、提高标定精度是机器视觉研究工作的重点所在。
相机标定技术大致上可以分为三类:传统的摄影测量学中的标定方法,主动视觉标定方法和自标定方法。传统标定方法需要使用经过精密加工的标定块,通过建立标定块上三维坐标已知的点与其图像点间的对应,来计算摄像机的内外参数。该方法的优点在于可以获得较高的精度,但标定过程费时费力,不适用于在线标定和不可能使用标定块的场合;基于主动视觉的标定方法需要控制摄像机做某些特殊运动,如绕光心旋转或纯平移等,利用这种运动的特殊性可以计算出内参数。该方法优点在于算法简单,往往能获得线性解,缺点是对摄像机做特殊运动时的精度要求高,且不能适用于摄像机运动未知或无法控制的场合;相机自标定方法利用摄像机内参数自身存在的约束实现标定,通过在某一静态场景对目标物体进行多次拍摄,或者多个相机对目标同时拍摄,利用图像序列之间的相互约束关系标定,灵活性强,有很广发的应用范围,但利用摄像机内参数自身约束的自标定方法标定过程复杂,不适用于实时性较强的场合,且由于自标定方法采用非线性标定,依赖良好的初值估计,鲁棒性不足,只适用于对精度要求不高的场合。
发明内容
针对现有技术中的上述不足,本发明提供的一种基于任意已知平面形状的单目相机自标定方法解决了传统标定技术中摄像机做特殊运动时操作复杂、精度要求高,以及自标定方法实时性差、标定结果精度差的问题。
为了达到上述发明目的,本发明采用的技术方案为:一种基于任意已知平面形状的单目相机自标定方法,包括以下步骤:
S1、基于目标中已知平面形状信息建立模板集合;
S2、通过单目相机获取目标图像,所述图像包含模板集合中的已知平面形状信息;
S3、在模板集合中搜索与目标图像中相匹配的模型,得到该模型的单应性矩阵;
S4、通过亚像素角点检测以及最小二乘拟合将单应性矩阵精确到亚像素级别;
S5、根据多幅图像的亚像素级别下单应性矩阵计算单目相机的内参矩阵及外参矩阵;
S6、对内参矩阵和外参矩阵进行极大似然估计,得到畸变系数和内参矩阵的最优解,完成单目相机的自标定。
进一步地:所述步骤S1具体包括:
S11、获取包含已知平面形状信息的ROI图像;
S12、对ROI图像生成ROI图像高斯金字塔;
S13、通过射影变换参数对ROI图像的每级金字塔进行转换得到图像集合,转换公式为:
S14、对图像集合中的图像通过canny算子检测边缘方法提取边缘特征,产生边缘像素及其方向向量,添加到模板集合中。
进一步地:所述步骤S3具体包括:
S31、将目标图像生成目标图像高斯金字塔;
S32、计算模型与目标图像金字塔最低分辨率层级像素梯度的匹配度,计算公式为:
S33、当匹配度大于阈值时,得到该模型的射影变换参数及单应性矩阵。
进一步地:所述步骤S4具体包括:
S41、利用Shi-Tomasi角点提取算法提取模型中的角点,并根据提取将其迭代至亚像素精度,得到亚像素坐标pj;
上式中,min为拟合后的最小值;
S44、根据亚像素级别射影变换参数A、b、C计算其单应性矩阵。
进一步地:所述步骤S5具体为:
S51、根据相机内参矩阵的约束条件计算内参矩阵B,约束条件为:
S52、根据内参矩阵B及单应性矩阵计算外参矩阵,计算公式为:
B[R1 R2 t]=H
上式中,[R1 R2 t]为外参矩阵。
进一步地:所述步骤S6中极大似然估计的公式为:
本发明的有益效果为:本发明充分利用目标中已知平面形状信息,对目标与相机间位姿关系无任何要求,适应尺度、旋转及射影变换,且对遮挡、杂波具有一定的鲁棒性;无需额外标定物及相机特殊运动,大大降低相机标定过程的复杂度;解决了自标定结果实时性、准确性差的技术问题。
附图说明
图1为本发明流程图。
具体实施方式
下面对本发明的具体实施方式进行描述,以便于本技术领域的技术人员理解本发明,但应该清楚,本发明不限于具体实施方式的范围,对本技术领域的普通技术人员来讲,只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内,这些变化是显而易见的,一切利用本发明构思的发明创造均在保护之列。
如图1所示,一种基于任意已知平面形状的单目相机自标定方法,包括以下步骤:
S1、基于目标中已知平面形状信息建立模板集合,具体包括:
S11、获取包含已知平面形状信息的ROI图像;
S12、对ROI图像生成ROI图像高斯金字塔,本发明实例中通过高斯平滑和亚采样获得lmax幅临时下采样图像,与原图像构成lmax+1层高斯图像金字塔;
S13、通过射影变换参数对ROI图像的每级金字塔进行转换得到图像集合,转换公式为:
S14、对图像集合中的图像通过canny算子检测边缘方法提取边缘特征,产生边缘像素及其方向向量,添加到模板集合中。
对各像素点的梯度值进行非极大值抑制,用梯度直方图的形式描述所有边缘标志不为0的像素,并计算其梯度累计直方图;选取梯度累计占比值大于0.95处的梯度值作为高阈值h2,低阈值h1=0.4×h2,由双阈值h1,h2分别得到强边缘像素和弱边缘像素,最后以强边缘像素点为中心点,比较领域边缘点和弱边缘像素中相应位置的边缘信息,进行强弱边缘的连接,得到边缘特征点集,其梯度作为该像素的方向向量。
S2、通过单目相机获取目标图像,所述图像包含模板集合中的已知平面形状信息。
S3、在模板集合中搜索与目标图像中相匹配的模型,得到该模型的单应性矩阵,具体包括:
S31、将目标图像生成目标图像高斯金字塔,并计算目标图像高斯金字塔的最低分辨率层级的像素梯度;
S32、计算模型与目标图像金字塔最低分辨率层级像素梯度的匹配度,计算公式为:
S33、当匹配度大于阈值时,得到该模型的射影变换参数及单应性矩阵。
S4、通过亚像素角点检测以及最小二乘拟合将单应性矩阵精确到亚像素级别,具体包括:
S41、利用Shi-Tomasi角点提取算法提取模型中的角点,并根据提取将其迭代至亚像素精度,得到亚像素坐标pj;
上式中,min为拟合后的最小值。
S44、根据亚像素级别射影变换参数A、b、C计算其单应性矩阵。
S5、根据多幅图像的亚像素级别下单应性矩阵计算单目相机的内参矩阵及外参矩阵,具体为:
S51、根据相机内参矩阵的约束条件计算内参矩阵B,约束条件为:
S52、根据内参矩阵B及单应性矩阵计算外参矩阵,计算公式为:
B[R1 R2 t]=H
上式中,[R1 R2 t]为外参矩阵。
S6、对内参矩阵和外参矩阵进行极大似然估计,得到畸变系数和内参矩阵的最优解,完成单目相机的自标定,极大似然估计的公式为:
Claims (6)
1.一种基于任意已知平面形状的单目相机自标定方法,其特征在于,包括以下步骤:
S1、基于目标中已知平面形状信息建立模板集合;
S2、通过单目相机获取目标图像,所述图像包含模板集合中的已知平面形状信息;
S3、在模板集合中搜索与目标图像中相匹配的模型,得到该模型的单应性矩阵;
S4、通过亚像素角点检测以及最小二乘拟合将单应性矩阵精确到亚像素级别;
S5、根据多幅图像的亚像素级别下单应性矩阵计算单目相机的内参矩阵及外参矩阵;
S6、对内参矩阵和外参矩阵进行极大似然估计,得到畸变系数和内参矩阵的最优解,作为单目相机的自标定结果。
2.根据权利要求1所述的基于任意已知平面形状的单目相机自标定方法,其特征在于,所述步骤S1具体包括:
S11、获取包含已知平面形状信息的ROI图像;
S12、对ROI图像生成ROI图像高斯金字塔;
S13、通过射影变换参数对ROI图像的每级金字塔进行转换得到图像集合,转换公式为:
S14、对图像集合中的图像通过canny算子检测边缘方法提取边缘特征,产生边缘像素及其方向向量,添加到模板集合中。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810952873.0A CN109242910B (zh) | 2018-08-21 | 2018-08-21 | 一种基于任意已知平面形状的单目相机自标定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810952873.0A CN109242910B (zh) | 2018-08-21 | 2018-08-21 | 一种基于任意已知平面形状的单目相机自标定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109242910A CN109242910A (zh) | 2019-01-18 |
CN109242910B true CN109242910B (zh) | 2021-10-08 |
Family
ID=65071643
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810952873.0A Expired - Fee Related CN109242910B (zh) | 2018-08-21 | 2018-08-21 | 一种基于任意已知平面形状的单目相机自标定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109242910B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110942484B (zh) * | 2019-11-26 | 2022-07-12 | 福州大学 | 基于遮挡感知和特征金字塔匹配的相机自运动估计方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2256576A1 (en) * | 2009-05-29 | 2010-12-01 | Datacard Corporation | System for the calibration of operating parameters of a laser engraver |
CN102651145A (zh) * | 2012-04-06 | 2012-08-29 | 哈尔滨工业大学 | 股骨三维模型可视化方法 |
CN103177439A (zh) * | 2012-11-26 | 2013-06-26 | 惠州华阳通用电子有限公司 | 一种基于黑白格角点匹配的自动标定方法 |
CN103440639A (zh) * | 2013-07-15 | 2013-12-11 | 江南大学 | 基于正三角形模板的摄像机内参数标定方法 |
CN103606149A (zh) * | 2013-11-14 | 2014-02-26 | 深圳先进技术研究院 | 一种双目摄像机的标定方法、装置及双目摄像机 |
CN107133986A (zh) * | 2017-04-26 | 2017-09-05 | 武汉科技大学 | 一种基于二维标定物的相机标定方法 |
CN107633536A (zh) * | 2017-08-09 | 2018-01-26 | 武汉科技大学 | 一种基于二维平面模板的相机标定方法及系统 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8798314B2 (en) * | 2008-07-14 | 2014-08-05 | National Ict Australia Limited | Detection of vehicles in images of a night time scene |
-
2018
- 2018-08-21 CN CN201810952873.0A patent/CN109242910B/zh not_active Expired - Fee Related
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2256576A1 (en) * | 2009-05-29 | 2010-12-01 | Datacard Corporation | System for the calibration of operating parameters of a laser engraver |
CN102651145A (zh) * | 2012-04-06 | 2012-08-29 | 哈尔滨工业大学 | 股骨三维模型可视化方法 |
CN103177439A (zh) * | 2012-11-26 | 2013-06-26 | 惠州华阳通用电子有限公司 | 一种基于黑白格角点匹配的自动标定方法 |
CN103440639A (zh) * | 2013-07-15 | 2013-12-11 | 江南大学 | 基于正三角形模板的摄像机内参数标定方法 |
CN103606149A (zh) * | 2013-11-14 | 2014-02-26 | 深圳先进技术研究院 | 一种双目摄像机的标定方法、装置及双目摄像机 |
CN107133986A (zh) * | 2017-04-26 | 2017-09-05 | 武汉科技大学 | 一种基于二维标定物的相机标定方法 |
CN107633536A (zh) * | 2017-08-09 | 2018-01-26 | 武汉科技大学 | 一种基于二维平面模板的相机标定方法及系统 |
Non-Patent Citations (3)
Title |
---|
View-Invariant Template Matching Using Homography Constraints;Sina Lotfian 等;《Computer Vision and Pattern Recognition》;20170512;1-8 * |
基于合作标定物的相机标定方法研究;游迪;《中国优秀博硕士学位论文全文数据库(硕士)信息科技辑》;20160815(第08期);I138-990 * |
多摄相机标定软件系统设计与实现;李现刚;《中国优秀博硕士学位论文全文数据库(硕士)信息科技辑》;20160215(第02期);I138-1697 * |
Also Published As
Publication number | Publication date |
---|---|
CN109242910A (zh) | 2019-01-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109345588B (zh) | 一种基于Tag的六自由度姿态估计方法 | |
CN109146980B (zh) | 基于单目视觉的优化的深度提取和被动测距方法 | |
CN111311679B (zh) | 一种基于深度相机的自由漂浮目标位姿估计方法 | |
CN110146099B (zh) | 一种基于深度学习的同步定位与地图构建方法 | |
CN105021124B (zh) | 一种基于深度图的平面零件三维位置和法向量计算方法 | |
CN108122256B (zh) | 一种逼近状态下旋转目标位姿测量的方法 | |
Koide et al. | General, single-shot, target-less, and automatic lidar-camera extrinsic calibration toolbox | |
CN112132874B (zh) | 无标定板异源图像配准方法、装置、电子设备及存储介质 | |
CN112200203B (zh) | 斜视场下弱相关散斑图像的匹配方法 | |
CN111897349A (zh) | 一种基于双目视觉的水下机器人自主避障方法 | |
Li et al. | Fast visual odometry using intensity-assisted iterative closest point | |
CN109493385A (zh) | 一种结合场景点线特征的移动机器人室内自主定位方法 | |
CN116912333B (zh) | 一种基于作业围栏标定杆的相机姿态自标定方法 | |
CN106952262B (zh) | 一种基于立体视觉的船板加工精度分析方法 | |
CN114022560A (zh) | 标定方法及相关装置、设备 | |
CN113642397B (zh) | 一种基于手机视频的物体长度测量方法 | |
Han et al. | Target positioning method in binocular vision manipulator control based on improved canny operator | |
CN108596947A (zh) | 一种适用于rgb-d相机的快速目标跟踪方法 | |
CN112767481B (zh) | 一种基于视觉边缘特征的高精度定位及建图方法 | |
CN109242910B (zh) | 一种基于任意已知平面形状的单目相机自标定方法 | |
CN110570473A (zh) | 一种基于点线融合的权重自适应位姿估计方法 | |
Alkaabi et al. | Iterative corner extraction and matching for mosaic construction | |
Chen et al. | High precision pose measurement for humanoid robot based on PnP and OI algorithms | |
Iida et al. | High-accuracy Range Image Generation by Fusing Binocular and Motion Stereo Using Fisheye Stereo Camera | |
Comport et al. | Efficient model-based tracking for robot vision |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20211008 |