CN113689507A - 一种基于共焦二次曲线标定针孔摄像机的方法及系统 - Google Patents
一种基于共焦二次曲线标定针孔摄像机的方法及系统 Download PDFInfo
- Publication number
- CN113689507A CN113689507A CN202111050387.8A CN202111050387A CN113689507A CN 113689507 A CN113689507 A CN 113689507A CN 202111050387 A CN202111050387 A CN 202111050387A CN 113689507 A CN113689507 A CN 113689507A
- Authority
- CN
- China
- Prior art keywords
- image
- confocal
- quadratic curve
- equation
- determining
- 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 45
- 239000011159 matrix material Substances 0.000 claims description 30
- 238000000354 decomposition reaction Methods 0.000 claims description 20
- 238000000605 extraction Methods 0.000 claims description 4
- 238000010586 diagram Methods 0.000 description 5
- 230000009977 dual effect Effects 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- PXFBZOLANLWPMH-UHFFFAOYSA-N 16-Epiaffinine Natural products C1C(C2=CC=CC=C2N2)=C2C(=O)CC2C(=CC)CN(C)C1C2CO PXFBZOLANLWPMH-UHFFFAOYSA-N 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 210000002230 centromere Anatomy 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 230000021615 conjugation Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000013519 translation 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)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明涉及一种基于共焦二次曲线标定针孔摄像机的方法级系统,方法包括:获取针孔摄像机拍摄的多张靶标图像;根据多张所述靶标图像确定共焦二次曲线的像的方程;根据所述共焦二次曲线的像的方程利用位置或者叉乘确定消失线;根据所述消失线和所述共焦二次曲线的像的方程确定圆环点的像或者正交消失点;根据所述圆环点的像求解绝对二次曲线的像,或者根据所述正交消失点求解绝对二次曲线的像;根据所述绝对二次曲线的像确定摄像机内参数。本发明能够提高摄像机内参数标定的准确性。
Description
技术领域
本发明涉及计算机视觉领域,特别是涉及一种基于共焦二次曲线标定针孔摄像机的方法及系统。
背景技术
计算机视觉可分为三个部分:图像特征的提取与处理,相关参数信息的计算,物体识别与三维重构,它的基本目的是根据若干幅场景的图像去获得三维场景中物体的形状、姿态、运动等结构信息。而摄像机标定是计算机视觉领域的关键技术之一,在二维图像平面上,像点的位置信息只与摄像机的内部结构和摄像机与场景的相对位置有关,一般称摄像机内部结构为内参数,称相对位置为外参数,而求解摄像机内参数和外参数的过程称为摄像机标定。在完成摄像机标定后,可以通过多视图几何原理来进行三维重构,从而恢复图像物体的三维信息,而摄像机标定的精度决定了恢复程度的好坏,因此研究摄像机标定具有重要意义。针孔摄像机模型通常称为线性模型,该模型成像简单,且具有固定的单一有效视点,因此被广泛研究。
在针孔摄像机标定发展过程中,最早出现的是摄像机自标定技术。有现有技术提出了一种不需要已知标定物形状的摄像机自标定方法,只需匹配图像序列之间的点对应,但算法比较复杂,难以实现。随着标定理论的成熟,后续出现了大量基于二次曲线的标定方法,其优点在于不用考虑点的匹配问题,与无穷远直线相交于绝对点,稳定性较好,在像平面易获取。有现有技术分析了同心圆的公共自极三角形性质,发现同心圆有无穷个公共自极三角形,且为直角三角形,该三角形两个点为无穷远点,另外一个点为圆心,由此提出了基于正交消失点和圆环点的像的两种标定方法。有现有技术以同心圆的投影的对偶为基础组成二次曲线束,讨论了同心圆的投影的代数约束和几何约束,进一步恢复了仿射结构和欧氏结构。有现有技术提出了一种利用同心圆的公切线性质求解圆心的像的算法,发现同心圆的公切线交于圆心,利用圆上的像点与圆心的像的连线可以构造正交消失点,进而标定摄像机,但是没有考虑二次曲线的情况。有现有技术分析了两条主轴对齐二次曲线的广义特征值分解,发现其中两个广义特征向量恰好是一对正交方向上的无穷远点,可以提供对绝对二次曲线的一个约束,并且通过求解绝对点的像可以获得绝对点的对偶的像,利用SVD分解可以获得相差一个平移变换的投影矩阵,通过张氏标定法可以进一步求解外参,但该方法在二次曲线离心率未知时不适用。有现有技术进一步分析了二次曲线主轴平行的情况,利用绝对点的对偶是圆环点的对偶的一般形式,同样可以恢复欧氏结构,并讨论了两条同心二次曲线的组合,分析了秩为1和2的情况。还有现有技术分析了有心共焦二次曲线的代数性质,通过分析共焦二次曲线的性质和广义特征值,讨论了在已知实焦点的像、已知半轴乘积、未知其它信息三种情况下如何求解圆环点的对偶的像,但是没有考虑抛物线的情况。
发明内容
本发明的目的是提供一种基于共焦二次曲线标定针孔摄像机的方法及系统,以提高摄像机内参数标定的准确性。
为实现上述目的,本发明提供了如下方案:
一种基于共焦二次曲线标定针孔摄像机的方法,包括:
获取针孔摄像机拍摄的多张靶标图像;
根据多张所述靶标图像确定共焦二次曲线的像的方程;
根据所述共焦二次曲线的像的方程利用位置或者叉乘确定消失线;
根据所述消失线和所述共焦二次曲线的像的方程确定圆环点的像或正交消失点;
根据所述圆环点的像求解绝对二次曲线的像,或者根据所述正交消失点求解绝对二次曲线的像;
根据所述绝对二次曲线的像确定摄像机内参数。
可选的,所述根据多张所述靶标图像确定共焦二次曲线的像的方程,具体包括:
提取多张所述靶标图像的边缘点的像素坐标;
利用最小二乘法拟合所述边缘点的像素坐标得到共焦二次曲线的像的方程。
可选的,所述根据所述共焦二次曲线的像的方程利用位置或者叉乘确定消失线,具体包括:
当所述共焦二次曲线的像的方程为共焦椭圆的像的方程时,通过位置确定消失线;
当所述共焦二次曲线的像的方程为共焦抛物线的像的方程时,通过叉乘确定消失线。
可选的,所述根据所述消失线和所述共焦二次曲线的像的方程确定圆环点的像或正交消失点,具体包括:
将所述消失线和所述共焦二次曲线的像的一对共轭复切线相交,确定圆环点的像;
当所述共焦二次曲线的像的方程为共焦椭圆的像的方程时,根据所述共焦椭圆的主轴的像与所述消失线的交点确定正交消失点;
当所述共焦二次曲线的像的方程为共焦抛物线的像的方程时,根据准线方向、所述共焦椭圆的像与消失线的切点确定正交消失点。
可选的,所述根据所述绝对二次曲线的像确定摄像机内参数,具体包括:
对所述绝对二次曲线的像进行Cholesky分解再求逆,得到摄像机内参数矩阵;
根据所述摄像机内参数矩阵得到摄像机内参数。
一种基于共焦二次曲线标定针孔摄像机的系统,包括:
获取模块,用于获取针孔摄像机拍摄的多张靶标图像;
共焦二次曲线的像的方程确定模块,用于根据多张所述靶标图像确定共焦二次曲线的像的方程;
消失线确定模块,用于根据所述共焦二次曲线的像的方程利用位置或者叉乘确定消失线;
圆环点的像或正交消失点确定模块,用于根据所述消失线和所述共焦二次曲线的像的方程确定圆环点的像或正交消失点;
绝对二次曲线的像求解模块,用于根据所述圆环点的像求解绝对二次曲线的像,或者根据所述正交消失点求解绝对二次曲线的像;
摄像机内参数确定模块,用于根据所述绝对二次曲线的像确定摄像机内参数。
可选的,所述共焦二次曲线的像的方程确定模块,具体包括:
提取单元,用于提取多张所述靶标图像的边缘点的像素坐标;
拟合单元,用于利用最小二乘法拟合所述边缘点的像素坐标得到共焦二次曲线的像的方程。
可选的,所述消失线确定模块,具体包括:
位置确定单元,用于当所述共焦二次曲线的像的方程为共焦椭圆的像的方程时,通过位置确定消失线;
叉乘确定单元,用于当所述共焦二次曲线的像的方程为共焦抛物线的像的方程时,通过叉乘确定消失线。
可选的,所述圆环点的像或正交消失点确定模块,具体包括:
圆环点的像确定单元,用于将所述消失线和所述共焦二次曲线的像的一对共轭复切线相交,确定圆环点的像;
第一正交消失点确定单元,用于当所述共焦二次曲线的像的方程为共焦椭圆的像的方程时,根据所述共焦椭圆的主轴的像与所述消失线的交点确定正交消失点;
第二正交消失点确定单元,用于当所述共焦二次曲线的像的方程为共焦抛物线的像的方程时,根据准线方向、所述共焦椭圆的像与消失线的切点确定正交消失点。
可选的,所述摄像机内参数确定模块,具体包括:
Cholesky分解和求逆单元,用于对所述绝对二次曲线的像进行Cholesky分解和求逆,得到摄像机内参数矩阵;
摄像机内参数确定单元,用于根据所述摄像机内参数矩阵得到摄像机内参数。
根据本发明提供的具体实施例,本发明公开了以下技术效果:
本发明提供的一种基于共焦二次曲线标定针孔摄像机的方法及系统,根据多张所述靶标图像确定共焦二次曲线的像的方程;根据所述共焦二次曲线的像的方程利用位置或者叉乘确定消失线;根据所述消失线和所述共焦二次曲线的像的方程确定圆环点的像或正交消失点;根据所述圆环点的像求解绝对二次曲线的像,或者根据所述正交消失点求解绝对二次曲线的像;根据所述绝对二次曲线的像确定摄像机内参数。本发明提供的方法可以提高曲线拟合的精确度,从而提高标定精度。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明提供的一种基于共焦二次曲线标定针孔摄像机的方法流程图;
图2为用靶标在针孔摄像机下的投影示意图;
图3为共焦椭圆在空间平面上的示意图;
图4为共焦抛物线在空间平面上的示意图;
图5为有心共焦二次曲线的像在图像平面上的示意图;
图6为无心共焦二次曲线的像在图像平面上的示意图;
图7为实验模板示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明的目的是提供一种基于共焦二次曲线标定针孔摄像机的方法及系统,以提高摄像机内参数标定的准确性。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
如图1所示,本发明提供的一种基于共焦二次曲线标定针孔摄像机的方法,包括:
步骤101:获取针孔摄像机拍摄的多张靶标图像。
步骤102:根据多张所述靶标图像确定共焦二次曲线的像的方程。其中,步骤102,具体包括:
提取多张所述靶标图像的边缘点的像素坐标。
利用最小二乘法拟合所述边缘点的像素坐标得到共焦二次曲线的像的方程。
步骤103:根据所述共焦二次曲线的像的方程利用位置或者叉乘确定消失线。其中,步骤103,具体包括:
当所述共焦二次曲线的像的方程为共焦椭圆的像的方程时,通过位置确定消失线。
当所述共焦二次曲线的像的方程为共焦抛物线的像的方程时,通过叉乘确定消失线。
步骤104:根据所述消失线和所述共焦二次曲线的像的方程确定圆环点的像或正交消失点。其中,步骤104,具体包括:
将所述消失线和所述共焦二次曲线的像的一对共轭复切线相交,确定圆环点的像。
或者,当所述共焦二次曲线的像的方程为共焦椭圆的像的方程时,根据所述共焦椭圆的主轴的像与所述消失线的交点确定正交消失点。当所述共焦二次曲线的像的方程为共焦抛物线的像的方程时,根据准线方向、所述共焦椭圆的像与消失线的切点确定正交消失点。
步骤105:根据所述圆环点的像求解绝对二次曲线的像,或者根据所述正交消失点求解绝对二次曲线的像。
步骤106:根据所述绝对二次曲线的像确定摄像机内参数。
其中,步骤106,具体包括:
对所述绝对二次曲线的像进行Cholesky分解再求逆,得到摄像机内参数矩阵。
根据所述摄像机内参数矩阵得到摄像机内参数。
本发明提供了一种制作简单,稳定性好的利用靶标求解针孔摄像机内参数的方法,该靶标由空间中两条共焦二次曲线构成。在射影空间中,根据是否存在中心,可将共焦二次曲线分为有心共焦二次曲线和无心共焦二次曲线,前者一般指椭圆和双曲线,后者为抛物线,本发明是利用空间中两条共焦二次曲线作为靶标来求解针孔摄像机内参数的方法。首先,用针孔摄像机从不同位置拍摄含有两个共焦椭圆和共焦抛物线的图像,从图像中提取靶标图像边缘点,使用最小二乘法拟合获得像的方程,利用共焦二次曲线中心的像的性质加以区分。在对偶空间中获得共焦二次曲线的像的对偶,从而得到它们的公切线。对于椭圆的像,通过判断位置关系区分出消失线,利用消失线和一组共轭复切线相交得到圆环点的像,也可以求解消失线和主轴的像的交点得到一对正交消失点。对于抛物线,通过判断叉乘是否为零区分消失线,利用消失线和一组共轭复切线相交得到圆环点的像,根据完全四点形的调和共轭性和拉盖尔定理获得两组正交消失点。利用圆环点的像或者正交消失点与绝对二次曲线的像的约束,可以求解针孔摄像机内参数。利用圆环点的像求解需要拍摄三幅图片,利用正交消失点求解至少需要拍摄五幅图片。具体分为以下步骤:
1.拟合靶标投影方程
利用Matlab平台上的Edge函数提取靶标图像边缘点的像素坐标,利用最小二乘法拟合获得靶标投影方程。
2.获得消失线
设平面π1上的两条共焦二次曲线Q1,Q2在像平面上的投影为C1,C2,取Q1上任意一点P,则P与其像点p之间满足λpp=K[R T]P,λp为P点投影的比例因子,K为摄像机内参数矩阵,记作其中fu,fv表示像平面上u轴和v轴上的尺度因子,s为倾斜因子,u0和v0为摄像机的主点坐标,R和T分别表示世界坐标系与摄像机坐标系的旋转和平移,其中fu,fv,s,u0,v0为针孔摄像机的5个内参数。
利用Matlab平台中的Edge函数提取出每幅图像中的靶标图像边缘点的像素坐标,利用最小二乘法拟合得到共焦二次曲线的像的方程,由于符号较多,这里以其中一幅图片为例。为简化表示,相同的符号既表示曲线,也表示其对应的系数矩阵,对于共焦椭圆的像C1,C2,通过位置判断区分出消失线,对于共焦抛物线像C3,C4,通过判断叉乘是否为零区分出消失线。
3.获得圆环点的像或者正交消失点
获得消失线后,将消失线与椭圆或者抛物线的像的一对共轭复切线相交,得到圆环点的像;对于椭圆像C1,C2,求解Q1,Q2的两条主轴的像与消失线的交点得到一组正交消失点,对于抛物线像C3,C4,求解C3,C4与消失线的切点以及准线方向上的消失点,可得两组正交消失点。
4.求解绝对二次曲线的像
利用圆环点的像或者正交消失点与绝对二次曲线的像的约束关系可以求解绝对二次曲线的像ω,利用圆环点的像或正交消失点求解至少需要拍摄三幅图片。
5.求解针孔摄像机内参数
绝对二次曲线的像ω与针孔摄像机内参数矩阵K满足ω=K-TK-1,对ω进行Cholesky分解再求逆得到内参数矩阵K,从而获得针孔摄像机5个内参数。其中,K-T为内参数矩阵K的逆的转置。
作为一个实施例,本发明的靶标由空间中两条共焦二次曲线构成,在求解针孔摄像机内参数过程中,使用针孔摄像机拍摄靶标至少三幅图像便可线性求解出针孔摄像机的5个内参数,具体步骤如下:
1.拟合靶标投影方程
利用Matlab程序中的Edge函数提取靶标图像边缘点的像素坐标,并用最小二乘法拟合获得靶标投影方程。
2.获得消失线
设Oc-xcyczc为摄像机坐标系,Ow-xwywzw为世界坐标系,平面共焦二次曲线在针孔摄像机下的投影示意图如图2所示,满足λC=H-TQH-1,其中λ为二次曲线Q投影的比例因子,C为共焦二次曲线的像,H为二次曲线所在支撑平面π1到像平面π的投影矩阵,H-T为H的逆的转置,Q为共焦二次曲线。如图3所示,设Q1,Q2为平面π1上的两条有心共焦二次曲线,与无穷远直线L∞相离,相交于一对共轭虚交点。平面π1与平面z=0重合,令二次曲线中心与原点重合,长轴与x轴重合,则Q1其矩阵表达式记作:
其中a和c分别为二次曲线的长轴和焦距,且c>0,为短轴,它的实焦点坐标为M=[-c 0 1]T,M'=[c 0 1]T,位于长轴LM上,虚焦点坐标为N=[0 ci 1]T,N'=[0 -ci1]T,其中i为虚数,位于短轴LN上,长轴和短轴相互垂直。ci为c与i的乘积。通过将实焦点M,M'代入到LN和L∞的直线方程中,发现小于零,而得出两个实焦点位于短轴LN的两侧,而位于无穷远直线L∞的一侧。
在图3和图4中,都存在无穷远直线L∞,从图中可以看出,图3中L∞与椭圆Q1,Q2相交于一对共轭虚点,也就是绝对点,图中没有标出;图4中L∞与抛物线Q3,Q4相切于两个实点V1∞和V2∞,这是图3和图4中最大的区别。图5和图6分别是图3和图4的投影,保持射影不变性,位置关系仍然成立。
将Q1,Q2的四条公切线记作Li(i=1~4),与无穷远直线L∞相交于圆环点I,J,将Q1,Q2投影记作C1,C2,如图5所示,在像平面π上,C1,C2公切线满足:
将C1,C2的四条公切线li(i=1~4)交点记作:
vj=li×li+1 (4)
通过排列组合vj有6种情况。两对共轭公切线(l1,l2)和(l3,l4)的交点分别记作v1和v2,即为(L1,L2)和(L3,L4)的交点M,M'的像m,m',位于长轴的像lm上。将l1与l3,l2与l4,l1与l4,l2与l3的交点分别记作v3,v4,v5,v6。由共焦二次曲线的性质可知,在(v3,v4)和(v5,v6)这两组交点中,其中一组为圆环点的像,另外一组为虚焦点的像,记l5=v3×v4,l6=v5×v6为(v3,v4)和(v5,v6)的连线。由于实焦点的像在短轴的像的两侧,在消失线的同一侧,通过将v1,v2代入到l5的直线方程中,若则消失线为l5,否则为l6,将消失线记作l∞,短轴记作ln。
如图4所示,设Q3,Q4为平面π1上两条抛物线,它们分别与无穷远直线L∞相切与一个实交点,记作V1∞和V2∞,它们的公共焦点为M1,对应的准线AA1和BB1分别位于直线LA和LB上,且LA,LB与无穷远直线的交点为V3∞,V4∞。四条公切线中两条为一对共轭复直线,记作LI和LJ,与无穷远直线相交于圆环点I,J,另外两条为重合的无穷远直线记作L∞,通过判断叉乘是否为零可以区分出L∞。将Q3,Q4投影记作C3,C4,如图6所示,在像平面π上,由射影不变性可知,C3,C4与消失线l∞相切于v1∞和v2∞,由(3)式可求解C4,C4的公切线。其四条公切线中,一对重合切线为消失线,记为l∞,另外一对公切线为共轭复直线,记为li和lj,相交于焦点的像m1,且与消失线l∞分别相交于圆环点的像mi和mj。准线的像aa1和bb1分别位于直线la和lb上,与消失线l∞相交于v3∞和v4∞。通过判断叉乘是否为0可以将l∞区分出来。
3.求解圆环点的像或者正交消失点
如图5,对于有心共焦二次曲线的像C1,C2,求解其任意一组共轭复切线与消失线l∞的交点,即为圆环点的像,求解消失线l∞与主轴的像lm和ln的交点获得一对正交消失点vm,vn,如下所示:
如图6,对于无心共焦二次曲线C3,C4,利用其共轭复切线li,lj和消失线l∞相交得到圆环点的像,记作:
C3,C4与消失线相切于两个消失点v1∞,v2∞,满足:
其焦点的像m1满足:
m1=li×lj (9)
由极点极线关系可以求解两条准线的像la,lb:
通过求解la,lb与消失线l∞的交点得到v3∞,v4∞:
v3∞=la×l∞ (11)
v4∞=lb×l∞ (12)
由完全四点形AA1JI和BB1JI的调和性和拉盖尔定理推论可知,其公共边L∞上存在两组正交方向的无穷远点,由射影不变性可知,v1∞,v3∞和v2∞,v4∞是两组正交消失点。
4.求解绝对二次曲线的像
因为Q1和Q2,Q3和Q4位于同一平面上,所求的圆环点的像mi,mj相同,圆环点的像位于绝对二次曲线的像ω上,则有:
因为mi,mj是一对共轭复点,因此只对ω提供两个约束,即:
其中,Re和Im表示取公式中的实部和虚部,三幅图片提供关于ω的六个约束,通过SVD分解求解ω。
由正交消失点与绝对二次曲线的像ω满足:
因此对于有心共焦二次曲线,五幅图片提供对ω的五个约束,对于无心共焦二次曲线,三幅图片提供对ω的六个约束,通过SVD分解可求解ω。
5.求解摄像机内参数
绝对二次曲线的像ω与针孔摄像机内参数矩阵K满足ω=K-TK-1,对ω进行Cholesky分解再求逆得到内参数矩阵K,从而获得针孔摄像机5个内参数。
作为一个实施例,本发明采用的实验模板示意图如图7所示,其中图7(a)为椭圆实验模板示意图,图7(b)为抛物线实验模板示意图。基于共焦二次曲线标定针孔摄像机采用的实验模板是共面的两条共焦椭圆和共焦抛物线。利用本发明中的方法来标定针孔摄像机,其具体步骤如下:
1.拟合靶标的投影方程
本发明利用针孔摄像机拍摄靶标的5幅图像,读入图像。利用Matlab中的Edge函数提取靶标图像边缘点的像素坐标,并利用最小二乘法拟合获得共焦二次曲线的像的方程,这里用Cn1和Cn2表示第n(n=1,2,…,5)幅图像中的两条共焦椭圆的系数矩阵,Cn3和Cn4表示第n(n=1,2,3)幅图像中的两条共焦抛物线的系数矩阵,结果如下:
2.获得消失线
因为Q1和Q2,Q3和Q4位于同一平面上,因此所求消失线相同,通过共焦二次曲线公切线性质可以判断出消失线,将五幅图片中的消失线的齐次线坐标矩阵分别记作ln∞(n=1,2,…,5),可得:
l1∞=[-0.001576096650242 0.000891226344602 1]T (34)
l2∞=[-0.000596434214352 -0.000034228805877 1]T (35)
l3∞=[-0.000277652722623 -0.0005450184685501]T (36)
l4∞=[-0.001073080263935 -0.000124870114277 1]T (37)
l5∞=[0.006502239658000 -0.0051308826738261]T (38)
3.获得圆环点的像或者正交消失点
利用消失线和共轭复切线相交可以得到圆环点的像m1i和m1j,用mni和mnj表示第n(n=1,2,3)幅图像中圆环点的像的齐次坐标矩阵,其中坐标中的i表示复数,如下所示:
求解正交消失点,将利用共焦椭圆求出的五对正交消失点记作如下:
v1m=[1428.135183943206 1403.548141369801 1]T (45)
v1n=[427.4237735895693 -366.1682850690036 1]T (46)
v2m=[1549.956902816710 2207.283323271339 1]T (47)
v2n=[1732.305092911781 -970.1193558313195 1]T (48)
v3m=[1229.420432816443 1208.487615715997 1]T (49)
v3n=[-3252.130774120768 3491.556843607697 1]T (50)
v4m=[812.3110955605783 1027.667796420715 1]T (51)
v4n=[1068.687076451496 -1175.517544075917 1]T (52)
v5m=[878.7094273063829 1308.464783328090 1]T (53)
v5n=[-203.9133729357064 -63.51609284726921 1]T (54)
利用前三幅图片中共焦抛物线求出六对正交消失点,用vn1和vn3,vn2和vn3表示第n(n=1,2,3)幅图像中的两对正交消失点,如下所示:
v11=[669.8074282172780 62.4771072557840 1]T (55)
v13=[199.5643030199445 -769.1282631552055 1]T (56)
v21=[1549.956908402741 2207.283225935183 1]T (57)
v22=[1878.035542712339 -3509.460828766771 1]T (58)
v23=[1584.783320405642 1600.435772478928 1]T (59)
v24=[1500.142285919463 3075.299052439041 1]T (60)
v31=[1229.420508234093 1208.487577295436 1]T (61)
v32=[968.134438246694 1341.596440392121 1]T (62)
vn3=[7689.338205002919 -2082.435281150562 1]T (63)
v34=[1496.640524812950 1072.3557406261321]T (64)
4.确定绝对二次曲线的像
对于圆环点的像,将(39)式到(44)式代入到(14)式中,建立关于ω的线性约束方程,对方程组进行SVD分解,可得ω的系数矩阵,结果如下:
对于共焦椭圆求出的五对正交消失点,建立关于ω的线性约束方程,对方程组进行SVD分解,可得ω的系数矩阵,结果如下:
对于共焦抛物线求出的六对正交消失点,建立关于ω的线性约束方程,对方程组进行SVD分解,可得ω的系数矩阵,结果如下:
5.求解针孔摄像机内参数
根据约束ω=K-TK-1,对(65)式中ω进行Cholesky分解再求逆得到内参数矩阵K,结果如下:
故针孔摄像机的5个内参数分别为:fu=799.9999999996533,fv=849.9999999995438,s=0.2000000003413,u0=399.9999999995597,v0=349.9999999996217。
对(66)式中的ω进行Cholesky分解再求逆得到内参数矩阵K,结果如下:
故针孔摄像机的5个内参数分别为:fu=800.0000000000920,fv=850.0000000001444,s=0.1999999999137,u0=400.0000000001802,v0=350.0000000001015。
对(67)式中的ω进行Cholesky分解再求逆得到内参数矩阵K,结果如下:
故针孔摄像机的5个内参数分别为:fu=799.9999999998992,fv=849.9999999999362,s=0.2000000000686519,u0=399.9999999999671,v0=349.9999999999783。
本发明还提供了一种基于共焦二次曲线标定针孔摄像机的系统,包括:
获取模块,用于获取针孔摄像机拍摄的多张靶标图像。
共焦二次曲线的像的方程确定模块,用于根据多张所述靶标图像确定共焦二次曲线的像的方程。
其中,所述共焦二次曲线的像的方程确定模块,具体包括:
提取单元,用于提取多张所述靶标图像的边缘点的像素坐标。
拟合单元,用于利用最小二乘法拟合所述边缘点的像素坐标得到共焦二次曲线的像的方程。
消失线确定模块,用于根据所述共焦二次曲线的像的方程利用位置或者叉乘确定消失线。
其中,所述消失线确定模块,具体包括:
位置确定单元,用于当所述共焦二次曲线的像的方程为共焦椭圆的像的方程时,通过位置确定消失线。
叉乘确定单元,用于当所述共焦二次曲线的像的方程为共焦抛物线的像的方程时,通过叉乘确定消失线。
圆环点的像或正交消失点确定模块,用于根据所述消失线和所述共焦二次曲线的像的方程确定圆环点的像或正交消失点。
其中,所述圆环点的像或正交消失点确定模块,具体包括:
圆环点的像确定单元,用于将所述消失线和所述共焦二次曲线的像的一对共轭复切线相交,确定圆环点的像。
第一正交消失点确定单元,用于当所述共焦二次曲线的像的方程为共焦椭圆的像的方程时,根据所述共焦椭圆的主轴的像与所述消失线的交点确定正交消失点。
第二正交消失点确定单元,用于当所述共焦二次曲线的像的方程为共焦抛物线的像的方程时,根据准线方向、所述共焦椭圆的像与消失线的切点确定正交消失点。
绝对二次曲线的像求解模块,用于根据所述圆环点的像求解绝对二次曲线的像,或者根据所述正交消失点求解绝对二次曲线的像。
摄像机内参数确定模块,用于根据所述绝对二次曲线的像确定摄像机内参数。其中,所述摄像机内参数确定模块,具体包括:
Cholesky分解和求逆单元,用于对所述绝对二次曲线的像进行Cholesky分解和求逆,得到摄像机内参数矩阵。
摄像机内参数确定单元,用于根据所述摄像机内参数矩阵得到摄像机内参数。
本发明具有以下优势:
(1)该靶标制作简单,只需两条共焦二次曲线。
(2)对该靶标的物理尺度没有要求,无需知道二次曲线中心在世界坐标系下的坐标。
(3)该靶标的图像边界点几乎可以全部提取,从而提高曲线拟合的精确度和标定精度。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处。综上所述,本说明书内容不应理解为对本发明的限制。
Claims (10)
1.一种基于共焦二次曲线标定针孔摄像机的方法,其特征在于,包括:
获取针孔摄像机拍摄的多张靶标图像;
根据多张所述靶标图像确定共焦二次曲线的像的方程;
根据所述共焦二次曲线的像的方程利用位置或者叉乘确定消失线;
根据所述消失线和所述共焦二次曲线的像的方程确定圆环点的像或正交消失点;
根据所述圆环点的像求解绝对二次曲线的像,或者根据所述正交消失点求解绝对二次曲线的像;
根据所述绝对二次曲线的像确定摄像机内参数。
2.根据权利要求1所述的基于共焦二次曲线标定针孔摄像机的方法,其特征在于,所述根据多张所述靶标图像确定共焦二次曲线的像的方程,具体包括:
提取多张所述靶标图像的边缘点的像素坐标;
利用最小二乘法拟合所述边缘点的像素坐标得到共焦二次曲线的像的方程。
3.根据权利要求1所述的基于共焦二次曲线标定针孔摄像机的方法,其特征在于,所述根据所述共焦二次曲线的像的方程利用位置或者叉乘确定消失线,具体包括:
当所述共焦二次曲线的像的方程为共焦椭圆的像的方程时,通过位置确定消失线;
当所述共焦二次曲线的像的方程为共焦抛物线的像的方程时,通过叉乘确定消失线。
4.根据权利要求1所述的基于共焦二次曲线标定针孔摄像机的方法,其特征在于,所述根据所述消失线和所述共焦二次曲线的像的方程确定圆环点的像或正交消失点,具体包括:
将所述消失线和所述共焦二次曲线的像的一对共轭复切线相交,确定圆环点的像;
当所述共焦二次曲线的像的方程为共焦椭圆的像的方程时,根据所述共焦椭圆的主轴的像与所述消失线的交点确定正交消失点;
当所述共焦二次曲线的像的方程为共焦抛物线的像的方程时,根据准线方向、所述共焦椭圆的像与消失线的切点确定正交消失点。
5.根据权利要求1所述的基于共焦二次曲线标定针孔摄像机的方法,其特征在于,所述根据所述绝对二次曲线的像确定摄像机内参数,具体包括:
对所述绝对二次曲线的像进行Cholesky分解再求逆,得到摄像机内参数矩阵;
根据所述摄像机内参数矩阵得到摄像机内参数。
6.一种基于共焦二次曲线标定针孔摄像机的系统,其特征在于,包括:
获取模块,用于获取针孔摄像机拍摄的多张靶标图像;
共焦二次曲线的像的方程确定模块,用于根据多张所述靶标图像确定共焦二次曲线的像的方程;
消失线确定模块,用于根据所述共焦二次曲线的像的方程利用位置或者叉乘确定消失线;
圆环点的像或正交消失点确定模块,用于根据所述消失线和所述共焦二次曲线的像的方程确定圆环点的像或正交消失点;
绝对二次曲线的像求解模块,用于根据所述圆环点的像求解绝对二次曲线的像,或者根据所述正交消失点求解绝对二次曲线的像;
摄像机内参数确定模块,用于根据所述绝对二次曲线的像确定摄像机内参数。
7.根据权利要求6所述的基于共焦二次曲线标定针孔摄像机的系统,其特征在于,所述共焦二次曲线的像的方程确定模块,具体包括:
提取单元,用于提取多张所述靶标图像的边缘点的像素坐标;
拟合单元,用于利用最小二乘法拟合所述边缘点的像素坐标得到共焦二次曲线的像的方程。
8.根据权利要求6所述的基于共焦二次曲线标定针孔摄像机的系统,其特征在于,所述消失线确定模块,具体包括:
位置确定单元,用于当所述共焦二次曲线的像的方程为共焦椭圆的像的方程时,通过位置确定消失线;
叉乘确定单元,用于当所述共焦二次曲线的像的方程为共焦抛物线的像的方程时,通过叉乘确定消失线。
9.根据权利要求6所述的基于共焦二次曲线标定针孔摄像机的系统,其特征在于,所述圆环点的像或正交消失点确定模块,具体包括:
圆环点的像确定单元,用于将所述消失线和所述共焦二次曲线的像的一对共轭复切线相交,确定圆环点的像;
第一正交消失点确定单元,用于当所述共焦二次曲线的像的方程为共焦椭圆的像的方程时,根据所述共焦椭圆的主轴的像与所述消失线的交点确定正交消失点;
第二正交消失点确定单元,用于当所述共焦二次曲线的像的方程为共焦抛物线的像的方程时,根据准线方向、所述共焦椭圆的像与消失线的切点确定正交消失点。
10.根据权利要求6所述的基于共焦二次曲线标定针孔摄像机的系统,其特征在于,所述摄像机内参数确定模块,具体包括:
Cholesky分解和求逆单元,用于对所述绝对二次曲线的像进行Cholesky分解和求逆,得到摄像机内参数矩阵;
摄像机内参数确定单元,用于根据所述摄像机内参数矩阵得到摄像机内参数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111050387.8A CN113689507B (zh) | 2021-09-08 | 2021-09-08 | 一种基于共焦二次曲线标定针孔摄像机的方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111050387.8A CN113689507B (zh) | 2021-09-08 | 2021-09-08 | 一种基于共焦二次曲线标定针孔摄像机的方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113689507A true CN113689507A (zh) | 2021-11-23 |
CN113689507B CN113689507B (zh) | 2024-04-05 |
Family
ID=78586042
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111050387.8A Active CN113689507B (zh) | 2021-09-08 | 2021-09-08 | 一种基于共焦二次曲线标定针孔摄像机的方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113689507B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH11355813A (ja) * | 1998-06-04 | 1999-12-24 | Honda Motor Co Ltd | カメラの内部パラメータ決定装置 |
US20130259403A1 (en) * | 2012-04-03 | 2013-10-03 | Oluwatosin Osinusi | Flexible easy-to-use system and method of automatically inserting a photorealistic view of a two or three dimensional object into an image using a cd,dvd or blu-ray disc |
CN111080715A (zh) * | 2019-12-18 | 2020-04-28 | 云南大学 | 三个球和无穷远点的极线性质标定摄像机内参数的方法 |
CN113345033A (zh) * | 2021-07-14 | 2021-09-03 | 云南大学 | 一种中心折反射摄像机内参数标定方法及系统 |
-
2021
- 2021-09-08 CN CN202111050387.8A patent/CN113689507B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH11355813A (ja) * | 1998-06-04 | 1999-12-24 | Honda Motor Co Ltd | カメラの内部パラメータ決定装置 |
US20130259403A1 (en) * | 2012-04-03 | 2013-10-03 | Oluwatosin Osinusi | Flexible easy-to-use system and method of automatically inserting a photorealistic view of a two or three dimensional object into an image using a cd,dvd or blu-ray disc |
CN111080715A (zh) * | 2019-12-18 | 2020-04-28 | 云南大学 | 三个球和无穷远点的极线性质标定摄像机内参数的方法 |
CN113345033A (zh) * | 2021-07-14 | 2021-09-03 | 云南大学 | 一种中心折反射摄像机内参数标定方法及系统 |
Non-Patent Citations (1)
Title |
---|
王秀平;白瑞林;刘子腾;: "由任意平行四边形确定摄像机内参数的方法", 上海交通大学学报, no. 03 * |
Also Published As
Publication number | Publication date |
---|---|
CN113689507B (zh) | 2024-04-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107886546B (zh) | 利用球像及公共自极三角形标定抛物折反射摄像机的方法 | |
CN111145232A (zh) | 一种基于特征信息变化度的三维点云自动配准方法 | |
Ying et al. | Geometric interpretations of the relation between the image of the absolute conic and sphere images | |
Ying et al. | Self-calibration of catadioptric camera with two planar mirrors from silhouettes | |
Ly et al. | Extrinsic calibration of heterogeneous cameras by line images | |
CN116129037B (zh) | 视触觉传感器及其三维重建方法、系统、设备及存储介质 | |
Cai et al. | Equivalent constraints for two-view geometry: Pose solution/pure rotation identification and 3D reconstruction | |
CN114202632A (zh) | 网格线性结构恢复方法、装置、电子设备及存储介质 | |
CN108921904B (zh) | 利用单个球及渐近线的性质标定针孔摄像机的方法 | |
CN113436269B (zh) | 图像稠密立体匹配方法、装置和计算机设备 | |
Yang et al. | Camera calibration using projective invariants of sphere images | |
Gårding et al. | Direct estimation of local surface shape in a fixating binocular vision system | |
CN113689507A (zh) | 一种基于共焦二次曲线标定针孔摄像机的方法及系统 | |
CN107958468B (zh) | 利用空间位置不同的三个球标定中心折反射摄像机的方法 | |
Lim et al. | Directions of egomotion from antipodal points | |
CN110120078B (zh) | 基于直径与共轭直径性质的摄像机标定方法和系统 | |
CN111210476B (zh) | 一种同时定位与建图的方法及装置 | |
CN111429522B (zh) | 利用共面圆的公共极点极线性质标定摄像机的方法和系统 | |
CN109272558B (zh) | 分离圆的公共自极三角形及圆环点标定针孔摄像机的方法 | |
Zhu et al. | Homography estimation based on order-preserving constraint and similarity measurement | |
CN115423854A (zh) | 一种基于多尺度特征提取的多视角点云配准及点云融合方法 | |
Tola | Multiview 3D Reconstruction of a scene containing independently moving objects | |
CN109191528B (zh) | 利用球像与圆环点极线的性质标定针孔摄像机的方法 | |
CN108090934B (zh) | 根据球截面的几何特性标定抛物折反射摄像机的方法 | |
CN115147498B (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 |