CN107886546B - 利用球像及公共自极三角形标定抛物折反射摄像机的方法 - Google Patents
利用球像及公共自极三角形标定抛物折反射摄像机的方法 Download PDFInfo
- Publication number
- CN107886546B CN107886546B CN201711075273.2A CN201711075273A CN107886546B CN 107886546 B CN107886546 B CN 107886546B CN 201711075273 A CN201711075273 A CN 201711075273A CN 107886546 B CN107886546 B CN 107886546B
- Authority
- CN
- China
- Prior art keywords
- circle
- image
- sphere
- center
- small circle
- 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
- G06T7/85—Stereo 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/20—Special algorithmic details
- G06T2207/20068—Projection on vertical or horizontal image axis
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
本发明涉及利用球像及公共自极三角形标定抛物折反射摄像机的方法,利用空间中的一个球作为标定物,在抛物折反射摄像机下,拟合出镜面轮廓投影方程、球像方程,并计算出对拓球像方程和圆心与单位视球球心重合且与单位视球模型上球的投影小圆平行的大圆的像的方程,由大圆的像分别与球像,对拓球像构成的两个公共自极三角形所对应的两个公共极点,连接两公共极点求解出消失线,从而利用圆的投影与消失线的关系求出圆环点的像,通过拍摄三幅图像获取三对圆环点的像,利用圆环点的像与摄像机内参数矩阵的约束完成该摄像机的标定。
Description
技术领域
本发明属于计算机视觉领域,涉及一种利用空间中的一个球及公共自极三角形求解抛物折反射摄像机内参数的方法。
背景技术
计算机视觉是通过计算机模拟人类视觉,从而对目标进行识别、检测、跟踪。计算机视觉的任务包括获取,处理,分析和理解数字图像,应用图像处理和分析技术,将从空间中获取的三维信息转化到二维平面上,通过分析二维平面图像的信息恢复三维空间信息,使计算机能够从二维图像认知三维世界,并由此识别、检测物体。
随着计算机视觉技术的发展,对摄像机也提出了更高的要求。传统摄像机由于可视范围小、获取信息量少的缺点,难以满足视觉技术的要求,在当今科学技术发展的浪潮下,中心折反射摄像机应运而生,原因是能提供更大的视野和信息而被广泛运用到计算机视觉的各个领域。中心折反射摄像机是由传统摄像机和反射镜面构成,文献“A Theory ofSingle-Viewpoint Catadioptric Image Formation”,(Baker S,Nayer K.,International Journal of Computer Vision,35(2):175-196,1999)将反射镜面分为四种类型,分别为平面镜,抛物面镜,椭球面镜以及双曲面镜。其类型由镜面参数决定:当镜面参数为0时为平面镜,当镜面参数在0到1之间时为双曲面镜或椭球面镜,当镜面参数为1时为抛物面镜。文献“Catadioptric projective geometry”,(Geyer C,Daniilidis K,International Journal of Computer Vision,45(3):223-243,2001)证明了中心折反射摄像机的成像可分为两步投影过程:首先,空间点到单位球上的投影,其次,空间点在单位球上的投影通过虚拟摄像机光心到图像平面上的成像,从而统一了针孔摄像机、抛物折反射摄像机、椭球折反射摄像机和双曲折反射摄像机的成像模型。
传统摄像机标定方法具有精度较高的优点,但同时操作的难度较大。文献“Aflexible new technique for camera calibration”,(Zhang Z.,IEEE Transactions onPattern Analysis and Machine Intelligence,22(11):1330-1334,2000)利用平面模板而非传统摄像机标定块进行标定,该方法简单、成本较低,同时在确定点的坐标时提供了较高的精度。文献“Stereo from uncalibrated cameras”,(Hartley R.,Proceedings ofIEEE Conference on Computer Vision and Pattern Recognition,62(2):761-7641992)利用平面模板和它的像之间的二次曲线对应的方法进行摄像机标定。由于二次曲线更简单,能提高摄像机标定的稳定性,因此,利用二次曲线进行摄像机标定更广泛。文献“Anew easy camera calibration technique based on circular points”,(Meng X,HuZ.,Pattern Recognition,36(5):1155-1164,2000)利用一个平面上的圆和经过这个圆心的直线作为平面模板,利用圆环点和绝对二次曲线求解摄像机的内参数,该方法首次将圆环点以及绝对二次曲线应用到摄像机标定中。文献“一种新的基于消失点的摄像机标定方法”,(胡培成,黎宁,赵亮亮.,第十三届全国图象图形学学术会议论文集,2006)中利用消失点(也就是无穷远点的像)对绝对二次曲线的约束来估计摄像机内参数,文献“Cameracalibration based on the common self-polar triangle of sphere images”,(HuangH.,Zhang H.,Cheung Y.M.,Lecture Notes in Computer Science,9004:19-29,2015)利用空间球标定针孔摄像机。
因为球有外轮廓,且具有从任意位置观察都是可见的和自身遮挡的优势,所以利用球作为参照物标定中心折反射摄像机成为了当前研究的一个热点,并由此提出了很多标定的方法。文献“基于平面和直线的仿射重构”,(胡占义,王光辉,吴福朝.,计算机学报,26(6):722-728,2003)中首次利用球来计算两幅图像的纵横比,文献“Camera Calibrationfrom Spheres Images”,(Daucher N.,Dhome M.,LaprestéJ.T.,Proceedings of ThirdEuropean Conference on Computer Vision,1994)提出了关于球的非线性方法进行摄像机标定,然而这种方法误差较大。文献“Camera Calibration by a single Image ofBalls:From Conics to the Absolute Conic”,(Teramoto H.,Xu G.,Proceedings ofFifth Asian Conference on Computer Vision,2002)提出利用非线性优化的方法把绝对二次曲线投影与球的像的代数关系进行摄像机的标定,该方法比较复杂,计算量大。文献“Camera calibration using spheres:A semi-definite programming approach”,(Agrawal M.,Davis L.S.,Proceedings of IEEE International Conference onComputer Vision,2003)利用球来计算摄像机的外参数。文献“Geometricinterpretations of the relation between the image of the absolute conic andsphere images”,(Ying X.,Zha H.,IEEE Transactions on Pattern Analysis andMachine Intelligence,28(12):2031-6,2006)提出利用双接触理论求解摄像机内参数,该方法精确度较高。文献“A stratified approach for camera calibration usingspheres”,(Wong K.Y.,Zhang G.,Chen Z.,IEEE Transactions on Image Processing,20(2):305-316,2011)利用单位球球心和球面的切线旋转形成圆锥求两幅图像的基础矩阵,通过基础矩阵求解摄像机的内参数。文献“Camera Calibration from Images ofSpheres”,(Zhang H,Wong K.Y.K.,Zhang G.,IEEE Transactions on Pattern Analysisand Machine Intelligence,29(3):499-502,2007)利用球像的对偶关系和绝对二次曲线的投影之间的代数关系求解折反射摄像机内参数。文献“秩1约束下基于圆球的摄像机标定算法”,(贾静,吴成柯,西安电子科技大学学报(自然科学版),40(2):1-4,2013)提出利用秩1约束下基于单位球的摄像机标定方法,将平面与单位球结合起来,这也是一个线性算法,但计算过程中需要优化,计算量较大。文献“Catadioptric camera calibration usinggeometric invariants”,(Ying X.,Hu Z.,IEEE Transactions on Pattern Analysisand Machine Intelligence,26(10):1260-71,2004)利用直线和球的几何不变性标定折反射摄像机,该方法较为简单,但只能求出其中一部分内参数。文献“Geometric propertiesof central catadioptric line images”,(Barreto J.P.,Araujo H.,IEEETransactions on Pattern Analysis and Machine Intelligence,27(8):1327-33,2005)研究了空间中任意直线成像模型下成像的性质来求解其内参数。文献“A ball-target-based extrinsic calibration technique for high-accuracy 3-D metrology usingoff-the-shelf laser-stripe sensors”,(Che C.,Ni J.,Precision Engineering,24(3):210-219,2000)提出利用球作为标定物来求解摄像机外参数。
摄像机标定的目的是求解摄像机内参数和外参数,是三维重建技术的基础。三维重建技术通过将三维空间信息映射到二维图像平面上,然后从获取的图像中找出这些点在像平面上的投影过程,并确定这些图像间的运动参数,求出摄像机内参数以及外参数,再从图像平面恢复三维立体信息。摄像机标定的结果对三维重建的结果有极为重要的影响,因此对摄像机标定的研究具有一定的理论意义和实用价值。
发明内容
本发明利用空间中的一个球作为标定物,在抛物折反射摄像机下,拟合出镜面轮廓投影方程、球像方程,并计算出对拓球像方程和圆心与单位视球球心重合且与单位视球模型上球的投影小圆平行的大圆的像的方程,由大圆的像分别与球像,对拓球像构成的两个公共自极三角形所对应的两个公共极点,连接两公共极点求解出消失线,从而利用圆的投影与消失线的关系求出圆环点的像,通过拍摄三幅图像获取三对圆环点的像,利用圆环点的像与摄像机内参数矩阵的约束完成该摄像机的标定。
本发明采用如下技术方案:
用抛物折反射摄像机从3个不同的位置拍摄一个空间球的3幅图像。本发明是利用空间中的一个球作为标定物用于求解抛物折反射摄像机内参数的方法,其特征在于仅利用了球元素。首先,利用Canny算子分别从3幅图像中提取空间球图像边缘点和其中一幅图像的镜面轮廓投影的边缘点,然后对获取的边缘点利用最小二乘法拟合出镜面轮廓投影方程、球像方程。其次,根据像点与其对拓像点的关系,获得对拓像点,从而拟合出球像的对拓球像方程。接着,根据抛物折反射摄像机的单位视球模型,球在其单位视球上的见到投影为小圆S+,不见到投影为小圆的对拓圆S-,圆心与单位视球球心重合且与小圆平行的单位视球球面上的圆称为大圆S,因为大圆与小圆平行,所以小圆所在平面的单位法向量就是大圆的单位法向量,可求出圆心与单位视球球心重合且与小圆平行的大圆的像的方程。因为小圆和大圆的圆心在同一条直线上,因此它们圆心的像也在同一条直线上。射影几何中,图像平面上两不相交的二次曲线有四个不同的交点时,它们有且仅有一个公共自极三角形。从而大圆的像分别与球像,对拓球像各有且仅有一个公共自极三角形。根据公共自极三角形顶点的性质,求出大圆的像分别与球像,对拓球像构成的公共自极三角形所对应的公共极点,连接两公共极点求解出消失线。消失线与球像的交点,即为圆环点的像。一幅图像所对应的圆环点的像对绝对二次曲线的像提供两个约束,三幅图像提供六个约束。最后,利用圆环点的像对绝对二次曲线的像的约束求解抛物折反射摄像机内参数。具体的步骤包括:拟合出镜面轮廓投影方程、球像投影方程,估计球像的对拓球像,求解圆心与单位视球球心重合且与对拓小圆平行的大圆的像的方程,最后确定圆环点的像,求解抛物折反射摄像机内参数。
1.拟合镜面轮廓投影方程及靶标投影方程
利用MATLAB中的Canny算子对拍摄的每幅图像进行边缘检测,提取其边缘点,然后对获取的边缘点利用最小二乘法拟合得到镜面轮廓投影方程、球像方程。
2.估计球像的对拓球像
一个空间球Q在拋物折反射摄像机的单位视球模型下的投影过程分为两步。第一步,球Q以单位视球球心O为中心透视投影到单位视球球面上形成两个平行小圆Si+和Si-(i=1,2,3表示拍摄三幅图像),并且Si+和Si-以单位视球球心O对称,单位视球直径的两个端点互为对拓点。第二步,以单位视球上的点Oc为摄像机的光心将单位视球上的两个平行小圆Si+和Si-分别投影为图像平面上的两条二次曲线Ci+和Ci-,其中可见的二次曲线Ci+称为空间球Q的像,不可见的二次曲线Ci-称为球像Ci+的对拓球像。令以Oc为光心的摄像机的内参数矩阵为其中r是纵横比,f是有效焦距,s是倾斜因子,[u0 v0 1]T是摄像机主点p的齐次坐标矩阵形式。利用MATLAB中的Canny算子提取第一幅图像中的镜面轮廓投影边缘点和3幅空间球在图像平面上的像的边缘点的像素坐标,通过最小二乘法拟合得到相应的二次曲线方程。这里用C0表示第1幅图像镜面轮廓的系数矩阵,Ci+分别表示第i幅图像中的球像的系数矩阵,Ci表示第i幅图像中的圆心与单位视球球心重合且与对拓小圆平行的大圆的像的系数矩阵。本文为了简化表述,用相同字母表示曲线和它的系数矩阵。通过C0可得到折反射摄像机内参数矩阵K的一个初始值从而得到绝对二次曲线的像ω的初始值这里:ω=K-TK-1,其中 是纵横比的初始值,是有效焦距的初始值,是倾斜因子的初始值,是折反射摄像机主点的初始齐次坐标矩阵形式,记为取Ci+(i=1,2,3)上的一组点则与它相对应的一组对拓像点为因为本文只讨论抛物折反射摄像机的情况,从而关于对拓像点的关系式如下:这里 用齐次坐标表示。根据对拓像点的定义,点在球像Ci+(i=1,2,3)的对拓球像Ci-(i=1,2,3)上,因此用最小二乘法拟合得到对拓球像Ci-的方程。
3.求解圆心与单位视球球心重合且与对拓小圆平行的大圆的像的方程
取单位视球上小圆S+上任意两点的坐标矩阵为[xs1 ys1 zs1]T和[xs2 ys2 zs2]T,小圆S+上的点与该小圆所在平面的法向量[n′x n′y n′z]T的向量积为0,同时这两点也位于单位视球上,因此小圆S+所在平面的法向量满足方程:
其中d0为单位视球球心O到小圆S+所在平面(基础平面)的距离,这两个方程组可求出法向量[n′x n′y n′z]T和d0,从而可以得到小圆S+所在平面的单位法向量[nx ny nz]T:
令圆心与单位视球球心重合并在单位视球球面上所形成的大圆S与小圆S+及对拓小圆S-平行,所以小圆S+所在平面的单位法向量就是大圆S所在平面的单位法向量。从而由:其中C'是圆心与单位视球球心重合且与小圆平行的大圆S的系数矩阵,求得圆心与单位视球球心重合且与小圆平行的大圆S的像的方程系数矩阵C。
4.确定圆环点的像
一个空间球Q在单位视球上的投影为两平行小圆,称为对拓小圆S+和S-,圆心与单位视球球心重合且与这两个小圆平行的圆,称为大圆S。由于大圆S与对拓小圆S+,S-平行,因此大圆S与对拓小圆S+,S-各有四个交点且这四个点均为虚点,由射影变换保持结合性知,大圆S和对拓小圆S+,S-在像平面上的投影为C和C+,C-,它们也各有四个交点,从而C与C+、C-分别有且仅有一个公共自极三角形。因为大圆S和对拓小圆S+,S-的圆心在同一条直线上,因此它们圆心的像也在同一条直线上,并且小圆S+的像C+位于C内,小圆S-的像C-位于C外。ΔA1B1C1为C与C+的公共自极三角形,A1,B1,C1为自极三角形的顶点,ΔA2B2C2为C和C-的公共自极三角形,A2,B2,C2为自极三角形顶点。根据公共自极三角形顶点的性质,可以求出C分别与C+、C-构成的公共自极三角形所对应的公共极点,连接两公共极点就求解出了消失线。消失线与球像的交点,即为圆环点的像,分别记为mI和mJ。5.求解抛物折反射摄像机的内参数
根据射影几何的性质,无穷远直线与其所在的平面上的所有圆相交于一对圆环点I,J。圆环点I,J在摄像机下的投影称为圆环点的像,分别记为mI和mJ,它们是一对共轭复点。由圆环点的定义可知,圆环点I,J在绝对二次曲线C∞上,因此C∞可以看作空间所有圆环点的集合。由于在射影空间中射影变换保持结合性不变,因此圆环点的像mI,mJ在绝对二次曲线的像ω上,于是可得到下述关于ω的约束方程
由于mI,mJ是两共轭复点,因此在上述方程中只能提供下述两个关于绝对二次曲线的像ω的实线性约束方程
其中Re,Im分别表示复数的实部和虚部。由实线性约束方程可以求出绝对二次曲线的像ω,然后利用Cholesky分解法对ω进行分解可唯一确定K-1,再对K-1求逆便得到摄像机的内参数矩阵K。
本发明优点:
(1)标定物制作简单,成本低廉,只需将一个球固定在一个支架上。
(2)对标定物的大小没有要求,无需知道任何物理度量。
(3)球的图像边界点几乎可以全部提取,这样可以提高曲线拟合的精确度,从而提高标定精度。
附图说明
图1是用于求解抛物折反射摄像机内参数的靶标在单位视球上的示意图。
图2是靶标在抛物折反射图像平面上的投影。
具体实施方式
本发明提供了一种利用空间中的一个球求解抛物折反射摄像机内参数的方法,如图1。用空间的一个球完成抛物折反射摄像机内参数的求解需要经过以下步骤:从折反射图像中提取镜面轮廓投影边缘点、球像的边缘点,然后使用最小二乘法拟合获得镜面轮廓投影方程、球像方程,根据像点及对拓像点之间的关系获得对拓像点,从而求解出球像的对拓球像方程。再将球像与对拓球像反投影到单位视球上,求出圆心与单位视球球心重合且与对拓小圆平行的大圆的像的方程。ΔA1B1C1为C和C+的公共自极三角形,A1,B1,C1为公共自极三角形ΔA1B1C1的顶点,ΔA2B2C2为C和C-的公共自极三角形,A2,B2,C2为公共自极三角形
ΔA2B2C2的顶点。根据求公共自极三角形顶点的性质,求出大圆的像分别与球像,对拓球像构成的公共自极三角形所对应的公共极点,连接两公共极点求解出了消失线。消失线与球像的交点,即为圆环点的像。一幅图像所对应的圆环点的像对绝对二次曲线的像提供两个约束,三幅图像提供六个约束。最后,利用圆环点的像对绝对二次曲线的像的约束求解抛物折反射摄像机内参数。利用本发明中的方法对实验中的拋物折反射摄像机进行标定,具体步骤如下:
1.拟合镜面轮廓投影方程、球像方程
利用MATLAB程序中的Canny算子对拍摄的每幅图像进行边缘检测,提取其边缘点,然后对获取的边缘点利用最小二乘法拟合得到镜面轮廓投影方程、球像方程。
2.拟合球像的对拓球像
一个空间球Q在拋物折反射摄像机的单位视球模型下的投影分为两步:第一步,将球Q以O为中心透视投影到单位视球球面上形成两个平行小圆Si+和Si-,称Si-为Si+的对拓圆,(i=1,2,3),这一过程如图1所示(i=1为例);第二步,以单位视球上的点Oc为摄像机的光心将单位视球上的两个平行小圆Si+和Si-分别投影为图像平面上的两条二次曲线Ci+和Ci-,其中可见的二次曲线Ci+称为空间球Q的像,不可见的二次曲线Ci-称为球像Ci+的对拓球像。利用MATLAB中的Canny算子分别提取第1幅图像中的镜面轮廓投影的像素坐标和空间球在3幅图像中投影的图像边缘点,通过最小二乘法拟合得到相应的二次曲线方程,这里用C0表示第1幅图像镜面轮廓投影曲线的系数矩阵,Ci+表示第i幅图像中的球像的系数矩阵,Ci-表示第i幅图像中的球像的对拓球像的系数矩阵,通过C0可获得摄像机内参数矩阵K的一个初始值具体如公式(1):
3.求解圆心与单位视球球心重合且与对拓小圆平行的大圆的像的方程
取小圆S+上任意两点的坐标矩阵为[xs1 ys1 zs1]T和[xs2 ys2 zs2]T,小圆S+上的点与该小圆所在平面的法向量[n′x n′y n′z]T的向量积为0,同时这两点也位于单位视球上,因此小圆S+所在平面的法向量满足方程:
其中d0为球心O到小圆S+所在平面(基础平面)的距离,这两个方程组可求出法向量[n′x n′y n′z]T和d0,从而可以得到小圆S+所在平面的单位法向量[nx ny nz]T:
因为圆心与单位视球球心重合且与对拓小圆平行的大圆所在平面与小圆S+所在平面平行,所以小圆S+所在平面的单位法向量就是大圆S所在平面的单位法向量。
取圆心与单位视球球心重合且与对拓小圆平行的大圆S上任一点x在世界坐标系下的非齐次坐标矩阵为x=[xs ys zs]T,则大圆S的方程为xTC′x=0,点x在抛物折反射像平面上的投影的齐次坐标矩阵为m=[x's y's z's]T,又所以有从而圆心与单位视球球心重合且与对拓小圆平行的大圆S的像的方程为:其中C′是圆心与单位视球球心重合且与对拓小圆平行的大圆S的系数矩阵。
4.确定圆环点的像
一个空间球Q在单位视球上的投影为两个平行小圆,可见的为S+,不可见的为S-,又称为小圆S+和对拓小圆S-,圆心与单位视球球心重合且与小圆S+,S-平行的圆,称为大圆S,如图1。由于大圆S和小圆S+,S-平行,因此大圆S和小圆S+,S-各有四个交点且这四个点均为虚点,由公共自极三角形的定义知,如果两条二次曲线有四个交点,那么它们有且仅有一个公共自极三角形。从而在图1中,大圆S分别与小圆S+,S-有且仅有一个公共自极三角形,这两个公共自极三角形的顶点(不与小圆S+,S-的圆心在一条直线)在无穷远直线上,无穷远直线与二次曲线的交点,就称为圆环点。
如图2所示,由射影变换保持结合性知,大圆S和小圆S+,S-在像平面上的投影分别为C和C+,C-,其中C-称为C+的对拓球像,它们也各有四个交点,从而C分别与C+,C-各有且仅有一个公共自极三角形。因为大圆S和对拓小圆S+,S-的圆心在同一条直线上,它们圆心的像也在同一条直线上,并且小圆S+的像C+位于C内,小圆S-的像C-位于C外。ΔA1B1C1为C与C+的公共自极三角形,A1,B1,C1为公共自极三角形ΔA1B1C1的顶点,ΔA2B2C2为C和C-的公共自极三角形,A2,B2,C2为公共自极三角形ΔA2B2C2的顶点。因为C与C+有且仅有一个公共自极三角形,所以它们有公共极点A1即三角形的其中一个顶点和极线l,那么极点和极线满足如下代数表达式:
其中λ为比例参数。整理方程(4)可得:
同理可得方程:
ax+by+c=0, (8)
其中a,b,c为直线的系数,则消失线V与球像的交点即为圆环点的像,分别记为mI和mJ。
5.求解抛物折反射摄像机内参数
根据圆环点的像对绝对二次曲线的像的线性约束有:
从而用SVD分解优化求解(9)获得ω,然后对ω=K-TK-1进行Cholesky分解得K-1,再求逆便得到摄像机内参数矩阵K,即求出了摄像机的5个内参数。
实施例
本发明提出了一种利用空间的一个球作为标定物线性求解抛物折反射摄像机内参数的方法。本发明采用的实验模板结构示意图如图1所示。下面用一实例对本发明的实施方案做出更为详细的描述。
抛物折反射摄像机标定采用的实验模板是空间中的一个球,如图1所示,球为Q。利用本发明中的方法对用于实验的抛物折反射摄像机进行标定,具体步骤如下:
1.拟合镜面轮廓投影方程、球像方程
本发明采用的图像大小为1590×1296。用抛物折反射摄像机拍摄空间中一个球的3幅不同位置的图像,读入图像,利用MATLAB中的Canny算子对拍摄的每幅图像进行边缘检测,提取其边缘点,然后对获取的边缘点利用最小二乘法拟合得到镜面轮廓投影方程、球像方程。镜面轮廓投影方程的系数矩阵为C0,3幅图像中球像方程的系数矩阵分别为Ci+(i=1,2,3),其结果如下:
2.拟合球像的对拓球像
先分别在球像C1+、C2+和C3+上取至少5个互异的点,再将所取的点和(14)分别带入(3)得对拓球像上点的坐标,通过最小二乘法拟合获得对拓球像C1-、C2-和C3-的估计,系数矩阵结果分别如下:
3.求解圆心与单位视球球心重合且与对拓小圆平行的大圆的像的方程
由具体实施方式中的第三步:C′是圆心与单位视球球心重合且与对拓小圆平行的大圆S的系数矩阵,求得3幅图像中圆心与单位视球球心重合且与对拓小圆平行的大圆的像的方程的系数矩阵为Ci(i=1,2,3),系数矩阵的结果分别如下:
4.确定圆环点的像
将(11)和(19)代入(6)可得极点A1的齐次坐标矩阵为:
A1=[0.790726975608435 0.612168971808550 1.000000000000000]; (22)将(15)和(19)代入(7)可得极点A2的齐次坐标矩阵为:
A2=[-0.790726975608436 -0.612168971808549 1.000000000000000]。 (23)连接点A1和点A2的直线就是消失线V,则由(8)得消失线V的齐次线坐标矩阵为:
V=[-0.001514985463743 0.001956877804991 1.000000000000000]。 (24)从而消失线V与二次曲线C1+的交点就是圆环点的像,齐次坐标矩阵如下:
m1I=[1009.933407145656+429.363513521747i 270.857193955645+332.406796166864i 1], (25)
m1J=[1009.933407145656-429.363513521747i 270.857193955645-332.406796166864i 1]; (26)同理由(12)(20)和(16)(20)可求得球像C2+上的圆环点的像,齐次坐标矩阵如下:
m2I=[1463.654060077871+359.573543039835i 254.859133006675+453.537133847140i 1], (27)
m2J=[1463.654060077871-359.573543039835i 254.859133006675-453.537133847140i 1]; (27)同理由(13)(21)和(17)(21)可求得球像C3+上的圆环点的像,齐次坐标矩阵如下:
m3I=[-1508.706260377841+506.076938380817i 1457.862648165625+441.283181925061i 1], (28)
m3J=[-1508.706260377841-506.076938380817i 1457.862648165625-441.283181925061i 1]。 (29)
5.求解抛物折反射摄像机内参数
将(25)、(28)和(29)带入(9)得到ω中元素的线性方程组,使用SVD分解求解该线性方程组得到ω的系数矩阵。结果如下:
最后,对(31)中的ω进行Cholesky分解再求逆便可获得矩阵K,结果如下:
其中纵横比r=K(1,1)/K(2,2)(K(1,1)表示矩阵K的第1行第1列的元素,K(2,2)表示矩阵K的第2行第2列的元素),故抛物折反射摄像机的5个内参数分别为:r=1.51709255,f=408.6713599724734,s=0.836521765234959,u0=647.917040048090,v0=794.8259823712548。
Claims (1)
1.一种利用空间中的一个球及公共自极三角形的性质标定抛物折反射摄像机内参数的方法,其特征是把空间中的一个球作为标定物;该方法的具体步骤为:首先,用抛物折反射摄像机从不同的位置拍摄3幅球的图像,提取镜面轮廓投影的边缘点、球像的边缘点,然后用最小二乘法拟合得到镜面轮廓投影方程、球像方程;其次,根据像点及对拓像点之间的关系获得对拓像点,从而拟合出球像的对拓球像方程;接着,根据抛物折反射摄像机的单位视球模型,球在单位视球上的见到投影为小圆,表示为S+,不见到投影为小圆S+的对拓圆,表示为S-,圆心与单位视球球心重合且与小圆S+平行的单位视球球面上的圆称为大圆,表示为S,因为S与S+及S-平行,所以小圆S+所在平面的单位法向量就是大圆S所在平面的单位法向量,求解出圆心与单位视球球心重合且与小圆S+平行的大圆S的像的方程;根据球公共自极三角形顶点的性质,求出大圆S的像分别与球像,对拓球像构成的公共自极三角形所对应的公共极点,连接两公共极点求解出了消失线;消失线与球像的交点,即为圆环点的像;圆环点的像mI,mJ在绝对二次曲线的像ω上,由于mI,mJ是两共轭复点,一幅图像只能提供两个关于ω的实线性约束,因此三幅图像的实线性约束求出ω,然后利用Cholesky分解法对ω进行分解以确定K-1,再对K-1求逆便得到摄像机的内参数矩阵K;
(1)求解圆心与单位视球球心重合且与小圆S+平行的大圆S的像的方程
取单位视球上小圆S+上任意两点的坐标矩阵分别为[xs1 ys1 zs1]T和[xs2 ys2 zs2]T,小圆S+上的点与该小圆S+所在平面的法向量[n′x n′y n′z]T的向量积为0,同时这两点也位于单位视球上,因此小圆S+所在平面的法向量满足方程:
其中d0为单位视球球心O到小圆S+所在平面的距离,这两个方程组求出法向量[n′x n′yn′z]T和d0,从而得到小圆S+所在平面的单位法向量[nx ny nz]T:
令圆心与单位视球球心重合并在单位视球球面上所形成的大圆S与小圆S+及对拓小圆S-平行,所以小圆S+所在平面的单位法向量就是大圆S所在平面的单位法向量;通过镜面轮廓投影拟合得到折反射摄像机内参数矩阵K的一个初始值从而由:其中C'是圆心与单位视球球心重合且与小圆S+平行的大圆S的系数矩阵,求得圆心与单位视球球心重合且与小圆S+平行的大圆S的像的方程系数矩阵C”;
(2)确定圆环点的像
一个空间球Q在单位视球上的投影为两平行小圆,小圆S+和对拓小圆S-,圆心与单位视球球心重合且与这两个小圆平行的圆,称为大圆S;由于大圆S与小圆S+及对拓小圆S-平行,因此大圆S与小圆S+及对拓小圆S-各有四个交点且这四个点均为虚点,由射影变换保持结合性知,大圆S和小圆S+及对拓小圆S-在像平面上的投影为C和C+,C-,它们也各有四个交点,从而C与C+、C-分别有且仅有一个公共自极三角形;因为大圆S和小圆S+及对拓小圆S-的圆心在同一条直线上,因此它们圆心的像也在同一条直线上,并且小圆S+的像C+位于C内,对拓小圆S-的像C-位于C外;ΔA1B1C1为C与C+的公共自极三角形,A1,B1,C1为自极三角形的顶点,ΔA2B2C2为C和C-的公共自极三角形,A2,B2,C2为自极三角形顶点;根据公共自极三角形顶点的性质,以求出C分别与C+、C-构成的公共自极三角形所对应的公共极点,连接两公共极点就求解出了消失线;消失线与球像的交点,即为圆环点的像,分别记为mI和mJ。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711075273.2A CN107886546B (zh) | 2017-11-06 | 2017-11-06 | 利用球像及公共自极三角形标定抛物折反射摄像机的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711075273.2A CN107886546B (zh) | 2017-11-06 | 2017-11-06 | 利用球像及公共自极三角形标定抛物折反射摄像机的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107886546A CN107886546A (zh) | 2018-04-06 |
CN107886546B true CN107886546B (zh) | 2021-07-30 |
Family
ID=61778596
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711075273.2A Active CN107886546B (zh) | 2017-11-06 | 2017-11-06 | 利用球像及公共自极三角形标定抛物折反射摄像机的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107886546B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109215088B (zh) * | 2018-09-04 | 2021-07-09 | 云南大学 | 利用公共自极三角形标定针孔摄像机的方法 |
CN109523598A (zh) * | 2018-11-28 | 2019-03-26 | 云南大学 | 单个球的公共自极三角形及圆环点标定抛物摄像机的方法 |
CN109360247B (zh) * | 2018-11-28 | 2021-06-04 | 云南大学 | 单个球的公共自极三角形及正交消失点标定抛物摄像机 |
CN109544646A (zh) * | 2018-11-28 | 2019-03-29 | 云南大学 | 利用球及三正交消失点标定抛物折反射摄像机的方法 |
CN109272558B (zh) * | 2018-11-28 | 2021-07-09 | 云南大学 | 分离圆的公共自极三角形及圆环点标定针孔摄像机的方法 |
CN109712195A (zh) * | 2018-12-17 | 2019-05-03 | 云南大学 | 利用球像的公共自极三角形进行单应性估计的方法 |
CN110033493B (zh) * | 2019-04-17 | 2021-05-11 | 深圳金三立视频科技股份有限公司 | 摄像机3d标定方法及终端 |
CN109934879B (zh) * | 2019-04-22 | 2022-05-20 | 云南大学 | 利用球和公共自极三角形标定抛物折反射摄像机的方法 |
CN110148183B (zh) * | 2019-05-08 | 2022-09-06 | 云南大学 | 利用球及极点极线标定摄像机的方法、存储介质和系统 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1998005157A2 (en) * | 1996-07-12 | 1998-02-05 | Real-Time Geometry Corporation | High accuracy calibration for 3d scanning and measuring systems |
CN101706960A (zh) * | 2009-10-30 | 2010-05-12 | 北京航空航天大学 | 一种共面圆圆心投影点的定位方法 |
CN102855620A (zh) * | 2012-07-13 | 2013-01-02 | 南开大学 | 基于球形投影模型的纯旋转摄像机自标定方法 |
CN102930548A (zh) * | 2012-11-19 | 2013-02-13 | 云南大学 | 利用两个相同的相交椭圆线性求解摄像机内参数 |
JP2013238567A (ja) * | 2012-05-17 | 2013-11-28 | Tetsuya Tagawa | レーダ試験装置 |
CN104835144A (zh) * | 2015-04-09 | 2015-08-12 | 云南大学 | 利用一个球的球心的像及正交性求解摄像机内参数 |
CN105787920A (zh) * | 2014-12-26 | 2016-07-20 | 秦永进 | 球幕标定方法、标定系统及控制设备 |
CN106340044A (zh) * | 2015-07-09 | 2017-01-18 | 上海振华重工电气有限公司 | 摄像机外参自动标定方法及标定装置 |
CN106558082A (zh) * | 2016-11-28 | 2017-04-05 | 云南大学 | 利用Veronese映射棋盘格的投影矩阵标定中心折反射摄像机 |
CN106651956A (zh) * | 2016-09-22 | 2017-05-10 | 云南大学 | 利用单个球及平行圆性质标定拋物折反射摄像机的方法 |
-
2017
- 2017-11-06 CN CN201711075273.2A patent/CN107886546B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1998005157A2 (en) * | 1996-07-12 | 1998-02-05 | Real-Time Geometry Corporation | High accuracy calibration for 3d scanning and measuring systems |
CN101706960A (zh) * | 2009-10-30 | 2010-05-12 | 北京航空航天大学 | 一种共面圆圆心投影点的定位方法 |
JP2013238567A (ja) * | 2012-05-17 | 2013-11-28 | Tetsuya Tagawa | レーダ試験装置 |
CN102855620A (zh) * | 2012-07-13 | 2013-01-02 | 南开大学 | 基于球形投影模型的纯旋转摄像机自标定方法 |
CN102930548A (zh) * | 2012-11-19 | 2013-02-13 | 云南大学 | 利用两个相同的相交椭圆线性求解摄像机内参数 |
CN105787920A (zh) * | 2014-12-26 | 2016-07-20 | 秦永进 | 球幕标定方法、标定系统及控制设备 |
CN104835144A (zh) * | 2015-04-09 | 2015-08-12 | 云南大学 | 利用一个球的球心的像及正交性求解摄像机内参数 |
CN106340044A (zh) * | 2015-07-09 | 2017-01-18 | 上海振华重工电气有限公司 | 摄像机外参自动标定方法及标定装置 |
CN106651956A (zh) * | 2016-09-22 | 2017-05-10 | 云南大学 | 利用单个球及平行圆性质标定拋物折反射摄像机的方法 |
CN106558082A (zh) * | 2016-11-28 | 2017-04-05 | 云南大学 | 利用Veronese映射棋盘格的投影矩阵标定中心折反射摄像机 |
Non-Patent Citations (3)
Title |
---|
Camera Calibration Based on the Common Self-polar Triangle of Sphere Images;Haifeng Huang;《Computer Vision》;20150416;第19-29页 * |
Homography Estimation from the Common Self-polar Triangle of Separate Ellipses;Hanfei Huang;《Computer Vision&Pattern Recognition》;20161212;第1737-1744页 * |
Paracatadioptric camera calibration using sphere images;Duan Huixian 等;《International Conference on Image Processing》;20110930;第1-5页 * |
Also Published As
Publication number | Publication date |
---|---|
CN107886546A (zh) | 2018-04-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107886546B (zh) | 利用球像及公共自极三角形标定抛物折反射摄像机的方法 | |
CN104835144A (zh) | 利用一个球的球心的像及正交性求解摄像机内参数 | |
CN103106661B (zh) | 空间二条相交直线线性求解抛物折反射摄像机内参数 | |
JP6483832B2 (ja) | Rgb−dセンサを使用して物体を走査する方法とシステム | |
Liu et al. | Near-light photometric stereo using circularly placed point light sources | |
CN111080715B (zh) | 三个球和无穷远点的极线性质标定摄像机内参数的方法 | |
CN105279758A (zh) | 使用双球的相切像与圆环点的像标定拋物折反射摄像机 | |
CN109360248A (zh) | 利用单个球及共轭直径的性质标定拋物折反射摄像机 | |
CN109325983A (zh) | 利用无穷远点关于圆极线性质标定拋物折反射摄像机 | |
CN104200476A (zh) | 利用双平面镜装置中的圆周运动求解摄像机内参数的方法 | |
CN105321181A (zh) | 使用双球的相离像与圆环点的像标定拋物折反射摄像机 | |
Cao et al. | Camera calibration using symmetric objects | |
Tao et al. | A convenient and high-accuracy multicamera calibration method based on imperfect spherical objects | |
CN108921904B (zh) | 利用单个球及渐近线的性质标定针孔摄像机的方法 | |
CN105303570A (zh) | 使用双球的相交像及圆环点的像标定拋物折反射摄像机 | |
CN110148183B (zh) | 利用球及极点极线标定摄像机的方法、存储介质和系统 | |
CN109360247B (zh) | 单个球的公共自极三角形及正交消失点标定抛物摄像机 | |
CN107958468B (zh) | 利用空间位置不同的三个球标定中心折反射摄像机的方法 | |
CN113345033B (zh) | 一种中心折反射摄像机内参数标定方法及系统 | |
CN109523598A (zh) | 单个球的公共自极三角形及圆环点标定抛物摄像机的方法 | |
CN109712195A (zh) | 利用球像的公共自极三角形进行单应性估计的方法 | |
CN111429522B (zh) | 利用共面圆的公共极点极线性质标定摄像机的方法和系统 | |
CN103810697A (zh) | 空间四条不平行直线的像标定抛物折反射摄像机内参数 | |
CN109191528B (zh) | 利用球像与圆环点极线的性质标定针孔摄像机的方法 | |
CN105354839A (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 |