CN113920206A - 透视移轴相机的标定方法 - Google Patents
透视移轴相机的标定方法 Download PDFInfo
- Publication number
- CN113920206A CN113920206A CN202111528806.4A CN202111528806A CN113920206A CN 113920206 A CN113920206 A CN 113920206A CN 202111528806 A CN202111528806 A CN 202111528806A CN 113920206 A CN113920206 A CN 113920206A
- Authority
- CN
- China
- Prior art keywords
- shift
- perspective
- tilt
- coordinate system
- camera
- 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
- 238000000034 method Methods 0.000 title claims abstract description 21
- 238000003384 imaging method Methods 0.000 claims abstract description 58
- 238000005457 optimization Methods 0.000 claims abstract description 18
- 239000011159 matrix material Substances 0.000 claims description 33
- 230000036544 posture Effects 0.000 claims description 12
- 230000003287 optical effect Effects 0.000 claims description 10
- 230000009466 transformation Effects 0.000 claims description 8
- 230000000694 effects Effects 0.000 claims description 4
- 238000012897 Levenberg–Marquardt algorithm Methods 0.000 claims description 3
- 238000013507 mapping Methods 0.000 claims description 3
- 230000004075 alteration Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000006467 substitution reaction Methods 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
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10004—Still image; Photographic image
Landscapes
- Engineering & Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Studio Devices (AREA)
Abstract
本发明实施例公开了一种透视移轴相机的标定方法,包括:步骤1:在待标定的透视移轴相机的视场范围内,拍摄多幅不同角度下标定板的图像,并提取出标靶点的像素坐标;步骤2:构建透视移轴相机的模型,获取模型所需的参数作为内参初始值,同时根据内参初始值与提取出标靶点的像素坐标代入模型中,计算出透视移轴相机的外参初始值;步骤3:在透视移轴相机的模型中加入畸变模型的影响,设定畸变参数初始值;步骤4:将内参初始值、外参初始值、畸变参数初始值代入模型进行优化,获得优化后的内参、外参、畸变参数。本发明的方向角度和倾斜角度与实际的角度更加接近,能够较好的反应实际成像方式,且标定精度更高、重投影误差更小。
Description
技术领域
本发明涉及计算机视觉技术领域,尤其涉及一种透视移轴相机的标定方法。
背景技术
对于在需要使用多个相机进行三维成像的场景中,相机与拍摄物体表面通常是呈一定角度。相机的景深在沿着相机方向上可以取得最大,但对于拍摄物体表面与相机呈一定倾斜角时会降低垂直于拍摄物体表面方向上的景深。因此,在三维成像中相机可以采用移轴镜头以增大拍摄物体表面方向的景深。
移轴镜头是利用了Scheimpflug原理,让镜头与成像面呈一定的角度,可以使得在相机倾斜的情况下,增大垂直于拍摄物体表面上的景深,使得多个相机进行三维成像时也能获得更大的三维成像深度范围。
透视移轴相机的标定都是通过在视场范围内放置包含标靶点的标定板,再拍摄进行标定。目前现有的基于Scheimpflug原理的透视成像相机标定的模型都是基于成像中心产生角度倾斜,但是实际成像面产生倾斜后,成像面中心位置并不一定在光轴上,而是可能产生一定的偏移。同时标定获得的结果大多与实际透视移轴相机的给出的实际物理参数相差较大,标定获得的Scheimpflug成像角度与实际相机内部的Scheimpflug成像角度不符,且标定重投影误差较大。因此,现有的标定模型不足以描述透视移轴相机的成像原理。
发明内容
本发明实施例所要解决的技术问题在于,提供一种透视移轴相机的标定方法,以提升标定精度,减小重投影误差。
为了解决上述技术问题,本发明实施例提出了一种透视移轴相机的标定方法,包括:
步骤1:在待标定的透视移轴相机的视场范围内,拍摄多幅不同角度下标定板的图像,并从标定板图像中提取出标靶点的像素坐标;
步骤2:构建透视移轴相机的模型,从所述透视移轴相机的出厂默认参数中获取模型所需的参数作为内参初始值,同时根据内参初始值与提取出标靶点的像素坐标代入构建的透视移轴相机的模型中,计算出透视移轴相机的外参初始值;
步骤3:在透视移轴相机的模型中加入畸变模型的影响,设定畸变参数初始值;
步骤4:将内参初始值、外参初始值、畸变参数初始值代入模型进行优化,获得优化后的内参、外参、畸变参数。
进一步地,所述内参包括焦距、旋转中心距、相机主点、移轴角度,外参为标定板至透视移轴相机垂直坐标系的旋转平移变换。
进一步地,所述透视移轴相机的模型为:
其中,为透视移轴相机成像的像素坐标,sx、sy分别为透视移轴相机成像
面在x方向和y方向的像素宽度,为透视移轴相机的主点坐标,为从垂直坐标系
下齐次坐标至倾斜坐标系下成像面的上点的变换矩阵,R表示外参的旋转矩阵,T表示外参
的平移向量,为在标定板坐标系的点坐标,z为经过变换后z轴的坐标值。
进一步地,所述倾斜坐标系设定为z轴垂直于相机成像平面,同时z轴经过透视移
轴相机镜头光轴上的旋转中心点,且倾斜坐标系原点在成像平面后的焦距处,倾斜坐标
系的x轴与成像平面的x轴平行,y轴与成像平面的y轴平行;所述垂直坐标系满足z轴平行于
镜头光轴方向,且坐标系原点距离旋转中心为。
进一步地,步骤2中,
进一步地,步骤3中,透视移轴相机的模型中加入畸变后,畸变发生在垂直坐标系齐次坐标至倾斜成像平面坐标的过程中,表达式有:
进一步地,步骤4中,将获得的内参初始值、外参初始值、畸变参数初始值代入优化,利用Levenberg–Marquardt算法优化,使得目标函数最小:
其中,表示实际成像点与标靶点经过透视移轴相机模型重投影点之间的误差,
n表示拍摄姿态总数,m表示标定板标靶点总数,表示第i个姿态下,第j个标靶点的实
际成像点,表示第i个姿态,参数下,第j个标靶点的重投影点;
将参数分为两步进行优化,第一步优化如下参数:
将第一步优化的参数结果作为第二步优化参数的初值,第二步优化如下参数:
其中,为焦距,为移轴旋转中心位置至坐标原点的距离,、分别表示
移轴的方向角度和倾斜角度,为透视移轴相机的主点坐标,为畸变参数,
为第i个姿态下外参旋转矩阵对应的旋转向量,为第i个姿态下外参平移向量,为
第一次优化的参数部分,为第二次优化的参数部分;
附图说明
图1是本发明实施例的透视移轴相机的标定方法的流程示意图。
图2是本发明实施例构建的透视移轴相机的模型示意图。
具体实施方式
需要说明的是,在不冲突的情况下,本申请中的实施例及实施例中的特征可以相互结合,下面结合附图和具体实施例对本发明作进一步详细说明。
请参照图1,本发明实施例的透视移轴相机的标定方法包括步骤1~步骤4。
步骤1,在透视移轴相机的视场范围内,拍摄多幅不同角度下标定板的图像,并从标定板图像中提取出标靶点的像素坐标。
步骤2,构建如图2所示的透视移轴相机的模型,从透视移轴相机的出厂默认参数中获取模型所需的参数作为内参初值,同时根据初始内参值与提取出标靶点的像素坐标代入构建的透视移轴相机的模型中,计算出透视移轴相机的外参初始值。透视移轴相机模型内参包括焦距、旋转中心距、相机主点、移轴角度,外参为标定板至透视移轴相机垂直坐标系的旋转平移变换。
步骤3:在透视移轴相机的模型中加入畸变模型的影响,设定初始畸变参数。
步骤4:将内参初始值、外参初始值、畸变参数初始值代入模型进行优化,获得优化后的内参、外参、畸变参数。
本发明实施例的步骤2中构建的透视移轴相机的模型可以表示为:
其中,为透视移轴相机成像的像素坐标,sx、sy分别为透视移轴相机成像
面在x方向和y方向的像素宽度,为透视移轴相机的主点坐标,为从垂直坐标
系下齐次坐标至倾斜坐标系下成像面的上点的变换矩阵,R表示外参的旋转矩阵,T表示外
参的平移向量,为在标定板坐标系的点坐标,z为经过变换后z轴的坐标值。
根据透视移轴相机出厂参数可以获得透视移轴相机的模型内参的初始值,包括:
倾斜坐标系和垂直坐标系的建立的方式为:所述透视移轴相机模型中成像面产生
倾斜可以看作位于透视移轴相机镜头的光轴上某一点进行旋转的,则倾斜坐标系可以设定
为z轴垂直于相机成像平面,同时z轴经过镜头光轴上的旋转中心点,且倾斜坐标系原点在
成像平面后的焦距处,倾斜坐标系的x轴与成像平面的x轴平行,y轴与成像平面的y轴平
行。由于倾斜坐标系是按照垂直坐标系在光轴中心距处按移轴角度、进行倾斜的,
可以计算得到垂直坐标系的位置,同时垂直坐标系满足z轴平行于镜头光轴方向,且坐标系
原点距离旋转中心为。
由于已经构建了透视移轴相机的模型,标定板标靶点实际坐标与拍摄图像标靶点成像像素坐标有如下关系:
由于:
尺度因子可以计算:
该角度下外参旋转矩阵和平移矩阵的初值可以获得:
考虑畸变对透视移轴相机成像的影响,将畸变模型加在在垂直坐标系齐次坐标至倾斜成像平面坐标的过程中,相应过程可以表示为:
将获得的内参初值、外参初值、畸变参数初值代入优化,利用Levenberg–Marquardt算法优化,使得目标函数最小:
其中,表示实际成像点与标靶点经过透视移轴相机模型重投影点之间的误差,n
表示不同拍摄角度姿态总数,m表示标定板标靶点总数,表示第i个姿态下,第j个标靶
点的实际成像点,表示第i个姿态,参数下,第j个标靶点的重投影点。
在实施例中,将参数分为两步进行优化,第一步对参数进行如下优化:
将第一步优化的参数结果作为第二步优化的参数初值,第二步优化如下参数:
其中,为焦距,为移轴旋转中心位置至坐标原点的距离,分别表示
移轴的方向角度和倾斜角度,为透视移轴相机的主点坐标,为畸变参数,
为第i个姿态下外参旋转矩阵对应的旋转向量,为第i个姿态下外参平移向量,为
第一次优化的参数部分,为第二次优化的参数部分。得到优化后的内参、外参、畸变系
数。
尽管已经示出和描述了本发明的实施例,对于本领域的普通技术人员而言,可以理解在不脱离本发明的原理和精神的情况下可以对这些实施例进行多种变化、修改、替换和变型,本发明的范围由所附权利要求及其等同范围限定。
Claims (8)
1.一种透视移轴相机的标定方法,其特征在于,包括:
步骤1:在待标定的透视移轴相机的视场范围内,拍摄多幅不同角度下标定板的图像,并从标定板图像中提取出标靶点的像素坐标;
步骤2:构建透视移轴相机的模型,从所述透视移轴相机的出厂默认参数中获取模型所需的参数作为内参初始值,同时根据内参初始值与提取出标靶点的像素坐标代入构建的透视移轴相机的模型中,计算出透视移轴相机的外参初始值;
步骤3:在透视移轴相机的模型中加入畸变模型的影响,设定畸变参数初始值;
步骤4:将内参初始值、外参初始值、畸变参数初始值代入模型进行优化,获得优化后的内参、外参、畸变参数。
2.如权利要求1所述的透视移轴相机的标定方法,其特征在于,所述内参包括焦距、旋转中心距、相机主点、移轴角度,外参为标定板至透视移轴相机垂直坐标系的旋转平移变换。
8.如权利要求7所述的透视移轴相机的标定方法,其特征在于,步骤4中,将获得的内参初始值、外参初始值、畸变参数初始值代入优化,利用Levenberg–Marquardt算法优化,使得目标函数最小:
其中,表示实际成像点与标靶点经过透视移轴相机模型重投影点之间的误差,n表
示拍摄姿态总数,m表示标定板标靶点总数,表示第i个姿态下,第j个标靶点的实际
成像点,表示第i个姿态,参数下,第j个标靶点的重投影点;
将参数分为两步进行优化,第一步优化如下参数:
将第一步优化的参数结果作为第二步优化参数的初值,第二步优化如下参数:
其中,为焦距,为移轴旋转中心位置至坐标原点的距离,、分别表示移
轴的方向角度和倾斜角度,为透视移轴相机的主点坐标,为畸变参数,
为第i个姿态下外参旋转矩阵对应的旋转向量,为第i个姿态下外参平移向量,为
第一次优化的参数部分,为第二次优化的参数部分;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111528806.4A CN113920206B (zh) | 2021-12-15 | 2021-12-15 | 透视移轴相机的标定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111528806.4A CN113920206B (zh) | 2021-12-15 | 2021-12-15 | 透视移轴相机的标定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113920206A true CN113920206A (zh) | 2022-01-11 |
CN113920206B CN113920206B (zh) | 2022-03-08 |
Family
ID=79248813
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111528806.4A Active CN113920206B (zh) | 2021-12-15 | 2021-12-15 | 透视移轴相机的标定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113920206B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115546311A (zh) * | 2022-09-28 | 2022-12-30 | 中国传媒大学 | 一种基于场景信息的投影仪标定方法 |
CN115690226A (zh) * | 2022-10-27 | 2023-02-03 | 合肥中科君达视界技术股份有限公司 | 一种基于Scheimpflug定律的大视场3D轮廓测量仪标定方法 |
CN116883294A (zh) * | 2023-09-08 | 2023-10-13 | 中国科学院光电技术研究所 | 一种镜头畸变参数的两步法估计方法 |
CN117523008A (zh) * | 2024-01-04 | 2024-02-06 | 北京友友天宇系统技术有限公司 | 一种相机主点标定方法、系统、装置及存储介质 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108510547A (zh) * | 2018-03-13 | 2018-09-07 | 深圳大学 | 一种远心移轴相机标定方法和系统 |
CN108510548A (zh) * | 2018-03-13 | 2018-09-07 | 深圳大学 | 一种远心移轴相机标定方法和系统 |
CN109934878A (zh) * | 2019-03-25 | 2019-06-25 | 合肥工业大学 | 一种基于相机坐标系的线性标定系统及其方法 |
CN111540022A (zh) * | 2020-05-14 | 2020-08-14 | 深圳市艾为智能有限公司 | 一种基于虚拟相机的图像一致化方法 |
CN111915685A (zh) * | 2020-08-17 | 2020-11-10 | 沈阳飞机工业(集团)有限公司 | 一种变焦摄像机标定方法 |
-
2021
- 2021-12-15 CN CN202111528806.4A patent/CN113920206B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108510547A (zh) * | 2018-03-13 | 2018-09-07 | 深圳大学 | 一种远心移轴相机标定方法和系统 |
CN108510548A (zh) * | 2018-03-13 | 2018-09-07 | 深圳大学 | 一种远心移轴相机标定方法和系统 |
CN109934878A (zh) * | 2019-03-25 | 2019-06-25 | 合肥工业大学 | 一种基于相机坐标系的线性标定系统及其方法 |
CN111540022A (zh) * | 2020-05-14 | 2020-08-14 | 深圳市艾为智能有限公司 | 一种基于虚拟相机的图像一致化方法 |
CN111915685A (zh) * | 2020-08-17 | 2020-11-10 | 沈阳飞机工业(集团)有限公司 | 一种变焦摄像机标定方法 |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115546311A (zh) * | 2022-09-28 | 2022-12-30 | 中国传媒大学 | 一种基于场景信息的投影仪标定方法 |
CN115690226A (zh) * | 2022-10-27 | 2023-02-03 | 合肥中科君达视界技术股份有限公司 | 一种基于Scheimpflug定律的大视场3D轮廓测量仪标定方法 |
CN115690226B (zh) * | 2022-10-27 | 2024-02-13 | 合肥中科君达视界技术股份有限公司 | 一种基于Scheimpflug定律的大视场3D轮廓测量仪标定方法 |
CN116883294A (zh) * | 2023-09-08 | 2023-10-13 | 中国科学院光电技术研究所 | 一种镜头畸变参数的两步法估计方法 |
CN116883294B (zh) * | 2023-09-08 | 2023-12-22 | 中国科学院光电技术研究所 | 一种镜头畸变参数的两步法估计方法 |
CN117523008A (zh) * | 2024-01-04 | 2024-02-06 | 北京友友天宇系统技术有限公司 | 一种相机主点标定方法、系统、装置及存储介质 |
CN117523008B (zh) * | 2024-01-04 | 2024-03-29 | 北京友友天宇系统技术有限公司 | 一种相机主点标定方法、系统、装置及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN113920206B (zh) | 2022-03-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113920206B (zh) | 透视移轴相机的标定方法 | |
CN110276808B (zh) | 一种单相机结合二维码测量玻璃板不平度的方法 | |
CN105096329B (zh) | 一种精确校正超广角摄像头图像畸变的方法 | |
CN109029299B (zh) | 舱段销孔对接转角的双相机测量装置及测量方法 | |
CN109859272B (zh) | 一种自动对焦双目摄像头标定方法及装置 | |
CN109272574B (zh) | 基于投影变换的线阵旋转扫描相机成像模型构建方法和标定方法 | |
CN107025670A (zh) | 一种远心相机标定方法 | |
CN109767476A (zh) | 一种自动对焦双目摄像头标定及深度计算方法 | |
EP1378790A2 (en) | Method and device for correcting lens aberrations in a stereo camera system with zoom | |
CN110099267A (zh) | 梯形校正系统、方法以及投影仪 | |
CN113920205B (zh) | 一种非同轴相机的标定方法 | |
CN110874854B (zh) | 一种基于小基线条件下的相机双目摄影测量方法 | |
CN111340737B (zh) | 图像矫正方法、装置和电子系统 | |
CN111915685B (zh) | 一种变焦摄像机标定方法 | |
CN109859137A (zh) | 一种广角相机非规则畸变全域校正方法 | |
CN111445537A (zh) | 一种摄像机的标定方法及系统 | |
CN110136068A (zh) | 基于双侧远心镜头相机之间位置标定的音膜球顶装配系统 | |
CN114998448A (zh) | 一种多约束双目鱼眼相机标定与空间点定位的方法 | |
CN111432117B (zh) | 图像矫正方法、装置和电子系统 | |
CN111754584A (zh) | 一种远距离大视场相机参数标定系统和方法 | |
CN115239816A (zh) | 一种相机标定方法、系统、电子设备及存储介质 | |
CN113012279B (zh) | 一种非接触三维成像测量方法、系统及计算机可读存储介质 | |
CN115797460A (zh) | 一种水下双目标定方法 | |
CN110232715B (zh) | 一种多深度相机自校准的方法、装置及系统 | |
JP2002350131A (ja) | 多眼カメラのキャリブレーション方法および装置並びにコンピュータプログラム |
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 |