专利名称：：采用两个相同圆进行摄像机标定的方法
技术领域：
：本发明涉及一种采用两个相同圆进行摄像机标定的方法，在所有摄像机参数均未知的情况下能够准确获得摄像机的内部参数和外部参数。本发明属于先进测量
技术领域：
，尤其适用于计算机立体视觉中的摄像机标定、非接触式工业检测和基于视觉的自主导航系统。
背景技术：
：摄像机标定是计算机视觉的基础，在机器视觉、虚拟现实等计算机图像应用领域有着广泛的应用，其主要目标是求出摄像机的内部参数和外部参数。由于圆形标志物容易获取、对阈值分割的不敏感性、可遮挡性以及亚像素定位精确性，基于圆形标志物的标定方法得到广泛应用。文献(HEIKKILAJ."Geometriccameracalibrationusingcircularcontrolpoints",Trara.PflWem爿w。/.Mac/zz"e/"/W/，2000,22(10):1066~1077)是基于圆形标志物的传统标定方法，其标定精度高，抗噪性能好，但是传统标定是利用标准参照物与其对应图像的约束关系来获得摄像机的参数的，而标准参照物的设置给标定带来了许多不便。相对于传统标定而言，自标定是一种灵活、简单的方法，它不需要标定物的实际坐标位置即可求出摄像机参数。近年来，不少学者对多个圆进行自标定进行了广泛研究。吴毅红利用了圆环点的特性，通过从不同角度拍摄三幅包含两个圆的模板图像，提出了基于两平《亍圆(WUYH，ZHUH."Cameracalibrationfromthequasi-affineinvarianceoftwoparallelcircles",ProceedingofECCV'04:EuropeanConferenceonComputerVision,2004:190—202)和两共面圆(WUYH,LIX,WUFC,""Z."Coplanarcircles,quasi-a伍neinvarianceandcalibration",imageam/FfswwCo附/wri"g,2006,24(4):319—326)的摄f象机内部参数的求取方法，但并没给出摄像机外部参数的获取方法。Chen等(CHENQ，WUH,WADAT."CameraCalibrationwithtwoarbitrarycoplanarcircles",/VoceeJ/"go/五CCK'O(Co"/^^wceo"Oww/w/"P7j/ow,2004:521-532)从解析几何的角度利用两个共面圆获取摄像机的内外参数，并给出了迭代优化算法，但是文中推导的解析几何方程比较复杂，并且标定过程需要已知主点坐标和纵横比等内参信息，而且只有当两圆中一圆的半径信息或者两圆圆心间距已知时才能完全恢复摄像机参数中的平移向量信息。Gurdjos等(GURDJOSP，STURMP，衡YH."EuclideanstructurefromN>=2parallelcircles:theoryandalgorithms",五CCT加￡^o/7ea7iCow，ewceCWp她rP7s/o"，2006:238-252)通过辨别圆环点像分析了平行圆的欧几里得结构，为基于平行圆或共面圆的标定奠定了坚实的基础；受到Gurdjos的启发，Zheng(ZHENGY，LIUY."Cameracalibrationusingoneperspectiveviewoftwoarbitrarycoplanarcircles",Qp"'ca/五"g/ween.wg,2008,47(6):067203)利用圆环点像及射影几何中二次包络线的知识，求取了摄像机的内外参数，避免了Chen釆用复杂的解析几何方法，但在求取旋转阵时要利用两圆圆心的真实投影位置，由于透视投影偏差(HEIKKILAJ,SILVENO."Afour-stepcameracalibrationprocedurewithimplicitimagecorrection",iVocgecz'wgo/Co，她rSoc/砂Co"，ececmCo附/WwPfe/朋a"t/i^股raiecogm'"o",1997:1106—1112)，圆心的投影点并非圆所成像的拟合圆心，而文中并没给出解决方法；此外，与Chen方法一样，标定过程仍然需要已知主点坐标和纵横比等内参信息，而且只有当两圆中一圆的半径信息或者两圆圆心间距已知时才能完全恢复摄像机参数中的平移向量信息。上述所有方法的核心思想都是利用圆环点和绝对二次曲线的不变性。但是由于圆的高度对称性，基于圆标志物的自标定方法往往只能实现摄像机内参的标定，然而计算机的三维重建中，摄像机的外部参数对于三维点云的恢复是至关重要的。目前已有方法不是假定部分内参信息己知，就是必须借助如直线等其他几何辅助物方能同时确定摄像机的内外参数。鉴于两相同圆容易绘制且自然界普遍存在，而透视投影变换保持相切不变，因此，两圆的公切线对于求解摄像机的参数很有意义。
发明内容技术问题针对现有技术所存在的缺点和限制，本发明的目的在于提供一种采用两个相同圆进行摄像机标定的方法，能够简便而精确地求取摄像机的内部参数和外部参数，而且整个标定过程无需人工干预。技术方案本发明设计一种采用两个相同圆进行摄像机标定的方法，该方法所采用的标定物是两个不相切的相同圆；将此标定物放置于摄像机视域内，用摄像机从三个不同角度拍摄标定物，得到三幅清晰图像；检测每幅图像中两个椭圆影像的边界，利用边界像素点拟合两个椭圆方程；根据每幅图中所拟合出的两个椭圆方程，求取该图像中圆环点像的坐标；根据三组圆环点像的坐标，求取摄像机的内部参数；根据每幅图中的两个椭圆方程，求取这两个椭圆的内公切线和外公切线；利用拉盖尔定理，恢复世界坐标系中两相同圆的内公切线间的夹角；利用射影变换的交比不变性，修正所求的两相同圆的内公切线间的夹角；根据所求夹角建立世界坐标系，获取两个相同圆的内公切线和外公切线的世界坐标；根据两个椭圆的内外公切线和世界坐标系下两个相同圆的内外公切线求取单映射阵；根据己经求取的单映射阵和内部参数矩阵，求取对应于该幅图的摄像机的外部参数，本发明采用如下技术方案步骤l:以两个不相切的相同圆作为标定物；步骤2:用摄像机从三个不同角度拍摄标定物，分别得到三幅具有两个椭圆影像的标定物图像，确保成像清晰；步骤3:对每幅标定物图像进行以下操作步骤3.1:利用经典的边缘检测算子Canny算子，分别检测出图像上每一个椭圆影像的边界，从而得到两组由椭圆影像边界像素点构成的点集，分别利用这两组点集，拟合出两个椭圆影像的方程，两个椭圆影像中左椭圆标记为C;，右椭圆标记为《，左椭圆《、右椭圆《的方程分别表示为C;:m2+4^wv+_Selv2++￡>elv+=0C;:++"Be2v2+Cc2m+_Dc2v+￡c2=0其中，w、v分别为图像中像素点的横纵坐标，4,、&pC￡l、_Dt)1、A,为左椭圆q'方程的系数，<2、A2、Ce2、^2、《2为右椭圆《方程的系数；同时，根据两椭圆方程系数，可分别求得左椭圆C;和右椭圆《的拟合圆心0:、《的坐标；步骤3.2:利用两个椭圆影像的方程，在复域内求解两个椭圆的交点坐标，然后提取圆环点像m,，m,的坐标，其中提取圆环点像m,，m,的方法为任意选取两个椭圆的交点中的一对共轭解，如果这对共轭解所组成的直线与两个椭圆拟合圆心《、《所组成的直线的交点在0:和点《之间，表明这对共轭解为非圆环点对；反之，il果这对共轭解所组成的直线与两个椭il拟合圆心O:、《所组成的直线的交点不在点《和点《之间，表明这对共轭解为圆环点对，这对共'轭解2为圆环点像m,，m,的坐标；'2步骤4:根据三幅图像的圆环点像坐标，求取摄像机的内部参数矩阵《，计算公式为/T/r'^T^i^0;其中，《是圆环点像的齐次坐标；步骤5:分别求取三副图像的摄像机的外部参数，具体操作如下-步骤5.1:根据图像中两个椭圆影像的方程，求取两个椭圆影像的内公切线和外公切线，分另廿表示为—、/z'we—/w;、//"e—cw/、//"e—<3"/2;步骤5.2:利用拉盖尔定理，恢复世界坐标系一两相同圆的内公切线间的初始夹角e，计算方法为首先求得两个椭圆影像的内公切线的交点坐标，交点标记为O'，然后由点O分另!j与圆环点f^，/w,组成直线0'附〃0附〃根据直线""e—、""e—z77;、、"求出世界坐标系中两相同圆的内公切线间的初始夹角P:6=丄ln(crcm(//we—/巧,/z'"e—i,0附,，O附j))2/一一其中，cra^()表示仿射变换中的交比，i表示虚部单位；步骤5.3:定义两内公切线的不夹两圆的夹角为^—/rwe，利用射影变换的交比不变性，根据步骤5.2所求的两相同圆的内公切线间的初始夹角6求取夹角0—frwe，具体方法如下分别求取两个椭圆影像的两条外公切线/z'加—、与椭圆《的交点尸3'、g'，分别求取两个椭圆影像的两条外公切线"'"e_0<、与椭圆C;的交点尸4、尸6，求取两个椭圆影像的两条外公切线/z'朋—0<、—0^2的交点/;以及两个椭圆影像的两条内公切线//朋—/";、"股—&2的交点0';由点/"和点0'组成直线0'i;，由点尸3和点尸5'组成直线尸3>5，由点A和点尸e组成直线A《；分别求取直线0'尸^分别与直线尸;g、gg的交点d;、《；由点5;和点《组成直线5;《，分别求取直线d、'《与椭圆7C;、C;的交点《'、尸2';如果《、尸2'、S'、《四点的交比值大于—0.4268，夹角9_/r"e大于冗/2，此时夹角^—^"e为;z"减去歩骤5.2所求的角度《；《'、尸2、5;、《四点的交比值不大于-0.4268，夹角^一fr狄不大于;r/2，此时夹角0jrae等于歩骤5.2所求的步骤5.4:根据歩骤5.3所求夹角^—fr"e，建立如下世界坐标系将世界坐标系的原点定为两个相同圆内公切线的交点，x轴定为两个相同圆圆心连线的方向，y轴方向垂直于x轴方向；根据夹角^一frwe和所建立的世界坐标系，两个相同圆的两外公切线分别表示为//"e—ow,-3，//we—oW2=7其中，为两圆的半径，将其作为世界坐标系的物理度量单位；两个相同圆的两内公切线表示为//we—=[-tan(Or—6>—&認)/2)10]r，//we一/"2=[tan(O—6>」rw)/2)10]r根据两相同圆内外公切线的齐次坐标//"e—;力,、//恥—!'"2、//"e_、//"e—oz#2和两椭圆影像的内外公切线的齐次坐标/z'"e—f";、//"e—/"2、//"e_o<、//"e—o^2，求得单映射矩阵/f为'=(//we—om^，//we—，//"e—/巧)—rx(//we一ow、,//"e—o<，/z.we—)r其中，<formula>formulaseeoriginaldocumentpage8</formula>其中，[ve"w]j表示向量vector的第y个分量；步骤5.5:根据已经求取的单映射阵/Z和内部参数矩阵K，求取摄像机的外部参数，具体方法如下求取矩阵B一/ew;，计算公式为令6,，z'-l，2,3为矩阵5—"，的列向量，摄像机的外部参数为<formula>formulaseeoriginaldocumentpage8</formula>其中，l卜ll表示向量的第二范数，det(，)表示矩阵的行列式。有益效果与现有技术相比，本发明具有如下优点1.本发明在不需要已知任何摄像机参数情况下，利用圆环点像和两相同圆的内外公切线，简单、方便地获得摄像机的内部参数和外部参数；2.本发明简单、灵活，不需要特制的标定物以及繁琐的图像坐标匹配过程，能够实现全自动标定，无人机交互。图1是采用两个相同圆进行摄像机标定的具体步骤流程图。图2是圆环点像的获取示意图。图3(a)是两个相同圆内外公切线示意图；图3(b)是图3(a)经过透视投影变换后的两个椭圆的内外公切线示意图。图4是两组相同圆示意图。图5是两椭圆内公切线的正确夹角获取示意图。图6是关键点获取示意图。具体实施例方式为了更好地理解本发明，下面结合附图对本发明的具体实施方式作进一步描述，具体步骤的流程图如图1所示，应用本方法进行摄像机标定的具体步骤如下步骤l:首先制定标定物。在一个平整的物体上任意绘制两个相同圆或任意放置两个相同圆形物体，注意两个相同圆不能相切。以这两个不相切的相同圆作为标定物。步骤2:将标定物平整地放置在摄像机的视域范围内，用摄像机从三个不同角度拍摄标定物，分别得到三幅具有两个椭圆影像的标定物图像，确保成像清晰。步骤3:对每幅标定物图像进行以下操作-步骤3丄利用经典的边缘检测算子Canny算子，分别检测出图像上每一个椭圆影像的边界，从而得到两组由椭圆影像边界像素点构成的点集，分别利用这两组点集，拟合出两个椭圆影像的方程，两个椭圆影像中左椭圆标记为C,'，右椭圆标记为C;，左椭圆《、右椭圆《的方程分别表示为fC;:m2+_4戸+AiV2+C>+Delv+￡el=0S.7，(1)LC2:w2+^e2wv+Se2v2+Ce2w+Z)e2v+￡e2=0其中，w、v分别为图像中像素点的横纵坐标，4,、A,、Cel、化,、A,为左椭圆C;方程的系数，<2、5e2、Ce2、De2、￡。2为右椭圆<^方程的系数；同时，根据两椭圆方程系数，可分别求得左椭圆《和右椭圆《的拟合圆心《、《的坐标。具体的拟合椭圆方法及拟合圆心的求法参见(FITZGIBBONAW,PILUMandFISHERRB,"DirectLeast-SquaresFittingofEllipses",压五五rra舰_P"e/-"她cA!'"e/"/e//，14(2):239-256)。步骤3.2:利用两个椭圆影像的方程，在复域内求解出两个椭圆的交点坐标，然后提取圆环点像m,，w,的坐标，其中提取圆环点像m,，m,的方法为根据文献(WUYH,LIX,WUFC，"a/."Coplanarcircles,quasi-affineinvarianceandcalibration",/mageflWComp^"g,2006,24(4):319—326)，在复域内两个椭圆的交点有两对共轭解，其中的一对共轭解是圆环点像w/、wj，并且由这对共轭解所组成的直线不在两个椭圆之间，而由另外一对共轭解所组成的直线在两个椭圆之间。因此，任意选取两个椭圆的交点中的一对共轭解，如果这对共轭解所组成的直线与两个椭圆拟合圆心O'、O'所组成的直线的交点在0:和点《之间，表明这对共轭解为非圆环点对；反之，如果这对共轭解所组成的直线^两个椭圆拟合圆心o:、《所组成的直线的交点不在点o:和点《之间，表明这对共轭解为圆环点对，这对共轭解为圆环点像m,，w,的坐标。如图2所示，4为圆环点像m/、mj所组成的直线，《为另外一对共轭解所组成的直线，《、《分别为两椭圆c;、《的拟合圆心，《为直线4和两个椭圆拟合圆心《、《所成直线的交点。步骤4:根据三幅图像的lll环点像坐标，求取摄像机的内部参数:具体方法为根据射影几何知识，圆环点的像在绝对二次曲线的投影上，即满足(2)其中，《为摄像机的内参阵，S是圆环点像的齐次坐标。由于内参阵《中包含5个未知参数，通过上述所求的三对圆环点，利用式(2)就可以求出内参阵f。步骤5:分别求取三副图像的摄像机的外部参数，具体操作如下步骤5丄根据图像中两个椭圆影像的方程，求取两个椭圆影像的两条内公切线和两条外公切线，分别表示为""e——、""e—/fwe—ow〈、//we—，如图3(b)所示；步骤5.2:利用拉盖尔定理，恢复世界坐标系中两相同圆的内公切线间的初始夹角0，计算方法为将世界坐标系的原点定于两圆的两内公切线的交点处，如图3(a)所示的0点。据射影几何知识，图3(b)中的/z'ne一—、Z/恥j":为两条非迷向直线的像，O'为两个椭圆影像的两条内公切线的交点，由点O'分别与li环点w,，组成的直线Om,、0'w^为两条迷向直线的像，其中叫，wj为圆环点的像。根据拉盖尔定理及透视投影变换的交比不变性，可知两内公切线/z'"ej",、//從_/"2的夹角为-^=丄ln(cro5\y(/iwe—〖"],〃we_z'"2,0m,,Om,))(3)2/一_'其中，CTOM()表示仿射变换中的交比，/表示虚部单位。因此，首先求得两个椭圆影像的内公切线的交点坐标O'，然后由点O'分别与圆环点w;，w^组成直线0'm,，Om"将直线/Z"e—/w,、/z'"e—/"2、、O'n^分别代入式(3)，就可求出世界坐标系中两相同圆的内公切线间的初始夹角^。步骤5.3:定义两内公切线的不夹两圆的夹角为0Jr欲，利用射影变换的交比不变性，根据步骤5.2所求的两相同圆的内公切线间的初始夹角e求取夹角0jn/e，如图5所示。由于步骤5.2中所求夹角为锐角或直角，对于如图4所示的两组相同圆(C,，CJ和{C,Crf}，这两组圆的内公切线相同，内公切线夹角相同，但是两组相同圆的几何形态不同，因此，为了区分这两种情况，确定夹角^—"M的方法如下如图5所示，在世界坐标系下，Q、^分别为圆C,、圆q的圆心，S为线段O,Q与圆C,的交点，A为线段0,(^与圆q的交点。不妨设两圆的半径为1，点(9与点P2的距离为J^，显然，随着tfc的增大，两内公切线间的夹角0—&從也逐渐增大，当P—rrae:r/2时，=—1。因此，当d&〉V^—l时，0—&恥>兀/2。此夕卜，根据.0,尸2^92四点的交比，有0謹(。2)=——^——(4)由式(4)可知，当cfo〉0时，CTOM(0,iy^92)差关于^i的递增函数。综上，可以得出结论tr,((92)>—0.42686>—她>"/2;ete汐—加e2;z72(5)对两个椭圆影像，分别求取两个椭圆影像的两条外公切线"恥_0<、与椭圆C;的交点g、户5，分别求取两个椭圆影像的两条外公切线/z'"e—OM/:、"恥—oW;与椭圆《的交点g、尸6，求取两个椭圆影像的两条外公切线/Z"e—o"〈、"恥—cn^的交点^p以及两个椭圆影像的两条内公切线Z/"e—"、—!、的交点0';由点^,p和点"组成直线由点g和点g组成直线gg，由点尸4和点《组成直线尸4尸6;分别求取直线"i^分别与直线尸3>5'、尸4'尸6的交点5,'、52，由点a和点《组成直线a《，分别求取线5;52与椭圆C;、《的交点《、P2'，如图6所示。根据透视投影不变性可知，如图6中所示的点if、尸2、《、52，分别是图5中所示的点《、尸2、Q、02经过透视投影变换后的像，根据射影变换中的交比不变性和式(5)，有.-z/cmw(W尸，;)>-0.42686」麼>"/2;,,w"W2(6)将点A、g、《的坐标代入式(6)，得到夹角0—/n^与;r/2的大小关系；然后做出判断，如果P—frwe>;r/2,夹角0—为;r减去步骤5.2所求的角度0;如果S兀/2时，夹角6—fr狄等于步骤5.2所求的角度P。步骤5.4:根据步骤5.3所求夹角0」n^，建立如下世界坐标系将世界坐标系的原点定为两个相同圆内公切线的交点，x轴定为两个相同圆圆心连线的方向，；;轴方向垂直于x轴方向，如图3(a)所示。根据夹角^—fr"e和所建立的世界坐标系，两个相同圆的两外公切线分别表示为-//"e—owf，=r，//we—oW2=1，/f"e—=[tan(O-^—f，)/2)10]r而透视投影变换的直线变换公式为/■='/(7)其中，/f为单映射阵，/为实际直线的齐次坐标，/'为其对应的图像直线的齐次坐标。将两相同圆内外公切线的齐次坐标""e—!"、""e—/"2、/f"e—o^、//加_0^2和两椭圆影像的内外公切线的齐次坐标/z'"e—!Vj;、/!'"e—f"2、""e_o<、""e—o^2代入式(7)，求得单映射矩阵W为=(//we一，〃we一cwf2，//we一/打,)—rx(//we—o",//we—owf,//we—)7(8)其中，乂.其中，[ve"or]j表示向量vector的第个分量。步骤5.5:根据已经求取的单映射阵i/和内部参数矩阵〖，求取摄像机的外部参数，具体方法为:令图像平面的世界坐标Z-O，由针孔透视模型可知r2。(9)求取矩阵S—Zewp，计算公式为5—femp=尺—17/(10)令&,^'=1,2,3为矩阵尺—W的列向量，根据式(9)、式(10)及旋转矩阵的特性，可得摄像机的外部参数为-及=6,6，det(6，116,iiii62ii^ii-ii^ii6,其中，H.ll表示向量的第二范数，det(O表示矩阵的行列式.权利要求1、一种采用两个相同圆进行摄像机标定的方法，其特征在于该方法包括以下步骤步骤1以两个不相切的相同圆作为标定物；步骤2用摄像机从三个不同角度拍摄标定物，分别得到三幅具有两个椭圆影像的标定物图像，确保成像清晰；步骤3对每幅标定物图像进行以下操作步骤3.1利用经典的边缘检测算子-Canny算子，分别检测出图像上每一个椭圆影像的边界，从而得到两组由椭圆影像边界像素点构成的点集，分别利用这两组点集，拟合出两个椭圆影像的方程，两个椭圆影像中左椭圆标记为C′1，右椭圆标记为C′2，左椭圆C′1、右椭圆C′2的方程分别表示为<mathsid="math0001"num="0001"><math><![CDATA[<mfencedopen='{'close=''><mtable><mtr><mtd><msubsup><mi>C</mi><mn>1</mn><mo>′</mo></msubsup><mo>:</mo><msup><mi>u</mi><mn>2</mn></msup><mo>+</mo><msub><mi>A</mi><mrow><mi>e</mi><mn>1</mn></mrow></msub><mi>uv</mi><mo>+</mo><msub><mi>B</mi><mrow><mi>e</mi><mn>1</mn></mrow></msub><msup><mi>v</mi><mn>2</mn></msup><mo>+</mo><msub><mi>C</mi><mrow><mi>e</mi><mn>1</mn></mrow></msub><mi>u</mi><mo>+</mo><msub><mi>D</mi><mrow><mi>e</mi><mn>1</mn></mrow></msub><mi>v</mi><mo>+</mo><msub><mi>E</mi><mrow><mi>e</mi><mn>1</mn></mrow></msub><mo>=</mo><mn>0</mn></mtd></mtr><mtr><mtd><msubsup><mi>C</mi><mn>2</mn><mo>′</mo></msubsup><mo>:</mo><msup><mi>u</mi><mn>2</mn></msup><mo>+</mo><msub><mi>A</mi><mrow><mi>e</mi><mn>2</mn></mrow></msub><mi>uv</mi><mo>+</mo><msub><mi>B</mi><mrow><mi>e</mi><mn>2</mn></mrow></msub><msup><mi>v</mi><mn>2</mn></msup><mo>+</mo><msub><mi>C</mi><mrow><mi>e</mi><mn>2</mn></mrow></msub><mi>u</mi><mo>+</mo><msub><mi>D</mi><mrow><mi>e</mi><mn>2</mn></mrow></msub><mi>v</mi><mo>+</mo><msub><mi>E</mi><mrow><mi>e</mi><mn>2</mn></mrow></msub><mo>=</mo><mn>0</mn></mtd></mtr></mtable></mfenced>]]></math></maths>其中，u、v分别为图像中像素点的横纵坐标，Ae1、Be1、Ce1、De1、Ee1为左椭圆C′1方程的系数，Ae2、Be2、Ce2、De2、Ee2为右椭圆C′2方程的系数；同时，根据两椭圆方程系数，分别求得左椭圆C′1和右椭圆C′2的拟合圆心O′1、O′2的坐标；步骤3.2利用两个椭圆影像的方程，在复域内求解两个椭圆的交点坐标，然后提取圆环点像mI，mJ的坐标，其中提取圆环点像mI，mJ的方法为任意选取两个椭圆的交点中的一对共轭解，如果这对共轭解所组成的直线与两个椭圆拟合圆心O′1、O′2所组成的直线的交点在O′1和点O′2之间，表明这对共轭解为非圆环点对；反之，如果这对共轭解所组成的直线与两个椭圆拟合圆心O′1、O′2所组成的直线的交点不在点O′1和点O′2之间，表明这对共轭解为圆环点对，这对共轭解为圆环点像mI，mJ的坐标；步骤4根据三幅图像的圆环点像坐标，求取摄像机的内部参数矩阵K，计算公式为<mathsid="math0002"num="0002"><math><![CDATA[<mrow><msup><mover><mi>m</mi><mo>~</mo></mover><mi>T</mi></msup><msup><mi>K</mi><mrow><mo>-</mo><mi>T</mi></mrow></msup><msup><mi>K</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mover><mi>m</mi><mo>~</mo></mover><mo>=</mo><mn>0</mn><mo>;</mo></mrow>]]></math>id="icf0002"file="A2009101832220002C2.tif"wi="30"he="4"top="174"left="34"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/></maths>其中，id="icf0003"file="A2009101832220002C3.tif"wi="3"he="3"top="174"left="82"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>是圆环点像的齐次坐标；步骤5分别求取三副图像的摄像机的外部参数，具体操作如下步骤5.1根据图像中两个椭圆影像的方程，求取两个椭圆影像的内公切线和外公切线，分别表示为line_in′1、line_in′2、line_out′1、line_out′2；步骤5.2利用拉盖尔定理，恢复世界坐标系中两相同圆的内公切线间的初始夹角θ，计算方法为首先求得两个椭圆影像的内公切线的交点坐标，交点标记为O′，然后由点O′分别与圆环点mI，mJ组成直线O′mI，O′mJ，根据直线line_in′1、line_in′2、O′mI、O′mJ求出世界坐标系中两相同圆的内公切线间的初始夹角θ<mathsid="math0003"num="0003"><math><![CDATA[<mrow><mi>θ</mi><mo>=</mo><mfrac><mn>1</mn><mrow><mn>2</mn><mi>i</mi></mrow></mfrac><mi>ln</mi><mrow><mo>(</mo><mi>cross</mi><mrow><mo>(</mo><mi>line</mi><mo>_</mo><msubsup><mi>in</mi><mn>1</mn><mo>′</mo></msubsup><mo>,</mo><mi>line</mi><mo>_</mo><msubsup><mi>in</mi><mn>2</mn><mo>′</mo></msubsup><mo>,</mo><msup><mi>O</mi><mo>′</mo></msup><msub><mi>m</mi><mi>I</mi></msub><mo>,</mo><msup><mi>O</mi><mo>′</mo></msup><msub><mi>m</mi><mi>J</mi></msub><mo>)</mo></mrow><mo>)</mo></mrow></mrow>]]></math></maths>其中，cross()表示仿射变换中的交比，i表示虚部单位；步骤5.3定义两内公切线的不夹两圆的夹角为θ_true，利用射影变换的交比不变性，根据步骤5.2所求的两相同圆的内公切线间的初始夹角θ求取夹角θ_true，具体方法如下分别求取两个椭圆影像的两条外公切线line_out′1、line_out′2与椭圆C′1的交点P′3、P′5，分别求取两个椭圆影像的两条外公切线line_out′1、line_out′2与椭圆C′2的交点P′4、P′6，求取两个椭圆影像的两条外公切线line_out′1、line_out′2的交点Pvsp以及两个椭圆影像的两条内公切线line_in′1、line_in′2的交点O′；由点Pvsp和点O′组成直线O′Pvsp，由点P′3和点P′5组成直线P′3P′5，由点P′4和点P′6组成直线P′4P′6；分别求取直线O′Pvsp分别与直线P′3P′5、P′4P′6的交点id="icf0005"file="A2009101832220003C1.tif"wi="15"he="5"top="73"left="57"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>由点id="icf0006"file="A2009101832220003C2.tif"wi="3"he="5"top="73"left="85"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>和点id="icf0007"file="A2009101832220003C3.tif"wi="4"he="5"top="73"left="99"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>组成直线id="icf0008"file="A2009101832220003C4.tif"wi="11"he="5"top="73"left="122"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>分别求取直线id="icf0009"file="A2009101832220003C5.tif"wi="9"he="5"top="73"left="163"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>与椭圆C′1、C′2的交点P′1、P′2；如果P′1、P′2、id="icf0010"file="A2009101832220003C6.tif"wi="13"he="5"top="80"left="95"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>四点的交比值大于-0.4268，夹角θ_true大于π/2，此时夹角θ_true为π减去步骤5.2所求的角度θ；P′1、P′2、id="icf0011"file="A2009101832220003C7.tif"wi="13"he="5"top="87"left="159"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>四点的交比值不大于-0.4268，夹角θ_true不大于π/2，此时夹角θ_true等于步骤5.2所求的角度θ；步骤5.4根据步骤5.3所求夹角θ_true，建立如下世界坐标系将世界坐标系的原点定为两个相同圆内公切线的交点，x轴定为两个相同圆圆心连线的方向，y轴方向垂直于x轴方向；根据夹角θ_true和所建立的世界坐标系，两个相同圆的两外公切线分别表示为line_out1＝T，line_out2＝T其中，r为两圆的半径，将其作为世界坐标系的物理度量单位；两个相同圆的两内公切线表示为line_in1＝[-tan((π-θ_true)/2)10]T，line_in2＝[tan((π-θ_true)/2)10]T根据两相同圆内外公切线的齐次坐标line_in1、line_in2、line_out1、line_out2和两椭圆影像的内外公切线的齐次坐标line_in′1、line_in′2、line_out′1、line_out′2，求得单映射矩阵H为H＝(line_out1，line_out2，line_in1)-T×diag(e1，e2，e3)×(line_out′1，line_out′2，line_in′1)T其中，<mathsid="math0004"num="0004"><math><![CDATA[<mrow><msub><mi>e</mi><mi>j</mi></msub><mo>=</mo><mfrac><msub><mrow><mo>[</mo><msup><mrow><mo>(</mo><mi>line</mi><mo>_</mo><msubsup><mi>out</mi><mn>1</mn><mo>′</mo></msubsup><mo>,</mo><mi>line</mi><mo>_</mo><msubsup><mi>out</mi><mn>2</mn><mo>′</mo></msubsup><mo>,</mo><msubsup><mrow><mi>line</mi><mo>_</mo><mi>in</mi></mrow><mn>1</mn><mo>′</mo></msubsup><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><msubsup><mrow><mi>line</mi><mo>_</mo><mi>in</mi></mrow><mn>2</mn><mo>′</mo></msubsup><mo>]</mo></mrow><mi>j</mi></msub><msub><mrow><mo>[</mo><msup><mrow><mo>(</mo><mrow><mi>line</mi><mo>_</mo></mrow><msub><mi>out</mi><mn>1</mn></msub><mo>,</mo><mi>line</mi><mo>_</mo><msub><mi>out</mi><mn>2</mn></msub><mo>,</mo><msub><mrow><mi>line</mi><mo>_</mo><mi>in</mi></mrow><mn>1</mn></msub><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><msub><mrow><mi>line</mi><mo>_</mo><mi>in</mi></mrow><mn>2</mn></msub><mo>]</mo></mrow><mi>j</mi></msub></mfrac></mrow>]]></math></maths>其中，[vector]j表示向量vector的第j个分量；步骤5.5根据已经求取的单映射阵H和内部参数矩阵K，求取摄像机的外部参数，具体方法如下求取矩阵B_temp，计算公式为B_temp＝K-1H令bi，i＝1，2，3为矩阵B_temp的列向量，摄像机的外部参数为<mathsid="math0005"num="0005"><math><![CDATA[<mfencedopen='{'close=''><mtable><mtr><mtd><mi>R</mi><mo>=</mo><mrow><mo>(</mo><mo>-</mo><mfrac><msub><mi>b</mi><mn>1</mn></msub><mrow><mo>|</mo><mo>|</mo><msub><mi>b</mi><mn>1</mn></msub><mo>|</mo><mo>|</mo></mrow></mfrac><mo>,</mo><mo>-</mo><mfrac><msub><mi>b</mi><mn>2</mn></msub><mrow><mo>|</mo><mo>|</mo><msub><mi>b</mi><mn>2</mn></msub><mo>|</mo><mo>|</mo></mrow></mfrac><mo>,</mo><mo>-</mo><mfrac><mrow><mi>det</mi><mrow><mo>(</mo><mfrac><msub><mi>b</mi><mn>1</mn></msub><mrow><mo>|</mo><mo>|</mo><msub><mi>b</mi><mn>1</mn></msub><mo>|</mo><mo>|</mo></mrow></mfrac><mo>,</mo><mfrac><msub><mi>b</mi><mn>2</mn></msub><mrow><mo>|</mo><mo>|</mo><msub><mi>b</mi><mn>2</mn></msub><mo>|</mo><mo>|</mo></mrow></mfrac><mo>,</mo><mfrac><mrow><mo>(</mo><msub><mi>b</mi><mn>1</mn></msub><mo>×</mo><msub><mi>b</mi><mn>2</mn></msub><mo>)</mo></mrow><mrow><mo>|</mo><mo>|</mo><msub><mi>b</mi><mn>1</mn></msub><mo>|</mo><mo>|</mo><mo>·</mo><mo>|</mo><mo>|</mo><msub><mi>b</mi><mn>2</mn></msub><mo>|</mo><mo>|</mo></mrow></mfrac><mo>)</mo></mrow><mrow><mo>(</mo><msub><mi>b</mi><mn>1</mn></msub><mo>×</mo><msub><mi>b</mi><mn>2</mn></msub><mo>)</mo></mrow></mrow><mrow><mo>|</mo><mo>|</mo><msub><mi>b</mi><mn>1</mn></msub><mo>|</mo><mo>|</mo><mo>·</mo><mo>|</mo><mo>|</mo><msub><mi>b</mi><mn>2</mn></msub><mo>|</mo><mo>|</mo></mrow></mfrac><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mi>t</mi><mo>=</mo><mo>-</mo><mfrac><msub><mi>b</mi><mn>3</mn></msub><mrow><mo>|</mo><mo>|</mo><msub><mi>b</mi><mn>1</mn></msub><mo>|</mo><mo>|</mo></mrow></mfrac></mtd></mtr></mtable></mfenced>]]></math></maths>其中，‖·‖表示向量的第二范数，det(·)表示矩阵的行列式。全文摘要本发明涉及一种采用两个相同圆进行摄像机标定的方法，利用两个不相切的相同圆作为标定物，用摄像机从三个不同角度拍摄标定物得到三幅清晰的图像，确定出每幅图中的两个椭圆影像的方程，求出每幅图中的一对圆环点像的坐标。利用三幅图像的圆环点像的坐标和绝对二次曲线的投影方程来确定摄像机的内部参数。对每幅图像，求取两个椭圆影像的内外公切线；恢复并修正世界坐标系中两相同圆的内公切线间的夹角；根据所建立的两个相同圆的内外公切线的世界坐标和其投影的图像坐标，求取摄像机的外部参数。本发明能简便而精确地求取摄像机内外参数，而且能实现全自动标定。本发明能够广泛应用于非接触式工业检测和基于视觉的自主导航系统中。文档编号H04N17/00GK101621711SQ20091018322公开日2010年1月6日申请日期2009年7月23日优先权日2009年7月23日发明者虎张,达飞鹏申请人:东南大学