2008年a题数学建模论文

2008 高教社杯全国大学生数学建模竞赛承 诺 书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。我们参赛选择的题号是（从 A/B/C/D 中选择一项填写）： A 我们的参赛报名号为（如果赛区设置报名号的话）： 所属学校（请填写完整的全名）： 吉林师范大学 参赛队员 (打印并签名) ：1. 宋关华 2. 李婷婷 3. 卢时伟 指导教师或指导教师组负责人 (打印并签名)： 滕飞 日期： 2008 年 9 月 22 日赛区评阅编号（由赛区组委会评阅前进行编号）：22008 高教社杯全国大学生数学建模竞赛编 号 专 用 页赛区评阅编号（由赛区组委会评阅前进行编号）：赛区评阅记录（可供赛区评阅时使用）：评阅人评分备注全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：1数码相机定位摘要数码相机在计算机视觉中的应用逐渐普及和深入，在应用中数码相机的标定相当重要.本文根据数码相机的成像原理，由摄影几何学中的相关知识，提出了一种标定方法.将数码相机的成像过程转化为图像在四个坐标系的变换，实现了三个变换步骤，分别将世界坐标系中的信息转换到摄像机坐标系，再由摄像机坐标系转换到图像坐标系，最后由图像坐标系转换到像素坐标系.本文通过一台相机拍摄的二维图像，实现空间坐标系与像平面坐标的透视变换，设计了一种方法检验了建立的模型算法，模型具有较好的稳定性.同时完成了对双目标定的重构，给出了两部固定相机相对位置的数学模型和方法.1、对问题中给定的标定图像及由相机拍摄得到的图像，引入四个坐标系即世界坐标系、摄像机坐标系、图像坐标系、像素坐标系，由坐标系的转换关系得到确定像平面坐标的数学模型.通过对内部参数矩阵及外部参数矩阵的确定实现了以圆为标定图像时，圆心坐标的转换，得到了圆心在像坐标系的坐标，建立确定靶标上圆的圆心在像平面坐标的数学模型.2、由给定的像平面建立坐标系，确定各圆心的坐标，根据（1）的模型利用Matlab 编 程 解得投影矩阵 （即内外参数） ，再根据坐标于内外参数 的关系得到MM靶标上圆的圆心在像平面的坐标。算的结果如下表:表（一）世界坐标(）毫米YX,（12,112） （42,112） （112,112） （12,12） (112,12)像素坐标（）像素vu,（326.5290,191.3318）（419.8290,195.1718）（637.5290,204.1318）（279.9890,496.8318）(590.9890,509.9318)图像坐标（ yx）毫米（-49.0664，50.9704）（-24.3839，49.9546）（33.2087，47.5842）（-61.3786，-29.8497）（20.8966，-33.3153）3、根据射影几何中直线与点结合性经过摄影变换的保持不变的性质可知两个共面的圆具有公共切线的性质经过透视投影之后仍成立，利用几何板测出像平面上靶标圆心坐标，即可检验模型的稳定性.4、进一步研究两部相机成像从而来确定实物的位置，即双目定位，给出一种简单的相机空间三维重建，简单地可以让两部相机公面且光轴互相平行，由中心摄影比例关系可知两部相机相对位置关系.最后，将简化模型推广建立了一般模型，实现了摄像机标定与三维重构.本文理论上提出了一种对相机直接标定的数学模型于算法，此算法简化了计算机过程，避免了对多余参数的计算.关键词 双目定位 三维重建 射影几何 摄像机标定2一、问题重述与分析数码相机在现实生活中有广泛的应用，如通过数码相机可以确定物体表面某些点的位置的特征，来给物体定位.常用的定位方法是双目定位,它的主要原理是对物体上一个特征点用两部固定在不同位置的相机摄得物体的像,分别获得该点在两部相机像平面上的坐标.只要知道两部相机精确的相对位置,就可用几何的方法得到该特征点在固定一部相机的坐标系中的坐标,即确定了特征点的位置.而对于双目定位的关键在于两部相机的位置,这一过程为系统标定.标定的过程由于无法直接得到没有几何尺寸的“点” ，利用平面画圆的标靶方法来把对应圆与像的圆心找到，标定实现.（1） 、建立数学模型和算法确定靶标上的圆心在该相机像平面的像坐标,以该相机的光心为原点建立坐标系,我们计算各个坐标系之间的转化,从而得到此坐标.（2） 、根据所给的图把靶标上的圆丰富的圆心在像平面上的图像坐标,通过坐标之间的转化求出相机的内外参数,利用计算机编程求得圆心的像坐标.（3） 、根据摄影几何中直线和点的结合性经过摄影变换仍然保持不变的性质,利用切线和切点来求出像点坐标进行检验.（4） 、确定两部固定相机的相对位置.二、条件假设1 、假设光轴与相机的成像平面垂直,两台相机的 轴重合, 轴相互平行.将其中一部相xy机沿着其 轴方向平移一段距离 （称为基线）后与令一部相机重合.xb2 、在问题四中假设两台相机的焦距及其他内部参数均相等.三、符号说明CCD-英文 Charge Coupled Device 即电荷耦合器件的缩写，它是一种特殊半导体器件（ ）-世界坐标系，其中 ， ， 分别为坐标系的坐标轴 wZYX, wXYwZ-摄像机坐标系，其中 ， ， 分别为坐标系的坐标轴c cc-图像坐标系，其中 ， 分别为坐标轴yx, xy-像素坐标系，其中 ， 分别为坐标轴vu uv-摄像机光心，即摄像机坐标系的原点O-坐标原点在 CCD 图像平面的中心1oU-像素坐标系坐标原点在 CCD 图像平面左上角 U 轴V -像素坐标系坐标原点在 CCD 图像平面左上角 V 轴-储存图像的 行 列数组NMMN- 在 ， 坐标系中的坐标0,vu1Ouv3-每一个像素在 轴方向上的物理尺寸dx x-每一个像素在 轴方向上的物理尺寸 y y-摄相机焦距1O-旋转矩阵R-三维平移向量t-假设空间中的某一点 P-任何点 在图像的投影位置为pPp-光心 与 点的连线OO- 空间某点 的图像点 的位置Tvu,-为空间第 个点的坐标TwiiZYX1,i-第 点的图像坐标iv,i-投影矩阵 的第 行 列元素ijmMij-为 矩阵K12n-为 维向量U- 矩阵的第 行的前三个元素组成的行向量3,21imT i， -分别为平移向量 的三个分量xtzy t-一条外公切线L、 -圆1C2-椭圆，圆 、 的像E、 1C2- 的齐次坐标mTwvu,- 的齐次坐标XWZYX- 表示这条外公切线与圆 、圆 的切点21,p21,P1C2-非零常数因子4-一个非奇异的 矩阵H3-空间平面上这两点之间的距离d四、模型的建立与求解问题一：为了定量描述光学成像过程,在线性模型摄像机标定中我们不妨引入四个坐标系：世界坐标系、摄像机坐标系、图像坐标系、像素坐标系.（1）世界坐标系-根据自然环境所选用的坐标系，坐标值用（ ）表wZYX,示.（2）摄像机坐标系-以摄像机的光心 点为坐标原点, 轴、 轴分别平行于OccCCD 平面的垂直边,摄像机的光轴为 轴,坐标值用 表示.cZZYX,（3）图像坐标系-坐标原点在 CCD 图像平面的中心 ,X 轴、Y 轴分别为平行于1oCCD 平面两条垂直边,坐标值用 表示.yx,（4）像素坐标系-坐标原点在 CCD 图像平面左上角 U 轴、V 轴分别为平行于图像坐标系的 X 轴、Y 轴,坐标值用 表示,且为离散的整数值.vu在解决此问题时,我们将图像坐标系与像素坐标系建立在同一个面上.在计算机内 N 可以用 的数组来储存图像, 行 列图像中的每一个元素（成为像素）的数NMMN值是该图像点的亮度,在下图 1，定义直角坐标系 , 是以像素为单位的图像vu坐标系的坐标即： 分别是该像素在数组中的列数与行数,因此并没有用物理单位vu表示出该像素在图像中的位置,需要在建立以物理单位（例如毫米）表示的图像坐标系.该坐标系以图像内某点 为原点. 轴与 轴分别与 轴、 轴平行.如下图 1 所1Oxyuv示,xyv图 1 图像坐标系图像坐标系 表示以像素为单位的图像坐标系的坐标, 表示以毫米为单位vu yx,的图像坐标,这两个坐标采用的单位坐标不同但都是图像坐标系,在 ， 坐标系中,原点 是摄像机光轴与坐标系平面的交点,这一点一般位于图像中心处,由于摄像镜头的1O畸变也许会是这个点稍微与理想值有点偏差,若 在 , 坐标系中的坐标为 假1Ouv0,vu5设每一个像素在 , 轴方向上的物理尺寸为 , ,这样可以通过下面的式子来联系xydxy两个坐标系 0uvdyv借助齐次坐标与矩阵将上面的形式转化如下式： 10110yxvdyuxvu 推出它的逆变关系式： 1010vudyyxdxy 2从而方便了两个坐标系的量化转换，摄像机成像关系可如下图 2 表示. 图 2 相机坐标系于世界坐标系其中 点称为摄像机光心，在图 2 中光轴与图像平面的交点，即为图像坐标系的O原点，摄像机坐标系即由点 与 , , 轴组成的直角坐标系。 为摄相机焦距.OCXYZ1O在环境中还可以选择一个基准坐标系来描述摄像机的位置，并用它来描述环境中任何物体的位置，该坐标系称为世界坐标系.相机坐标系与世界坐标系之间的关系可以用旋转矩阵 与平移向量 来描述.因此，现在在这里我们假设空间中某一点 在世界坐标Rt P系与摄像机坐标系下的齐次坐标分别为 ,这样可以用下面,1,TWZYTCZYX1,这个矩阵来描述这两个坐标系之间的关系 111WWTC ZMZOtRZYX 36其中， 是 正交单位矩阵；t 为三维平移向量； ； 为 矩阵.R3 TO,1M4通过上面的介绍，非常清楚地知道了在研究相机成像经常用到的三个坐标系以及它们之间的转换关系、在这基础上来探讨一下线性模型摄像机模型的.其实，空间上的任何一点 在图像上的成像位置可以近似于真空模型，也就是任何P点 在图像的投影位置为 ，从图 2，可以知道位置 是光心 与 点的连线 与图PppP像平面的交点.可以得到如下的比例关系式：CZfYXx 4其中， 为点 点图像坐标； 为空间点 在摄影坐标系下的坐标.这样yx,pX,P利用其次坐标与矩阵可以表示透视摄影关系如下： 1001CC ZYXfYZ 5然后将式 与 代入 式中就可以得到世界坐标系表示的点 与其投影点 之间的345 Pp关系：WWTyxWTC MXZYXOtRvu ZYXOtRfvdyuxvuZ 2100110 10101 6其中， ； 为 矩阵，称为投影矩阵； 1完全由 决定，dfxfy,M43 0,vuyx只与摄像机内部结构有关，因此，把这些参数称为摄像机外部单数；确定0,vuyx一个摄像机的内外参数的过程，称之为摄像机标定.由上面的 式可以知道如果确定了摄像机的内外参数，就知道投影矩阵 ，在知6 M道 的情况下，只要得到任何空间点 ,假设点 的坐标为 ，就MPTWWZYX1,可以求出图像点 的坐标（位置） ，这是因为在已知 与 时，由式 可以pTvu 67给出三个方程，在这三个方程中消去 就可以求出 。反过来，在已知空间某点CZTvu,的图像点 的位置 以及摄像机内外参数， 也不能唯一确.其实，在式 中，PTvu, WX6是 不可逆矩阵，当已知 与 时，由式 给出的三个方程中消去 ,只可M43MTvu,6CZ得到关于 ， , 的两个线性方程，由这两个线性方程组成的方程组即为射线wXYZ的方程，也就是说，投影点位 的所有点均在该射线上，该空间点并不都是 点，OPP p因此，该空间点并不能唯一地确定下来。这样来提一下由 演变过来的规一化投影关6系： 100101 CC ZYXfyxZ把式 代入得到：1 1010 10100110Cyx CC ZYXvu ZYXfvdyuxvuZ定义 为规一化的图像坐标，则得到：Tnyx1,1011vuyxyxn7这样由以上的两个式子得到 10011CnC ZYXyxZ由此进一步得到：8CnCnZYyXx 8上面使用规一化坐标，其实相当于假设摄像机焦距等于 1 的情况推导出开的，将规一化坐标转换成实际的图像坐标。7问题二：首先给出如何由参照物图像来求投影矩阵图 .将上面的式 写成： M8 113433231 222 1413121 wiiici ZYXmmvuZ 9其中， 为空间第 个点的坐标； 为第 点的图像坐标； 为TwiiYX, iTivu,iijm投影矩阵 的第 行 列元素。由此得出如下三个方程：Mj 3433231 222 1413121 mZYmXZvuwiwiwici iiiic 10将 中的 消去后得到两个关于 的线性方程如下：10ci ij 3433231243221 11 vZvYXvZYmXuuiwiwiwiwiiwi i 1我们观察上面的式子了解到，如果知道了标定块上 个已知点的话，就可以通过它n们的空间坐标 以及它们的图像点坐标 ，得Tii,n,Tiu,n,到由 个关于 中组成元素 的线性方程，这些线性方程也可以写成矩阵的形式n2Mijm如下：9 Tnnwnwnnwnwnwnwnwwwwww mvumvumZYXZYXuZYX vvZuYXu 341341 3231243212 111 111111100 000 0 2由上面的式子 可以知道， 矩阵乘以不为零的常数并不影响 与9MTwiiZYX,的关系，因此，在式 中我们指定 ，从而得到关于 矩阵其他元素的Tivu, 12134M个线性方程，这些位置元素个数为 11 维向量 ，将 简写成：n2 m2UK13其中， 为式 左边 矩阵， 为未知的 11 维向量； 为 右边的 维向K1312n Un2量； ， 为已知向量，当 时，采用最小二乘法求出上述线性方程的解为: UmTT114向量 与 构成了所求解的 矩阵.m134M求出 矩阵后，可由式 表示的关系算出摄像机的全部内外参数。但要注意一M12点，求出 矩阵与式 所表示的矩阵 相差一个常数因子 ，这一点可以从解方34m程 时指定 并不影响投影关系，但在分解 矩阵时是必须考虑的。 投影矩1234m M阵与摄像机内外参数的关系可以写成 1010132132424 zTyxyxT trvu15其中， 为由式 求得的 矩阵的第 行的前三个元素组成的行向量；,im9Mi为 矩阵第 行第四列元素； 为旋转矩阵 的第 行； ，4i3,21Mi3,2irTRixt分别为平移向量 的三个分量。由式 得到zyt, t1510 zTyy zxXT trvvtutum3002313242114 16比较 与 可知， ，由于 是正交单位矩阵的第三行， 因此，1564m3r 3r可以从 推导出 。在通过如下式来求 3431 ;,03yxvur4tz013umtxx0243vtyy343mr 3234my31423010 mruuTTx3011urx3243020vrvTTy 302342vy3134mx综上所说，根据公式 我们利用 MATLAB 得到 的相应取值范)2( 2321312,mm围即：表（二） 1m1213m21223（1.5337，1.62290）（0.1597，0.2517）（117.9308，121.6136）（0.0510，0.0867）（1.5583，1.5217）（258.9572，266.9674）任意取得三组数据利用 MATLAB 编程得到三组不同的 ，把它们2321312,mm对应相加并求得平均数，再把得到的平均数当作新的 用于计算，再利用 MATLAB 编程（程序见附录（一） ）把摄影圆心的坐标得到，如下表：11表（三） 世界坐标(）毫米YX,（12,112） （42,112） （112,112） （12,12） (112,12)像素坐标（）像素vu,（326.5290,191.3318） （419.8290,195.1718） （637.5290,204.1318） （279.9890,496.8318） (590.9890,509.9318)图像坐标（ yx）毫米（-49.0664，50.9704）（-24.3839，49.9546）（33.2087，47.5842）（-61.3786，-29.8497）（20.8966，-33.3153）问题三：用基于摄影不变性的测量对我们的模型进行检验.在以下模型中，我们设计了一种新型的平面模板，并基于这一模板实现了相机的简易标定和空间点与点之间距离的测量.（一）平面模板平面模板是由大小相同且不相交的共面圆组成，如图 3 为两个圆的平面模板基本结构与图像，其中 表示一条外公切线， 表示这条外公切线与圆 、圆 的切点.L21,P1C2经过透视投影后，平面模板的图像如图 4 所示，直线 是外公切线 像，椭圆 是lL1,E圆 、圆 的像，直线 是椭圆 的公切线，并且切点 是 的图像点.1C2l21,E21,p21,P图 3 两个共面圆和它们的公切线图 4：图 3 的透视投影由摄影几何中的性质（线和点的结合性经过摄影变换仍然保持不变）可知，两个共面的圆具有公切线的性质在经过透视投影之后仍然是成立的.根据椭圆的公切线，很容易计算出特征点来，将形如图 3 中的模板扩展可以得到12一个方格和圆组成的模板，它的图像如图 4 所示，图中的黑点表示公切线与椭圆的切点，并且由这个模板可以计算出椭圆的中心来.图 5 应用公切线计算特征点方法 1、利用几何画板做圆的切线，得到切点，两对切点连线得到相对应的圆心，得出标靶圆心坐标 ， ，82.19,.31O128.95,8360.4O, ,0641536.0,7图 6 几何画板确定圆心坐标（二）平面测量根据摄影几何中直线和点的结合性保持不变的性质，本模型取出图 5 中的特征点，并由这些特征点计算出空间平面与图像平面间的单应矩阵，从而能够测量出空间平面上任意两点之间的距离.文中假设相像机模型为常见的针孔模型，则三维空间点和它们对应的图像点满足下面的关系(17)XpPXm4321其中, 和 是齐次坐标且分别为 和 , 为非零 TwvumTWZYX常数因子.对于任何一个空间平面上的三维点，模型假设其中一个坐标分量为零， （也就是说把模板放在世界坐标系的 平面上） ，即 ，于是YO0Z13(18) ZYXhZYXpWYXPwvuH323113210因此，空间平面上的三维点和它相对应的图像点之间的关系可以由一个非奇异的矩阵 来表示，这个矩阵称为单应矩阵，在相差一个非零的常数因3321pH子的意义下，这个矩阵具有 个自由度。根据式（18）可得，已知一对空间平面和图像8平面上的对应点，可以求得两个关于单应矩阵的 个元素为未知量的线性约束.所以，9若已知空间平面和它的图像平面之间的四对对应点（其中任意三点不共线） ，则在相差一个常数因子的意义下可以唯一确定单应矩阵 .H若已知空间平面和它的图像平面之间的单应矩阵 ，则由图像点可以求出空间平面上的三维点坐标，因而，可以测量出在空间平面上任意两点之间的距离.假如已知空间平面的图像平面上两点 的坐标，则由式（17）可得21,x（19）211XH在计算出空间平面上两点 的坐标之后，从,11TWYTWYX22而得到空间平面上这两点之间的距离2121XWd上述提出了一种基于摄影不变性的单幅图像测量方法.在该发法中，设计了一种新型模板，该方法与基于其他模板方法相比较，具有一定的精度和鲁棒性.问题四：以上我们研究了一部相机定位问题，在此我们进一步研究两部相机成像来确定实物的位置，即双目定位.（1）我们给出一种简单的相机配置下的空间三维重建.如图 7图 7 双相机配置假如 与 相机的内部构造均相同，而且两部相机均在同一平面，其中，它们1C2的光轴相互平行， 轴重合， 轴互相平行。只允许两相机有位移上的偏差，但没有角xy度上的偏差.当用此方法配置时，两相机坐标只相差 轴方向上的一个平移，此平移距x离记为 为 与 坐标系的 轴.若 平面与两个图像平面 及 的交线分别21O12x21OP1I214为 与 ，由于两个图像平面位于同一个平面，则 与 为该平面上的同一直线，1E2 1E2又由于图像平面平行于 轴，故 与 与 轴平行，若 与 分别为 与 及x1E2xpPO1I与 的交点.则 与 分别在 与 上，事实上， 与 就是巡航面所说的极PO2I1p2 12线.在进一步假设 坐标系为 ， 坐标系为 ，则在上述相机配置下，C1zyOx2C2zyOx若任何空间点 的坐标在 坐标系下为 ，在 坐标系下为 .由1Tz1, Tzybx11,中心摄影的比例关系可得 10xu101zyv102bxu102zyv其中 为相机内部参数， 分别为 与 的图像坐标.过上述yxvu,0 TTvu21, 1p2式子可以解出： 210ub210vyyx21ubzx其中， 为基线长度， 为视差.有上面的式子可知，有 与 的图像坐标b 1p2可求出空间点 的三维坐标 . TTvu21, P1,zyx(2)可以用立体视觉方法进行三维重建，在计算机视觉中是指两幅或多幅两维图像恢复物体三位几何形状的方法。由图 6 看出，对于空间物体表面任意一点 ，用 相P1C机观察，它在 相机的图像点位于 ，但我们仍无法由 得知 的确切三维位置.事1C1p1p实上，在 ( 为 相机的光心)连线上的任意点 的图像点都是 ，因此，由POP115点的位置，这样只知道空间点位于 连线上的某一位置，但是如果用 与 两1p 1pO1C2个相机同时观察 点，并且如果能确定在 相机图像上的点 与在 相机图像上的点PC1p2是空间同一点 的图像点，这样情况就不同了，因为知道了空间点 既位于 上，2 P1pO有位于 上，因此 点是 与 两条直线的交点，即它的三维位置是唯一确2pO1pO2定的.以上介绍了立体视觉的基本原理.在就要从实际计算的角度出发.在什么条件下，能得到 点的三维坐标.P假如物体表面上所有点的三维坐标已知，则三维物体的形状与位置就是唯一确定的.图 8 双摄像机看空间点最简单的方法是用立体视觉的方法获取三维点的坐标.假定，空间任意点 在两个P相机 与 上的图像点 与 为空间同一点的对应点，相机 与 已标定，它们的1C21p2 1C2投影矩阵分别为 与 ，于是有：1M 11343132131 222 1431211 ZYXmmvuZC 17 113423232312ZYXmmvuZC 8其中， 与 为 与 点在各自图像中的齐次坐标，Tvu,1Tv,21p2为点 在世界坐标系下的齐次坐标； 分别为 的第TZYX,Pij4,32;,jKM行和第 列元素.将 式展开后共产生六个式子，按将一式除以二式，二式除以ij7816三式.依次下去消去 和 得到四个线性方程：1CZ2 3422423322322132 111 34243312312131 1111 mvZmvYmvXmv uuuu 9只需要其中三个就可以解出 ,就是说，式子 中，只有三个独立方程，因为Z9假设 与 是空间同一点 的对应点，所以就相当于已经假设直线 与直线1p2P1pO一定相交，或者可以说，四个方程必定有解，而且解是唯一的，在现实应用当2O中，由于数据总是有噪声的，因此必须采用最小二乘法求出 .ZYX,(1)这种相机的几何关系最简单，但一般情况下，这样的配置很难做到，因为在实际相机安装时，可以用相机标定方法求出任意配置的双相机的投影矩阵，因此用上面给出的普遍方法求解出 点的三维坐标.P五、模型的优缺点与推广模型的优点：（1）本论文比较全面的对每个问均有建立模型,且稳定性好.（2）同传统标定方法相比，该方法简单、快捷、鲁棒性好，可用于现场标定.（3）算法向并行化发展，提高速度，减少运算量，增强系统的实用性.模型的不足：本文的第四个问只给出了两部相机的相对位置模型，并没有对其进行检验.本模型有着广阔的科学应用前景，不仅能应用于机器人导航、工业检测、生物医学、反求工程、虚拟现实等领域，还有可能应用于航天遥测、军事侦察等领域.六、参考文献1 赵继广，李建峰.一种基于平面标定物体的摄像机标定方法.装备指挥技术学学报，第 15 卷，第 2 期，2004 年 2 月.2 周立丽，车仁生.一种基于立体标定件的双摄像机外参数标定法.航天计测技术，第1725 卷，第 4 期，2005 年 8 月.3 鲁建波，徐培全，姚舜，唐新华.圆形扫描结构光传感器的标定和焊缝检测.焊接学报，第 27 卷，第 4 期，2006 年 4 月.4 陈西.摄像机标定与三维重建研究.北京化工大学硕士研究生学位论文.2007 年 6 月.5 苏小华.基于双目立体视觉的计算机三维重建 .南京: 南京航空航天大学，2003.D6 姜启源，谢金星，叶俊.数学模型.高等教育出版社，2003 年 8 月.7 姜树民，刘德朋 .高等几何学.陕西人民教育出版社，2000 年 10 月.附录（一）：A=12,112,1,0,0