CN101661617B - 摄像机标定的方法及装置 - Google Patents

摄像机标定的方法及装置 Download PDF

Info

Publication number
CN101661617B
CN101661617B CN2008101309229A CN200810130922A CN101661617B CN 101661617 B CN101661617 B CN 101661617B CN 2008101309229 A CN2008101309229 A CN 2008101309229A CN 200810130922 A CN200810130922 A CN 200810130922A CN 101661617 B CN101661617 B CN 101661617B
Authority
CN
China
Prior art keywords
matching characteristic
calibrated
image
point
initial alignment
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.)
Expired - Fee Related
Application number
CN2008101309229A
Other languages
English (en)
Other versions
CN101661617A (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.)
Global Innovation Polymerization LLC
Original Assignee
Huawei Device Co Ltd
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 Huawei Device Co Ltd filed Critical Huawei Device Co Ltd
Priority to CN2008101309229A priority Critical patent/CN101661617B/zh
Publication of CN101661617A publication Critical patent/CN101661617A/zh
Application granted granted Critical
Publication of CN101661617B publication Critical patent/CN101661617B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Image Analysis (AREA)
  • Studio Devices (AREA)

Abstract

本发明实施例提供一种摄像机的标定方法及装置,为实现摄像机的在线标定,并提高摄像机标定的精确度而发明。其中,所述方法包括:利用包含有标定参照物的初始图像,对摄像机进行初始标定,获取摄像机的初始标定参数;获取所述已进行初始标定的初始图像与待标定图像的匹配特征参数,其中所述待标定图像是在摄像机参数发生变化时获取的,且所述待标定图像不包括所述标定参照物;利用所述摄像机的初始标定参数及匹配特征参数,对所述待标定图像进行标定,获取摄像机新的标定参数。本发明实施例主要应用于计算机视觉相关技术领域。

Description

摄像机标定的方法及装置
技术领域
本发明涉及计算机视觉检测技术领域,尤其涉及一种摄像机标定的方法及装置。
背景技术
在计算机视觉检测过程中,为了获取计算机图像象素点与实际物理空间点的对应关系,摄像机标定是必不可少的过程。所述摄像机标定,就是在一定的摄像机模型下,经过对图像的处理,利用一系列数学变换和计算方法,获得摄像机内部的几何和光学特性(也即内部参数),以及摄像机坐标系相对于空间坐标系的位置关系(也即外部参数)。
在现有技术中,有多种摄像机的标定方法,如传统标定方法,自标定方法等。下面分别介绍一下传统标定法和自标定法的原理。
(一)传统标定法
利用传统标定法进行摄像机标定的时候,需要使用标定参照物。传统标定法可包括直接线性变换标定法(DLT),径向排列约束(RAC)标定法,主动视觉标定法和平面标定法等。下面以直接线性变换标定法和径向排列约束标定法为例,描述一下传统标定法的基本原理。
直接线性变换标定法是由Abdal-Aziz和Karara于70年代初提出的。该方法首先需建立摄像机成像模型线性方程组,并测出场景中一组点的世界坐标和其在成像平面上的对应坐标,然后将这些坐标值代入该线性方程组中求出该线性方程组的未知系数,该线性方程组的系数包含了摄像机的内部参数。
但是,直接线性变换标定法忽略了镜头畸变,只适合视野较小的摄像机标定。当镜头畸变明显时(尤其是采用广角镜头),所述摄像机成像模型线性模型无法准确地描述摄像机的成像过程,这时在摄像机的标定过程中需要引入畸变因子进行校正。
Tsai于80年代中期提出了基于RAC(径向排列约束)的标定方法。该方法的核心是利用径向一致约束,来求解除摄像机光轴方向的平移外的其它摄像机外部参数,然后再求解摄像机的其它参数。由于RAC方法考虑了径向畸变,因此相比DLT方法,它的精度更高。
(二)自标定法
1992年,O.D.Faugeras,Q.T.Luong,和S.J.Maybank首先提出了自标定的概念,使得在场景未知和摄像机任意运动的一般情形下的摄像机标定成为可能。也就是说,利用自标定方法进行摄像机标定时,不需要标定参照物,仅仅是通过图像点之间的对应关系就可以对摄像机进行标定。自标定依据的是多幅图像中成像点间存在的特殊约束关系(如极线约束),因此可以不需要场景的结构信息。常用的自标定方法包括:基于Kruppa方程的自标定方法和基于绝对二次曲面、无穷远平面的自标定方法等。
但是,在实现本发明的过程中,发明人发现现有技术中至少存在如下的问题:
由于传统标定法在拍摄和标定过程中一直需要使用标定参照物,因而给拍摄操作及标定方法的使用带来了很大的不便。特别是在一些计算机视觉和摄影测量任务的执行过程中,摄像机的焦距、位置等参数由于受视觉或测量系统的机械震动或热效应影响而变化,或者由于执行某一具体的任务而需要有意改变,从而使得初始的摄像机标定参数不再有效,因而需要进行在线重标定。而由于标定参照物的使用和位置调整将导致在线任务的中断,所以传统标定方法不适合在线重标定。而自标定方法的标定精度不高,鲁棒性不足;并且利用自标定方法进行摄像机标定时一般无法获取摄像机的外部参数。
发明内容
本发明的实施例提供一种摄像机的标定方法及装置,以实现在线摄像机标定,并提高摄像机标定的精确度。
本发明实施例摄像机的标定方法采用以下技术方案:
一种摄像机的标定方法包括如下步骤:
利用包含有标定参照物的初始图像,对摄像机进行初始标定,获取摄像机的初始标定参数;
获取所述已进行初始标定的初始图像与待标定图像的匹配特征参数,其中所述待标定图像是在摄像机参数发生变化时获取的,且所述待标定图像不包括所述标定参照物;
利用所述摄像机的初始标定参数及所述匹配特征参数,对所述待标定图像进行标定,获取摄像机新的标定参数。
本发明实施例摄像机的标定装置采用以下技术方案:
一种摄像机的标定装置包括:
初始标定参数获取单元,用于利用包含有标定参照物的初始图像,对摄像机进行初始标定,获取摄像机的初始标定参数;
匹配特征参数获取单元,用于获取所述已进行初始标定的初始图像与待标定图像的匹配特征参数,其中所述待标定图像是在摄像机参数发生变化时获取的,且所述待标定图像不包括所述标定参照物;
重标定参数获取单元,用于利用由所述初始标定参数获取单元获取的摄像机的初始标定参数及由所述匹配特征参数获取单元获取的匹配特征参数,对所述待标定图像进行标定,获取摄像机新的标定参数。
本发明实施例摄像机的标定方法及装置,可先对包含有标定参照物初始图像进行初始标定,并获取摄像机的初始标定参数。当摄像机的参数发生变化时,可获得初始图像与待标定图像之间的匹配特征参数。继而通过已经获得的摄像机的初始标定参数以及匹配特征参数,计算出待标定图像的摄像机参数。与现有技术中的传统标定法相比,由于本发明实施例计算出了包含标定参照物的初始图像以及不包含标定参照物的待标定图像之间的匹配特征参数,因此,当摄像机的参数发生变化时,能够准确的获得摄像机的新的标定参数,实现对摄像机的在线标定。而相较于现有技术中的自标定法而言,在进行摄像机的初始标定时,利用的是包含标定参照物的初始图像先进行摄像机的初始标定,并以获得的摄像机的初始标定参数作为进行对摄像机进行新的标定的基础,因此,本发明实施例提高了摄像机标定的精确度。
附图说明
图1为本发明实施例一摄像机的标定方法的流程图;
图2为本发明实施例二摄像机的标定方法的流程图;
图3为本发明实施例三摄像机的标定装置的示意图;
图4为本发明实施例三摄像机的标定装置的结构图;
图5为本发明实施例三中匹配特征参数获取单元的示意图。
具体实施方式
为了更清楚地说明本发明实施例的技术方案,下面将对实施例描述中所需要使用的附图作一简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
为解决现有技术中的传统标定法以及自标定法无法实现摄像机的在线标定以及摄像机标定精确度不高的问题,本发明实施例提供了一种摄像机标定的方法。下面结合附图对本发明实施例做详细说明。
实施例一
如图1所示,本发明实施例一摄像机的标定方法包括如下步骤:
步骤11、利用包含有标定参照物的初始图像,对摄像机进行初始标定,获取摄像机的初始标定参数。
在此步骤中,为了保证摄像机标定的精确度,可利用传统的两步法,即径向排列约束标定法进行摄像机标定。当然,还可利用其他的传统标定法,如DLT方法等对所述初始图像进行摄像机标定。所述的标定参照物可任意选定。
步骤12、获取所述已进行初始标定的初始图像与待标定图像的匹配特征参数;其中所述待标定图像是在摄像机参数发生变化时获取的,且所述待标定图像不包括所述标定参照物。
所述的匹配特征参数包括:获取的匹配特征点,以及所述匹配特征点的像素坐标以及世界坐标等。
在本发明实施例中,可以利用尺度不变特征变换(SIFT,Scale InvariantFeature Transform),获取已进行初始标定的初始图像与待标定图像的匹配特征参数。下面对尺度不变特征变换算法作一简单介绍。尺度不变特征变换算法首先在尺度空间进行特征检测,并确定关键点的位置和关键点所处的尺度,然后使用关键点邻域梯度的主方向作为该关键点的方向特征,以实现算子对尺度和方向的无关性。尺度不变特征变换算法的主要计算步骤如下:
1)尺度空间极值检测:搜索整个尺度和图像位置。通过使用高斯差分函数确定对尺度和方向具有不变性的兴趣点,以初步确定关键点位置和所在尺度。
2)关键点定位:在每一个候选位置,详细地模型拟合以确定所述关键点的位置和尺度,并根据各关键点的稳定性的高低,选择关键点。然后通过拟合三维二次函数以精确地确定关键点的位置和尺度,同时消除低对比度的关键点和不稳定的边缘响应点,以增强匹配稳定性、提高抗噪声能力。
3)方向分配:根据局部图像梯度,为每一个关键点方向分配一个或多个方向。其后对该关键点的所有操作都根据为其分配的方向、尺度、位置进行。利用关键点邻域像素的梯度方向分布特性,为每个关键点指定方向参数,使算子具备旋转不变性。
4)关键点描述符:在每个关键点周围,在选定的尺度下测量局部图像梯度。为了增强匹配的稳健性,对每个关键点使用4×4共16个种子点来描述,这样对于一个关键点就可以产生128个数据,最终形成128维的SIFT特征向量。
利用尺度不变特征变换算法提取的SIFT特征向量具有如下特性:
a)SIFT特征是图像的局部特征,其对旋转、尺度缩放、亮度变化保持不变性,对视角变化、仿射变换、噪声也保持一定程度的稳定;
b)独特性(Distinctiveness)好,信息量丰富,适用于在海量特征数据库中进行快速、准确的匹配等。
当然,除了利用上述尺度不变特征变换方法获取所述匹配特征参数外,还可利用其他的方法获得,利用Harris角点提取法等。
步骤13、利用所述摄像机的初始标定参数及匹配特征参数,对所述待标定图像进行标定,获取摄像机新的标定参数。
通过上述步骤可以看出,与现有技术中的传统标定法相比,由于本发明实施例计算出了包含标定参照物的初始图像以及不包含标定参照物的待标定图像之间的匹配特征参数,因此,当摄像机的参数发生变化时,能够准确的获得摄像机的新的标定参数,实现对摄像机的在线标定。而相较于现有技术中的自标定法而言,在进行摄像机的初始标定时,利用的是包含标定参照物的初始图像先进行摄像机的初始标定,并以获得的摄像机的初始标定参数作为进行对摄像机进行新的标定的基础,因此,本发明实施例提高了摄像机标定的精确度。
实施例二
下面结合实施例二详细描述一下,本发明实施例所述的摄像机的标定方法的实现过程。在实施例二中,以选取两组包含有标定参照物初始图像为例进行描述。如图2所示,实施例二所述的方法包括如下步骤:
步骤21、获取进行初始标定的初始图像。
在标定开始时,利用摄像机拍摄两组包含有标定参照物的初始图像,这两组图像分别用来提供非共面的特征点。在此,为方便描述,分别将该两组初始图像称为第一组图像和第二组图像。在拍摄第二组图像时,可使摄像机的内部参数和外部参数相对第一组图像稍微有所变化,以保证两组图像中有一定的相同的内容,并可保证摄像机标定的客观性和准确性。
步骤22、在获得了初始图像后,利用角点提取法提取所述初始图像的特征点。所述的角点提取法包括Harris提取法,SIFT提取法等。但由于Harris角点提取法操作简便,在此利用其进行特征点提取时能够保证摄像机标定的效率。其中Harris角点提取法的原理与现有技术中的相同,在此不再赘述。
步骤23、在提取出特征点后,将这些角点分别在两组图像中的像素坐标(xf,yf)记录下来,并获取这些角点对应的世界坐标(xw,yw,zw)。
步骤24、根据传统的两步法(RAC标定法),获取每组图像所对应的摄像机参数,包括摄像机的内部参数和外部参数。其中,传统的两步法的计算原理与现有技术中的相同,在此做简单描述。
计算摄像机参数的具体过程如下:首先根据所述特征点的图像坐标以及所述世界坐标,利用最小二乘法求解超定线性方程,计算摄像机的外部参数,然后根据所述特征点的图像坐标以及所述世界坐标,计算摄像机的内部参数,分别用以下参数表示:
(f1
Figure G2008101309229D00081
Figure G2008101309229D00082
Figure G2008101309229D00084
Figure G2008101309229D00085
Figure G2008101309229D00086
Figure G2008101309229D00089
Figure G2008101309229D000810
Figure G2008101309229D000812
Figure G2008101309229D000813
(f2
Figure G2008101309229D000816
Figure G2008101309229D000817
Figure G2008101309229D000818
Figure G2008101309229D000819
Figure G2008101309229D000820
Figure G2008101309229D000822
Figure G2008101309229D000823
Figure G2008101309229D000824
Figure G2008101309229D000825
Figure G2008101309229D000827
Figure G2008101309229D000828
Figure G2008101309229D000829
Figure G2008101309229D000830
)。其中上标1、2分别表示第一图像组和第二组图像中的相应参数。其中f是内参焦距(mm),(Cx,Cy)分别是光心的象素坐标(pix),k1是摄像机镜头径向畸变的一阶系数,Sx是非确定性标度尺寸因子。外部参数r、t分别是三维空间世界坐标系与摄像机坐标系之间的旋转矩阵和平移向量。其中,Tx,Ty,Tz是从世界坐标系到摄像机坐标系变换的沿三个坐标轴的平移量(旋转之后)。其中r,t可以表示成如公示(1)的形式:
r = r 11 r 12 r 13 r 21 r 22 r 23 r 31 r 32 r 33 , t = t x t y t z     (1)
如果摄像机坐标系在三维空间世界坐标系下的方向绕X轴逆时针旋转角度α,绕Y轴逆时针旋转角度β,绕Z轴逆时针旋转角度γ,则旋转矩阵为r=rαrβrγ,其中,rα,rβ,rγ可表示成如公示(2),(3),(4)的形式;
r α = 1 0 0 0 cos α sin α 0 - sin α cos α , - - - ( 2 )
r β = cos β 0 - sin β 0 1 0 sin β 0 cos β , - - - ( 3 )
r γ = cos γ sin γ 0 - sin γ cos γ 0 0 0 1 . - - - ( 4 )
步骤25、在对两组初始图像进行完摄像机的初始标定后,移开标定参照物,相对于上述的第二组图像,获取待标定的图像。
步骤26、利用尺度不变特征变换提取法提取所述已进行初始标定的初始图像和待标定图像的特征点,以便保证提取的特征点对图像旋转、尺度缩放、平移甚至亮度变化保持不变性,对视角变化、仿射变换、噪声也保持一定程度的稳定性,且不需要假设场景中存在角点、直线、圆等固定的几何特征,从而对拍摄场景没有限制或要求。其中所述特征点的提取法与实施例一中所述的相同。当然,在此步骤中也可利用Harris角点提取法等进行特征点的提取。
步骤27、将所述已进行初始标定的初始图像与待标定图像进行特征点匹配,获取匹配特征点。
匹配特征点的提取是基于SIFT特征的提取方法,具体过程是:在第一组图像或第二组图像分别选取一幅图像,从待标定的图像中选取一幅图像,在所述选取的两幅图像之间进行SIFT特征点匹配,采用两幅图像中特征点的特征向量的欧式距离作为两幅图像中SIFT特征点的相似性判定度量。
首先,利用尺度不变特征变换提取法,为所述已进行初始标定的初始图像中的各个特征点,在所述待标定图像的各特征点中,分别选取与所述已进行初始标定的初始图像中的相应特征点间的欧式距离最近的特征点和次近的特征点。若所述最近的特征点所对应的欧式距离与所述次近的特征点所对应的欧式距离之间的比例阈值小于预设值,则所述最近的特征点为所述已进行初始标定的初始图像中相应特征点的匹配特征点。
例如若假设已进行初始标定的初始图像中的特征点分别为P1,P2,P3......Pn,在待标定的图像中的特征点分别为Q1,Q2,Q4......Qm。以P1点为例,首先在待标定图像的特征点中选取与P1点的欧式距离最近的特征点和次近的特征点,假设为分别Q2点和Q3点,那么若P1点与Q2点之间的欧式距离与P1点与Q3点之间的欧式距离之间的比例阈值小于预设值,则可确定Q2点为P1点的匹配特征点。对所述已进行初始标定的初始图像中的其他特征点,确定其匹配特征点的过程相同。
该预设值可以设为0.6,当然,还可根据需要设置为其他的数值。预设值设置的越小,则两幅图像中匹配特征点的个数就越少,但是摄像机标定的稳定性就越高。
步骤28、分别获取所述匹配特征点在所述已进行初始标定的包含有标定参照物的初始图像中的像素坐标以及在待标定图像中的像素坐标。
在具体的查找匹配特征点的过程中,为了保证特征点匹配的正确性,可从第一组图像、第二组图像以及待标定图像中各选取一幅图像,可采用“两两匹配”的方式,以确定在三幅图像中对应同一匹配特征点的像素坐标。
首先利用由所述第一组图像中选取的图像,由第二组图像中选取的图像,选取这两幅图像中的对应某三维空间特征点A的匹配特征向量V1、V2。然后,将第一轮匹配的特征向量V1与由待标定图像中所选取的图像进行匹配。
若能在由待标定图像中所选取的图像中找到匹配的特征向量V3,则再以V3为基础,在由第二组图像中选取的图像中寻找匹配的特征向量;若找到的特征向量正好是V2,则表明该特征点A是三幅图像共有的匹配特征点,特征向量分别为V1、V2、V3。然后将从SIFT特征向量V1、V2、V3中抽取出的特征点的象素坐标记录下来,分别为
Figure G2008101309229D00101
Figure G2008101309229D00102
Figure G2008101309229D00103
这种三幅图像之间循环的两两匹配方法能最大可能保证特征点匹配准确性。
根据所述匹配特征点的像素坐标以及摄像机的初始标定参数,按照公式(5)-(10)所示的计算所述匹配特征点在理想坐标系下的理想图像坐标(xu,yu)。
x d = d x ′ ( x f - C x ) / s x - - - ( 5 )
yd=dy(yf-Cy)     (6)
xu=xd(1+k1r2)    (7)
yu=yd(1+k1r2)    (8)
d x ′ = d x N cx N fx - - - ( 9 )
r 2 = x d 2 + y d 2 - - - ( 10 )
其中(xd,yd)是匹配特征点的像素坐标,dx,dy分别是x方向(扫描线方向)和y方向相邻CCD感光单元中心之间的距离(mm);Ncx是X方向感光单元的数目(由摄像机厂家提供),Nfx是计算机每行采样的象素数目,即图像的X方向尺寸(pix个数)。
步骤29、计算所述匹配特征点的世界坐标(xw,yw,zw)。
在此步骤中,首先利用所述摄像机初始标定的摄像机参数,以及所述匹配特征点的像素坐标,形成超定线性方程组(11);
( f 1 r 11 1 - x u 1 r 31 1 ) x w + ( f 1 r 12 1 - x u 1 r 32 1 ) y w + ( f 1 r 13 1 - x u 1 r 33 1 ) z w = x u 1 T z 1 - f 1 T x 1 ( f 1 r 21 1 - y u 1 r 31 1 ) x w + ( f 1 r 22 1 - y u 1 r 32 1 ) y w + ( f 1 r 23 1 - y u 1 r 33 1 ) z w = y u 1 T z 1 - f 1 T y 1 ( f 2 r 11 2 - x u 2 r 31 2 ) x w + ( f 2 r 12 2 - x u 2 r 32 2 ) y w + ( f 2 r 13 2 - x u 2 r 33 2 ) z w = x u 2 T z 2 - f 2 T x 2 ( f 2 r 21 2 - y u 2 r 31 2 ) x w + ( f 2 r 22 2 - y u 2 r 32 2 ) y w + ( f 2 r 23 2 - y u 2 r 33 2 ) z w = y u 2 T z 2 - f 2 T y 2 - - - ( 11 )
通过最小二乘法求解所述超定线性方程组,计算出所述匹配特征点的世界坐标(xw,yw,zw)。
步骤30、利用上述步骤21-29求出匹配特征点的像素坐标,以及所述匹配特征点的世界坐标后,即可再次利用传统的两步法求解出待标定图像的摄像机参数。其中传统两步法的计算原理与步骤24中描述的相同。
此后若摄像机参数再次发生改变,需要对其他待标定图像进行标定时,即可执行上述的步骤25-步骤30即可。
利用实施利一和实施例二的摄像机的标定方法,可先对包含有标定参照物初始图像进行初始标定,并获取摄像机的初始标定参数。当摄像机的参数发生变化时,可获得初始图像与待标定图像之间的匹配特征参数。继而通过已经获得的摄像机的初始标定参数以及匹配特征参数,计算出待标定图像的摄像机参数。与现有技术中的传统标定法相比,由于本发明实施例计算出了包含标定参照物的初始图像以及不包含标定参照物的待标定图像之间的匹配特征参数,因此,当摄像机的参数发生变化时,能够准确的获得摄像机的新的标定参数,实现对摄像机的在线标定。而相较于现有技术中的自标定法而言,在进行摄像机的初始标定时,利用的是包含标定参照物的初始图像先进行摄像机的初始标定,并以获得的摄像机的初始标定参数作为进行对摄像机进行新的标定的基础,因此,本发明实施例提高了摄像机标定的精确度。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程,是可以通过计算机程序来指令相关的硬件来完成,所述的程序可存储于一计算机可读取存储介质中,该程序在执行时,可包括如上述各方法的实施例的流程。其中,所述的存储介质可为磁碟、光盘、只读存储记忆体(Read-Only Memory,ROM)或随机存储记忆体(Random Access Memory,RAM)等。
实施例三
与本发明实施例摄像机的标定方法相对应,本发明实施例三还提供了一种摄像机的标定装置。
如图3所示,本发明实施例三所述的摄像机的标定装置,包括:初始标定参数获取单元42,匹配特征参数获取单元43,以及重标定参数获取单元44。
其中,初始标定参数获取单元42,用于利用获取的包含有标定参照物的初始图像,对摄像机进行初始标定,获取摄像机的初始标定参数;匹配特征参数获取单元43,用于获取所述已进行初始标定的初始图像与待标定图像的匹配特征参数;重标定参数获取单元44,用于利用由所述初始标定参数获取单元42获取的摄像机的初始标定参数及由所述匹配特征参数获取单元43获取的匹配特征参数,对所述待标定图像进行标定,获取摄像机新的标定参数。
为进一步提高摄像机标定的效率,如图4所示,所述装置还包括:图像获取单元41,用于获取包含有标定参照物的初始图像,以及用于当摄像机参数发生变化时,获取不包含所述标定参照物的待标定图像。
其中,如图5所示,所述匹配特征参数获取单元43包括:特征点提取模块431,用于提取所述已进行初始标定的初始图像和待标定图像的特征点;匹配特征点提取模块432,用于将由所述特征点提取模块431提取的所述已进行初始标定的初始图像的特征点与待标定图像中的特征点进行匹配,获取匹配特征点;坐标获取模块433,用于根据所述匹配特征点在已进行初始标定的初始图像中的像素坐标、所述匹配特征点在所述待标定图像中的像素坐标以及所述摄像机的初始标定参数,计算所述匹配特征点的世界坐标。所述特征点提取模块431在进行特征点提取的过程中,可以采用Harris角点提取法,SIFT提取法等。
为了提高摄像机标定的准确性,所述匹配特征点提取模块432可包括:特征点选取子模块4321,用于利用尺度不变特征变换提取法,为所述已进行初始标定的初始图像中的各个特征点,在所述待标定图像的各特征点中,分别选取与所述已进行初始标定的初始图像中的相应特征点间的欧式距离最近的特征点和次近的特征点;匹配特征点确定子模块4322,用于在所述最近的特征点所对应的欧式距离与所述次近的特征点所对应的欧式距离之间的比例阈值小于预设值时,确定所述最近的特征点为所述已进行初始标定的初始图像中相应特征点的匹配特征点。其中所述的预设值可以为0.6,或其他任意的数值。当该预设值的值越小时,匹配特征点的个数越小,但是摄像机标定越稳定。
其中所述坐标获取模块433包括:图像坐标获取子模块4331,用于分别获取由所述匹配特征点提取模块431获取的匹配特征点,在所述已进行初始标定的初始图像中的像素坐标及待标定图像中的像素坐标;世界坐标获取子模块4332,用于利用所述摄像机的初始标定参数,以及所述匹配特征点在已进行初始标定的初始图像中的像素坐标,形成超定线性方程组,并根据所述超定线性方程组,计算出所述匹配特征点的世界坐标。
此时,所述重标定参数获取单元44,用于根据所述世界坐标获取子模块4332计算得到的匹配特征点的世界坐标,以及由所述图像坐标获取子模块4331获得的所述匹配特征点在所述待标定图像中的像素坐标,计算摄像机新的标定参数。
本发明实施例三所述的设摄像机的标定装置的工作原理,与实施例一和实施例二所描述的实现过程相同,在此不再赘述。
综上所述,本发明实施例摄像机的标定方法及装置,可先对包含有标定参照物初始图像进行初始标定,并获取摄像机的初始标定参数。当摄像机的参数发生变化时,可获得初始图像与待标定图像之间的匹配特征参数。继而通过已经获得的摄像机的初始标定参数以及匹配特征参数,计算出待标定图像的摄像机参数。与现有技术中的传统标定法相比,由于本发明实施例计算出了包含标定参照物的初始图像以及不包含标定参照物的待标定图像之间的匹配特征参数,因此,当摄像机的参数发生变化时,能够准确的获得摄像机的新的标定参数,实现对摄像机的在线标定。而相较于现有技术中的自标定法而言,在进行摄像机的初始标定时,利用的是包含标定参照物的初始图像先进行摄像机的初始标定,并以获得的摄像机的初始标定参数作为进行对摄像机进行新的标定的基础,因此,本发明实施例提高了摄像机标定的精确度。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应所述以权利要求的保护范围为准。

Claims (9)

1.一种摄像机的标定方法,其特征在于,所述方法包括如下步骤:
利用包含有标定参照物的初始图像,对摄像机进行初始标定,获取摄像机的初始标定参数;
获取所述已进行初始标定的初始图像与待标定图像的匹配特征参数,其中所述待标定图像是在摄像机参数发生变化时获取的,且所述待标定图像不包括所述标定参照物,所述匹配特征参数包括:匹配特征点,以及所述匹配特征点的像素坐标以及世界坐标;所述获取所述已进行初始标定的初始图像与待标定图像的匹配特征参数包括:
提取所述已进行初始标定的初始图像和待标定图像的特征点;将所述已进行初始标定的初始图像与待标定图像的特征点进行匹配,获取匹配特征点;分别获取所述匹配特征点在所述已进行初始标定的初始图像中的像素坐标及待标定图像中的像素坐标;根据所述匹配特征点在已进行初始标定的初始图像中的像素坐标、所述匹配特征点在所述待标定图像中的像素坐标以及所述摄像机的初始标定参数,计算所述匹配特征点的世界坐标;
利用所述摄像机的初始标定参数及所述匹配特征参数,对所述待标定图像进行标定,获取摄像机新的标定参数。
2.根据权利要求1所述的摄像机的标定方法,其特征在于,当利用尺度不变特征变换提取法时,将所述已进行初始标定的初始图像与待标定图像的特征点进行匹配,获取匹配特征点的步骤包括:
为所述已进行初始标定的初始图像中的各个特征点,在所述待标定图像的各特征点中,分别选取与所述已进行初始标定的初始图像中的相应特征点间的欧式距离最近的特征点和次近的特征点;
若所述最近的特征点所对应的欧式距离与所述次近的特征点所对应的欧式距离之间的比例阈值小于预设值,则所述最近的特征点为所述已进行初始标定的初始图像中相应特征点的匹配特征点。
3.根据权利要求1所述的摄像机的标定方法,其特征在于,根据所述匹配特征点在已进行初始标定的初始图像中的像素坐标、所述匹配特征点在所述待标定图像中的像素坐标以及所述摄像机的初始标定参数,计算所述匹配特征点的世界坐标的步骤具体为:
利用所述摄像机的初始标定参数,以及所述匹配特征点在已进行初始标定的初始图像中的像素坐标,形成超定线性方程组;
根据所述超定线性方程组,计算出所述匹配特征点的世界坐标。
4.根据权利要求3所述的摄像机的标定方法,其特征在于,利用所述初始标定参数及匹配特征参数,对所述待标定图像进行标定,获取摄像机新的标定参数的步骤具体为:
根据所述匹配特征点的世界坐标,以及所述匹配特征点在所述待标定图像中的像素坐标,计算摄像机新的标定参数。
5.一种摄像机的标定装置,其特征在于,所述装置包括:
初始标定参数获取单元,用于利用包含有标定参照物的初始图像,对摄像机进行初始标定,获取摄像机的初始标定参数;
匹配特征参数获取单元,用于获取所述已进行初始标定的初始图像与待标定图像的匹配特征参数,其中所述待标定图像是在摄像机参数发生变化时获取的,且所述待标定图像不包括所述标定参照物,所述匹配特征参数包括:匹配特征点,以及所述匹配特征点的像素坐标以及世界坐标;
重标定参数获取单元,用于利用由所述初始标定参数获取单元获取的摄像机的初始标定参数及由所述匹配特征参数获取单元获取的匹配特征参数,对所述待标定图像进行标定,获取摄像机新的标定参数;
所述匹配特征参数获取单元包括:
特征点提取模块,用于分别提取所述已进行初始标定的初始图像的特征点和所述待标定图像的特征点;
匹配特征点提取模块,用于将由所述特征点提取模块提取的所述已进行初始标定的初始图像的特征点与所述待标定图像中的特征点进行匹配,获取匹配特征点;
坐标获取模块,用于根据所述匹配特征点在已进行初始标定的初始图像中的像素坐标、所述匹配特征点在所述待标定图像中的像素坐标以及所述摄像机的初始标定参数,计算所述匹配特征点的世界坐标。
6.根据权利要求5所述的摄像机的标定装置,其特征在于,所述装置还包括:
图像获取单元,用于获取包含有标定参照物的初始图像,以及用于当摄像机参数发生变化时,获取不包含所述标定参照物的待标定图像。
7.根据权利要求5所述的摄像机的标定装置,其特征在于,所述匹配特征点提取模块包括:
特征点选取子模块,用于利用尺度不变特征变换提取法,为所述已进行初始标定的初始图像中的各个特征点,在所述待标定图像的各特征点中,分别选取与所述已进行初始标定的初始图像中的相应特征点间的欧式距离最近的特征点和次近的特征点;
匹配特征点确定子模块,用于在所述最近的特征点所对应的欧式距离与所述次近的特征点所对应的欧式距离之间的比例阈值小于预设值时,确定所述最近的特征点为所述已进行初始标定的初始图像中相应特征点的匹配特征点。
8.根据权利要求7所述的摄像机的标定装置,其特征在于,所述坐标获取模块包括:
图像坐标获取子模块,用于分别获取由所述匹配特征点提取模块获取的匹配特征点,在所述已进行初始标定的初始图像中的像素坐标及待标定图像中的像素坐标;
世界坐标获取子模块,用于利用所述摄像机的初始标定参数以及所述匹配特征点在已进行初始标定的初始图像中的像素坐标,形成超定线性方程组,并根据所述超定线性方程组,获取所述匹配特征点的世界坐标。
9.根据权利要求8所述的摄像机的标定装置,其特征在于,
所述重标定参数获取单元,用于根据所述匹配特征点的世界坐标,以及所述匹配特征点在所述待标定图像中的像素坐标,计算摄像机新的标定参数。
CN2008101309229A 2008-08-30 2008-08-30 摄像机标定的方法及装置 Expired - Fee Related CN101661617B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2008101309229A CN101661617B (zh) 2008-08-30 2008-08-30 摄像机标定的方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2008101309229A CN101661617B (zh) 2008-08-30 2008-08-30 摄像机标定的方法及装置

Publications (2)

Publication Number Publication Date
CN101661617A CN101661617A (zh) 2010-03-03
CN101661617B true CN101661617B (zh) 2011-11-02

Family

ID=41789621

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2008101309229A Expired - Fee Related CN101661617B (zh) 2008-08-30 2008-08-30 摄像机标定的方法及装置

Country Status (1)

Country Link
CN (1) CN101661617B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106022281A (zh) * 2016-05-27 2016-10-12 广州帕克西软件开发有限公司 一种面部数据测量方法及系统

Families Citing this family (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101957990B (zh) * 2010-08-13 2012-07-04 武汉大学 一种摄像机定标方法、图像处理设备和机动车辆
CN102646277B (zh) * 2012-03-05 2015-06-10 上海海事大学 一种具有实时参数监控和在线标定功能的视觉系统
CN102629381B (zh) * 2012-03-05 2014-12-10 上海海事大学 一种无人水面艇视觉系统的标定参数实时监测方法
CN102693539A (zh) * 2012-03-13 2012-09-26 夏东 一种用于智能监控系统的宽基线快速三维标定方法
CN105157681B (zh) * 2015-08-23 2018-07-24 西北工业大学 室内定位方法、装置和摄像机以及服务器
CN105405135B (zh) * 2015-11-06 2018-08-10 中国人民解放军信息工程大学 基于基本配置点的两步法摄影物点、像点自动匹配方法
CN105513056B (zh) * 2015-11-30 2018-06-19 天津津航技术物理研究所 车载单目红外相机外参自动标定方法
CN106441160A (zh) * 2016-08-31 2017-02-22 上海衣得体信息科技有限公司 基于人体扫描仪的标定方法及人体扫描系统
CN106447734A (zh) * 2016-10-10 2017-02-22 贵州大学 一种采用人脸标定物的智能手机相机标定算法
CN106815870B (zh) * 2016-12-16 2020-05-19 珠海研果科技有限公司 一种机内标定全景摄像装置的方法和系统
CN107566720B (zh) * 2017-08-25 2019-09-27 维沃移动通信有限公司 一种更新移动终端的标定值的方法及移动终端
WO2019084726A1 (zh) * 2017-10-30 2019-05-09 深圳市柔宇科技有限公司 基于标志物的摄像头图像处理方法、增强现实设备
CN108881898A (zh) * 2018-06-07 2018-11-23 歌尔股份有限公司 景深模组非线性标定的测试方法
CN110174093B (zh) * 2019-05-05 2022-10-28 腾讯科技(深圳)有限公司 定位方法、装置、设备和计算机可读存储介质
CN113538588A (zh) * 2020-04-17 2021-10-22 虹软科技股份有限公司 标定方法、标定装置及应用其的电子设备

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1590040A (zh) * 2003-09-03 2005-03-09 中国科学院自动化研究所 一种基于机器人运动的摄像机自标定方法
CN1620153A (zh) * 2004-12-02 2005-05-25 武汉大学 利用平面控制点场进行非量测数码相机标定的方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1590040A (zh) * 2003-09-03 2005-03-09 中国科学院自动化研究所 一种基于机器人运动的摄像机自标定方法
CN1620153A (zh) * 2004-12-02 2005-05-25 武汉大学 利用平面控制点场进行非量测数码相机标定的方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
方晓莹 等.摄像机内参数异构情况下交通参数的获取.《信号处理》.2005,第21卷(第4A期), *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106022281A (zh) * 2016-05-27 2016-10-12 广州帕克西软件开发有限公司 一种面部数据测量方法及系统

Also Published As

Publication number Publication date
CN101661617A (zh) 2010-03-03

Similar Documents

Publication Publication Date Title
CN101661617B (zh) 摄像机标定的方法及装置
US11570423B2 (en) System and methods for calibration of an array camera
CN101630406B (zh) 摄像机的标定方法及摄像机标定装置
Thirthala et al. The radial trifocal tensor: A tool for calibrating the radial distortion of wide-angle cameras
Devernay et al. Automatic calibration and removal of distortion from scenes of structured environments
CN107833181B (zh) 一种基于变焦立体视觉的三维全景图像生成方法
Moulon et al. Adaptive structure from motion with a contrario model estimation
Gonzalez-Aguilera et al. An automatic approach for radial lens distortion correction from a single image
CN105701827A (zh) 可见光相机与红外相机的参数联合标定方法及装置
Sun et al. Review of calibration methods for Scheimpflug camera
Lin et al. Cylindrical panoramic image stitching method based on multi-cameras
Xu et al. Hand-held 3-D reconstruction of large-scale scene with kinect sensors based on surfel and video sequences
Camposeco et al. Non-parametric structure-based calibration of radially symmetric cameras
Ventura et al. P1ac: Revisiting absolute pose from a single affine correspondence
CN107256563B (zh) 基于差异液位图像序列的水下三维重建系统及其方法
CN110310337B (zh) 一种基于光场基本矩阵的多视光场成像系统全参数估计方法
CN110322514B (zh) 一种基于多中心投影模型的光场相机参数估计方法
CN111753644A (zh) 在三维人脸扫描上检测关键点的方法及装置
Omidalizarandi et al. Robust external calibration of terrestrial laser scanner and digital camera for structural monitoring
Qian et al. Image distortion correction for single-lens stereo vision system employing a biprism
Liu et al. Image Mosaic Technology Based on Harris Corner Feature
Reichel et al. Camera calibration as machine learning problem using dense phase shifting pattern, checkerboards, and different cameras
Dhome Visual Perception Through Video Imagery
Liu et al. A fast epipolar line matching method based on 3D spherical panorama
Al Assaad et al. Homography-based model with light calibration for plenoptic cameras

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20180206

Address after: 518129 California

Patentee after: Tanous Co.

Address before: 518129 Longgang District, Guangdong, Bantian HUAWEI base B District, building 2, building No.

Patentee before: HUAWEI DEVICE Co.,Ltd.

Effective date of registration: 20180206

Address after: California, USA

Patentee after: Global innovation polymerization LLC

Address before: California, USA

Patentee before: Tanous Co.

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

Granted publication date: 20111102

Termination date: 20210830