CN109345598A - 球与公共自极三角形的性质标定中心折反射摄像机的方法 - Google Patents

球与公共自极三角形的性质标定中心折反射摄像机的方法 Download PDF

Info

Publication number
CN109345598A
CN109345598A CN201811215376.9A CN201811215376A CN109345598A CN 109345598 A CN109345598 A CN 109345598A CN 201811215376 A CN201811215376 A CN 201811215376A CN 109345598 A CN109345598 A CN 109345598A
Authority
CN
China
Prior art keywords
ball
picture
polar
video camera
width
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
Application number
CN201811215376.9A
Other languages
English (en)
Other versions
CN109345598B (zh
Inventor
杨丰澧
赵越
汪雪纯
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Yunnan University YNU
Original Assignee
Yunnan University YNU
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Yunnan University YNU filed Critical Yunnan University YNU
Priority to CN201811215376.9A priority Critical patent/CN109345598B/zh
Publication of CN109345598A publication Critical patent/CN109345598A/zh
Application granted granted Critical
Publication of CN109345598B publication Critical patent/CN109345598B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/80Analysis of captured images to determine intrinsic or extrinsic camera parameters, i.e. camera calibration
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20068Projection on vertical or horizontal image axis
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30204Marker
    • G06T2207/30208Marker 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

球与公共自极三角形的性质标定中心折反射摄像机的方法
技术领域
本发明属于计算机视觉领域,涉及一种利用空间中一个球与公共自极三角形的性质求解中心折反射摄像机内参数的方法。
背景技术
计算机视觉的中心任务就是对图像进行理解,而它的最终目标是使计算机具有通过二维图像认知三维环境信息的能力。这种能力将不仅使机器能感知包括形状、姿态、运动等在内的三维环境中物体的几何信息,而且能对它们进行描述、存储、识别与理解。摄像机标定就是确定从三维空间点到它的二维图像点之间的映射关系,它是许多计算机视觉应用必不可少的步骤。为了确定这一映射过程,需要建立摄像机的几何成像模型,几何模型的参数称为摄像机参数,摄像机参数可分为内参数和外参数两类。内参数描述成像系统的成像几何特性,外参数描述成像系统关于世界坐标系的方向和位置。摄像机标定可分为传统标定、自标定和基于几何实体的标定。无论哪种标定方法,都旨在建立二维图像与摄像机内参数之间的约束关系,特别是线性约束关系,这是目前摄像机标定所追求的目标,也是目前计算机视觉领域研究的热点之一。
中心折反射摄像机由一个反射镜面和一个摄像机组成,它保持了投影视点的唯一性,并且它的成像视野大,是全景视觉领域研究的热点之一。文献“Catadioptric self-calibration”,(Kang S.B.,Proceedings of IEEE Conference on Computer Vision andPattern Recognition,vol.1,pp.201-207,2000)提出了一种折反射摄像机自标定方法,这类方法的优点是不需要使用标定块,缺点是必须获得图像之间的对应点。而在计算机视觉中,实现一个十分有效的寻找对应点的方法是很困难的。文献“Geometric properties ofcentral catadioptric line images and their application in calibration”,(Barreto J.P.,Araujo H.,IEEE Transactions on Pattern Analysis and MachineIntelligence,vol.27,no.8,pp.1327-1333,2005)研究了中心折反射摄像机下直线的像的几何性质,并将这些性质应用于中心折反射摄像机的标定。文献“A new linear algorithmfor calibrating central catadioptric cameras”,(Wu F.,Duan F.,Hu Z.et al.,Pattern Recognition,vol.41,no.10,pp.3166-3172,2008)介绍了对拓点和对拓像点,导出了空间中的一个点在视球上的投影和它的折反射图像点之间的关系,使用这个关系建立了中心折反射摄像机内参数的线性约束,通过此线性约束即可获得中心折反射摄像机内参数。文献“Calibration of central catadioptric cameras using a DLT-likeapproach”,(Puig L.,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)首次提出了利用球标定中心折反射摄像机。他们证明了球在中心折反射摄像机的单位球投影模型下的像为椭圆,并且在非退化情况下一个球的投影二次曲线提供两个不变量。为了降低求解的复杂度,他们提出了一种分步标定方法,该方法至少需要3个球的投影才能完成摄像机的标定。但是该文献提出的标定方法是非线性的,计算的复杂度较高。文献“Identical projective geometric properties of central catadioptricline images and sphere images with applications to calibration”,(Ying X.,ZhaH.,International Journal of Computer Vision,vol.78,no.1,pp.89-105,2008)介绍了修正的绝对二次曲线的像(MIAC)在中心折反射摄像机标定中的作用。他们通过研究球在中心折反射摄像机下的像与MIAC的几何与代数关系提出了两种线性标定算法。它们得出的结论对于对偶形式也是成立的。但是这篇文献中的理论和标定方法对于抛物折反射摄像机的情况是退化的。文献“A calibration method for paracatadioptric camera fromsphere images”,(Duan H.,Wu Y.,Pattern Recognition Letters,vol.33,no.6,pp.677-684,2012)基于圆环点理论提出了一种利用对拓球像标定抛物折反射摄像机的线性方法,但是这篇文献中关于圆环点的像的选取比较复杂。
发明内容
本发明提供了一种制作简单,适用广泛,稳定性好的利用靶标求解中心折反射摄像机内参数的方法,该靶标由空间中一个球构成。在求解中心折反射摄像机内参数的过程中,需使用中心折反射摄像机拍摄靶标的3幅图像求解出中心折反射摄像机的5个内参数。
本发明采用如下技术方案:
用中心折反射摄像机拍摄3幅含有一个球的图像。本发明是利用空间中一个球作为靶标用于求解中心折反射摄像机内参数的方法,其特征在于仅利用球元素。首先,从该幅图像中提取靶标图像边缘点,使用最小二乘法拟合获得球像的投影。其次,根据任意两幅球像的公共自极三角形的关系,从而得到两幅球像的公共自极三角形的顶点。然后由极点极线的关系,确定顶点对应的公共自极三角形的边。三幅图像提供三组公共自极三角形的顶点和边。接下来,利用公共自极三角形的顶点和边关于修正的绝对二次曲线的像的极点极线关系,从而确定修正的绝对二次曲线的像,通过分解修正的绝对二次曲线的像,得到修正的内参数矩阵。最后,利用修正的内参数矩阵和镜面参数确定摄像机内参数。具体的步骤包括:拟合靶标投影方程,估计两幅球像的公共自极三角形的顶点和边,确定修正的绝对二次曲线的像,求解中心折反射摄像机内参数。
1.拟合靶标投影方程
利用Matlab程序中的Edge函数提取靶标图像边缘点的像素坐标,并用最小二乘法拟合获得球像的方程。
2.估计两幅球像的公共自极三角形的顶点和边
空间中的球Q,在中心折反射摄像机的单位球模型下的投影分为两步。第一步,将球Q投影为以O为中心的单位视球上的小圆Sn(n=1,2,3表示拍摄的三幅图像),这里假设以单位视球中心O为原点的坐标系为世界坐标系;第二步,以单位视球内上的一点Oc为投影中心,这里Oc可看作一个摄像机的光心,其中单位视球中心O到投影中心Oc的距离为l=|OOc|,(0≤l<1)。然后将小圆Sn投影到垂直于光轴OcO的中心折反射图像平面上的二次曲线Cn,这里假设以投影中心Oc为原点的坐标系为摄像机坐标系。利用Matlab中的Edge函数提取每幅图像中的靶标图像边缘点的像素坐标,通过最小二乘法拟合得到相应的二次曲线方程。这里用Cn分别表示第n幅图像中的球像系数矩阵。本文为了简化表述,用相同字母表示曲线和它的系数矩阵。任意取两幅球像方程Ci,Cj(i,j=1,2,3,i≠j),通过矩阵Ci -1Cj的特征向量可以确定两幅球像Ci,Cj的公共自极三角形的顶点xij,则该顶点xij所对应的极线lij可由关系式lij=Ci·xij确定,其中·表示内积,lij,xij用齐次坐标矩阵表示。
3.确定修正的绝对二次曲线的像
若令单位视球上的小圆Sn的所在平面π1的单位法向量矩阵为[mnx mny mnz]T,该平面π1到单位视球的中心O的距离为d0,则在摄像机坐标系下投影中心Oc与小圆Sn形成的锥可以表示为若令以Oc为光心的摄像机的内参数矩阵为修正的内参数矩阵为其中l(0≤l<1)是镜面参数,表示单位视球中心O到投影中心Oc的距离,fx,fy分别表示在u轴和v轴方向上的尺度因子,s是畸变因子,[u0 v0 1]T是摄像机主点p的齐次坐标矩阵形式,fx,fy,s,u0,v0为摄像机的5个内参数。则球的成像方程Cn可以通过锥Qns与摄像机内参数矩阵Kc确定,即有关系式Cn=Kc -TQnsKc -1。通过化简,球像方程Cn满足关系式其中是修正的绝对二次曲线的像,λ是非零比例因子,根据射影几何的对偶理论,则在对偶空间中,球像的对偶Cn *满足关系式其中是修正的绝对二次曲线的像的对偶,其中λ*是非零比例因子,是点[mnx mny l·d0-mnz]T的投影。任意取两幅球像方程Ci,Cj(i,j=1,2,3,i≠j),则存在由点和点构成的直线使两幅球像的对偶Ci *,Cj *满足关系式则由上述关系式可知,lij *是Ci,Cj的公共自极三角形的边,xij *是Ci,Cj的公共自极三角形的顶点,即有lij *=lij,xij *=xij。因此,通过可以利用xij和lij关于修正的绝对二次曲线的像的极点极线的约束确定
4.求解中心折反射摄像机的内参数
根据进行Cholesky分解再求逆得到修正的内参数矩阵最后,利用镜面参数l(0≤l<1),通过关系式得到的内参数的矩阵Kc,即获得摄像机的5个内参数。
本发明优点:
(1)该靶标制作简单,只需将一个球固定在一个支架上。
(2)对该靶标的物理尺度没有要求,可线性的确摄像机内参数。
(3)该靶标的图像边界点几乎可以全部提取,这样可以提高曲线拟合的精确度,从而提高标定精度。
附图说明
图1是用于求解中心折反射摄像机内参数的靶标在单位视球上的示意图。
图2是靶标在中心折反射图像平面上的投影。
具体实施方式
本发明提供了一种利用靶标求解中心折反射摄像机内参数的方法,靶标是由空间中的一个球构成,如图1。用此靶标完成中心折反射摄像机内参数的求解需要经过以下步骤:从中心折反射图像中提取靶标图像边缘点,使用最小二乘法拟合获得球像方程。根据任意两幅球像的公共自极三角形的关系,从而得到两幅球像的公共自极三角形的顶点。然后由极点极线的关系,确定顶点对应的公共自极三角形的边。三幅图像提供三组公共自极三角形的顶点和边。接下来,利用公共自极三角形的顶点和边关于修正的绝对二次曲线的像的极点极线的约束,从而确定修正的绝对二次曲线的像,通过分解修正的绝对二次曲线的像,得到修正的内参数矩阵。最后,利用修正的内参数矩阵和镜面参数得到摄像机内参数矩阵。利用本发明中的方法对用于实验的中心折反射摄像机进行标定,具体步骤如下:
1.拟合靶标投影方程
利用Matlab程序中的Edge函数提取靶标图像边缘点的像素坐标,并用最小二乘法拟合获得球像的方程。
2.估计两幅球像的公共自极三角形的顶点和边
空间中的球Q(如图1),在中心折反射摄像机的单位球模型下的投影分为两步。第一步,将球Q投影到以O为中心的单位视球上的小圆Sn(n=1,2,3),在这里以O为原点建立世界坐标系O-xwywzw,其中zw轴与光轴OcO重合,Oc是单位视球内的一点,点Oc到世界坐标系原点O的距离为l=|OOc|(0≤l<1),xw,yw轴分别与中心折反射图像平面的水平和垂直的轴方向一致,这一过程如图1(以n=1为例)所示。第二步,以点Oc为原点建立投影坐标系Oc-xcyczc,其中zc轴与zw轴重合,xc,yc轴分别与xw,yw轴平行。又以Oc为投影中心,这里Oc可看作一个摄像机的光心,通过摄像机的光心Oc将这小圆Sn投影为中心折反射图像平面上的二次曲线Cn,其中中心折反射图像平面是垂直于光轴OcO的一个平面π,如图1所示(以n=1为例)。利用Matlab中的Edge函数提取每幅图像中的靶标图像边缘点的像素坐标,通过最小二乘法拟合得到相应的二次曲线方程,这里用Cn表示第n幅图像中的球像的系数矩阵。本文为了简化表述,用相同字母表示曲线和它的系数矩阵。若任意取两幅球像方程Ci,Cj(i,j=1,2,3,i≠j),根据公共自极三角形的定义可知,两幅球像Ci,Cj的公共自极三角形的顶点xij和边lij满足如下关系式:
其中,λij是非零比例因子,“·”表示点积。若消去(1)中的lij,即有:
其中I是3×3阶单位矩阵。则由(2)式可知,两幅球像Ci,Cj的公共自极三角形的顶点xij可以由矩阵Ci -1Cj的特征值所对应的特征向量来确定。通过极点极线的关系,即利用(1)式则可以确定公共自极三角形的顶点xij对应的边lij
3.确定修正的绝对二次曲线的像
如图1所示,若令单位视球上的小圆Sn(以n=1为例)的所在平面π1的单位法向量矩阵为[mnx mny mnz]T,其中下标n与小圆Sn的下标一致,且该平面π1到单位视球的中心O的距离为d0,则在世界坐标系O-xwywzw下平面π1的方程表示为:
mnxXW+mnyYW+mnzZW+d0=0。 (3)
因此,在世界坐标系O-xwywzw下小圆Sn上的点齐次坐标矩阵Mnw=[XnW YnW ZnW 1]T应该满足:
因为世界坐标系O-xwywzw和摄像机坐标系Oc-xcyczc之间的旋转矩阵R和平移向量T分别为R=I和T=[0 0 l]T,其中I为3×3阶单位矩阵,则在摄像机坐标Oc-xcyczc下的小圆Sn上点Mnc=[XnC YnC ZnC 1]T有:
其中λnC是非零比例因子。联立(4)和(5)式,通过消去XnW,YnW,ZnW和λnC,则可得方程:
其中是Mnc的非齐次坐标,而表示在摄像机坐标系下投影中心Oc与小圆Sn形成的锥。若令以Oc为光心的摄像机的内参数矩阵为修正的内参数矩阵为 其中l(0≤l<1)是镜面参数,即单位视球中心O到投影中心Oc的距离,fx,fy分别表示在u轴和v轴方向上的尺度因子,s是畸变因子,[u0 v0 1]T是摄像机主点p的齐次坐标矩阵形式。若点Mnc的中心折反射投影点为mnc,则有下式成立:
其中λm是非零比例因子,联立(6)和(7)式,则有:
mnc TKc -TQnsKc -1mnc=0。 (8)
则由(8)式可知,球像方程Cn可以通过锥Qns与摄像机内参数Kc确定,即:
Cn=Kc -TQnsKc -1。 (9)
通过化简(9)式,则有:
其中是修正的绝对二次曲线的像,λ是非零比例因子,根据射影几何的对偶理论,则在对偶空间中,由(10)式可知,球像的对偶Cn *满足:
其中是修正的绝对二次曲线的像的对偶,其中λ*是非零比例因子,是点[mnx mny l·d0-mnz]T的投影。若任意取两幅球像方程Ci,Cj(i,j=1,2,3,i≠j),如图2所示(以i=1,j=2为例)。利用(11)式可知,两幅球像的对偶Ci *,Cj *满足下式:
其中λi *j *是非零比例因子。因为分别是[mix miy l·d0-miz]T和[mjx mjy l·d0-mjz]T的投影点,则在中心折反射图像平面π上,存在一条由点和点构成的直线(如图2所示),其中“×”表示叉积。若同时在(12)式两端乘以直线则有:
若令根据对偶空间的定义可知Cn *∝Cn -1其中∝表示相差一个非零比例因子。则由上述关系式可知,lij *是球像Ci,Cj的公共自极三角形的边,xij *是球像Ci,Cj的公共自极三角形的顶点,即有lij *=lij,xij *=xij。由于已知lij *和xij *,因此,通过可以利用xij和lij关于修正的绝对二次曲线的像的极点极线的约束确定即有关系式:
4.求解中心折反射摄像机的内参数
首先,根据进行Cholesky分解再求逆得到修正的内参数矩阵最后,利用镜面参数l(0≤l<1),内参数矩阵Kc可以通过如下关系式确定:
即获得摄像机的5个内参数。
实施例
本发明提出了一种利用空间的一个球作为靶标线性确定中心折反射摄像机内参数的方法。本发明采用的实验模板结构示意图如图1所示。下面以一实例对本发明的实施方案做出更为详细的描述。
基于空间中球的中心折反射摄像机标定采用的实验模板是空间中的一个球,如图1所示,球为Q。利用本发明中的方法对用于实验的中心折反射摄像机进行标定,具体步骤如下:
1.拟合靶标曲线方程
用抛物折反射摄像机拍摄靶标的3幅实验图像,读入图像,利用Matlab中的Edge函数提取靶标图像边缘点的像素坐标,并用最小二乘法拟合获得球像的方程。3幅球像方程的系数矩阵分别为Cn(n=1,2,3),结果如下:
2.估计两幅球像的公共自极三角形的顶点和边
分别计算矩阵C2 -1C1,C1 -1C3,,C3 -1C2结果如下:
将(19)式带入(2)式,则可得球像C2,C1之间的公共自极三角形的顶点x21的齐次坐标矩阵,结果如下:
x21=[346.1379 -431.17891]T, (22)
将(20)式带入(2)式,则可得球像C1,C3之间的公共自极三角形的顶点x13的齐次坐标矩阵,结果如下:
x13=[-335.1375 -1319.27421]T, (23)
将(21)式带入(2)式,则可得球像C3,C2之间的公共自极三角形的顶点x32的齐次坐标矩阵,结果如下:
x32=[677.3112 272.21271]T。 (24)
将(22)式带入(1)式,即可确定顶点x21对应的边l21的齐次线坐标矩阵,结果如下:
l21=[0.00003283 -0.000742301]T, (25)
将(23)式带入(1)式,即可确定顶点x13对应的边l13的齐次线坐标矩阵,结果如下:
l13=[-0.00054218 -0.001143601]T, (26)
将(24)式带入(1)式,即可确定顶点x32对应的边l32的齐次线坐标矩阵,结果如下:
l32=[0.00065500 0.000052151]T。 (27)
3.确定修正的绝对二次曲线的像
将(22-27)式带入(14)式,即可确定修正的绝对二次曲线的像矩阵,结果如下:
4.求解中心折反射摄像机的内参数
根据对(28)式进行Cholesky分解再求逆得到修正的内参数矩阵有:
因为已知镜面参数l=0.2,则将(28)式带入(15)式,则可获得Kc,有:
故中心折反射摄像机的5个内参数分别为:fx=800.0000,fy=850.0000,s=0.2000,u0=320.0000,v0=240.0000。

Claims (1)

1.一种利用球与公共自极三角形的性质标定中心折反射摄像机的方法;本发明是利用空间中一个球作为靶标用于求解中心折反射摄像机内参数的方法,其特征在于仅利用球元素;首先,从该幅图像中提取靶标图像边缘点,使用最小二乘法拟合获得球像的投影;其次,根据任意两幅球像的公共自极三角形的关系,从而得到两幅球像的公共自极三角形的顶点;然后由极点极线的关系,确定顶点对应的公共自极三角形的边;三幅图像提供三组公共自极三角形的顶点和边;接下来,利用公共自极三角形的顶点和边关于修正的绝对二次曲线的像的极点极线关系,从而确定修正的绝对二次曲线的像,通过分解修正的绝对二次曲线的像,得到修正的内参数矩阵;最后,利用修正的内参数矩阵和镜面参数确定摄像机内参数;具体的步骤包括:拟合靶标投影方程,估计两幅球像的公共自极三角形的顶点和边,确定修正的绝对二次曲线的像,求解中心折反射摄像机内参数;
1)估计两幅球像的公共自极三角形的顶点和边
空间中的球Q,在中心折反射摄像机的单位球模型下的投影分为两步;第一步,将球Q投影为以O为中心的单位视球上的小圆Sn,其中n=1,2,3表示拍摄的三幅图像,这里假设以单位视球中心O为原点的坐标系为世界坐标系;第二步,以单位视球内上的一点Oc为投影中心,这里Oc看作一个摄像机的光心,其中单位视球中心O到投影中心Oc的距离为l=|OOc|,其中0≤l<1;然后将小圆Sn投影到垂直于光轴OcO的中心折反射图像平面上的二次曲线Cn,这里假设以投影中心Oc为原点的坐标系为摄像机坐标系;这里用Cn分别表示第n幅图像中的球像系数矩阵;本文为了简化表述,用相同字母表示曲线和它的系数矩阵;任意取两幅球像方程Ci,Cj,其中i,j=1,2,3,i≠j,通过矩阵Ci -1Cj的特征向量以确定两幅球像Ci,Cj的公共自极三角形的顶点xij,则该顶点xij所对应的极线lij由关系式lij=Ci·xij确定,其中·表示内积,lij,xij用齐次坐标矩阵表示;
2)确定修正的绝对二次曲线的像
若令单位视球上的小圆Sn的所在平面π1的单位法向量矩阵为[mnx mny mnz]T,该平面π1到单位视球的中心O的距离为d0,则在摄像机坐标系下投影中心Oc与小圆Sn形成的锥以表示为若令以Oc为光心的摄像机的内参数矩阵为修正的内参数矩阵为其中0≤l<1是镜面参数,表示单位视球中心O到投影中心Oc的距离,fx,fy分别表示在u轴和v轴方向上的尺度因子,s是畸变因子,[u0 v0 1]T是摄像机主点p的齐次坐标矩阵形式,fx,fy,s,u0,v0为摄像机的5个内参数;则球的成像方程Cn以通过锥Qns与摄像机内参数矩阵Kc确定,即有关系式Cn=Kc -TQnsKc -1;通过化简,球像方程Cn满足关系式其中是修正的绝对二次曲线的像,λ是非零比例因子,根据射影几何的对偶理论,则在对偶空间中,球像的对偶Cn *满足关系式其中是修正的绝对二次曲线的像的对偶,其中λ*是非零比例因子,是点[mnx mny l·d0-mnz]T的投影;任意取两幅球像方程Ci,Cj,其中i,j=1,2,3,i≠j,则存在由点和点构成的直线其中“×”表示叉积,使两幅球像的对偶Ci *,Cj *满足关系式则由上述关系式知,lij *是Ci,Cj的公共自极三角形的边,xij *是Ci,Cj的公共自极三角形的顶点,即有lij *=lij,xij *=xij;因此,通过利用xij和lij关于修正的绝对二次曲线的像的极点极线的约束确定其中“·”表示点积;
3)求解中心折反射摄像机的内参数
根据进行Cholesky分解再求逆得到修正的内参数矩阵最后,利用镜面参数l,其中0≤l<1,通过关系式得到的内参数的矩阵Kc,即获得摄像机的5个内参数。
CN201811215376.9A 2018-10-18 2018-10-18 球与公共自极三角形的性质标定中心折反射摄像机的方法 Expired - Fee Related CN109345598B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811215376.9A CN109345598B (zh) 2018-10-18 2018-10-18 球与公共自极三角形的性质标定中心折反射摄像机的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811215376.9A CN109345598B (zh) 2018-10-18 2018-10-18 球与公共自极三角形的性质标定中心折反射摄像机的方法

Publications (2)

Publication Number Publication Date
CN109345598A true CN109345598A (zh) 2019-02-15
CN109345598B CN109345598B (zh) 2021-06-08

Family

ID=65310184

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811215376.9A Expired - Fee Related CN109345598B (zh) 2018-10-18 2018-10-18 球与公共自极三角形的性质标定中心折反射摄像机的方法

Country Status (1)

Country Link
CN (1) CN109345598B (zh)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102999894A (zh) * 2011-09-16 2013-03-27 云南大学 利用两个相交圆线性求解摄像机内参数
CN106327504A (zh) * 2016-09-22 2017-01-11 云南大学 利用单个球及圆环点的像标定拋物折反射摄像机的方法
CN106558081A (zh) * 2016-11-28 2017-04-05 云南大学 标定非中心轴对称系统的圆锥折反射摄像机的方法
CN109035341A (zh) * 2018-07-12 2018-12-18 北京师范大学-香港浸会大学联合国际学院 一种基于平面同心圆的相机标定的计算方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102999894A (zh) * 2011-09-16 2013-03-27 云南大学 利用两个相交圆线性求解摄像机内参数
CN106327504A (zh) * 2016-09-22 2017-01-11 云南大学 利用单个球及圆环点的像标定拋物折反射摄像机的方法
CN106558081A (zh) * 2016-11-28 2017-04-05 云南大学 标定非中心轴对称系统的圆锥折反射摄像机的方法
CN109035341A (zh) * 2018-07-12 2018-12-18 北京师范大学-香港浸会大学联合国际学院 一种基于平面同心圆的相机标定的计算方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
HAIFEI HUANG等: "The common self-polar triangle of concentric circles and its application to camera calibration", 《2015 IEEE CONFERENCE ON COMPUTER VISION AND PATTERM RECOGNITION(CVPR)》 *
YUE ZHAO等: "Method of Camera Calibration Using Concentric Circles and Lines through Their Centres", 《ADVANCES IN MULTIMEDIA》 *
顾万家: "基于完全四点形利用空间球标定拋物折反射摄像机", 《中国优秀硕士学位论文全文数据库信息科技辑》 *

Also Published As

Publication number Publication date
CN109345598B (zh) 2021-06-08

Similar Documents

Publication Publication Date Title
WO2015096806A1 (zh) 智能机姿态测定、全景影像生成及目标识别方法
CN106327504A (zh) 利用单个球及圆环点的像标定拋物折反射摄像机的方法
CN107886546B (zh) 利用球像及公共自极三角形标定抛物折反射摄像机的方法
CN110288657A (zh) 一种基于Kinect的增强现实三维注册方法
CN106846417A (zh) 基于视觉里程计的单目红外视频三维重建方法
CN106447731A (zh) 利用单个球及正交消失点标定拋物折反射摄像机的方法
CN107657645B (zh) 利用直线及圆的共轭直径的性质标定拋物折反射摄像机
CN104217435A (zh) 两个相互遮挡的球线性确定拋物折反射摄像机内参数方法
CN102982551B (zh) 空间三条不平行直线线性求解抛物折反射摄像机内参数
CN109360248A (zh) 利用单个球及共轭直径的性质标定拋物折反射摄像机
CN103106661B (zh) 空间二条相交直线线性求解抛物折反射摄像机内参数
CN107145224B (zh) 基于三维球面泰勒展开的人眼视线跟踪方法和装置
CN109163657A (zh) 一种基于双目视觉三维重建的圆形目标位姿检测方法
CN109523595A (zh) 一种建筑工程直线棱角间距视觉测量方法
CN109325983A (zh) 利用无穷远点关于圆极线性质标定拋物折反射摄像机
CN107644445A (zh) 利用单个球及圆切线的性质标定拋物折反射摄像机的方法
CN105279758A (zh) 使用双球的相切像与圆环点的像标定拋物折反射摄像机
CN105321181A (zh) 使用双球的相离像与圆环点的像标定拋物折反射摄像机
CN103106662A (zh) 空间一条直线线性求解抛物折反射摄像机内参数
Liu et al. Dense stereo matching strategy for oblique images that considers the plane directions in urban areas
CN108921904B (zh) 利用单个球及渐近线的性质标定针孔摄像机的方法
CN105303570A (zh) 使用双球的相交像及圆环点的像标定拋物折反射摄像机
CN107958468B (zh) 利用空间位置不同的三个球标定中心折反射摄像机的方法
CN109360247A (zh) 单个球的公共自极三角形及正交消失点标定抛物摄像机
CN109523598A (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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20210608

Termination date: 20211018

CF01 Termination of patent right due to non-payment of annual fee