CN110930461B - 利用靶标求解针孔摄像机内参数的方法、存储介质和系统 - Google Patents
利用靶标求解针孔摄像机内参数的方法、存储介质和系统 Download PDFInfo
- Publication number
- CN110930461B CN110930461B CN201911138630.4A CN201911138630A CN110930461B CN 110930461 B CN110930461 B CN 110930461B CN 201911138630 A CN201911138630 A CN 201911138630A CN 110930461 B CN110930461 B CN 110930461B
- Authority
- CN
- China
- Prior art keywords
- image
- pinhole camera
- solving
- camera
- parameters
- 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
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/30—Subject of image; Context of image processing
- G06T2207/30204—Marker
- G06T2207/30208—Marker matrix
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)
- Image Analysis (AREA)
Abstract
本发明公开了一种利用靶标求解针孔摄像机内参数的方法、存储介质和系统。方法包括:从图像中提取靶标图像边缘点,使用最小二乘法拟合获得球像。根据球像的射影不变量的性质获得消失点,从而估计消失线。在获得消失线的基础上,根据球像和消失线的交点是圆环点的像,从而确定该平面上圆环点的像。利用圆环点的像对摄像机内参数的约束求解摄像机内参数。存储介质存储有运行该方法的程序。系统包括针孔摄像机和处理器,处理器包括运算器和上述存储介质,运算器运算存储介质中存储的程序。本发明为一个线性求解的过程,计算量小,场景布置简单,对靶标无特殊尺寸要求,适用广泛。
Description
技术领域
本发明涉及,尤其是一种利用靶标求解针孔摄像机内参数的方法、存储有运行该方法的程序的存储介质,以及对应的系统。
背景技术
计算机视觉的中心任务就是对图像进行理解,而它的最终目标是使计算机具有通过二维图像认知三维环境信息的能力。这种能力将不仅使机器能感知包括形状、姿态、运动等在内的三维环境中物体的几何信息,而且能对它们进行描述、存储、识别与理解。摄像机标定就是确定从三维空间点到它的二维图像点之间的映射关系,它是许多计算机视觉应用必不可少的步骤。为了确定这一映射过程,需要建立摄像机的几何成像模型,几何模型的参数称为摄像机参数,摄像机参数可分为内参数和外参数两类。内参数描述成像系统的成像几何特性,外参数描述成像系统关于世界坐标系的方向和位置。摄像机标定可分为传统标定、自标定和基于几何实体的标定。无论哪种标定方法,都旨在建立二维图像与摄像机内参数之间的约束关系,特别是线性约束关系,这是目前摄像机标定所追求的目标,也是目前计算机视觉领域研究的热点之一。
针孔摄像机成像模型简单,几何原理清晰,不需要一些特殊的镜面,在视觉领域中有重要的应用。文献“An algorithm for self calibration from several views”,(R.Hartley,In Proc.IEEE Conference on Computer Vision and PatternRecognition,pages 908–912,June 1994.)提出了一种针孔摄像机自标定方法,这类方法的优点是不需要使用标定块,缺点是必须获得图像之间的对应点。而在计算机视觉中,实现一个十分有效的寻找对应点的方法是很困难的。文献“Camera calibration by a singleimage of balls:From conics to the absolute conic”,(Teramoto H.and Xu G.,InProc.of 5th ACCV,2002,pp.499–506.)研究了在针孔摄像机下球像和绝对二次曲线的关系,通过最小化重投影误差非线性的标定内参数。而该方法需要一个好的初始化步骤,不然在最小化过程时会得到一个局部最小值。文献“Camera Calibration from Images ofSpheres”,(Hui Zhang and Kwan-Yee K.,IEEE Transactions on Pattern Analysis&Machine Intelligence.2007,29(3):499-502)研究了球像的对偶和绝对二次曲线的像的关系,并将该关系应用于摄像机的标定。文献“Interpreting Sphere Images Using theDouble-Contact Theorem”,(X.Ying,H.Zha,Springer Berlin Heidelberg,2006,3851(91):724-733)介绍了双接触原理,利用双接触原理可以确定三个球像与绝对二次曲线的像的关系,使用这个关系建立了针孔摄像机内参数的线性约束,通过此线性约束即可获得针孔摄像机的内参数。
球作为一种常见的几何体,其最重要的优点在于无自身遮挡,从任何一个方向看空间中一个球的封闭轮廓线是一个圆,并且它的投影轮廓线可全部提取。由于球具的视觉几何特性,因此利用球进行摄像机标定已成为近年来的一个研究方向。文献“Cameracalibration:a quick and easy way to determine the scale factor”,(M.A.Penna,IEEE Trans on Pattern Analysis and Machine Intelligence,vol.13,no.12,pp.1240-1245,Dec.1991.)首次提出了利用球像计算,通过在像平面上最小二乘拟合球像估计纵横比,但该方法的误差较大,不能得到精确解。文献“Camera calibration using spheres:Asemi-definite programming approach”,(Agrawal M.and Davis L.S.,In Proc.of IEEEInternational Conf.on Computer Vision,2003,pp.782–789.)研究了球图像的对偶形式与绝对二次曲线投影的对偶之间的代数关系,通过半正定规划算法确定摄像机内参数,但该方法计算复杂,且在边界拟合或椭圆提取不够好的情况下,该方法存在退化的情况。文献“Camera calibration from spheres images”,(Daucher D.,Dhome M.,in Proc.ECCV,1994,pp.449–454.)提出了一种非线性的估计摄像机的内参数的方法,且该方法需要经过多个步骤,首先,利用椭圆的两个主要轴的交点确定主点和纵横比,接下来,根据椭圆的系数方程确定焦距。然而,该方法会造成误差累计且只能估计部分摄像机内参数。
发明内容
本发明的发明目的在于:针对上述存在的问题,提供一种利用靶标求解针孔摄像机内参数的方法、存储介质和系统。以提供一种制作简单、适用广泛、稳定性好、计算量小(线性计算)的标定针孔摄像机的方案。
本发明采用的技术方案如下:
一种利用靶标求解针孔摄像机内参数的方法,其包括以下步骤:
A.在空间中设置3个在水平面上投影分离的球S1,S2和S3,利用针孔摄像机在同一位姿采集至少3幅不同视角的场景图像,针孔摄像机与其中一个球S1的球心在同一竖轴上;对于每一幅场景图像,提取场景图像边缘点的像素坐标,基于所提取的像素坐标,计算出对应的球像方程,所述球像为球朝针孔摄像机方向在像平面上的投影;
B.针对于每一幅场景图像,基于其球像方程,计算出在像平面上的消失点对,并基于消失点对,计算出对应的消失线;
C.针对于每一幅场景图像,基于其消失线和球像方程,计算出圆环点的像对;
D.根据圆环点的像与绝对二次曲线的像的线性约束关系,基于各场景图像对应的圆环点的像对,计算出绝对二次曲线的像,再根据该绝对二次曲线的像,计算出针孔摄像机的内参数。
上述方法的试验场景布置简单,方法操作过程简单;试验场景对标靶和环境无特殊要求,可以适用于各种摄像机的标定;参数求解过程为线性求解过程,因此,计算量小,计算过程简单。
进一步的,所述步骤B中,计算消失点对的方法为:
分别提取S1,S2和S3的球像方程的系数矩阵c1,c2和c3,将c1分别与c2、c3构成两个矩阵对,分别计算两个矩阵对的广义特征向量,根据计算结果即可得到对应的消失点对。
进一步的,所述步骤C中,基于其消失线和球像方程,计算出圆环点的像对的方法为:
求解消失线与S1的球像之间的两个交点,取求解结果的系数平均值作为圆环点的像的系数。取系数平均值是为了减小图像采集时遭受的噪声的影响。
进一步的,步骤A中,所述基于所提取的像素坐标,计算出对应的球像方程的方法为:基于所提取的像素坐标,利用最小二乘法拟合出对应的球像方程。
进一步的,所述D具体为:将各场景图像对应的圆环点的像对带入圆环点的像与绝对二次曲线的像的线性约束关系,通过SVD分解方法对该线性约束关系求解得到绝对二次曲线的像,对绝对二次曲线的像进行Cholesky分解再求逆得到针孔摄像机内参数矩阵,对应提取该内参数矩阵中的参数即得到所求内参数。
本发明提供了一种存储介质,该存储介质存储有程序,运行该程序可执行上述的利用靶标求解针孔摄像机内参数的方法。
本发明提供了一种利用靶标求解针孔摄像机内参数的系统,该系统包括相互连接的针孔摄像机和处理器,所述处理器包括运算器和上述的存储介质,所述运算器用于运行所述存储介质中存储的程序,以对所述针孔摄像机所采集的场景图像进行计算。
综上所述,由于采用了上述技术方案,本发明的有益效果是:
1、本发明涉及的结构非常简单,操作过程仅需将物体标靶至于环境中,采集(拍摄)几幅场景图像进行导入即可完成对摄像机内参数的计算。
2、本发明对于物体的尺寸没有要求,无需对物体位置进行严格限制,无需知道球心在世界坐标系下的坐标。
3、本发明的应用的标靶的图像边界点几乎可以全部提取,因此,可以提高曲线拟合的精确度,从而提高标定精度。
4、该方法是一个线性的算法,计算简单,仅仅只需要对球像方程特征值分解即可完成标定,较传统方式中的3次及以上多项式求解,可以大幅减小计算量。
附图说明
本发明将通过例子并参照附图的方式说明,其中:
图1是本发明求解过程所用到的试验场景(靶标)示意图。
图2是图1中靶标在像平面上的投影示意图。
具体实施方式
本说明书中公开的所有特征,或公开的所有方法或过程中的步骤,除了互相排斥的特征和/或步骤以外,均可以以任何方式组合。
本说明书(包括任何附加权利要求、摘要)中公开的任一特征,除非特别叙述,均可被其他等效或具有类似目的的替代特征加以替换。即,除非特别叙述,每个特征只是一系列等效或类似特征中的一个例子而已。
本发明实施例所用到的试验场景:
如图1所示,以空间中任意一点Ow为原点建立世界坐标系Ow-XwYwZw,在该世界坐标系中,布置三个球S1,S2和S3作为标靶,其中Zw经过其中一个球S1的球心O1。以垂直于球心与原点Ow的方向,过球心O1、O2和O3分别得到三个支撑平面π1、π2和π3,三个球S1,S2和S3分别于对应的支撑平面π1、π2和π3(球Sn对应πn,n=1、2、3)的相交线C1,C2和C3称为轮廓圆,那么三个球S1,S2和S3的轮廓圆C1,C2和C3与原点Ow形成三个正圆锥Q1,Q2和Q3。在三个球与原点Ow之间,存在像平面π,像平面π与支撑平面π1平行。三个正圆锥Q1,Q2和Q3与像平面π之间的交线cn为对应球Sn朝原点Ow方向的投影,即在像平面π上的球像,像平面π与Zw垂直。
实施例一
本实施例公开了一种利用靶标求解针孔摄像机内参数的方法,该方法的求解过程为:从图像中提取靶标图像边缘点,使用最小二乘法拟合获得球像。根据球像的射影不变量的性质获得消失点,从而估计消失线。在获得消失线的基础上,根据球像和消失线的交点是圆环点的像,从而确定该平面上圆环点的像。利用圆环点的像对摄像机内参数的约束求解摄像机内参数。具体的,方法包括场景图像采集、拟合球像方程、估计消失线、确定圆环点的像和求解针孔摄像机的内参数的步骤。
1、场景图像采集
利用针孔摄像机(下称摄像机)在原点Ow处采集至少三幅不同视角下的场景图像。摄像机的光心Oc为三个正圆锥Q1,Q2和Q3的顶点作为原点所构建的相机坐标系Oc-XcYcZc的原点。摄像机的光心Oc与原点Ow重合,这样,相机坐标系和世界坐标系之间仅仅只有旋转变换。
2、拟合球像方程
利用Matlab程序中的Edge函数分别提取各场景图像边缘点的像素坐标,并用最小二乘法拟合获得各场景图像中,各球的球像方程。
3、估计消失线
令以Oc为光心的摄像机的内参数矩阵为其中rc是纵横比,fc是有效焦距,s是倾斜因子,[u0 v0 1]T是摄像机主点p的齐次坐标矩阵形式,其中rc,fc,u0,v0,s为摄像机的5个内参数。轮廓圆C1,C2和C3与点Oc形成三个正圆锥Q1,Q2和Q3。则通过特征值分解不难计算出Q1和Q2的其中一个广义特征值λ1所对应的广义特征向量V1是平面π1上的无穷远点。同理的,另外两个正圆锥对(Q1,Q3)的广义特征向量中也封装了平面π1上的另外一个无穷远点V2。用cni表示第n(n=1,2,3)幅图像中的第i(i=1,2,3)个球像的系数矩阵。本文为了简化表述,用相同字母表示曲线和它的系数矩阵。
一般来说,球像方程cni可以通过正圆锥Qi与摄像机内参数K和旋转矩阵Rn确定,即有关系式其中λcni是非零比例因子,是正圆锥Qi的非齐次坐标表示形式。在第n幅透视图像平面上,任意取两幅球像方程cni,cnj(i,j=1,2,3,i≠j),那么矩阵对(cni,cnj)的广义特征向量等价于矩阵cnj -1cni的特征向量,通过化简,矩阵cnj -1cni满足关系式这里∝表示相差一个非零比例因子。因为矩阵cnj - 1cni和矩阵由一个非奇异的单应Hn=KRn所关联,即而两个正圆锥和的广义特征分解是射影不变量,那么通过矩阵cnj -1cni的特征向量可以确定无穷远点V1在像平面π上的像vn1和无穷远点V2在像平面π上的像vn2。连接两个消失点vn1和vn2可以获得消失线ln∞,即λnlln∞=vn1×vn2,其中λnl是非零比例因子,×表示叉积。
4、确定圆环点的像
5、求解针孔摄像机的内参数
由圆环点的像In,Jn(n=1,2,3)对绝对二次曲线的像ω的线性约束求解ω,即:其中Re,Im分别表示复数的实部和虚部。最后,对ω进行Cholesky分解再求逆便得到内参数矩阵K,即获得摄像机的5个内参数。
实施例二
本实施例公开了一种利用靶标求解针孔摄像机内参数的方法,包括以下步骤;
1、场景图像采集
利用摄像机在原点Ow处采集至少三幅不同视角下的场景图像。
2、拟合靶标投影方程
利用Matlab程序中的Edge函数提取靶标图像边缘点的像素坐标,并用最小二乘法拟合获得球像的方程。
3、估计消失线
以空间中的三个球为标定物体,由于球的投影过程同构于球面轮廓的成像过程,而球面轮廓从不同的方向观察都是一个圆。如图1所示,若以空间中任意一点Ow为原点建立世界坐标系Ow-XwYwZw,其中Zw经过其中一个球S1的中心O1。那么三个球S1,S2和S3的轮廓圆C1,C2和C3与点Ow形成三个正圆锥Q1,Q2和Q3。若假设在世界坐标系下三个球的球心的齐次坐标分别是O1=[0 0 d1 1]T,O2=[nx ny nz 1]T和O3=[mx my mz 1]T,由于Q1,Q2和Q3是正圆锥,则三个轮廓圆C1,C2和C3的支撑平面π1,π2和π3的方程为:
Z1-d1=0, (1)
假设三个球S1,S2和S3的半径分别是r1,r2和r3,则三个球的方程可以写成:
X1 2+Y1 2+(Z1-d1)2=r1 2, (4)
因为轮廓圆C1,C2和C3分别是三个球S1,S2和S3与平面π1,π2和π3相截的截口圆,则联立(1-6)式,三个正圆锥Q1,Q2和Q3在世界坐标系Ow-XwYwZw下用矩阵的形式分别表示为:
一般来说,两个正圆锥的广义特征值分解满足下式:
Q1u1=λQ2u1, (10)
联立(7)式和(8)式,则通过Maple不难计算出Q1和Q2的其中一个广义特征值及其对应的广义特征向量为:
又因为在世界坐标系下,平面OwXwYw的单位法向量为[0 0 1]T,则点V1=[-ny/nx 10]T是平面OwXwYw上的无穷远点。而平面OwXwYw平行于平面π1。因此,根据射影空间的性质,V1=[-ny/nx 1/0]T也是平面π1上的无穷远点。同理的,两个正圆锥Q1和Q3的广义特征向量中也封装了平面π1上的另外一个无穷远点V2=[-my/mx 1 0]T。
用Matlab中的Edge函数提取3幅图像中的靶标图像边缘点的像素坐标,通过最小二乘法拟合得到相应的二次曲线方程,这里用cni表示第n(n=1,2,3)幅图像中的第i(i=1,2,3)个球像的系数矩阵。
若假设球S1在垂直于Zc轴的图像平面π上的投影为cn1,根据射影变换的同素性,则图像点mn在球像cn1上,即有:
mn Tcn1mn=0。 (14)
因为KRn是3×3阶可逆矩阵,则联立(12-14)式,我们可得:
其中λcn1是非零比例因子。
同理的,若空间中球S2和S3的在图像平面π上的投影分别为cn2和cn3,则有下面的式子成立:
其中λcn2,λcn3是非零比例因子。首先考虑矩阵对(cn1,cn2),从代数上,因为矩阵对(cn1,cn2)的广义特征向量等价于矩阵cn2 -1cn1的特征向量,从(15)和(16)式可知,下面的等式满足:
这里∝表示相差一个非零比例因子。由于二次曲线对(cn1,cn2)和正圆锥对由一个非奇异的单应Hn=KRn所关联,即从文献“Euclidean Structure form Confocal Conics:Theory and Aplication to CameraCalibration”,(P.Gurdjos,J.-S.Kim,and I.-S.Kweon,In Proceeding of IEEEConference on Computer Vision and Pattern Recognition,2006,pp.1214-1221.)可知,两个正圆锥Q1和Q2的广义特征分解是射影不变量,即若wk(k=1,2,3)是矩阵Q1和Q2的广义特征向量,那么一定是矩阵cn1和cn2的广义特征向量。则矩阵cn2 -1cn1的其中一个特征向量对应于无穷远点V1的像vn1。因为V1是平面π1上的无穷远点,则其像点vn1就是平面π1上消失点。同理的,对于二次曲线对(cn1,cn3),其广义特征向量对应于平面π1上的另外一个消失点vn2。
如图2所示,在像平面上,连接两个消失点vn1和vn2可以获得消失线ln∞,即:
λnlln∞=vn1×vn2, (19)
其中λnl是非零比例因子,×表示叉积。
4、确定圆环点的像
如图2所示,已知平面π1上消失线ln∞,而平面π1与锥Q1相截于轮廓圆C1,则ln∞与球像cn1相交于cn1上的圆环点的像In和Jn,则可得到两个方程,一个是:
另一个是:
理论上,方程(20)得到的解是[a+bi c+di 1]T,方程(21)得到的解为[a-bi c-di1]T。但是由于噪声的影响,它们不能得到理想的解。若记方程(20)的解为[a1+b1i c1+d1i 1]T,方程(21)的解为[a2-b2i c2-d2i 1]T。取方程(20)和方程(21)解的系数平均值作为圆环点的像In和Jn的系数,则有:
于是可以得到平面π1上的圆环点的像In和Jn。
5、求解摄像机的内参数
3幅球像可以估计3组圆环点的像分别是I1,J1,I2,J2,I3,J3。其次,由圆环点的像对绝对二次曲线的像ω的线性约束有
其中Re,Im分别表示复数的实部和虚部。可通过SVD分解方法求解方程组(23)获得ω。最后,对ω进行Cholesky分解再求逆便可获得K,即获得了针孔摄像机的内参数。
实施例三
本实施例公开了一种利用靶标求解针孔摄像机内参数的方法,采用的试验场景为图1中的试验场景。基于空间中球的针孔摄像机标定采用的实验模板是空间中的球,如图1所示,三个球分别记为S1,S2和S3。利用本发明中的方法对用于实验的针孔摄像机进行标定,具体步骤如下:
1.拟合图像边界及靶标曲线方程
本发明采用的图像大小为1038×1048。用针孔摄像机拍摄靶标的3幅实验图像,读入图像,利用Matlab中的Edge函数提取靶标图像边缘点的像素坐标,并用最小二乘法拟合获得球像的方程。第n(n=1,2,3)幅图像中第i(i=1,2,3)个球像的方程的系数矩阵分别为cni,结果如下:
2.估计消失线
因为任意两幅球像方程cni,cnj(i,j=1,2,3,i≠j),矩阵对(cni,cnj)的广义特征向量等价于矩阵的特征向量。因此,我们计算矩阵c12 -1c11,c13 -1c11,c22 -1c21,c23 -1c21,c32 - 1c31,c33 -1c31,结果如下:
将(33)式特征值分解,则可得消失点v11的齐次坐标矩阵,结果如下:
v11=[0.3239260420991590 -0.9460823109545044 -0.0004244395816469]T,(39)
将(34)式特征值分解,则可得消失点v12的齐次坐标矩阵,结果如下:
v12=[0.5023356802841776 -0.8646726238535217 -0.0003433244800152]T,(40)
将(35)式特征值分解,则可得消失点v21的齐次坐标矩阵,结果如下:
v21=[0.1292022079882120 0.9916181939742900 0.0003831842091463]T,(41)
将(36)式特征值分解,则可得消失点v22的齐次坐标矩阵,结果如下:
v22=[-0.0360674662367494 0.9993493176587721 0.0002813791165151]T,(42)
将(37)式特征值分解,则可得消失点v31的齐次坐标矩阵,结果如下:
v31=[0.3490742252669690 0.9370910695294886 0.0027409236297956]T,(43)
将(38)式特征值分解,则可得消失点v32的齐次坐标矩阵,结果如下:
v32=[-0.0737611552888697 0.9972696041528682 0.0035536746098802]T。(44)
将(39)和(40)式带入(19)式,则可得消失线l1∞的齐次线坐标矩阵,结果如下:
l1∞=[-0.0000421880692698 -0.0001019994060190 0.1951609205217916]T,(45)
将(41)和(42)式带入(19)式,则可得消失线l2∞的齐次线坐标矩阵,结果如下:
l2∞=[-0.0001039142266072 -0.0000501752866613 0.1648832941239406]T,(46)
将(43)和(44)式带入(19)式,则可得消失线l3∞的齐次线坐标矩阵,结果如下:
l3∞=[0.0005966769176329 -0.0014426699047871 0.4172420343513371]T。(47)
3.确定圆环点的像
已知平面π1上消失线l1∞,l2∞,l3∞,而平面π1与锥Q1相截于轮廓圆C1,则将(45)式带入(20),(21)和(22)式可得c11上的圆环点的像I1和J1,结果如下:
I1=[781.6671954234+1437.7185025445i 1590.0474038550-594.6560881397i1]T, (48)
J1=[781.6671954234-1437.7185025445i 1590.0474038550+594.6560881397i1]T。 (49)
将(46)式带入(20),(21)和(22)式可得c21上的圆环点的像I2和J2,结果如下:
I2=[1217.5506060557+606.6892960676i 764.5689161555-1256.4681377356i1]T, (50)
J2=[1217.5506060557-606.6892960676i 764.5689161555+1256.4681377356i1]T。 (51)
将(47)式带入(20),(21)和(22)式可得c31上的圆环点的像I3和J3,结果如下:
I3=[265.6352244416-762.6964891388i 399.0798167865-315.4452648653i 1]T,(52)
J3=[265.6352244416+762.6964891388i 399.0798167865+315.4452648653i 1]T。(53)
4.求解针孔摄像机的内参数
将(48-53)代入(23)得到ω中元素的线性方程组,使用SVD分解求解该线性方程组得到ω的系数矩阵。结果如下:
最后,对(54)中的ω进行Cholesky分解再求逆便可获得K,结果如下:
其中纵横比rc=K(1,1)K(2,2)(K(1,1)表示矩阵K的第1行第1列的元素,K(2,2)表示矩阵K的第2行第2列的元素),故=针孔摄像机的5个内参数分别为:rc=0.9091231667311,fc=880.0078101213656509,s=0.1302122543477692,u0=319.9324197701554908,v0=239.9799653447098819。
实施例四
本实施例公开了一种存储介质,该存储介质存储有程序,运行该程序可执行上述实施例所述的利用靶标求解针孔摄像机内参数的方法。
实施例五
本实施例公开了一种利用靶标求解针孔摄像机内参数的系统,该系统包括一处理器,该处理器搭载有上述实施例四中的存储介质。
实施例六
本实施例公开了一种利用靶标求解针孔摄像机内参数的系统,该系统包括一处理器,该处理器用于执行如上述实施例一-三的利用靶标求解针孔摄像机内参数的方法。
实施例七
本实施例公开了一种利用靶标求解针孔摄像机内参数的系统,该系统包括处理器和针孔摄像机,该针孔摄像机用于采集三幅小球无重叠的场景图像,处理器包括运算器和实施例四的存储介质,运算器运行存储介质中的程序,对针孔摄像机采集的场景图像进行运算,并输出运算结果。
本发明并不局限于前述的具体实施方式。本发明扩展到任何在本说明书中披露的新特征或任何新的组合,以及披露的任一新的方法或过程的步骤或任何新的组合。
Claims (4)
1.一种利用靶标求解针孔摄像机内参数的方法,其特征在于,包括以下步骤:
A.在空间中设置3个在水平面上投影分离的球S1,S2和S3,利用针孔摄像机在同一位姿采集至少3幅不同视角的场景图像,针孔摄像机与其中一个球S1的球心在同一竖轴上;对于每一幅场景图像,提取场景图像边缘点的像素坐标,基于所提取的像素坐标,计算出对应的球像方程,所述球像为球朝针孔摄像机方向在像平面上的投影;
B.针对于每一幅场景图像,基于其球像方程,计算出在像平面上的消失点对,计算消失点对的方法为:分别提取S1,S2和S3的球像方程的系数矩阵c1,c2和c3,将c1分别与c2、c3构成两个矩阵对,分别计算两个矩阵对的广义特征向量,根据计算结果即可得到对应的消失点对;
并基于消失点对,计算出对应的消失线;
C.针对于每一幅场景图像,基于其消失线和球像方程,计算出圆环点的像对,方法为:求解消失线与S1的球像之间的两个交点,取求解结果的系数平均值作为圆环点的像的系数;
D.将各场景图像对应的圆环点的像对带入圆环点的像与绝对二次曲线的像的线性约束关系,通过SVD分解方法对该线性约束关系求解得到绝对二次曲线的像,对绝对二次曲线的像进行Cholesky分解再求逆得到针孔摄像机内参数矩阵,对应提取该内参数矩阵中的参数即得到所求内参数。
2.如权利要求1所述的利用靶标求解针孔摄像机内参数的方法,其特征在于,步骤A中,所述基于所提取的像素坐标,计算出对应的球像方程的方法为:基于所提取的像素坐标,利用最小二乘法拟合出对应的球像方程。
3.一种存储介质,其特征在于,该存储介质存储有程序,运行该程序可执行如权利要求1或2所述的利用靶标求解针孔摄像机内参数的方法。
4.一种利用靶标求解针孔摄像机内参数的系统,该系统包括相互连接的针孔摄像机和处理器,其特征在于,所述处理器包括运算器和如权利要求3所述的存储介质,所述运算器用于运行所述存储介质中存储的程序,以对所述针孔摄像机所采集的场景图像进行计算。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911138630.4A CN110930461B (zh) | 2019-11-20 | 2019-11-20 | 利用靶标求解针孔摄像机内参数的方法、存储介质和系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911138630.4A CN110930461B (zh) | 2019-11-20 | 2019-11-20 | 利用靶标求解针孔摄像机内参数的方法、存储介质和系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110930461A CN110930461A (zh) | 2020-03-27 |
CN110930461B true CN110930461B (zh) | 2022-09-20 |
Family
ID=69850353
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911138630.4A Active CN110930461B (zh) | 2019-11-20 | 2019-11-20 | 利用靶标求解针孔摄像机内参数的方法、存储介质和系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110930461B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105321181A (zh) * | 2015-10-28 | 2016-02-10 | 云南大学 | 使用双球的相离像与圆环点的像标定拋物折反射摄像机 |
CN109191528A (zh) * | 2018-07-13 | 2019-01-11 | 云南大学 | 利用球像与圆环点极线的性质标定针孔摄像机的方法 |
CN110148183A (zh) * | 2019-05-08 | 2019-08-20 | 云南大学 | 利用球及极点极线标定摄像机的方法、存储介质和系统 |
-
2019
- 2019-11-20 CN CN201911138630.4A patent/CN110930461B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105321181A (zh) * | 2015-10-28 | 2016-02-10 | 云南大学 | 使用双球的相离像与圆环点的像标定拋物折反射摄像机 |
CN109191528A (zh) * | 2018-07-13 | 2019-01-11 | 云南大学 | 利用球像与圆环点极线的性质标定针孔摄像机的方法 |
CN110148183A (zh) * | 2019-05-08 | 2019-08-20 | 云南大学 | 利用球及极点极线标定摄像机的方法、存储介质和系统 |
Non-Patent Citations (6)
Title |
---|
Flexible and accurate camera calibration using grid spherical;Z. Liu等;《Optics Express》;20171231;第25卷(第13期);15269–15285 * |
Spherical-projection-model-based self-calibration of rotation cameras;Li Baoquan等;《Proceedings of the 31st Chinese Control Conference》;20121224;4991-4996 * |
Unwrapping and stereo rectification for omnidirectional images;Jie Lei等;《Journal of Zhejiang University-SCIENCE A》;20090801;第10卷;1125–1139 * |
利用双球相对位置的成像关系标定抛物折反射摄像机;王亚林;《中国优秀硕士学位论文全文数据库 (信息科技辑)》;20150815(第08期);I138-1455 * |
基于球几何特征的摄像机内参分步标定方法;赵云涛等;《北京航空航天大学学报》;20140430;第40卷(第4期);558-563 * |
球像和圆像的射影几何性质及摄像机标定应用;杨丰澧;《中国优秀硕士学位论文全文数据库 (信息科技辑)》;20210815(第08期);I138-226 * |
Also Published As
Publication number | Publication date |
---|---|
CN110930461A (zh) | 2020-03-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Usenko et al. | The double sphere camera model | |
Puig et al. | Calibration of central catadioptric cameras using a DLT-like approach | |
CN111750820B (zh) | 影像定位方法及其系统 | |
Zhou | A new minimal solution for the extrinsic calibration of a 2D LIDAR and a camera using three plane-line correspondences | |
CN107886546B (zh) | 利用球像及公共自极三角形标定抛物折反射摄像机的方法 | |
Orghidan et al. | Camera calibration using two or three vanishing points | |
CN109064516B (zh) | 一种基于绝对二次曲线像的相机自标定方法 | |
CN111080715B (zh) | 三个球和无穷远点的极线性质标定摄像机内参数的方法 | |
CN107657645B (zh) | 利用直线及圆的共轭直径的性质标定拋物折反射摄像机 | |
CN109523595A (zh) | 一种建筑工程直线棱角间距视觉测量方法 | |
Hansen et al. | Online continuous stereo extrinsic parameter estimation | |
US20200294269A1 (en) | Calibrating cameras and computing point projections using non-central camera model involving axial viewpoint shift | |
CN108921904B (zh) | 利用单个球及渐近线的性质标定针孔摄像机的方法 | |
Takahashi et al. | Mirror-based camera pose estimation using an orthogonality constraint | |
Tao et al. | A convenient and high-accuracy multicamera calibration method based on imperfect spherical objects | |
Cao et al. | Camera calibration using symmetric objects | |
CN111223148B (zh) | 一种基于相同圆和正交性质标定摄像机内参数的方法 | |
CN110148183B (zh) | 利用球及极点极线标定摄像机的方法、存储介质和系统 | |
CN109035342B (zh) | 利用一条直线及圆环点极线标定拋物折反射摄像机的方法 | |
CN110930461B (zh) | 利用靶标求解针孔摄像机内参数的方法、存储介质和系统 | |
CN107958468B (zh) | 利用空间位置不同的三个球标定中心折反射摄像机的方法 | |
CN111429522B (zh) | 利用共面圆的公共极点极线性质标定摄像机的方法和系统 | |
CN109191528B (zh) | 利用球像与圆环点极线的性质标定针孔摄像机的方法 | |
CN113345033B (zh) | 一种中心折反射摄像机内参数标定方法及系统 | |
CN111223149A (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 |