CN112308926B - 一种无公共视场的相机外参标定方法 - Google Patents
一种无公共视场的相机外参标定方法 Download PDFInfo
- Publication number
- CN112308926B CN112308926B CN202011109791.3A CN202011109791A CN112308926B CN 112308926 B CN112308926 B CN 112308926B CN 202011109791 A CN202011109791 A CN 202011109791A CN 112308926 B CN112308926 B CN 112308926B
- Authority
- CN
- China
- Prior art keywords
- camera
- points
- dimensional
- marking
- coordinates
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 39
- 238000006243 chemical reaction Methods 0.000 claims description 17
- 239000003550 marker Substances 0.000 claims description 8
- 239000011159 matrix material Substances 0.000 claims description 8
- 238000012216 screening Methods 0.000 claims 1
- 238000005259 measurement Methods 0.000 description 4
- 238000001514 detection method Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 1
- 230000001131 transforming effect 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
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10004—Still image; Photographic image
- G06T2207/10012—Stereo images
Landscapes
- Engineering & Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明公开了一种无公共视场的相机外参标定方法,两部相机的相对位置固定,其分别采集标定板不同区域的图像,两区域之间无交集;计算区域一中标记点和编码点的二维图像坐标、解算初始RT,利用初始RT建立区域一中标记点二维坐标与三维坐标之间的对应关系;并建立共线性约束方程;采用相同的方法,对相机二获取的区域二图像进行处理,建立共线性约束方程;结合,1RTw=1RT2×2RTw;变换拍摄位姿关系,每变换一次,两个相机分别采集标定板的图像,完成所有位姿的拍摄后,联立共线性约束方程,求解出相机一和相机二之间的外参1RT2;本方法可以实现自动化、高精度的无公共视场相机间外参标定。
Description
技术领域
本发明涉及相机标定领域,具体涉及一种无公共视场的相机外参标定方法。
背景技术
相机标定是视觉测量与检测领域的关键技术之一,精确的相机标定结果是保证视觉测量和检测结果准确可靠的前提;其中双目相机被广泛应用于三维测量领域,其主要利用两相机的公共视场,在公共视场中设置控制点,利用外部三维坐标测量装置测量控制点三维坐标或者给定基准距离长度,带入双目立体视觉模型求解双目相机外参;但在大尺寸工件测量或检测场景下,两个相机的结构不构成双目立体模型,没有公共视场,此时,两个相机无法拍到同一区域,无法在公共视场设置控制点,现有标定方法无法进行两相机外参标定。
发明内容
针对上述问题,本发明提出一种无公共视场的相机外参自动化标定方法,本方法利用标记点二维-三维点对、相机内参、两相机与世界坐标系的转换关系,构建共线性约束方程;通过变换不同位姿,构建一系列共线性约束方程,结合两相机之间的外参与单相机与世界坐标系之间的转换关系,通过联立共线性约束方程(一系列的方程组),能够同时求得1RTw、2RTw和1RT2,准确得到两相机之间的外参1RT2;本方法可以实现自动化、高精度的无公共视场相机间外参标定。
一种无公共视场的相机外参标定方法,预先定义世界坐标系并记录标定板中各个标记点和编码点的空间三维坐标、标定相机一和相机二的内参;
两部相机的相对位置固定,其分别采集标定板不同区域的图像,记为区域一和区域二,所述区域一和区域二内均包含有多个标记点和编码点,两区域之间无交集;
利用以下方法,求解相机一和相机二之间的外参1RT2:
S1、对相机一获取的区域一图像进行以下处理:
计算区域一中标记点和编码点的二维图像坐标;
利用编码点的二维图像坐标和空间三维坐标解算相机一与世界坐标系之间的初始RT,利用初始RT建立区域一中标记点二维图像坐标与空间三维坐标之间的对应关系;
并利用对应关系,建立共线性约束方程:
其中,(xp,yp)为第p个标记点二维图像坐标;p=1,2……n,n表示区域一中标记点的总数;(Xp,Yp,Zp)为第p个标记点的空间三维坐标;
(x01,y01,c1,Δx1,Δy1)为相机一的内参;(x01,y01)表示图像主点坐标,c1表示相机一的有效焦距,(Δx1,Δy1)表示图像主点坐标水平和竖直方向的偏移量;
式中的未知参数利用矩阵1RTw表示,其代表了当前拍摄位姿下,相机一与世界坐标系之间的转换关系:
S2、采用与步骤S1相同的方法,对相机二获取的区域二图像进行处理,建立共线性约束方程:
其中,q=1,2……t,t表示区域二中标记点的总数;(xq,yq)为第q个标记点二维图像坐标;(Xq,Yq,Zq)为第q个标记点的空间三维坐标;
(x02,y02,c2,Δx2,Δy2)为相机二的内参;其中,(x02,y02)表示图像主点坐标,c2表示相机一的有效焦距,(Δx2,Δy2)表示图像主点坐标水平和竖直方向的偏移量。
此外,1RTw=1RT2×2RTw;
S3、变换N次标定板与相机一和相机二之间的拍摄位姿关系,每变换一次,两个相机分别采集标定板的图像,并对采集到的图像重复进行步骤S1、S2;
待完成所有位姿的拍摄后,联立共线性约束方程,求解出相机一和相机二之间的外参1RT2。(以及各个位姿下的1RTw和2RTw)
进一步,在单个拍摄位姿处,区域一/区域二中的编码点数量不少于4个,标记点数量不少于18个。
进一步,步骤S1中,所述利用初始RT建立区域一中标记点二维图像坐标与空间三维坐标之间的对应关系,具体为:
利用初始RT,将标记点空间三维坐标转换到相机一坐标系,再反投到图像坐标系,计算反投得出的图像坐标点与各个标记点二维图像坐标之间的距离,将距离值小于阈值T的两个点标记为匹配点对;据此为每个标记点空间三维坐标匹配对应的二维图像坐标,完成图像上所有点的空间三维坐标与二维图像坐标之间的匹配。
优选,所述标记点包括平面标记点和立体标记点,二者均设有多个并均布在标定板上,所述立体标记点下表面与标定板固定,上表面高出标定板形成凸起并设有圆形标记点。
进一步,预先定义世界坐标系并记录标定板中各个标记点和编码点的空间三维坐标的方法包括:
方法一:预先定义标定板坐标系为世界坐标系,根据标记点和编码点的加工尺寸,分别记录其在世界坐标系中的空间三维坐标;
方法二:设置摄影测量系统,利用摄影测量系统获取标定板中各个标记点和编码点的空间三维坐标,预先定义摄影测量坐标系为世界坐标系。
进一步,利用对齐误差对外参标定结果进行评价,所述对齐误差error,计算如下:
其中,、ave表示求取平均值;i=1,2……N,表示第i个拍摄位姿下相机一与世界坐标系之间的转换关系、/>表示第i个拍摄位姿下相机二与世界坐标系之间的转换关系;1RT2为相机一和相机二之间的外参;j=1,2……S,/>表示标定板中第j个标记点的齐次世界坐标,S表示标定板中标记点的总数。
附图说明
图1为具体实施方式中标定板及两相机采集区域示意图;
图2为具体实施方式标定板与相机一和相机二之间的多个拍摄位姿示意图。
具体实施方式
一种无公共视场的相机外参标定方法,预先定义世界坐标系并记录标定板(如图1所示)中各个标记点1和编码点2的空间三维坐标、标定相机一和相机二的内参;
具体的,预先定义世界坐标系并记录标定板中各个标记点和编码点的空间三维坐标的方法包括:
方法一:预先定义标定板坐标系为世界坐标系,根据标记点和编码点的加工尺寸,分别记录其在世界坐标系中的空间三维坐标;
方法二:设置摄影测量系统,利用摄影测量系统获取标定板中各个标记点和编码点的空间三维坐标,预先定义摄影测量坐标系为世界坐标系。
两部相机的相对位置固定,其分别采集标定板不同区域的图像(如图1所示),记为区域一3和区域二4,区域一3和区域二4内均包含有多个标记点和编码点,两区域之间无交集;
利用以下方法,求解相机一和相机二之间的外参1RT2:
S1、对相机一获取的区域一图像进行以下处理:
计算区域一中标记点和编码点的二维图像坐标;
利用编码点2的二维图像坐标和空间三维坐标解算相机一与世界坐标系之间的初始RT,利用初始RT建立区域一中标记点1二维图像坐标与空间三维坐标之间的对应关系;
并利用对应关系,建立共线性约束方程:
其中,(xp,yp)为第p个标记点二维图像坐标;p=1,2……n,n表示区域一中标记点的总数;(Xp,Yp,Zp)为第p个标记点的空间三维坐标;
(x01,y01,c1,Δx1,Δy1)为相机一的内参;(x01,y01)表示图像主点坐标,c1表示相机一的有效焦距,(Δx1,Δy1)表示图像主点坐标水平和竖直方向的偏移量;
式中的未知参数利用矩阵1RTw表示,其代表了当前拍摄位姿下,相机一与世界坐标系之间的转换关系:
S2、采用与步骤S1相同的方法,对相机二获取的区域二图像进行处理,建立共线性约束方程:
其中,q=1,2……t,t表示区域二中标记点的总数;(xq,yq)为第q个标记点二维图像坐标;(Xq,Yq,Zq)为第q个标记点的空间三维坐标;(x02,y02,c2,Δx2,Δy2)为相机二的内参;
(x02,y02,c2,Δx2,Δy2)为相机二的内参;其中,(x02,y02)表示图像主点坐标,c2表示相机一的有效焦距,(Δx2,Δy2)表示图像主点坐标水平和竖直方向的偏移量;
此外,1RTw=1RT2×2RTw;
S3、变换N次标定板与相机一和相机二之间的拍摄位姿关系,如图2所示,每变换一次,两个相机分别采集标定板的图像,并对采集到的图像重复进行步骤S1、S2;
其中,为了满足解算结果精度要求,在单个拍摄位姿处,区域一/区域二中的编码点数量不少于4个,标记点数量不少于18个。
在步骤S1中,利用初始RT建立区域一中标记点二维图像坐标与空间三维坐标之间的对应关系,具体为:
利用初始RT,将标记点空间三维坐标转换到相机一坐标系,再反投到图像坐标系,计算反投得出的图像坐标点与各个标记点二维图像坐标之间的距离,将距离值小于阈值T的两个点标记为匹配点对;若匹配点对为多对,则将得到距离值最小的一对匹配点对记为最终的匹配点对,据此为每个标记点空间三维坐标匹配对应的二维图像坐标,完成图像上所有点的空间三维坐标与二维图像坐标之间的匹配。
其中,初始RT根据相机一/相机二的内参、利用PNP方法解算编码点的二维图像坐标和空间三维坐标之间的转换关系得到。
本实施例中,阈值T小于10像素。
作为本发明的一种优选实施方式,本实施例的标定板中:如图1所示,标记点包括平面标记点和立体标记点,二者均设有多个并均布在标定板上,立体标记点下表面与标定板固定,上表面高出标定板形成凸起并设有圆形标记点;立体标记点本身蕴含有空间三维信息,在解算转换关系矩阵时,其空间三维信息会比平面标记点的二维坐标解算更加准确,采用两者交叉分布的标定板,更有利于得出接近理想情况的外参转换关系。
进一步的,本实施例还涉及一种利用对齐误差对外参标定结果进行评价的方法,计算对齐误差error:
其中,、ave表示求取平均值;i=1,2……N,表示第i个拍摄位姿下相机一与世界坐标系之间的转换关系、/>表示第i个拍摄位姿下相机二与世界坐标系之间的转换关系;1RT2为相机一和相机二之间的外参;j=1,2……S,/>表示标定板中第j个标记点的齐次世界坐标,S表示标定板中标记点的总数。
当对其误差结果数值满足预设要求时,证明本实施例中得出的外参标定结果1RT2达到预期,否则需要调整标定板或增加拍摄位姿,获取更多共线性约束方程,再次求解外参,计算对其误差,判断外参是否达到预期。
Claims (6)
1.一种无公共视场的相机外参标定方法,预先定义世界坐标系并记录标定板中各个标记点和编码点的空间三维坐标、标定相机一和相机二的内参;
两部相机的相对位置固定,其分别采集标定板不同区域的图像,记为区域一和区域二,所述区域一和区域二内均包含有多个标记点和编码点,两区域之间无交集;
其特征在于,利用以下方法,求解相机一和相机二之间的外参1RT2:
S1、对相机一获取的区域一图像进行以下处理:
计算区域一中标记点和编码点的二维图像坐标;
利用编码点的二维图像坐标和空间三维坐标解算相机一与世界坐标系之间的初始RT,利用初始RT建立区域一中标记点二维图像坐标与空间三维坐标之间的对应关系;
并利用对应关系,建立共线性约束方程:
其中,(xp,yp)为第p个标记点二维图像坐标;p=1,2……n,n表示区域一中标记点的总数;(Xp,Yp,Zp)为第p个标记点的空间三维坐标;(x01,y01,c1,Δx1,Δy1)为相机一的内参;
式中的未知参数利用矩阵1RTw表示,其代表了当前拍摄位姿下,相机一与世界坐标系之间的转换关系:
S2、采用与步骤S1相同的方法,对相机二获取的区域二图像进行处理,建立共线性约束方程:
其中,q=1,2……t,t表示区域二中标记点的总数;(xq,yq)为第q个标记点二维图像坐标;(Xq,Yq,Zq)为第q个标记点的空间三维坐标;(x02,y02,c2,Δx2,Δy2)为相机二的内参;
此外,1RTw=1RT2×2RTw;
S3、变换N次标定板与相机一和相机二之间的拍摄位姿关系,每变换一次,两个相机分别采集标定板的图像,并对采集到的图像重复进行步骤S1、S2;
待完成所有位姿的拍摄后,联立共线性约束方程,求解出相机一和相机二之间的外参1RT2。
2.如权利要求1所述无公共视场的相机外参标定方法,其特征在于:在单个拍摄位姿处,区域一/区域二中的编码点数量不少于4个,标记点数量不少于18个。
3.如权利要求1所述无公共视场的相机外参标定方法,其特征在于:步骤S1中,所述利用初始RT建立区域一中标记点二维图像坐标与空间三维坐标之间的对应关系,具体为:
利用初始RT,将标记点空间三维坐标转换到相机一坐标系,再反投到图像坐标系,计算反投得出的图像坐标点与各个标记点二维图像坐标之间的距离,将距离值小于阈值T的两个点标记为匹配点对;据此为每个标记点空间三维坐标匹配对应的二维图像坐标,完成图像上所有点的空间三维坐标与二维图像坐标之间的匹配。
4.如权利要求1所述无公共视场的相机外参标定方法,其特征在于:所述标记点包括平面标记点和立体标记点,二者均设有多个并均布在标定板上,所述立体标记点下表面与标定板固定,上表面高出标定板形成凸起并设有圆形标记点。
5.如权利要求1或4所述无公共视场的相机外参标定方法,其特征在于:预先定义世界坐标系并记录标定板中各个标记点和编码点的空间三维坐标的方法包括:
方法一:预先定义标定板坐标系为世界坐标系,根据标记点和编码点的加工尺寸,分别记录其在世界坐标系中的空间三维坐标;
方法二:设置摄影测量系统,利用摄影测量系统获取标定板中各个标记点和编码点的空间三维坐标,预先定义摄影测量坐标系为世界坐标系。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011109791.3A CN112308926B (zh) | 2020-10-16 | 2020-10-16 | 一种无公共视场的相机外参标定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011109791.3A CN112308926B (zh) | 2020-10-16 | 2020-10-16 | 一种无公共视场的相机外参标定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112308926A CN112308926A (zh) | 2021-02-02 |
CN112308926B true CN112308926B (zh) | 2023-03-31 |
Family
ID=74327964
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011109791.3A Active CN112308926B (zh) | 2020-10-16 | 2020-10-16 | 一种无公共视场的相机外参标定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112308926B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113436277A (zh) * | 2021-07-15 | 2021-09-24 | 无锡先导智能装备股份有限公司 | 3d相机标定方法、装置及标定系统 |
CN114373019B (zh) * | 2021-12-29 | 2024-05-28 | 易思维(杭州)科技股份有限公司 | 一种利用最优化方法对无公共视场相机进行标定的方法 |
CN117848403A (zh) * | 2024-03-08 | 2024-04-09 | 宁德时代新能源科技股份有限公司 | 标定件和使用标定件标定检测系统的方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE2631226A1 (de) * | 1976-07-12 | 1978-01-19 | Helmut Prof Dipl Ing D Kellner | Verfahren und vorrichtung zur photogrammetrischen vermessung von raeumlichen objekten, insbesondere nach verkehrsunfaellen |
CN109522935A (zh) * | 2018-10-22 | 2019-03-26 | 易思维(杭州)科技有限公司 | 一种对双目视觉测量系统的标定结果进行评价的方法 |
CN109903341A (zh) * | 2019-01-25 | 2019-06-18 | 东南大学 | 一种车载摄像机外参动态自标定方法 |
CN111210478A (zh) * | 2019-12-31 | 2020-05-29 | 重庆邮电大学 | 一种无共同视野多相机系统外参标定方法、介质及系统 |
CN111595289A (zh) * | 2020-05-25 | 2020-08-28 | 湖北三江航天万峰科技发展有限公司 | 一种基于图像处理的三维角度测量系统及方法 |
-
2020
- 2020-10-16 CN CN202011109791.3A patent/CN112308926B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE2631226A1 (de) * | 1976-07-12 | 1978-01-19 | Helmut Prof Dipl Ing D Kellner | Verfahren und vorrichtung zur photogrammetrischen vermessung von raeumlichen objekten, insbesondere nach verkehrsunfaellen |
CN109522935A (zh) * | 2018-10-22 | 2019-03-26 | 易思维(杭州)科技有限公司 | 一种对双目视觉测量系统的标定结果进行评价的方法 |
CN109903341A (zh) * | 2019-01-25 | 2019-06-18 | 东南大学 | 一种车载摄像机外参动态自标定方法 |
CN111210478A (zh) * | 2019-12-31 | 2020-05-29 | 重庆邮电大学 | 一种无共同视野多相机系统外参标定方法、介质及系统 |
CN111595289A (zh) * | 2020-05-25 | 2020-08-28 | 湖北三江航天万峰科技发展有限公司 | 一种基于图像处理的三维角度测量系统及方法 |
Non-Patent Citations (2)
Title |
---|
《An optimal on site calibration method in measurement system of dynnamic envelope curve of high-speed vehicle》;LIU Meng 等;;《Journal of Measurement Science and Instrumentation》;20161231;第7卷(第3期);第398-404页; * |
《非参数模型摄像机标定方法及误差分析》;李华晋;《中国优秀硕士学位论文全文数据库 (工程科技Ⅱ辑)》;20140815;第C028-84页; * |
Also Published As
Publication number | Publication date |
---|---|
CN112308926A (zh) | 2021-02-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112308926B (zh) | 一种无公共视场的相机外参标定方法 | |
CN110285793B (zh) | 一种基于双目立体视觉系统的车辆智能测轨迹方法 | |
CN102376089B (zh) | 一种标靶校正方法及系统 | |
CN107886547B (zh) | 一种鱼眼相机标定方法及系统 | |
CN104240262B (zh) | 一种用于摄影测量的相机外参数标定装置及标定方法 | |
CN103033132B (zh) | 基于单目视觉的平面测量方法及装置 | |
CN111192235B (zh) | 一种基于单目视觉模型和透视变换的图像测量方法 | |
CN109522935A (zh) | 一种对双目视觉测量系统的标定结果进行评价的方法 | |
CN106097300A (zh) | 一种基于高精度运动平台的多相机标定方法 | |
CN105913439A (zh) | 一种基于激光跟踪仪的大视场摄像机标定方法 | |
CN111996883B (zh) | 一种检测公路路面宽度的方法 | |
CN107589069B (zh) | 一种物体碰撞恢复系数的非接触式测量方法 | |
CN106709955B (zh) | 基于双目立体视觉的空间坐标系标定系统和方法 | |
US20190082173A1 (en) | Apparatus and method for generating a camera model for an imaging system | |
CN106447733A (zh) | 颈椎活动度及活动轴线位置的确定方法、系统及装置 | |
CN112212788A (zh) | 基于多台手机的视觉空间点三维坐标测量方法 | |
CN109373912A (zh) | 一种基于双目视觉的非接触式六自由度位移测量方法 | |
CN107966100A (zh) | 基于相机阵列的测量方法及测量系统 | |
CN107123147A (zh) | 双目摄像机的标定方法、装置及双目摄像机系统 | |
CN108362205A (zh) | 基于条纹投影的空间测距方法 | |
CN104123726B (zh) | 基于消隐点的大锻件测量系统标定方法 | |
CN110033046A (zh) | 一种计算特征匹配点分布可信度的量化方法 | |
CN108257182A (zh) | 一种立体摄像头模组的标定方法及装置 | |
CN106500729A (zh) | 一种无需控制信息的智能手机自检校方法 | |
CN110378967B (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 | ||
CP01 | Change in the name or title of a patent holder | ||
CP01 | Change in the name or title of a patent holder |
Address after: Room 495, building 3, 1197 Bin'an Road, Binjiang District, Hangzhou City, Zhejiang Province 310051 Patentee after: Yi Si Si (Hangzhou) Technology Co.,Ltd. Address before: Room 495, building 3, 1197 Bin'an Road, Binjiang District, Hangzhou City, Zhejiang Province 310051 Patentee before: ISVISION (HANGZHOU) TECHNOLOGY Co.,Ltd. |