CN109360248A - 利用单个球及共轭直径的性质标定拋物折反射摄像机 - Google Patents
利用单个球及共轭直径的性质标定拋物折反射摄像机 Download PDFInfo
- Publication number
- CN109360248A CN109360248A CN201811478846.0A CN201811478846A CN109360248A CN 109360248 A CN109360248 A CN 109360248A CN 201811478846 A CN201811478846 A CN 201811478846A CN 109360248 A CN109360248 A CN 109360248A
- Authority
- CN
- China
- Prior art keywords
- point
- picture
- ball
- opening
- straight line
- 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.)
- Pending
Links
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
Landscapes
- Engineering & Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Image Analysis (AREA)
Abstract
一种利用单个球及共轭直径的性质标定拋物折反射摄像机的方法,其特征在于仅利用球元素。首先,分别从3幅图像中提取靶标图像边缘点和其中一幅图像的镜面轮廓投影的边缘点,使用最小二乘法拟合获得镜面轮廓投影和球像的投影;其次,根据像点和其对拓像点的关系获得对拓像点,从而拟合出球像的对拓球像;根据圆的共轭直径的性质,消失点关于球像和对拓球像的极线的交点与该消失点为一组正交消失点;在球像上取两组互异的两个点,可获得两组正交消失点;一个平面所对应的所有正交消失点对绝对二次曲线的像只提供两个约束,三幅图像提供六个约束;最后,利用正交消失点对绝对二次曲线像的约束求解摄像机内参数。
Description
技术领域
本发明属于计算机视觉领域,涉及一种利用空间中一个球及共轭直径的性质求解抛物折反射摄像机内参数的方法。
背景技术
计算机视觉的中心任务就是对图像进行理解,而它的最终目标是使计算机具有通过二维图像认知三维环境信息的能力。这种能力将不仅使机器能感知包括形状、姿态、运动等在内的三维环境中物体的几何信息,而且能对它们进行描述、存储、识别与理解。摄像机标定就是确定从三维空间点到它的二维图像点之间的映射关系,它是许多计算机视觉应用必不可少的步骤。为了确定这一映射过程,需要建立摄像机的几何成像模型,几何模型的参数称为摄像机参数,摄像机参数可分为内参数和外参数两类。内参数描述成像系统的成像几何特性,外参数描述成像系统关于世界坐标系的方向和位置。摄像机标定可分为传统标定、自标定和基于几何实体的标定。无论哪种标定方法,都旨在建立二维图像与摄像机内参数之间的约束关系,特别是线性约束关系,这是目前摄像机标定所追求的目标,也是目前计算机视觉领域研究的热点之一。
抛物折反射摄像机由一个抛物镜面和一个正交摄像机组成,它的成像视野大,是全景视觉领域研究的热点之一。文献“Catadioptric self-calibration”,(Kang S.B.,Proceedings of IEEE Conference on Computer Vision and Pattern Recognition,vol.1,pp.201-207,2000.)提出了一种折反射摄像机自标定方法,这类方法的优点是不需要使用标定块,缺点是必须获得图像之间的对应点。而在计算机视觉中,实现一个十分有效的寻找对应点的方法是很困难的。文献“Geometric properties of centralcatadioptric line images and their application in calibration”,(Barreto J.P.,Araujo H.,IEEE Transactions on Pattern Analysis and Machine Intelligence,vol.27,no.8,pp.1327-1333,2005)研究了中心折反射摄像机下直线的像的几何性质,并将这些性质应用于中心折反射摄像机的标定。文献“A new linear algorithm forcalibrating central catadioptric cameras”,(Wu F.,Duan F.,Hu Z.et al.,PatternRecognition,vol.41,no.10,pp.3166-3172,2008)介绍了对拓点和对拓像点,导出了空间中的一个点在视球上的投影和它的折反射图像点之间的关系,使用这个关系建立了中心折反射摄像机内参数的线性约束,通过此线性约束即可获得中心折反射摄像机内参数。文献“Calibration of central catadioptric cameras using a DLT-like approach”,(PuigL.,Bastanlar Y.,Sturm P.,et al.International Journal ofComputer Vision,vol.93,no.1,pp.101-114,2011)提出了一种基于三维控制点的标定方法,通过使用Veronese映射对三维点和其图像点的坐标进行了扩展,在扩展坐标的基础上基于DLT(直接线性变换)——相似方法实现了中心折反射摄像机的标定,但是这类方法需要已知三维点的位置,并且容易从图像中提取其图像点。
球作为一种常见的几何体,其最重要的优点在于无自身遮挡,从任何一个方向看空间中一个球的封闭轮廓线总是一个圆,并且它的投影轮廓线可全部提取。由于球具有丰富的视觉几何特性,因此利用球进行摄像机标定已成为近年来的一个热点。文献“Catadioptric camera calibration using geometric invariants”,(Ying X.,Hu Z.,IEEE Transactions on Pattern Analysis and Machine Intelligence,vol.26,no.10,pp.1260-1271,2004)首次提出了利用球标定中心折反射摄像机。他们证明了球在中心折反射摄像机的单位球投影模型下的像为椭圆,并且在非退化情况下一个球的投影二次曲线提供两个不变量。为了降低求解的复杂度,他们提出了一种分步标定方法,该方法至少需要4个球的投影才能完成摄像机的标定。但是该文献提出的标定方法是非线性的,计算的复杂度较高,并且该标定方法只能标定抛物折反射摄像机的部分内参数。文献“Identicalprojective geometric properties of central catadioptric line images andsphere images with applications to calibration”,(Ying X.,Zha H.,InternationalJournal of Computer Vision,vol.78,no.1,pp.89-105,2008)介绍了修正绝对二次曲线的像(MIAC)在中心折反射摄像机标定中的作用。他们通过研究球在中心折反射摄像机下的像与MIAC的几何与代数关系提出了两种线性标定算法。它们得出的结论对于对偶形式也是成立的。但是这篇文献中的理论和标定方法对于抛物折反射摄像机的情况是退化的。文献“A calibration method for paracatadioptric camera from sphere images”,(DuanH.,Wu Y.,Pattern Recognition Letters,vol.33,no.6,pp.677-684,2012)基于圆环点理论提出了一种利用对拓球像标定抛物折反射摄像机的线性方法。但是这篇文献中关于圆环点的像的选取比较复杂。
发明内容
本发明提供了一种制作简单,适用广泛,稳定性好的利用靶标求解抛物折反射摄像机内参数的方法,该靶标由空间中一个球构成。在求解抛物折反射摄像机内参数的过程中,需使用抛物折反射摄像机拍摄靶标的3幅图像线性求解出抛物折反射摄像机的5个内参数。
本发明采用如下技术方案:
用抛物折反射摄像机从不同的位置拍摄3幅含有一个球的图像。本发明是利用空间中一个球作为靶标用于求解抛物折反射摄像机内参数的方法,其特征在于仅利用球元素。首先,分别从3幅图像中提取靶标图像边缘点和其中一幅图像的镜面轮廓投影的边缘点,使用最小二乘法拟合获得镜面轮廓投影和球像的投影。其次,根据像点和其对拓像点的关系获得对拓像点,从而拟合出球像的对拓球像。在球像上取互异的两个点,并求出两个点的对拓像点。由对拓像点的定义,上述两组对拓像点提供一个消失点。根据圆的共轭直径的性质,消失点关于球像和对拓球像的极线的交点与该消失点为一组正交消失点。在球像上取两组互异的两个点,可获得两组正交消失点。一个平面所对应的所有正交消失点对绝对二次曲线的像只提供两个约束,三幅图像提供六个约束。最后,利用正交消失点对绝对二次曲线像的约束求解摄像机内参数。具体的步骤包括:拟合出镜面轮廓投影方程及靶标投影方程,估计球像的对拓球像,确定正交消失点,求解抛物折反射摄像机内参数。
1.拟合镜面轮廓投影方程及靶标投影方程
利用Matlab程序中的Edge函数提取镜面轮廓投影边缘点和靶标图像边缘点的像素坐标,并用最小二乘法拟合获得镜面轮廓投影方程和球像的方程。
2.估计球像的对拓球像
空间中的球Q,在拋物折反射摄像机的单位球模型下的投影分为两步。第一步,将球Q投影为以O为中心的单位视球上的平行小圆Sn+和Sn-(n=1,2,3表示拍摄的第n幅图像),称Sn-为Sn+的对拓圆,并且Sn-以单位视球中心O对称于Sn+,单位视球直径的两个端点互为对拓。第二步,以单位视球表面上的一点Oc为投影中心,这里Oc可看作一个摄像机的光心,将平行小圆Sn+和Sn-分别投影为抛物折反射图像平面上的二次曲线Cn+,Cn-。其中称可见的二次曲线Cn+为球Q的像,不可见的二次曲线Cn-为球像Cn+的对拓球像。对应于单位视球上的Sn-称为对拓圆,抛物折反射图像平面与直线OcO垂直。对应于单位视球直径的一个端点,该端点在抛物折反射图像平面不可见称为对拓像点,则单位视球直径的这个端点称为对拓点。令以Oc为光心的摄像机的内参数矩阵为其中rc是纵横比,fc是有效焦距,s是倾斜因子,[u0 v0 1]T是摄像机主点p的齐次坐标矩阵形式,rc,fc,s,u0,v0为抛物折反射摄像机的5个内参数。利用Matlab中的函数提取第一幅图像中的镜面轮廓投影边缘点和3幅图像靶标图像边缘点的像素坐标,通过最小二乘法拟合得到相应的二次曲线方程。这里用C0表示第一幅图像镜面轮廓投影曲线的系数矩阵,Cn+分别表示第n幅图像中的球像的系数矩阵。本文为了简化表述,用相同字母表示曲线和它的系数矩阵。通过C0可获得摄像机内参数矩阵Kc的一个初始矩阵值从而得到绝对二次曲线的像ω的初始值这里:其中 是纵横比的初始值,是有效焦距的初始值,是倾斜因子的初始值,是摄像机主点的初始齐次坐标矩阵表示,记取Cn+上的一组点则与它相对应的一组对拓像点可由关系确定,用齐次坐标矩阵表示。根据对拓像点的定义,点在球像Cn+的对拓球像Cn-上,因此可用最小二乘法拟合得到对拓球像Cn-的方程。
3.确定正交消失点
在球Q的投影小圆S1+上取两组互异的两个点A1i+和B1i+,其中下标i=1,2,用A1i-和B1i-分别表示A1i+和B1i+关于单位视球的中心O对称的点,即对拓点,则点A1i-和B1i-是S1+关于单位视球球心O的对称圆,即对拓圆S1-上。L1i+为A1i+,B1i+两点所在的直线,L1i-为A1i-,B1i-两点所在的直线,根据对拓点的定义,L1i+//L1i-,于是L1i+和L1i-具有相同的无穷远点,这里用表示L1i+,L1i-上的无穷远点。点关于小圆S1+,S1-的极线记为H1i+,H1i-,根据对极几何关系,H1i+//H1i-,于是H1i+和H1i-具有相同的无穷远点,这里用表示。点是圆S1+或S1-的一组共轭直径方向上的无穷远点。根据射影几何中圆的共轭直径性质,无穷远点为一组正交方向上的无穷远点。
用a1i+,b1i+,a1i-,b1i-分别表示A1i+,B1i+,A1i-,B1i-的像,则{a1i+,a1i-},{b1i+,b1i-}为两对对拓像点。记过a1i+,b1i+两点的直线为l1i+,过a1i-,b1i-两点的直线为l1i-,则根据射影变换的性质,直线l1i±为直线L1i±的像,于是可通过直线{l11+,l11-},{l12+,l12-}确定小圆S1+所在平面上的两个消失点d11,d12(的像)。消失点d1i关于球像C1+和对拓球像C1-的极线分别为h1i+,h1i-(H1i+,H1i-的像),极线h1i+,h1i-的交点为d'1i(的像)。消失点{d11,d'11}和{d12,d'12}为两组正交消失点。对于小圆S2+,S3+所在平面上的正交消失点{d2i,d'2i},{d3i,d'3i}(i=1,2)可用类似的方法可获得。
4.求解抛物折反射摄像机内参数
由正交消失点dij,d'ij(i=1,2,3;j=1,2)对绝对二次曲线的像ω的线性约束dij Tωd'ij=0获得ω。用SVD分解优化求解dij Tωd'ij=0。最后,对ω进行Cholesky分解再求逆便得到内参数矩Kc阵,即获得摄像机5个内参数。
本发明优点:
(1)该靶标制作简单,只需将一个球固定在一个支架上。
(2)对该靶标的物理尺度没有要求,无需知道球心在世界坐标系下的坐标。
(3)该靶标的图像边界点几乎可以全部提取,这样可以提高曲线拟合的精确度,从而提高标定精度。
附图说明
图1是用于求解抛物折反射摄像机内参数的靶标在单位视球上的示意图。
图2是靶标在抛物折反射图像平面上的投影。
具体实施方式
本发明提供了一种利用靶标求解抛物折反射摄像机内参数的方法,靶标是由空间中的一个球构成,如图1。用此靶标完成抛物折反射摄像机内参数的求解需要经过以下步骤:从折反射图像中提取镜面轮廓投影边缘点和靶标图像边缘点,使用最小二乘法拟合获得镜面轮廓投影和球的像。根据像点和其对拓像点的关系获得对拓像点,从而拟合出球像的对拓球像。在球像上取互异的四个点,并求得对拓像点,由对拓像点的定义及圆的共轭直径性质,上述四组对拓像点提供同一平面上两组正交消失点。从三个不同的方位对球拍摄图片,得到六组正交消失点。利用正交消失点对绝对二次曲线像的约束求解摄像机内参数。利用本发明中的方法对用于实验的拋物折反射摄像机进行标定,具体步骤如下:
1.拟合镜面轮廓投影方程及靶标投影方程
利用Matlab程序中的Edge函数提取镜面轮廓投影边缘点和靶标图像边缘点的像素坐标,并用最小二乘法拟合获得镜面轮廓投影方程和球像的方程。
2.估计球像的对拓球像
空间中的球Q(如图1),在拋物折反射摄像机的单位球模型下的投影分为两步。第一步,将球Q投影为以O中心的单位视球上的平行小圆Sn+,Sn-,称Sn-为Sn+的对拓圆(n=1,2,3),这一过程如图1所示(n=1为例)。第二步,通过摄像机的光心Oc将平行小圆Sn+,Sn-别投影为抛物折反射图像平面上的二次曲线Cn+,Cn-,这里称可见的二次曲线Cn+为球Q的像,不可见的二次曲线Cn-为球像Cn+的对拓球像,p是光轴OOc在像平面的投影,如图1所示。利用Matlab中的Edge函数分别提取3幅图像中的靶标图像边缘点和第1幅图像的镜面轮廓投影的像素坐标,通过最小二乘法拟合得到相应的二次曲线方程,这里用C0表示第1幅图像镜面轮廓投影曲线的系数矩阵,Cn+表示第n幅图像中的球像的系数矩阵。通过C0可获得摄像机内参数矩阵Kc的一个初始矩阵值具体如公式(1):
这里,C0(p,q)(p=1,2;q=1,2,3)表示矩阵C0的第p行第q列元素,φ为摄像机视场角的一半,ρ为抛物面镜轮廓投影椭圆C0的长半轴长。在获得的基础上可得到ω的初始值如公式(2):
取C+上的一组点则与它相对应的一组对拓像点可由关系式(3)确定:
根据对拓像点的定义,点在球像C+的对拓球像C_上,因此可用最小二乘法拟合得到对拓球像C_的方程。
3.确定正交消失点
在球Q的投影小圆S1+上分别取两组互异的两个点A1i+和B1i+,其中下标i=1,2,用A1i-和B1i-分别表示A1i+和B1i+的对拓点,如图1所示(图上以i=1为例),则点A1i-和B1i-在S1+的对拓圆S1-上。L1i+为A1i+,B1i+两点所在的直线,L1i-为A1i-,B1i-两点所在的直线。根据对拓点的定义,L1i+//L1i-,于是L1i+和L1i-具有相同的无穷远点,这里用表示L1i+,L1i-上的无穷远点。点关于小圆S1+,S1-的极线记为H1i+,H1i-,根据对极几何关系,H1i+//H1i-,于是H1i+和H1i-具有相同的无穷远点,这里用表示。点是圆S1+(S1-)的一组共轭直径方向上的无穷远点,根据射影几何中圆的共轭直径性质,无穷远点为一组正交方向上的无穷远点。
如图2所示,用a1i+,b1i+,a1i-,b1i-分别表示A1i+,B1i+,A1i-,B1i-的像,则{a1i+,a1i-},{b1i+,b1i-}为两组对拓像点,其中i=1,2。记过a1i+,b1i+两点的直线为l1i+,过a1i-,b1i-两点的直线为l1i-,设a1i±,b1i±的齐次坐标矩阵分别为[u1ia± v1ia± 1]T,[u1ib± v1ib± 1]T,l1i±的齐次线坐标矩阵为[u1il± v1il± 1]T,±中的+表示可见,-表示不可见,下标1ia,1ib,1il对应于投影的像点a1i,b1i,直线l1i,通过等式(4-7)可获得直线l1i±:
λ11+[u11l+ v11l+ 1]T=[u11a+ v11a+ 1]T×[u11b+ v11b+ 1]T, (4)
λ11-[u11l- v11l- 1]T=[u11a- v11a- 1]T×[u11b- v11b- 1]T; (5)
λ12+[u12l+ v12l+ 1]T=[u12a+ v12a+ 1]T×[u12b+ v12b+ 1]T, (6)
λ11-[u12l- v12l- 1]T=[u12a- v12a- 1]T×[u12b- v12b- 1]T, (7)
其中λ1i±(i=1,2)是非零比例因子,×表示两向量的差积。根据射影变换的性质,直线l1i±为直线L1i±的像。记过直线{l11+,l11-},{l12+,l12-}的交点分别为d11,d12,设d11和d12的齐次坐标矩阵分别为[u11d v11d 1]T和[u12d v12d 1]T,通过联立l11+和l11-的方程可得d11的坐标:
λ1d[u11d v11d 1]T=[u11l+ v11l+ 1]T×[u11l- v11l- 1]T, (8)
同理,通过联立l12+和l12-的方程可得d12的坐标:
λ2d[u12d v12d 1]T=[u12l+ v12l+ 1]T×[u12l- v12l- 1]T, (9)
于是可通过直线{l11+,l11-},{l12+,l12-}确定小圆S1+所在平面上的两个消失点d11,d12。下标11d,12d与d11,d12对应,11l±,12l±与l11±,l12±对应。
消失点d11关于球像和对拓球像C1+,C1-的极线分别为h11+,h11-,则
其中·表示矩阵与向量的乘积。记过直线{h11+,h11-}的交点为d'11,齐次坐标矩阵为[u'11d v'11d 1]T,h11±的齐次坐标为齐次坐标矩阵为[u11h± v11h± 1]T,通过联立h1i+和h1i-的方程可得d'11的坐标:
λ11d[u'11d v'11d 1]T=[u11h+ v11h+ 1]T×[u11h- v11h- 1]T, (11)
消失点d12关于球像和对拓球像C1+,C1-的极线分别为h12+,h12-,则
记过直线{h12+,h12-}的交点为d'12,齐次坐标矩阵为[u'12d v'12d 1]T,h12±的齐次坐标为齐次坐标矩阵为[u12h± v12h± 1]T,通过联立h1i+和h1i-的方程可得d'12的坐标:
λ12d[u'12d v'12d 1]T=[u12h+ v12h+ 1]T×[u12h- v12h- 1]T, (13)
点{d1i,d'1i}为一组正交消失点。则对于小圆S2+,S3+所在平面上的正交消失点{d2i,d'2i},{d3i,d'3i}(i=1,2)点可用类似的方法可获得。
4.求解抛物折反射摄像机内参数
由正交消失点对绝对二次曲线的像的线性约束有:
可用SVD分解优化求解(14)获得ω,最后对进行Cholesky分解得再求逆便得到内参数矩阵Kc阵,即获得摄像机5个内参数。
实施例
本发明提出了一种利用空间的一个球作为靶标线性确定抛物折反射摄像机内参数的方法。本发明采用的实验模板结构示意图如图1所示。下面以一实例对本发明的实施方案做出更为详细的描述。
基于空间中球的抛物折反射摄像机标定采用的实验模板是空间中的一个球,如图1所示,球为Q。利用本发明中的方法对用于实验的抛物折反射摄像机进行标定,具体步骤如下:
1.拟合图像边界及靶标曲线方程
本发明采用的图像大小为1300×1200。用抛物折反射摄像机拍摄靶标的3幅实验图像,读入图像,利用Matlab中的Edge函数提取第1幅图像镜面轮廓投影边缘点和3幅图像靶标图像边缘点的像素坐标,并用最小二乘法拟合获得镜面轮廓投影方程和球像的方程。第1幅图像镜面轮廓投影方程的系数矩阵为C0,3幅球像的方程的系数矩阵分别为Cn+(n=1,2,3),结果如下:
2.估计球像的对拓球像
将(15)代入(1)和(2)可得结果如下:
先分别在球像C1+、C2+和C3+上取至少5个互异的点,再将所取的点和(19)分别带入(3)得对拓球像上点的坐标,通过最小二乘法拟合获得对拓球像C1-、C2-和C3-的估计,结果如下:
3.确定正交消失点
在球像C1+上任取两组两个互异的点,记为{a11+,b11+},{a12+,b12+},它们的齐次坐标矩阵分别为:
a11+=[-2116.45842297415 -5239.77849311033 1]T, (23)
b11+=[-2124.27492846507 -1121.07919300893 1]T; (24)
a12+=[-5120.98102017421 -6914.58536653278 1]T, (25)
b12+=[-1045.71415540725 -3034.61929577788 1]T。 (26)
根据对拓像点的性质,将(23-26)和(19)分别带入(3)可获得与点a11+,b11+,a12+,b12+对应的对拓像点a11-,b11-,a12-,b12-,结果如下:
a11-=[761.83169955416 1000.88071125125 1]T, (27)
b11-=[1244.43965971482 1068.95686097796 1]T; (28)
a12-=[758.629642708597 859.748550065155 1]T, (29)
b12-=[85.337170385310 121.727850699579 1]T。 (30)
将(23)和(24)带入(4)可得过点a11+,b11+的直线l11+的齐次线坐标矩阵,结果如下:
l11+=[0.00047027784685 0.00000089249763 1]T, (31)
将(27)和(28)带入(5)可得过点a11-,b11-的直线l11-的齐次线坐标矩阵,结果如下:
l11-=[0.00015788687038 -0.00111929744495 1]T。 (32)
将(31)和(32)代入(8)可得消失点d11,结果如下:
d11=[-2127.52851345658 593.31072751158 1]T。 (33)
将(33)带入(10)可的消失点d11关于球像C1+,C1-的极线h11+,h11-,结果如下:
h11+=[-0.00002845876767 0.00051560236084 1]T, (34)
h11-=[-0.00111037084247 0.00010657542809 1]T。 (35)
将(34)和(35)带入(11)可获得消失点d'11,结果如下:
d'11=[718.25036004216 -1899.83513319264 1]T。 (36)
将(25)和(26)带入(6)可得过点a12+,b12+的直线l12+的齐次线坐标矩阵,结果如下:
l12+=[-0.00046692865378 0.00049043183272 1]T, (37)
将(29)和(30)带入(7)可得过点a12-,b12-的直线l12-的齐次线坐标矩阵,结果如下:
l12-=[-0.00188393314479 0.00049922448661 1]T。 (38)
将(37)和(38)代入(9)可得消失点d12,结果如下:
d12=[-12.72750872269 -2051.13691936468 1]T。 (39)
将(39)带入(12)可的消失点d12关于球像C1+,C1-的极线h12+,h12-,结果如下:
h12+=[0.00051090157822 0.00003133237666 1]T, (40)
h12-=[-0.00008031044905 -0.00086234026043 1]T。 (41)
将(40)和(41)带入(13)可获得消失点d'12,结果如下:
d'12=[-2040.09378790757 1349.63065231659 1]T。 (42)
在球像C2+上任取两组两个互异的点{a21+,b21+},{a22+,b22+},它们的齐次坐标矩阵分别为:
a21+=[2621.64034587542 2109.44584911265 1]T, (43)
b21+=[2932.67137635217 4885.30468734724 1]T; (44)
a22+=[2304.00696058136 2229.61718605529 1]T, (45)
b22+=[3431.77424997081 2221.28998124386 1]T。 (46)
利用公式(3-13)类似的方法获得正交消失点,结果如下:
d21=[2463.25229220234 695.87988437505 1]T, (47)
d'21=[512.55500622276 2369.01889414111 1]T; (48)
d22=[541.61643588448 2242.63032005460 1]T, (49)
d'22=[2597.93592726816 664.50894920648 1]T。 (50)
在球像C3+上任取两组两个互异的点{a31+,b31+},{a32+,b32+},它们的齐次坐标矩阵分别为:
a31+=[2385.73437099712 -397.36353457557 1]T, (51)
b31+=[3689.23701587339 -2081.46555485263 1]T; (52)
a32+=[2677.54946331018 -405.33897556082]T, (53)
b32+=[2190.23420069052 -547.17464901812 1]T。 (54)
利用公式(3-13)类似的方法获得正交消失点,结果如下:
d31=[1676.02017663872 519.57446569350 1]T, (55)
d'31=[4478.02259987014 -8228.44029300166 1]T; (56)
d32=[36.43216257061 -1174.05007444243 1]T, (57)
d'32=[1700.99656795318 312.31135034773 1]T。 (58)
4.求解抛物折反射摄像机内参数
将(55-58)(47-50)和(33,36,39,42)带入(14)得到ω中元素的线性方程组,使用SVD分解求解该线性方程组得到ω的系数矩阵。结果如下:
最后,对(59)中的ω进行Cholesky分解再求逆便可获得Kc,结果如下:
其中纵横比rc=Kc(1,1)/Kc(2,2)(Kc(1,1)表示矩阵Kc的第1行第1列的元素,Kc(2,2)表示矩阵Kc的第2行第2列的元素),故抛物折反射摄像机的5个内参数分别为:
rc=1.00470476660717,fc=399.710536550706,s=3.597610238596,u0=696.392566712461,v0=580.066537232092。
Claims (1)
1.一种利用单个球及共轭直径的性质标定拋物折反射摄像机的方法,其特征在于仅利用球元素;首先,分别从3幅图像中提取靶标图像边缘点和其中一幅图像的镜面轮廓投影的边缘点,使用最小二乘法拟合获得镜面轮廓投影和球像的投影;其次,根据像点和其对拓像点的关系获得对拓像点,从而拟合出球像的对拓球像;在球像上取互异的两个点,并求出两个点的对拓像点;由对拓像点的定义,上述两组对拓像点提供一个消失点;根据圆的共轭直径的性质,消失点关于球像和对拓球像的极线的交点与该消失点为一组正交消失点;在球像上取两组互异的两个点,可获得两组正交消失点;一个平面所对应的所有正交消失点对绝对二次曲线的像只提供两个约束,三幅图像提供六个约束;最后,利用正交消失点对绝对二次曲线像的约束求解摄像机内参数;具体的步骤包括:拟合出镜面轮廓投影方程及靶标投影方程,估计球像的对拓球像,确定正交消失点,求解抛物折反射摄像机内参数;
(1)确定正交消失点
在球Q的投影小圆S1+上取两组互异的两个点A1i+和B1i+,其中下标i=1,2,用A1i-和B1i-分别表示A1i+和B1i+关于单位视球的中心O对称的点,即对拓点,则点A1i-和B1i-是S1+关于单位视球球心O的对称圆,即对拓圆S1-上;L1i+为A1i+,B1i+两点所在的直线,L1i-为A1i-,B1i-两点所在的直线,根据对拓点的定义,L1i+//L1i-,于是L1i+和L1i-具有相同的无穷远点,这里用表示L1i+,L1i-上的无穷远点;点关于小圆S1+,S1-的极线记为H1i+,H1i-,根据对极几何关系,H1i+//H1i-,于是H1i+和H1i-具有相同的无穷远点,这里用表示;点是圆S1+或S1-的一组共轭直径方向上的无穷远点;根据射影几何中圆的共轭直径性质,无穷远点为一组正交方向上的无穷远点;用a1i+,b1i+,a1i-,b1i-分别表示A1i+,B1i+,A1i-,B1i-的像,则{a1i+,a1i-},{b1i+,b1i-}为两对对拓像点;记过a1i+,b1i+两点的直线为l1i+,过a1i-,b1i-两点的直线为l1i-,则根据射影变换的性质,直线l1i±为直线L1i±的像,于是通过直线{l11+,l11-},{l12+,l12-}确定小圆S1+所在平面上的两个消失点d11,d12(的像);消失点d1i关于球像C1+和对拓球像C1-的极线分别为h1i+,h1i-(H1i+,H1i-的像),极线h1i+,h1i-的交点为d'1i(的像);消失点{d11,d'11}和{d12,d'12}为两组正交消失点。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811478846.0A CN109360248A (zh) | 2018-12-05 | 2018-12-05 | 利用单个球及共轭直径的性质标定拋物折反射摄像机 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811478846.0A CN109360248A (zh) | 2018-12-05 | 2018-12-05 | 利用单个球及共轭直径的性质标定拋物折反射摄像机 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN109360248A true CN109360248A (zh) | 2019-02-19 |
Family
ID=65331186
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811478846.0A Pending CN109360248A (zh) | 2018-12-05 | 2018-12-05 | 利用单个球及共轭直径的性质标定拋物折反射摄像机 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109360248A (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109934879A (zh) * | 2019-04-22 | 2019-06-25 | 云南大学 | 利用球和公共自极三角形标定抛物折反射摄像机的方法 |
CN110060307A (zh) * | 2019-04-22 | 2019-07-26 | 云南大学 | 利用球投影的对拓关系标定抛物折反射摄像机的方法 |
CN110120078A (zh) * | 2019-05-08 | 2019-08-13 | 云南大学 | 基于直径与共轭直径性质的摄像机标定方法和系统 |
CN111223148A (zh) * | 2020-01-07 | 2020-06-02 | 云南大学 | 一种基于相同圆和正交性质标定摄像机内参数的方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106447731A (zh) * | 2016-09-22 | 2017-02-22 | 云南大学 | 利用单个球及正交消失点标定拋物折反射摄像机的方法 |
CN107644445A (zh) * | 2017-10-09 | 2018-01-30 | 云南大学 | 利用单个球及圆切线的性质标定拋物折反射摄像机的方法 |
CN107657645A (zh) * | 2017-10-09 | 2018-02-02 | 云南大学 | 利用直线及圆的共轭直径的性质标定拋物折反射摄像机 |
-
2018
- 2018-12-05 CN CN201811478846.0A patent/CN109360248A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106447731A (zh) * | 2016-09-22 | 2017-02-22 | 云南大学 | 利用单个球及正交消失点标定拋物折反射摄像机的方法 |
CN107644445A (zh) * | 2017-10-09 | 2018-01-30 | 云南大学 | 利用单个球及圆切线的性质标定拋物折反射摄像机的方法 |
CN107657645A (zh) * | 2017-10-09 | 2018-02-02 | 云南大学 | 利用直线及圆的共轭直径的性质标定拋物折反射摄像机 |
Non-Patent Citations (1)
Title |
---|
YUE ZHAO 等: "Calibrating a paracatadioptric camera by theproperty of the polar of a point at infinity with respect to a circle", 《APPLIED OPTICS》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109934879A (zh) * | 2019-04-22 | 2019-06-25 | 云南大学 | 利用球和公共自极三角形标定抛物折反射摄像机的方法 |
CN110060307A (zh) * | 2019-04-22 | 2019-07-26 | 云南大学 | 利用球投影的对拓关系标定抛物折反射摄像机的方法 |
CN110060307B (zh) * | 2019-04-22 | 2022-04-29 | 云南大学 | 利用球投影的对拓关系标定抛物折反射摄像机的方法 |
CN109934879B (zh) * | 2019-04-22 | 2022-05-20 | 云南大学 | 利用球和公共自极三角形标定抛物折反射摄像机的方法 |
CN110120078A (zh) * | 2019-05-08 | 2019-08-13 | 云南大学 | 基于直径与共轭直径性质的摄像机标定方法和系统 |
CN110120078B (zh) * | 2019-05-08 | 2022-02-22 | 云南大学 | 基于直径与共轭直径性质的摄像机标定方法和系统 |
CN111223148A (zh) * | 2020-01-07 | 2020-06-02 | 云南大学 | 一种基于相同圆和正交性质标定摄像机内参数的方法 |
CN111223148B (zh) * | 2020-01-07 | 2023-03-28 | 云南大学 | 一种基于相同圆和正交性质标定摄像机内参数的方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108470370B (zh) | 三维激光扫描仪外置相机联合获取三维彩色点云的方法 | |
CN109360248A (zh) | 利用单个球及共轭直径的性质标定拋物折反射摄像机 | |
CN106327504A (zh) | 利用单个球及圆环点的像标定拋物折反射摄像机的方法 | |
CN107886546B (zh) | 利用球像及公共自极三角形标定抛物折反射摄像机的方法 | |
CN109325983A (zh) | 利用无穷远点关于圆极线性质标定拋物折反射摄像机 | |
CN104217435B (zh) | 两个相互遮挡的球线性确定拋物折反射摄像机内参数方法 | |
CN107644445A (zh) | 利用单个球及圆切线的性质标定拋物折反射摄像机的方法 | |
CN107657645B (zh) | 利用直线及圆的共轭直径的性质标定拋物折反射摄像机 | |
CN106447731A (zh) | 利用单个球及正交消失点标定拋物折反射摄像机的方法 | |
CN106651956A (zh) | 利用单个球及平行圆性质标定拋物折反射摄像机的方法 | |
CN105279758A (zh) | 使用双球的相切像与圆环点的像标定拋物折反射摄像机 | |
CN105321181A (zh) | 使用双球的相离像与圆环点的像标定拋物折反射摄像机 | |
CN108921904B (zh) | 利用单个球及渐近线的性质标定针孔摄像机的方法 | |
CN105303570A (zh) | 使用双球的相交像及圆环点的像标定拋物折反射摄像机 | |
CN109215089A (zh) | 利用单个球及弦中点的性质标定拋物折反射摄像机的方法 | |
CN109360247B (zh) | 单个球的公共自极三角形及正交消失点标定抛物摄像机 | |
CN109035342B (zh) | 利用一条直线及圆环点极线标定拋物折反射摄像机的方法 | |
CN109325982A (zh) | 利用单个球及平行圆切线性质标定拋物折反射摄像机 | |
CN109523598A (zh) | 单个球的公共自极三角形及圆环点标定抛物摄像机的方法 | |
CN107958468B (zh) | 利用空间位置不同的三个球标定中心折反射摄像机的方法 | |
CN107610184A (zh) | 利用一条直线及圆环点的像标定拋物折反射摄像机的方法 | |
CN109191528B (zh) | 利用球像与圆环点极线的性质标定针孔摄像机的方法 | |
CN105354840A (zh) | 使用相离的双球像与正交消失点标定拋物折反射摄像机 | |
CN105354839A (zh) | 使用相切的双球像与正交消失点标定拋物折反射摄像机 | |
CN110120078B (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 | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20190219 |