专利名称：一种共面圆圆心投影点的定位方法
技术领域：
本发明涉及图像处理中的定位技术，尤其涉及一种共面圆圆心投影点的定位方法。
背景技术：
由于圆投影在图像平面上的椭圆是由大量边缘点定义，而曲线拟合过程对图像点 的噪声有平滑作用，因此，圆特征具备较强的抗噪声能力，利用现有的拟合算法便可以精确 的获得其参数。圆心投影点的定位是视觉任务中图像处理过程的基本任务，在完成圆心投 影点定位的基础上，才可以进一步完成摄像机标定、相对位置姿态测量、图像的度量矫正以 及一些涉及到圆特征的3D重构等图像处理过程。目前，圆心投影点的定位多应用于摄像机 标定、摄像机视觉姿态测量、以及无人机的自主着陆等方面。 圆心投影点的定位过程中，圆在透视投影中其圆心的投影点存在透视畸变，且该 圆心投影点与所投影的椭圆的圆心点在通常条件下并不重合。因此，直接使用椭圆圆心点 定位圆心的真实投影点来进行标定或者姿态测量，不仅理论上会存在误差，而且在实际应 用中也会影响后续处理过程的精度。Heikkila(海奇拉)等人在1997年的Proceedings of the IEEE ComputerSociety
Conference on Computer Vision and Pattern Recognition(计算机视觉与模式视觉
IEEE会议)上发表的论文"A Four-st印Camera Calibration Procedurewith Implicit
Image Correction.(含有图像校正的4步法摄像机标定)"中提到，利用已知各圆尺寸以
及圆之间相对距离的阵列圆进行摄像机标定，具体是将椭圆中心作为特征点，采用非线性
优化方法估计摄像机内外参数，并且由摄像机外部参数在非线性优化过程中对所得到的特
征点的位置进行矫正以获得圆心的真实投影点的位置。这种处理方法不仅处理过程较为复
杂、效率不高，而且需要已知阵列圆的排列方式和相对位置关系才能实现。 此夕卜，YihongWu(吴毅红)等人在2006年4月发表于Image and
VisionComputing(图像与视觉计算)上的论文"Coplanar circles, quasi-affine
invarianceand calibration.(共面圆，准仿射不变量以及摄像机标定)"中提出了一种利
用共面双圆进行摄像机自标定的方法，通过比较双圆在图像平面上的交点的位置关系，来
确定真实的虚圆点位置。该方法在计算过程中需要确定两个构造的二次曲线之间的交点，
以便通过极点-极线约束来获得圆心的真实投影点位置。显然，这种定位方法需要较为复
杂的处理过程，才可以确定两个二次曲线的交点，从而增加了定位处理过程的复杂度，降低
了定位处理过程中算法的稳定性，而且无法获得双圆的结构参数。
发明内容
有鉴于此，本发明的主要目的在于提供一种共面圆圆心投影点的定位方法，能解
决现有技术中圆心真实投影点定位中处理过程复杂的问题。 为达到上述目的，本发明的技术方案是这样实现的
本发明提供了一种共面圆圆心投影点的定位方法，所述方法包括 A、获取摄像机拍摄处于同一平面的两个共面圆得到的测量图像，对该测量图像进
行边缘提取处理，获取测量图像中的椭圆边缘点，并根据所获取的椭圆边缘点，得到所述共
面圆投影在所述测量图像平面上的椭圆表达式； B、根据得到的椭圆表达式构造二次曲线，得到所述二次曲线的特征值及对应的特 征向量； C、根据所述特征向量得到所述共面圆的投影不变点，确定所述共面圆圆心连线的 投影直线，并得到该投影直线与所述椭圆的交点；并根据所述特征值确定所述共面圆的结 构参数； D、根据所述共面圆圆心连线的投影直线与所述椭圆的交点、以及所述共面圆的结 构参数，由交比不变原理，确定所述共面圆圆心的真实投影点的位置。 上述方案中，步骤A中所述对测量图像进行处理为通过Hessian矩阵法对所述测 量图像进行处理。 上述方案中，步骤B中在得到所述二次曲线的特征值及对应的特征向量后，还包 括B1、根据计算的所述特征向量，得到所述各特征向量之间的欧式距离；并根据所述欧式 距离中的最小值，对所述特征向量进行排序；B2、根据所述特征向量的排序结果，对所述特 征向量对应的特征值进行排序。 上述方案中，步骤C中所述根据特征向量得到投影直线与所述椭圆的交点具体包 括C11、根据排序后的所述特征向量，确定所述共面圆的投影不变点的齐次坐标；C12、根 据所述共面圆的投影不变点的齐次坐标，确定所述共面圆圆心连线的投影直线的表达式； C13、根据所述共面圆圆心连线的投影直线的表达式及所述椭圆的表达式，得到该投影直线 与所述椭圆的交点的坐标。
上述方案中，步骤C中所述根据所述特征值确定所述共面圆的结构参数包括
C21、根据排序后的所述特征值，确定所述共面圆中两圆之间的圆心距与各圆半径的比值；
C22、根据所述圆心距与各圆半径的比值关系，得到所述共面圆中两圆半径的比值。 本发明所提供的共面圆圆心投影点的定位方法，通过共面双圆在摄像机图像平面
上的椭圆来构造二次曲线，获取所构造二次曲线的特征值和特征向量，并根据特征值和特
征向量得到共面双圆的透视投影不变点和双圆的结构参数，再利用交比不变原理来确定图
像平面上的双圆圆心的真实投影点。本发明无需已知共面双圆的任何结构参数，便可以确
定双圆圆心的真实投影点，并且可得到双圆的结构参数，方法简单易行，有效降低圆心投影
点定位处理过程的复杂度，并提高处理效率，且具有较好的实用性、稳定性和可靠性。
图1为本发明共面圆圆心投影点的定位方法的方法流程图；
图2为本发明共面双圆投影的示意图； 图3为由摄像机拍摄获得的共面双圆和多边形的特征图像示意图； 图4为对图3进行图像处理后得到的椭圆边缘点和多边形边缘点的示意图； 图5为从图4中分离出的椭圆边缘点和多边形边缘点的示意图。
具体实施例方式
下面结合附图和具体实施例对本发明的技术方案进一步详细阐述。
本发明中共面圆圆心投影点的定位方法，参照图l，该方法主要包括以下步骤
步骤101 :获取摄像机拍摄共面圆得到的测量图像，对该测量图像进行边缘提取处理，获取测量图像中的椭圆边缘点； 这里，所述共面圆包括两个处于同一平面的圆；所述测量图像可以通过将被测物体置于单目视觉传感器的测量空间内，再由摄像机拍摄共面双圆来得到。 所述获取测量图像的椭圆边缘点可以通过采用Hessian矩阵法对测量图像进行边缘提取处理来实现。其中，Hessian矩阵法的具体步骤可以参照Steger C(斯特格)于1998年发表的论文Unbiased Extraction of Curvilinear Structuresfrom 2D and3D images (二维图像和三维图像中曲线结构的无偏提取)。另外，还可以通过简单常用的Ca皿y算子，Robert算子等图像处理方法来实现。 步骤102 :根据所获取的椭圆边缘点，得到所述共面圆投影在所述测量图像平面上的椭圆表达式； 这里，由获取的椭圆边缘点得到所述共面圆投影在所述测量图像平面上的椭圆表达式的方法为现有技术，在此不再赘述。 步骤103 :根据得到的椭圆表达式构造二次曲线，并得到所述二次曲线的特征值及对应的特征向量； 这里，在得到所述二次曲线的特征值及对应的特征向量后，还可以包括根据计算的所述特征向量，得到所述各特征向量之间的欧式距离；根据所述欧式距离中的最小值，对所述特征向量进行排序；根据所述特征向量的排序结果，对所述特征向量对应的特征值进行排序。 步骤104 :根据所述特征向量得到所述共面圆的投影不变点，确定所述共面圆圆心连线的投影直线，并得到该投影直线与所述椭圆的交点； 具体可以为根据排序后的所述特征向量，通过对其进行归一化处理来确定所述
共面圆的至少两个投影不变点的齐次坐标；根据所述共面圆的至少两个投影不变点的齐次
坐标，由两点确定一条直线的原理即可确定所述共面圆圆心连线的投影直线的表达式；根
据所述共面圆圆心连线的投影直线的表达式及所述椭圆的表达式，计算投影直线与椭圆的
交点，从而得到该投影直线与所述椭圆的交点的坐标。
步骤105 :根据所述特征值确定所述共面圆的结构参数； 具体可以为根据排序后的所述特征值，分别计算出与两共面圆的圆心距与比例因子的比值、及两圆的圆半径与该比例因子的比值，从而确定所述共面圆中两圆之间的圆心距与各圆半径的比值；根据所述圆心距与各圆半径的比值，计算圆心距与各圆的圆半径的比值之间的比值，即可得到所述共面圆中两圆半径的比值关系。
步骤106 :根据所得到的共面圆圆心连线的投影直线与所述椭圆的交点，以及所得到的共面圆的结构参数，由交比不变原理，确定所述共面圆圆心的真实投影点的位置。
下面以分离的共面双圆为例，对上述的定位方法进一步详细阐述。如图2所示，空间平面上的共面双圆为圆Ci巧、圆Ciiv它们的圆心分别为圆心c(^、圆心c(v两圆的圆心距dist，两圆半径分别为半径rac^、半径rac^，两圆圆心的连线cclCc2与两圆的交点分别为点V Q、 D工；在摄像机图像平面上，双圆经过透视投影为椭圆e"椭圆^，两圆圆心连线
CClCc2投影为直线px，直线px与椭圆^、椭圆e2的交点依次为点^、bpCpC^。 步骤1，获取摄像机拍摄空间共面双圆得到的测量图像，并通过Hessian矩阵法对
测量图像进行边缘提取处理，获取测量图像中的椭圆边缘点。 步骤2，根据获取的测量图像椭圆边缘点，拟合得到椭圆参数，并将椭圆用矩阵表达式表示为e2 ; 步骤3，根据椭圆的矩阵表达式ei、e2，构造二次曲线Ql和二次曲线q2，并分别得到二次曲线Ql、q2的特征值以及对应的特征向量； 具体为根据椭圆的矩阵表达式ei、 e2，由公式(1)构造二次曲线Ql :
<formula>formula see original document page 6</formula> 并由公式(1)得到的二次曲线Ql的表达式，计算得到二次曲线Ql的特征值A u(i
=1，2，3)以及对应的特征向量Vu(i = 1，2，3); 由公式(2)计算特征向量Vli之间的欧式距离:<formula>formula see original document page 6</formula>
v2i、v22、v23和特征值入<formula>formula see original document page 6</formula>选取其中最小的欧式距离cLn，欧式距离cLn所对应的特征向量记为Vll、 V『第三
个特征向量则记为^3，并将二次曲线A的三个特征向量^i(i = 1，2，3)排序为vn、^、v『再将二次曲线A的特征值= 1，2，3)按照该特征向量的排序结果对应排序为入n、
入12、入13，其中入u、 A12、 Au分别对应特征向量Vn、V『V『 同理，由以下公式(3)构造二次曲线q2 :
<formula>formula see original document page 6</formula> 计算二次曲线q2的特征值A2i(i = 1，2，3)及对应的特征向并依照上述方法对其特征值和特征向量进行排序，记为特征向
入22、入23。 步骤4，根据排序后的特征向量，得到共面双圆的投影不变点PpP2、P3，确定该投影
不变点所在的直线px，以及直线px与椭圆ei、 e2的交点ai、 bp Cl、 ; 其中，投影不变点所在的直线px即为所述共面圆圆心连线的投影直线； 本步骤具体为二次曲线Ql中欧式距离cLn对应的特征向量v『^可视为摄像机
图像坐标系下两个点的齐次坐标，将该两点的齐次坐标归一化后即为共面圆的投影不变点
Pi、P2的齐次坐标，将投影不变点Pl、P2的齐次坐标分别记为<formula>formula see original document page 6</formula>;同样，再由第三个特征向量v^采用同样的方式得到共面圆的投影不变点P3的坐标，记
为<formula>formula see original document page 6</formula>
根据得到的投影不变点Pp^、P3的齐次坐标，计算其中任意两个不同点的叉积，便可以得到过投影不变点Pl、 P2、 P3的直线px的表达式，该直线px即为空间平面上共面圆圆心的连线cc^q在摄像机图像平面上的投影直线。这里，由公式(4)计算投影不变点P"P2的叉积来得到直线PX的表达式
<formula>formula see original document page 6</formula>
根据直线px的表达式，计算直线px与椭圆ei、 e2的四个交点，分别记为点a" VCpdp分别为空间平面上共面圆圆心的连线cclCC2与共面双圆Ciri、Cir2的四个交点A、BpQ、DJ勺投影点。
步骤5，根据排序后的特征值，得到圆半径的比值rir2和两圆圆心距dist与两圆的半径rad" rad2的比值、、k2 ; 具体可以这样实现根据下述公式(5)及得到的特征值A u、 A 12、 A 13，得到与共面双圆的圆心距dist相差一比例因子m的值dt :^ =#UHVI不3 I (5) 并根据公式(6)，得到与圆Cir2的半径rad2相差上述比例因子m的值r :
^ I (6) 再由得到的dt和r的值，由公式(7)得到圆心距dist与圆Cir2的半径rad2的比值k2: k2 = dist/rad2 = dt乂巧 (7) 同理，由二次曲线化的特征值A21、 A22、入23并结合上述公式(5)-(7)，可以得到圆心距dist与圆Ci巧的半径rac^的比值、; 再由、、k2的比值根据公式(8)便可以确定两圆Ciri、Cir2半径rac^、ra4的比值r^,艮卩: = 1^/、 (8) 步骤6，根据上述步骤4中得到投影直线px与椭圆ei、 e2的交点&1、 1^、 Cl、 和
步骤5中得到的两圆的结构参数，圆心距与半径的比值、、k2和两圆半径rac^、ra4的比值
巧iv由交比不变原理建立如下恒等式(9)， <j 、" 在式(9)中，crossratio表示求取4个点的交比 crossratio (x" x2， x3， x4) = (| | | | / | | x2x31 |) / (| | x!x41 | / | | x2x41 |) (10)
在式(10)中，I卜I I表示两点或者两个坐标值之间的距离。 其中，平面上共面圆圆心的连线cclCc2与共面双圆Ci巧、Cir2的交点A" Q、 D工可以根据两圆的圆心距dist与任一圆的半径的比值、或1^2及两圆半径的比值1^2来确定。 例如，可以根据得到的双圆结构参数，该结构参数包括圆心距与两圆半径的比值、、k2及两圆半径的比值rir2，以圆Ci巧的圆心CCl为原点，圆心cCl、Cc2的连线CclCc2为X轴，建立平面0-XY坐标系，并且令圆Ci巧的半径radi为1，则圆Cir2的半径rad2为两圆半径的比值巧iv两圆的圆心距dist为圆心距与圆Ci巧的半径radi的比值、，圆心cCl和圆心cq的坐标为(0，0)t、(kpO)t，而圆心连线cc^C2与两圆的交点Ai、BpQ、Di的坐标分别为(-1 ， 0)t、 (1 ， 0)t、 (VracVrad" 0)T、 (VracVrad" 0)T。其中，圆心距dist与圆Ci巧的半径racU勺比值、和两圆半径的比值i^在步骤5中已计算获得。因此，根据上述交点VB"Q、的坐标和圆心cCl、cc2的坐标，由式(9)和式(10)，可以得到交比= crossratio (Apcc丄，B， C工)禾口 CR2 = crossratio (B丄，C丄，cc2， D工)，在交比crossratio (a" pc" b" c》中交点A、V^的坐标均已知，那么可以通过下式(11)获得圆心cq在图像平面上的真实投影点的坐标值
7
<formula>formula see original document page 8</formula>(11)
<formula>formula see original document page 8</formula>
同理，可以计算得到圆心cc2在图像平面上的真实投影点pc2的坐标值。最后，由圆心的真实投影点的坐标值，便可以确定其位置了。 需要指出的是，共面双圆在空间中的分布按照其圆心距的不同可分为相离、相交、内嵌以及同心圆四种情况。其中，共面双圆为相交、内嵌的情况与上述实施方式基本相同，所不同的是步骤6中建立如式(9)所示的交比不变恒等式时，圆心的真实投影点与交点ai、h、 Cl、 的排列方式不同，这对于本领域人员来说，根据本领域公知常识即可确定，因此本发明不再详细赘述。 另夕卜，当共面双圆为同心圆时，步骤3-4中由计算得到的特征向量Vi(i = 1，2，3)所确定的图像平面上的三个投影不变点Pl、 P2、 P3，其中一个投影不变点必为两圆圆心的真实投影点。由于其余两个投影不变点的位置不会在同心圆投影在图像平面上的椭圆之内，而两圆圆心的真实投影点必然位于该椭圆之内。因此，也可通过判断三个投影不变点中哪个位于椭圆之内，来确定圆心的真实投影点的位置。 本发明通过共面双圆在摄像机图像平面上的像构造图像平面上的二次曲线，并计算该二次曲线对应的特征值及对应的特征向量，由特征向量确定投影不变点并得到双圆圆心的投影直线与椭圆的交点，再由特征值得到包含有空间共面双圆的半径比和圆心距与半径的比值的双圆的结构参数，最后由投影直线与椭圆的交点和双圆的结构参数利用交比不变原理来确定圆心的无畸变投影坐标即圆心的真实投影点。本发明的计算过程简单易行，且具有较好的实用性和可靠性，无需其它附加条件即可实现。 下面再以具体的实施例进一步详细说明本发明的效果，为了验证本发明的准确
性，本实施例将本发明的计算结果应用于图像的度量矫正中。 首先，在圆心投影点定位阶段 采用一个Mintron-468P的CCD摄像机，8mm精工镜头组成一个单目视觉传感器，图像分辨率为768X576pixe12 ;通过摄像机拍摄获得共面圆及多边形的特征图像，如图3所示。图4为经图像处理后得到的图3所示图像的边缘点图，图5为分离出的椭圆边缘点和多边形边缘点的示意图；提取图5中由图3的图像投影的椭圆和多边形的特征点图像坐标，并存储到计算机中，其中多边形的编号为1 4，多边形各边的编号次序为以某一边为起始边，在图像中逆时针排列编号；以计算机屏幕的左上角为原点建立靶标坐标系，水平向右，竖直向下分别为靶标坐标系的x、y轴正向。通过椭圆拟合获得投影的椭圆ei、e2的矩阵表达式分别为
—_1.685e313.193 6,525e5 _
q = 13.193-1.04k3 2.083e56.525e52.083e5 — 2.921e8
e，=
一1.618e3 —66.110_ 66.110 -774.2089.998e5 2.806e5
9.998e52.806e5-6.890e8 构造二次曲线仏=61^ 仏=
33.242 8.730 -2,262e417.107 5.327 -1.190e40.083 0.022 —56.647 _
计算二次曲线&的特征值分别为入u = -18. 783入12 = -0. 028入13
将得到的特征向量排序后为
0. 732
—0.885——0.889—_ 0.325 _ Vn =0.4650.456、3 =—0.946
0.0020.001_1.016e-4 对特征值排序后得到的序列为 入u = -18. 783入12 = -0. 028入13 = 0. 732 根据排序后的特征向量vn、 v12、 v13，归一化后得到的投影不变点的齐次坐标为
—399.557——599.1928——3196.358 — A =209.946307.4958A =-9311.672
111 两圆圆心连线在图像上的投影直线px的表达式可根据Pl、p2计算得到
px = (-97. 549， 199. 635，-2. 936e3)T贝U，由直线px的表达式与椭圆ei、e2的矩阵表达式，可以得到两者的交点按照图像坐标X的大小从小到大排列依次为 ^ -
—341.588—、A =—436.229——569.915—、、《=—640.622—
181.621227.866293.190327.740_
根据排序后的特征值序列，可以得到圆心距与圆的半径的比值为
k2 = 5. 303
同理，构造二次曲线^ =^^，并得到圆心距与另一圆的半径的比值
、=5. 261
根据得到的圆心距与圆的半径的比值、、k2得到，两圆半径的比值I^:
= 1. 008
根据计算得到双圆结构参数^A、巧iv可以按照下式计算得到空间平面上VBp
Q、DpCCpCC2之间的交比 其中，令半径rac^ = 1，则根据、、k2、 可求得rad2 = 1/1^2、 dist =、，结合
9
上式，可知 crossratio (a丄，pc丄，b丄，c》=1.6204， 式中^、bp(^均已知，通过式(11)便可获得一圆的圆心在图像平面上的真实投影 点的坐标为pCl = (390. 278,205. 412)T
同理，得到 crossratio (bp pc2， d》=1.6223 式中V^、c^均已知，通过式(11)可获得另一圆圆心在图像平面上的真实投影点 pc2的坐标为pc2 = (606. 149,310. 895) t
其次，在度量矫正阶段 下面通过对测量图像的度量矫正来验证圆心点计算的准确性，一旦圆心的真实投 影点在图像平面上确定下来，可根据极点-极线关系确定图像上的无穷远直线的像(IOL， Image of Line at Infinite) ， I0L与椭圆e!、 e2的交点则为图像平面上虚圆点的像(ICP， Imaged Circular Points)，确定ICP之后则可在图像平面上计算得到与虚圆点对偶的二次 曲线(CDCP :The conic dual to the circularpoints)，然后便可确定空间平面上与圆平 面共面的直线间的夹角。 由圆心投影点定位阶段得到的圆心的真实投影点的坐标，依据上述方法经计算得 到如下结果 多边形1的邻边夹角为60. 162° 、59. 924° 、59. 914° ; 多边形2的邻边夹角为89. 87° 、89.663° 、89.768° 、89.975° ; 多边形3的邻边夹角为89. 251° 、89.201° 、88.673° 、88.722° ; 多边形4的邻边夹角为60. 236° 、60.643 ° 、59.172 ° 、60.377 ° 、60.414° 、
59.158° 。 在图3中所示的真实场景中，多边形1 4分别为正三角形、长方形、菱形、正六边 形，从上述计算结果可以看出，采用本发明定位方法得到的投影点坐标与真实值非常接近， 如此，足以说明根据圆心投影点可以得到较高精度的图像处理结果，也就是说圆心投影点 的定位结果具有较高的准确性。 以上所述，仅为本发明的较佳实施例而已，并非用于限定本发明的保护范围，凡在 本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护 范围之内。
权利要求
一种共面圆圆心投影点的定位方法，其特征在于，所述方法包括A、获取摄像机拍摄处于同一平面的两个共面圆得到的测量图像，对该测量图像进行边缘提取处理，获取测量图像中的椭圆边缘点，并根据所获取的椭圆边缘点，得到所述共面圆投影在所述测量图像平面上的椭圆表达式；B、根据得到的椭圆表达式构造二次曲线，得到所述二次曲线的特征值及对应的特征向量；C、根据所述特征向量得到所述共面圆的投影不变点，确定所述共面圆圆心连线的投影直线，并得到该投影直线与所述椭圆的交点；并根据所述特征值确定所述共面圆的结构参数；D、根据所述共面圆圆心连线的投影直线与所述椭圆的交点、以及所述共面圆的结构参数，由交比不变原理，确定所述共面圆圆心的真实投影点的位置。
2. 根据权利要求1所述共面圆圆心投影点的定位方法，其特征在于，步骤A中所述对测量图像进行处理为通过Hessian矩阵法对所述测量图像进行处理。
3. 根据权利要求1所述共面圆圆心投影点的定位方法，其特征在于，步骤B中在得到所述二次曲线的特征值及对应的特征向量后，还包括Bl、根据计算的所述特征向量，得到所述各特征向量之间的欧式距离；并根据所述欧式距离中的最小值，对所述特征向量进行排序；B2、根据所述特征向量的排序结果，对所述特征向量对应的特征值进行排序。
4. 根据权利要求3所述共面圆圆心投影点的定位方法，其特征在于，步骤C中所述根据特征向量得到投影直线与所述椭圆的交点具体包括Cl 1 、根据排序后的所述特征向量，确定所述共面圆的投影不变点的齐次坐标；C12、根据所述共面圆的投影不变点的齐次坐标，确定所述共面圆圆心连线的投影直线的表达式；C13、根据所述共面圆圆心连线的投影直线的表达式及所述椭圆的表达式，得到该投影直线与所述椭圆的交点的坐标。
5. 根据权利要求3所述共面圆圆心投影点的定位方法，其特征在于，步骤C中所述根据所述特征值确定所述共面圆的结构参数包括C21、根据排序后的所述特征值，确定所述共面圆中两圆之间的圆心距与各圆半径的比值；C22、根据所述圆心距与各圆半径的比值关系，得到所述共面圆中两圆半径的比值。
全文摘要
本发明公开了一种共面圆圆心投影点的定位方法，通过共面双圆在摄像机图像平面上的椭圆来构造二次曲线，获取该二次曲线的特征值和特征向量来得到共面双圆的透视投影不变点和双圆的结构参数，再利用交比不变原理来确定图像平面上的双圆圆心的真实投影点。本发明无需已知共面双圆的任何结构参数，便可以确定双圆圆心的真实投影点，并且可得到双圆的结构参数，方法简单易行，且具有较好的实用性和可靠性。
文档编号G06T7/00GK101706960SQ200910236998
公开日2010年5月12日 申请日期2009年10月30日 优先权日2009年10月30日
发明者张广军, 王巍, 赵征, 魏振忠 申请人:北京航空航天大学