CN102982550B - 利用正五棱台求解摄像机内参数 - Google Patents
利用正五棱台求解摄像机内参数 Download PDFInfo
- Publication number
- CN102982550B CN102982550B CN201210542006.2A CN201210542006A CN102982550B CN 102982550 B CN102982550 B CN 102982550B CN 201210542006 A CN201210542006 A CN 201210542006A CN 102982550 B CN102982550 B CN 102982550B
- Authority
- CN
- China
- Prior art keywords
- point
- little
- positive
- terrace
- hidden
- 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
Links
Landscapes
- Image Analysis (AREA)
Abstract
本发明是由底边是正五边形的正棱台构成的用于摄像机自标定的靶标。具体的步骤包括:从图像上提取正五棱台的上底面和两个侧面上的特征点,根据两平行直线的交点是无穷远点,而无穷远点的像称为隐消点,以及共线四点交比的性质,得出图像平面上隐消点坐标,根据正交隐消点对绝对二次曲线的约束,实现对摄像机内参数的线性求解。对靶标从不同方向拍摄两幅图像,提取每幅图像上的特征点坐标,计算正五棱台上底面和两个侧面上的正交隐消点,建立正交隐消点关于摄像机内参数的约束方程,线性求解摄像机内参数矩阵。利用本发明中的靶标可实现全自动标定,减少了标定过程中由测量引起的误差。隐消点是一种更简洁更全局化的基元,在标定过程中提高了精度。
Description
技术领域
本发明属于计算机研究领域,涉及一种用于求解摄像机内参数的一个底边是正五边形的正棱台靶标,利用平面上任意一个正五棱台作为标定模板,通过求解两幅图像上的正交隐消点的坐标,线性确定摄像机的内参数。
背景技术
计算机视觉的基本任务之一,就是从摄像机获得的二维图像信息出发恢复物体在三维空间中的几何信息,从而识别和重建三维空间中物体的几何形状。在此过程中必须确定空间物体点的三维几何位置与其图像中的对应点之间的相互关系,而这种关系又由摄像机成像的几何模型决定的,这些几何模型的参数就是摄像机参数。在大多数条件下,这些参数都是通过实验得到的,这就是摄像机标定。它一般分为传统标定和自标定两种方法,无论哪种标定方法,标定物体都是采用一些特殊的几何模型,例如:平面正方形、三角形、圆、空间立方体及圆柱等等,如何建立这些几何模型与摄像机参数之间的关系尤其是某种线性的关系,是目前摄像机标定所追求的目标,也是目前计算机视觉领域研究的热点之一。
传统的摄像机标定方法虽然可以获得较高的精度,但是标定块制作困难,不便于操作,针对这一问题文献“Aflexiblenewtechniqueforcameracalibration”,(Zhengyou Zhang,IEEETransactionsonPatternAnalysisandMachineIntelligence,vol.22,no.11,pp.1330-1334,2000.)提出了用平面模板代替传统标定块的方法,这种方法简单方便,成本低,并且能获得较高的精度,但需要精确定位模板上点阵的物理坐标,由于隐消点是一种更简洁更全局化的基元,因而可以进一步提高方法的稳定性,于是用隐消点解决标定问题被广泛使用,文献”CriminisiA,ReidI.andA.Zisserman.SingleViewMetrology,InternationalJournalofComputerVision,2000,40(2),pp.123-148.”利用隐消点的属性,从单一投影图像中计算摄像机的内参数,该方法首次将射影几何中的隐消点融入到摄像机标定中,于是隐消点成了摄像机自标定方法的理论基础。
隐消点是像平面上一个特殊的点,平面上某个方向上的所有平行直线都交于一个无穷远点,无穷远点在像平面上的投影称为隐消点。文献(王年,范益政,梁栋,韦穗,“一种基于等腰梯形的摄像机自标定方法”,东南大学学报(自然科学版),vol.35,no.2,pp.195-198,2005.)根据等腰梯形的性质和交比不变性计算图像中的隐消点,文献(胡培成,黎宁,赵亮亮,“一种新的基于消失点的摄像机标定方法”,第十三届全国图像图形学学术会议,pp.1-5,2006.)利用隐消点标定摄像机的方法,根据两条平行线的交点是无穷远点,直接求解两个平行线的像的交点来完成标定。
发明内容
本发明提供了一种制作简单,适用广泛,稳定性好的用于求解摄像机内参数的靶标。该靶标是一个底边的是正五边形的正棱台,在求解摄像机内参数的过程中,只需摄像机从不同方位拍摄2幅图像就可以线性求解出摄像机的5个内参数。
本发明采用如下技术方案:
本发明是由底边是正五边形的正棱台构成的用于摄像机自标定的靶标。具体的步骤包括:从图像上提取正五棱台的上底面和两个侧面上的特征点,根据两平行直线的交点是无穷远点,而无穷远点的像称为隐消点,以及共线四点交比的性质,得出图像平面上隐消点坐标,根据两幅图像上正交隐消点对绝对二次曲线的约束,线性求解摄像机的内参数。
1.提取图像中的特征点
在Matlab软件中利用Harris角点检测提取出每个图像平面特征点的坐标。
2.求隐消点的坐标
(1)计算正五棱台上底面的隐消点
在世界坐标系的一个平面上,存在一个正五棱台(如图1),正五棱台的上、下底面是一个正五边形(如图2),它上底面的角点为,下底面的部分角点为,是五边形内切圆的圆心。是线段的中点,垂直于,方向上的无穷远点设为;方向上的无穷远点设为.则与正交。上底面的投影如图3所示,分别是的像,分别是,的像,根据共线四点交比的性质有:,,根据射影变换中对应点的射影不变性,有:,。从而可以求得的坐标。
(2)计算正五棱台侧面的隐消点
在正五棱台的一幅图片中,我们至少可以看到它的两个侧面是两个全等的等腰梯形,如图4所示,平行于,分别是和的中点,并连接,是的中点,设,上的无穷远点为,方向上的无穷远点是,如图5,分别是的像,是无穷远点的像点,即:,方向上的隐消点,是无穷远点的像点,即:方向上的隐消点。因为是和的交点,所以,根据射影变换中对应点的射影不变性,有:,根据射影几何中调和共轭理论,得:,根据射影变换中的交比不变性,得,从而解得的坐标。在另一个等腰梯形中也存在着一组正交隐消点,记为,,计算公式为,。
3.求解摄像机内参数
利用在两幅图像中所求的六组正交隐消点对绝对二次曲线的约束可以列出六个方程,线性求解出摄像机的5个内参数,即矩阵,为图像的畸变因子,为图像坐标系中轴轴的尺度因子,是主点坐标,为摄像机的5个内参数。
本发明优点:
(1)该靶标制作简单,将一个正五棱锥平行于底面切割,就可以得出一个正五棱台。
(2)只需用摄像机从不同方位拍摄两幅图像便可线性求解出摄像机的5个内参数。
(3)方法容易移植到底边是奇数正多边形的正棱台上。
附图说明
图1用于求解摄像机内参数的靶标结构示意图。
图2是靶标上底面结构示意图。
图3是靶标顶面的平面隐消点求解原理图。
图4是靶标侧面结构示意图。
图5是靶标侧面隐消点的求解原理图。
具体实施方式
一种用于求解摄像机内参数的靶标,它是由空间中的一个底边是正五边形的正棱台构成的,如图1。用此新型靶标完成摄像机内参数的求解需要经过以下步骤:正交隐消点的摄像机标定方法采用的实验模板是空间中的任意的一个正五棱台,如图1所示。其上底面是一个正五边形,侧面是五个全等的等腰梯形,利用本发明中的方法对用于实验的摄像机进行标定,具体步骤如下:
1.拟合图像中的特征点
在Matlab软件中利用Harris角点检测法提取出图像特征点的坐标。
2.计算隐消点
(1)计算五棱台上底面的隐消点
在世界坐标系的一个平面上,存在一个正五棱台(如图1),正五棱台的上、下底面是一个正五边形(如图2),它上底面的角点为,下底面的部分角点为,是五边形内切圆的圆心。是线段的中点,垂直于,方向上的无穷远点设为;方向上的无穷远点设为.则与正交。上底面的投影如图3所示,分别是的像,分别是,的像,根据共线四点交比的性质有:,,根据射影变换中对应应点的射影不变性,有,,从而求得的坐标。
(2)计算正五棱台侧面上的隐消点
在正五棱台的一幅图片中,至少可以看到它的两个侧面是两个全等的等腰梯形,如图4所示,平行于,分别是和的中点,并连接,是的中点,设,上的无穷远点为,方向上的无穷远点是,如图5,分别是的像,是无穷远点的像点,即:,方向上的隐消点,是无穷远点的像点,即:方向上的隐消点。因为是和的交点,所以,根据射影变换中相应点的对应性,有:,根据射影几何中调和共轭理论,得:,根据射影变换中的交比不变性,得,从而解得的坐标。同理,在另一个等腰梯形中也能找到一组正交隐消点,记为,。计算公式为,得到,。
3.求解摄像机内参数
获取两幅图像,利用正交隐消点对绝对二次曲线的约束,线性求解出摄像机的5个内参数,即矩阵,为图像的畸变因子,为图像坐标系中、轴的尺度因子,是主点坐标。
实施例
本发明提出了利用平面上任意一个底边为正五边形的正棱台作为靶标线性确定摄像机的内参数,本发明采用的实验模块是正五棱台,结构示意图如图1所示,下面以一实例对本发明的实施方案做出更为详细的描述。
利用正五棱台求解摄像机内参数采用的实验模板是平面上任意一个正五棱台,如图1所示,利用本发明中的方法对用于实验的摄像机进行标定,具体步骤如下:
(1)拟合图像中的特征点
本发明采用的图像分辨率为640×480个像素,用摄像机从不同方向拍摄多幅实验图片,选取两幅较为清晰的图片,读入图像在Matlab中利用Harris角点检测提取出图像特征点的坐标。
第一幅图像上的特征点如下,为 的像:
第二副图像上的特征点如下,为 的像:
(2)计算正五棱台上底面图像的隐消点
在世界坐标系的一个平面上,存在一个正五棱台(如图1),正五棱台的上、下底面是一个正五边形(如图2),它上底面的角点为,下底面的部分角点为,是五边形内切圆的圆心。是线段的中点,垂直于,方向上的无穷远点设为;方向上的无穷远点设为.则与正交。上底面的投影如图3所示,分别是的像,分别是,,的像,根据共线四点交比的性质有:,,根据射影变换中对应点的射影不变性,有:,。从而可以求得的坐标。第二副图像上底面的一组正交隐消点记为:。计算公式为,,可得。
由以上方法可以计算出两幅图像上底面的隐消点,计算结果分别如下:
第一幅图像,上底面中的一组正交隐消点:
;
第二幅图像,上底面的一组正交隐消点:
。
(3)计算正五棱台侧面图像的隐消点
在正五棱台的一幅图片中,至少可以看到它的两个侧面是两个全等的等腰梯形,如图4所示,平行于,分别是和的中点,并连接,是的中点,设,上的无穷远点为,方向上的无穷远点是,如图5,分别是的像,是无穷远点的像点,即:,方向上的隐消点,是无穷远点的像点,即:方向上的隐消点。因为是和的交点,所以,根据射影变换中相应点的对应性,有:,根据射影几何中调和共轭理论,得:,根据射影变换中的交比不变性,得,从而解得的坐标。同理,在另一个等腰梯形中也能找到一组正交隐消点,记为,,计算公式为,。第二副图像两个侧面的两组正交隐消点分别记为:和,计算公式分别为,;,。
由以上方法计算出两幅图像侧面的隐消点,计算结果分别如下:
第一幅图像两侧面的两组正交隐消点:
第二副图像两侧面上的两组正交隐消点:
(4)求解出摄像机内参数
经过以上步骤得到的两幅图像上的共六组正交隐消点的坐标,可线性地求解出摄像机内参数矩阵
五个内参数(单位:像素)分别为:,,,和。
Claims (1)
1.一种用靶标求解摄像机内参数的方法,其特征在于所述靶标是将空间中任意一个正五棱锥平行于底面切割得到,所述方法具体的步骤包括:
(1)计算正五棱台上底面的隐消点
在世界坐标系的一个平面上,存在一个正五棱台,正五棱台的上、下底面是一个正五边形,它上底面的角点为A,B,C,D,E,下底面的对应上底面C,D,E的角点为C1,D1,E1,O是上底面正五边形内切圆的圆心;F是线段CD的中点,AF垂直于CD,AF方向上的无穷远点设为P1∞;CD方向上的无穷远点设为P2∞,则P1∞与P2∞正交;在上底面的投影中,a,b,c,d,e,o特征点分别是A,B,C,D,E,O的像,p1,p2分别是P1∞,P2∞的像,根据共线四点交比的性质有:(AF|OP1∞)=-cos(π/5),(CD|FP2∞)=-1,根据射影变换中对应点的射影不变性,有:(af|op1)=-cos(π/5),(cd|fp2)=-1;从而可以求得p1,p2的坐标;
(2)计算正五棱台侧面的隐消点
在摄像机拍摄的正五棱台的一幅图片中包含正五棱台的两个侧面,并且这两个侧面是两个全等的等腰梯形;一个等腰梯形中,CD平行于C1D1,F,F1分别是CD和C1D1的中点,并连接F,F1,G是FF1的中点,设CD,C1D1上的无穷远点为P3∞,FF1方向上的无穷远点是P4∞;c,d,c1,d1,f,f1,g特征点分别是C,D,C1,D1,F,F1,G的像,p3是无穷远点P3∞的像点,即:CD,C1D1方向上的隐消点,p4是无穷远点P4∞的像点,即:FF1方向上的隐消点;因为P3∞是CD和C1D1的交点,用叉积计算,表示为:P3∞=(CD)×(C1D1),根据射影变换中对应点的射影不变性,用对应点构成直线的叉积计算,表示为:p3=(cd)×(c1d1),根据射影几何中调和共轭理论,得:(FF1|GP4∞)=-1,根据射影变换中的交比不变性,得(ff1|gP4)=-1,从而解得p4的坐标;在另一个等腰梯形DD1EE1中也存在着一组正交隐消点,记为p5,p6,计算公式表示为:p5=(de)×(d1e1),(ii1|g1p6)=-1,求出p5,p6;其中,d,e,d1,e1,i,i1,g1特征点分别是D,E,D1,E1,I,I1,G1的像;
(3)求解摄像机内参数
利用摄像机不同角度拍摄的两幅图像中所求的六组正交隐消点对绝对二次曲线的约束可以列出六个方程,线性求解出摄像机的5个内参数,即矩阵 s为图像的畸变因子,fu,fv为图像坐标系中u轴v轴的尺度因子,(u0,v0)是主点坐标,即s,fu,Fv,u0,v0为摄像机的5个内参数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210542006.2A CN102982550B (zh) | 2012-12-14 | 2012-12-14 | 利用正五棱台求解摄像机内参数 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210542006.2A CN102982550B (zh) | 2012-12-14 | 2012-12-14 | 利用正五棱台求解摄像机内参数 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102982550A CN102982550A (zh) | 2013-03-20 |
CN102982550B true CN102982550B (zh) | 2016-01-06 |
Family
ID=47856491
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210542006.2A Expired - Fee Related CN102982550B (zh) | 2012-12-14 | 2012-12-14 | 利用正五棱台求解摄像机内参数 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102982550B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103983961A (zh) * | 2014-05-20 | 2014-08-13 | 南京理工大学 | 一种3d激光雷达和摄像机联合标定立体标定靶 |
GB2557633A (en) * | 2016-12-14 | 2018-06-27 | Fuel 3D Tech Limited | Method of obtaining data characterizing 3D-imaging equipment |
CN108335333A (zh) * | 2018-04-03 | 2018-07-27 | 哈尔滨理工大学 | 一种线性摄像机标定方法 |
CN110132132B (zh) * | 2018-12-11 | 2021-08-06 | 中国航空工业集团公司北京长城计量测试技术研究所 | 一种多杆梯形六面体标准器 |
CN110163922B (zh) * | 2019-05-22 | 2020-06-23 | 四川深瑞视科技有限公司 | 鱼眼相机标定系统、方法、装置、电子设备及存储介质 |
WO2020232971A1 (zh) | 2019-05-22 | 2020-11-26 | 四川深瑞视科技有限公司 | 鱼眼相机标定系统、方法、装置、电子设备及存储介质 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102110290A (zh) * | 2009-12-28 | 2011-06-29 | 云南大学 | 一种利用正三棱柱靶标求解摄像机内参数方法 |
CN102402785A (zh) * | 2010-09-07 | 2012-04-04 | 云南大学 | 一种基于二次曲线的摄像机自标定方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR100386090B1 (ko) * | 2001-04-02 | 2003-06-02 | 한국과학기술원 | 동심원 패턴을 이용한 카메라 내부변수 보정시스템 및카메라 보정방법 |
-
2012
- 2012-12-14 CN CN201210542006.2A patent/CN102982550B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102110290A (zh) * | 2009-12-28 | 2011-06-29 | 云南大学 | 一种利用正三棱柱靶标求解摄像机内参数方法 |
CN102402785A (zh) * | 2010-09-07 | 2012-04-04 | 云南大学 | 一种基于二次曲线的摄像机自标定方法 |
Non-Patent Citations (2)
Title |
---|
A flexible new technique for camera calibration;zhengyou zhang等;《IEEE Transactions on Pattern Analysis and Machine Intelligence》;20001130;第22卷(第11期);第1330-1334页 * |
基于消隐点几何特性的摄像机自标定方法;霍炬等;《光学学报》;20100215;第30卷(第2期);第465-472页 * |
Also Published As
Publication number | Publication date |
---|---|
CN102982550A (zh) | 2013-03-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102982550B (zh) | 利用正五棱台求解摄像机内参数 | |
CN103035007B (zh) | 利用正六棱台求解摄像机内参数 | |
CN102930548B (zh) | 利用两个相同的相交椭圆线性求解摄像机内参数 | |
CN102982551B (zh) | 空间三条不平行直线线性求解抛物折反射摄像机内参数 | |
CN102136140B (zh) | 一种基于矩形图样的视频图像距离检测方法 | |
CN104835144A (zh) | 利用一个球的球心的像及正交性求解摄像机内参数 | |
CN103106661B (zh) | 空间二条相交直线线性求解抛物折反射摄像机内参数 | |
CN102930551B (zh) | 利用圆心的投影坐标和极线求解摄像机内参数 | |
CN104034269A (zh) | 一种单目视觉测量方法与装置 | |
CN109544628A (zh) | 一种指针式仪表的准确读数识别系统及方法 | |
CN103440638A (zh) | 利用双平面镜装置和圆环点性质求解摄像机内参数 | |
CN105513063A (zh) | Veronese映射和棋盘格标定拋物折反射摄像机 | |
CN103116892A (zh) | 两个相交相同圆及公切线求解摄像机内参数 | |
CN104200477A (zh) | 基于空间平行圆求解平面折反射摄像机内参数的方法 | |
CN102999895B (zh) | 利用两个同心圆线性求解摄像机内参数 | |
CN103400374B (zh) | 利用双平面镜装置和拉盖尔定理的推论求解摄像机内参数的方法 | |
CN103942784A (zh) | 棋盘格中三条互不平行直线求解抛物折反射摄像机内参数 | |
CN102012213A (zh) | 单副图像测量前景高度的方法 | |
CN103116888A (zh) | 利用平面三角形求解摄像机的内参数 | |
CN102999894B (zh) | 利用两个相交圆线性求解摄像机内参数 | |
CN102110290A (zh) | 一种利用正三棱柱靶标求解摄像机内参数方法 | |
CN202257686U (zh) | 求解摄像机内参数的两个相交圆标定物 | |
CN102982549B (zh) | 两个共中心和主轴正交的相同相交椭圆求解摄像机内参数 | |
CN104182986A (zh) | 由正交隐消点的性质求解平面折反射摄像机内参数的方法 | |
CN103810697A (zh) | 空间四条不平行直线的像标定抛物折反射摄像机内参数 |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160106 Termination date: 20161214 |
|
CF01 | Termination of patent right due to non-payment of annual fee |