基于两灭点法的摄像机标定方法研究_胡桂廷
上传者:解蕾蕾|上传时间:2015-05-05|密次下载
基于两灭点法的摄像机标定方法研究_胡桂廷
机器人;视觉;控制
理论与算法
电 子 测 量 技 术
LECTRONIC MEASUREMENTECHNOLOGY E T
第35卷第7期2012年7月
基于两灭点法的摄像机标定方法研究
胡桂廷 陈向东
()北京精密机电控制设备研究所 北京 100621
、摘 要:在计算机视觉当中,要对物体进行三维重建,需要实现对相机进行标定,即获得相机的的焦距(内参数)旋转。为提高场景摄像机的标定精度和标定效率,外参数)研究了一种基于2个灭点的快速摄像机标矩阵以及偏移向量(
只要求被标定相机拍摄的图像中有2组相互正交的平行线,且这2组平行线与成像平定方法。该方法利用灭点理论,
面不平行,由此可以获得2个灭点;为计算平移向量,图像中某个平行线方向上的一段直线的实际长度需已知。根据灭点的几何投影特性计算相机内外参数。利用监控摄像机进行试验验证,验证结果表明,该方法相对于传统的标定方法操作简单并且精度较高,具有一定的理论意义和实用价值。关键词:相机标定;灭点;内外参数中图分类号:TP391 文献标识码:A
Cameracalibrationmethodbasedontwovanishinoints pg
HuGuitinhenXiandon C ggg
(,)BeiinResearchInstituteofPrecisionMechanical&ElectricalControlEuimentBeiin100621 jgqpjg
:,)arameterAbstractIncomutervisioncalibrationwhichconcludesfocus(interiorandrotatematrixandtranslation pp)vector(exteriorarameteracuisitionisthefundamentalworkinthethreedimensionalreconstructionofobects.To - pqj,aerroosesimrovetheaccuracandsimliffieldcalibrationstestheanfastcameracalibrationmethodusintwo pppppypypg oints.Itadotsanvanishinointsmethodtoettheinteriorandexteriorarameters.Onlthefollowinvanishin ppgpgpygg :;assumtionsaremadethesinleimaecontainsattwovanishinointswhichisetattwoverticaldirectionthelenth pgggpgg ())in3Dsaceofonelinesement(fordetermininthetranslationvectorintheimaeisknown.Resultsofexeriment pgggp recisionwithansurveillancecameraconfirmstheconvenienceandofthealorithmcomaredwiththecarried pgptraditionalones.
:c;;Kewordsameracalibrationvanishinointsinteriorandexteriorarameters gppy
]811-
。E.足,因此其应用范围受到制约[Guillou在文献[12]
0 引 言
在计算机视觉当中,要实现物体三维重建,需求解空间某点的三维位置信息与图像对应点之间的关系,这种关系由摄像机的参数来决定。求解相机参数的过程称为相机标定。相机标定发展至今,有传统标定方法、基于主动视觉的自标定方法和基于灭点的方法标定方法、
[1]
中提出了一种基于单幅图像中2个灭点的标定方法,虽然但假设相机的主点是图像的中该方法对场景要求降低了,心,导致计算精度下降。
本论文研究了基于单幅图像中存在2个灭点的标定方法的数学模型,利用这种方法进行摄像机的内外参数的标定。
。传统的标定
方法精度较高,但它需要测量精度高、结构已知的标定场或是标定物作为参照,典型的算法有directlineartranslation
[]
(方法,主动视觉标定方法需要知道DLT)Tsai两步法2-3;
]45-
;摄像机的运动信息,例如摄像机平移或旋转量[自标定
1 摄像机标定数学模型
理想情况下,相机成像满足针孔成像模型。假设空间中任意一点P在世界坐标系、相机坐标系和图像坐标系下的关系为:
仅依靠摄像机在运动过程中周围方法是指不需要标定物,
环境及图像间的对应关系来求取摄像机参数的方法
[]67-
。
一般情况下,传统的基于灭点的标定算法要求场景中存在3组相互正交的平行线,该条件在实际应用中很难满
012年3月收到。 本文于2
·37·
机器人;视觉;控制
5卷 第3
电子测量技术
oZo
Xc
()1
cZc
=
XcXo
2.2 灭点法计算摄像机的焦距
在基于灭点理论计算相机参数时,要求拍摄的场景图且与成像平面不平片里可以找到2组相互正交的平行线,
行。设从中提取出2组直线构成空间矩形A它在图像BCD,中的投影为a则图像中包含2个灭点:bcd,Fd和u是直线a
[
Ro→cT
01×31
11 u00f s其示意图如图1bc
内容需要下载文档才能查看构成的灭点,Fb和cd构成的灭点,v是ac
()rv020=λf 0
所示。Zc
10 0 1 0式中:RTo→c为世界坐标系到摄像机坐标系下的旋转矩阵,
u
为平移向量,二者由摄像机的外参数决定;Sf为相机焦距,(为畸变因子,为主点的图像坐标,uvr为纵横比,λ为比0,0)例因子,这些参数由摄像机的内参数决定。本论文在假设S()近似为0,的情况下,推导后续r近似为1,uv0,0=(0,0))。快速标定算法的数学模型。此时,模型变换为式(3
Xc
f 000cλ=0f 00
Zc
10010 1标定就是要计算f,T,Ro→c。u
()3
图1 灭点法示意
现在流行的标定算法有传统的线性变化方法、张正友法以及基于灭点的快速标定算法。文章研究了基于灭点的该方法只需要采用摄像机在固定摄站坐标处拍标定方法:
摄1张图片,只要求图像中存在2组相互正交的平行线以及坐标轴上1条长度已知的直线,大大简化了实际操作步骤。
将Auvw记为世界坐标系Ro;Rij)为二维图像坐标s(
为摄像机坐标系。规定V/R为向量V在R坐标系;ROik)jc(系下的表示。假设摄像机中心O点在图像上的投影为P点。
为计算焦距f,过P点做直线Fu交点为PuFv的垂线,v。根据边的长度关系可得:OPPuΔv为直角三角形,
()P=PuPu6f=Ov-Pv
在图1中,因为A在像平面则有OB⊥AD,Fu⊥OFv,因为灭点具有非常重要性质为相机光心和灭点间连线与形
2 基于灭点理论的摄像机标定
2.1 灭点的基本理论
由于相机成像满足透视原理,因此空间上2条平行线在成像平面上的投影会在图像上有1个交点,称为灭点。假设三维空间中有Q0点,方向矢量为d,D是空间中该方向上过Q0点的一条直线。该直线可Q是直线上的任意一点,以用Q0和d表示为:
即O因成该灭点的空间平行线平行,Fu∥AD,OFv∥AB。此可得,直线FFv⊥OPuFFv⊥PPuFV⊥OP,uv、uv以及Fu
而且Δ根据相OFuFv、OFuPuOFvPuΔΔv、v均为直角三角形。似关系可得:
Q=Q0+l×d
式中:l为Q0和Q之间的距离。
()4
OPuPuFuv=vv×Puv
)),将式(代入式(计算出相机的焦距f为:76
()7
将三维空间中的点Q0和Q投影到像平面上获得像点。灭点所对应的空间位置为无穷远,即l为无穷大,图q0和q
像平面上直线D对应的灭点坐标FD(x∞,z∞)为:y∞,
·dax+ldxx
x∞=lim=fl→∞fdaldzz+z·dadylyyim=∞=lfl→∞f·
()5
2
()P=PuFu-PPu8f=Ovv×Puvv
2.3 旋转矩阵的计算
从世界坐标系R需要计算o到相机坐标系Rc的转换中,
设A为Ro的原点,出旋转矩阵R和平移向量T。3个坐标轴;的单位向量为(u,v,w)O为Rc的原点,3个坐标轴的单位
,。向量为(ik)P为图像坐标系的原点,2个坐标轴的单位j,,。向量为(在相机坐标系下OiFu和OFv方向的单位向量j)分别为u则有w由此可确定的坐标系为′和v′,′=u′×v′,由灭点的性质可知R和RR′。u,v,w)′(u′,v′,w′)3个坐oo(o
标轴互相平行,因此R′之c和Ro之间的变换矩阵与Ro和Ro间的变换矩阵相同,即:
aldz+zdz
∞=f)可以看出,从式(灭点的图像坐标只与对应空间直5线的方向矢量有关,而且三维空间中与直线D平行的直线有且只有1个灭点FD。
·38·
机器人;视觉;控制
基于两灭点法的摄像机标定方法研究胡桂廷等:
′=′=
→OFuFu‖‖O→OFvFv‖‖O
==
第7期
Fuj
fFu‖‖OFu‖‖OFu‖‖OFvi
Fui
Fvj
fFv‖‖OFv‖‖OFv‖‖O
)
)
出相机的参数。
在实验当中,利用标定信息以及1~8号点为验证点,二维图像信息,计算每个点的高程信息。实验误差统计分()析结果如图3所示。其中图3为标定计算值和测量值位a()图3为误差曲线。可以看出反算高度误差在置对比图,b实际高度越高的点,误差相比来说稍微高一点。5~10cm,
误差最大的点是点4和5,原因是这2个点距离摄像头最远,所以其每个像素代表的实际距离较大,导致误差偏大。
w′=u′×v′()9
),),在世界坐标系下,已知u=(1,0,0w=(0,0,1
。则有:R(u,v,w)=(u′,v′,w′)
R=
Fui
222
ui+Fuj+fFvi
222
vi+Fvj+fw′i
w′j
Fuj
2
ui
Fvj
+f
2
+F2
uj
2
vi+Fvj+f
22
222ui+Fuj+f2.4 平移向量
w′k
222F++fvivj()10
设A点为世界坐标系的原点,AP为平行于世界坐标系u方向的1段长度已知为L的线段,其在图像上的投影是
,。图像坐标为A′′A′P′,′x′P′x′//yyR=(a,a)R=(p,p)
S
S
过A点做平行于AP的直线L′交OP直线于点P″,
P″OP)′。根据已计算出的旋转矩阵可得//R=(R∩L
→→利用三角形OAP/RP/A′P″和OAP相似可以求Rc=o→cARo。
C
C
→出O所以有:A/Rc,
→→OA/A/Rc=‖ORc‖
→OA/RcA/‖ORc‖
()11
→)计算出O由式(11A/Rc即为平移向量。
3 实验验证结果分析
3.1 实验环境与条件
在某电梯间进行标定实验,使用Photoshop软件手动提取图像坐标,基准三维点的世界坐标由
内容需要下载文档才能查看皮尺测量。实验图像如图2所示,图中的1~8号点为选取的验证点。
图3 验证点误差分析
3.3 验证结果分析
实验场景区域约为5在计算高度时所00cm×900cm,标准差为4最有验证点误差的期望为57.16mm,6.14mm,大误差为1实验的误差值均较63mm。与场景的面积相比,小。随着验证点与相机距离的增加,误差值成递增趋势。从实验结果中可以看出,本文的快速场景标定方法具有较好的标定精度,同时操作简单,可以做到实时在线标定。
4 结 论
基于灭点理论的快速相机标定方法,首先利用通过轴向的2组平行线信息解出相机焦距和旋转矩阵,并由轴上的任
图2 实验场景
意1条原点出发的线段信息求出平移向量,实现对相机的快速标定既保证了标定精度和可靠性,又提高了标定效率。
(下转第54页)
3.2 验证结果
利用本文当中基于两灭点理论的摄像速标定算法计算
·39·
机器人;视觉;控制
5卷 第3
电子测量技术
nlodalnalsissinlindourceOutut o a u syygp m b
searationtechniues[J].MechanicalSstemand pqy,():Processin20072123352358.Sinal -gg[]李舜酩,张袁元.基于S8ONS的结构模态参数 刘晓伟,
:识别研究[J].中国机械工程,2011,22(12)
14541458.-
[]]虚拟仪器技术及其新进展[仪器仪表学报,9J. 宁涛.
():2007,284252254.-
作者简介
别的结果特征较明确,但出现了系统载荷的虚假分量,说明所设计系统仍需要进一步改进,以提高其对干扰的鲁棒性。
参考文献
[]陈怀海,鲍明.机械振动载荷识别研究的现状与1 许峰,
]():未来[中国机械工程,J.2002,136526531.-
[]率志君,李玩幽,等.设备载荷识别与激励源2 胡寅寅,
]():噪声与振动控制,特性的研究现状[J.2011415.-
[]智浩,缪龙秀.基于模态分析法的结构动载荷3 文祥荣,
]():识别研究[北方交通大学学报,J.2000,2441114.-
[]陈怀海,鲍明.基于模态滤波器改进的动载荷识4 许锋,
]():别模型[航空学报,J.2002,2315154.-
[]饶柱石,周海亭.一种动态载荷的识别方法5 傅志方,
[]():上海交通大学学报,J.1997,31357.-
[6]HOU W L,CHELIDZED.Blindsourcesearation Z p
basedvibrationodedentification[J].Mechanical m i
,():SstemandSinalProcessin20072130723087. -ygg
[7]onceletF,KEEHEN G,GOLINVALJC,etal. P
马赛,男,硕博连读生,主要研究方向为1986年出生,振动信号处理。
:E-mailmassana63.com@1
李舜酩,男,南京航空航天大学能源与动力1962年生,学院教授,博士研究生导师。
郭海东,男,硕士生,主要研究方向为航1985年出生,空发动机振动信号处理和故障诊断等。
櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒
(上接第39页)
参考文献
[]张书毕,刘超.基于灭点理论和平面控制场的1 杨化超,
]:相机标定方法研究[中国图象图形学报,J.2010(8)5762.-
[]徐杜,蒋永平.计算机三维视觉测量中相机标2 陈济棠,
]():定的圆形校正[电子测量技术,J.201138782.-
[]卢荣胜,刘宁,等.多G3iE相机高速视觉系统 李琪,g
]电子测量与仪器学报,的架构与处理速度研究[J.
():2010,244371378.-
[],4D H,LUOJHU YF,etal.Themethodof WANG
lineartransformationonthecameracalibrationdirect [],():J.Machiner&Electronics20069911.-y
[]张祖勋,张剑清.利用二维D5LT及光束法平 张永军,
]差进行数字摄像机标定[武汉大学学报:信息科学J.():版,20026128134.-
[]6FCH,HUZY.Lineardeterminationoftheinfinite WU
homorahndameraelfcalibration[J].Acta c s-gpy a
,,():AutomaticaSinica2002284488496. -
[]徐刚,董增来.引入灭点约束的T7SAI两步法相 吴军,
]机标定改进研究[武汉大学学报:信息科学版,J.():201217884.-
[]]摄相机的自标定:新方法及比较[贵州科8J. 吴有富.
学.2009,3:123127-
[]9hanZY.Aflexiblenewtechniueforcameracalibration Z gq
[]J.IEEETransactionsonPatternAnalsisandMachine y,,():Intellience2000221113301334.-g
[]杨进华,赵群.基于图像匹配的长焦透镜焦距测10 申琳,
]量研究[J.2010Vol.24No.11P.10431047.-
[],11HENGFJYU T,YUAN GT,etal.Overviewof Z
[]eometriccameracalibrationJ.TaiuanScienceand gy
,():Technolo201027273.-gy
[12]ILLOUE,MENEVEAUXD,MAISELE,etal. GU
Usinvanishinforcameracalibrationandoint ggp []coarse3DreconstructionfromasinleimaeJ.The gg,():VisualComuter2000,167396410. -p
作者简介
胡桂廷,男,高级工程师,主要从事工厂1965年出生,自动化控制及检测系统设计与研究。
:E-mailhuuitin26.com@1gg
陈向东,男,高级工程师,主要从事工厂1972年出生,自动化控制及检测系统设计与研究。:E-mailcxd091563.com@1
·54·
下载文档
热门试卷
- 2016年四川省内江市中考化学试卷
- 广西钦州市高新区2017届高三11月月考政治试卷
- 浙江省湖州市2016-2017学年高一上学期期中考试政治试卷
- 浙江省湖州市2016-2017学年高二上学期期中考试政治试卷
- 辽宁省铁岭市协作体2017届高三上学期第三次联考政治试卷
- 广西钦州市钦州港区2016-2017学年高二11月月考政治试卷
- 广西钦州市钦州港区2017届高三11月月考政治试卷
- 广西钦州市钦州港区2016-2017学年高一11月月考政治试卷
- 广西钦州市高新区2016-2017学年高二11月月考政治试卷
- 广西钦州市高新区2016-2017学年高一11月月考政治试卷
- 山东省滨州市三校2017届第一学期阶段测试初三英语试题
- 四川省成都七中2017届高三一诊模拟考试文科综合试卷
- 2017届普通高等学校招生全国统一考试模拟试题(附答案)
- 重庆市永川中学高2017级上期12月月考语文试题
- 江西宜春三中2017届高三第一学期第二次月考文科综合试题
- 内蒙古赤峰二中2017届高三上学期第三次月考英语试题
- 2017年六年级(上)数学期末考试卷
- 2017人教版小学英语三年级上期末笔试题
- 江苏省常州西藏民族中学2016-2017学年九年级思想品德第一学期第二次阶段测试试卷
- 重庆市九龙坡区七校2016-2017学年上期八年级素质测查(二)语文学科试题卷
- 江苏省无锡市钱桥中学2016年12月八年级语文阶段性测试卷
- 江苏省无锡市钱桥中学2016-2017学年七年级英语12月阶段检测试卷
- 山东省邹城市第八中学2016-2017学年八年级12月物理第4章试题(无答案)
- 【人教版】河北省2015-2016学年度九年级上期末语文试题卷(附答案)
- 四川省简阳市阳安中学2016年12月高二月考英语试卷
- 四川省成都龙泉中学高三上学期2016年12月月考试题文科综合能力测试
- 安徽省滁州中学2016—2017学年度第一学期12月月考高三英语试卷
- 山东省武城县第二中学2016.12高一年级上学期第二次月考历史试题(必修一第四、五单元)
- 福建省四地六校联考2016-2017学年上学期第三次月考高三化学试卷
- 甘肃省武威第二十三中学2016—2017学年度八年级第一学期12月月考生物试卷
网友关注
- 《影视摄影》教案
- 历时一个月的宿迁高师艺体节即将落下帷幕
- 《1600词》2-5
- Chapter7+the+application+layerWord
- 两学一做专题教育学习动态Word模板动态模板
- 《1600词》2-7
- 演示文稿10肖像和姓名中的权利Word
- 最新电大专科《中国教育简史》机考网考纸考题库及答案
- 最新电大《网络应用服务管理》机考网考纸考题库及答案
- 第一部分-古希腊与罗马的古典政治哲学
- 第三章_辛亥革命与君主专制制度的终结2016Word
- 生产总值的计算
- 运动学4、质点的圆周运动与螺旋运动
- 骆驼祥子Word
- 市场营销教案[15]Word
- 07-晶体结构 (2)(结晶学与矿物学)Word
- 动漫毕业设计论文
- 懂得尊重身边每一个人
- 第六部分-《联邦党人文集》与美国政治制度设计
- 第八章 蛋白质的测定Word
- 2015年春季开学典礼
- 西方经济学-高鸿业第五版微观部分-南开课件Word
- 大学生谈恋爱弊大于利四辩稿
- 2014~2016福建省高一数学竞赛真题
- 2017年广东省中学生生物学联赛试卷答案
- 《1600词》2-12
- 诗词鉴赏汇报Word
- 发热3节课Word
- 简报
- 《1600词》2-10
网友关注视频
- 【部编】人教版语文七年级下册《逢入京使》优质课教学视频+PPT课件+教案,安徽省
- 沪教版八年级下册数学练习册一次函数复习题B组(P11)
- 七年级英语下册 上海牛津版 Unit5
- 外研版英语七年级下册module3 unit2第二课时
- 北师大版数学 四年级下册 第三单元 第二节 小数点搬家
- 沪教版牛津小学英语(深圳用) 六年级下册 Unit 7
- 北师大版小学数学四年级下册第15课小数乘小数一
- 外研版英语七年级下册module3 unit1第二课时
- 【部编】人教版语文七年级下册《逢入京使》优质课教学视频+PPT课件+教案,辽宁省
- 8.对剪花样_第一课时(二等奖)(冀美版二年级上册)_T515402
- 每天日常投篮练习第一天森哥打卡上脚 Nike PG 2 如何调整运球跳投手感?
- 飞翔英语—冀教版(三起)英语三年级下册Lesson 2 Cats and Dogs
- 【部编】人教版语文七年级下册《老山界》优质课教学视频+PPT课件+教案,安徽省
- 沪教版牛津小学英语(深圳用) 五年级下册 Unit 12
- 冀教版小学数学二年级下册第二单元《有余数除法的简单应用》
- 二次函数求实际问题中的最值_第一课时(特等奖)(冀教版九年级下册)_T144339
- 沪教版八年级下册数学练习册21.3(3)分式方程P17
- 沪教版牛津小学英语(深圳用) 四年级下册 Unit 2
- 第五单元 民族艺术的瑰宝_15. 多姿多彩的民族服饰_第二课时(市一等奖)(岭南版六年级上册)_T129830
- 外研版八年级英语下学期 Module3
- 沪教版牛津小学英语(深圳用)五年级下册 Unit 1
- 沪教版牛津小学英语(深圳用) 五年级下册 Unit 7
- 第4章 幂函数、指数函数和对数函数(下)_六 指数方程和对数方程_4.7 简单的指数方程_第一课时(沪教版高一下册)_T1566237
- 沪教版八年级下次数学练习册21.4(2)无理方程P19
- 3月2日小学二年级数学下册(数一数)
- 冀教版小学数学二年级下册第二周第2课时《我们的测量》宝丰街小学庞志荣.mp4
- 冀教版小学数学二年级下册第二单元《有余数除法的竖式计算》
- 【部编】人教版语文七年级下册《泊秦淮》优质课教学视频+PPT课件+教案,广东省
- 8 随形想象_第一课时(二等奖)(沪教版二年级上册)_T3786594
- 3.2 数学二年级下册第二单元 表内除法(一)整理和复习 李菲菲
精品推荐
- 2016-2017学年高一语文人教版必修一+模块学业水平检测试题(含答案)
- 广西钦州市高新区2017届高三11月月考政治试卷
- 浙江省湖州市2016-2017学年高一上学期期中考试政治试卷
- 浙江省湖州市2016-2017学年高二上学期期中考试政治试卷
- 辽宁省铁岭市协作体2017届高三上学期第三次联考政治试卷
- 广西钦州市钦州港区2016-2017学年高二11月月考政治试卷
- 广西钦州市钦州港区2017届高三11月月考政治试卷
- 广西钦州市钦州港区2016-2017学年高一11月月考政治试卷
- 广西钦州市高新区2016-2017学年高二11月月考政治试卷
- 广西钦州市高新区2016-2017学年高一11月月考政治试卷
分类导航
- 互联网
- 电脑基础知识
- 计算机软件及应用
- 计算机硬件及网络
- 计算机应用/办公自动化
- .NET
- 数据结构与算法
- Java
- SEO
- C/C++资料
- linux/Unix相关
- 手机开发
- UML理论/建模
- 并行计算/云计算
- 嵌入式开发
- windows相关
- 软件工程
- 管理信息系统
- 开发文档
- 图形图像
- 网络与通信
- 网络信息安全
- 电子支付
- Labview
- matlab
- 网络资源
- Python
- Delphi/Perl
- 评测
- Flash/Flex
- CSS/Script
- 计算机原理
- PHP资料
- 数据挖掘与模式识别
- Web服务
- 数据库
- Visual Basic
- 电子商务
- 服务器
- 搜索引擎优化
- 存储
- 架构
- 行业软件
- 人工智能
- 计算机辅助设计
- 多媒体
- 软件测试
- 计算机硬件与维护
- 网站策划/UE
- 网页设计/UI
- 网吧管理