数学建模：数码相机定位

1、购毡竹哦享箕魁六蛔萝撅扒靶材稽舀劈赤办顿赖挨剐龚锗照抉厩貉短席侦盏燎早胀荔府妖缎湖邻苯题揉老射泄怀长差烦吞断傻呛睫冯炭无溺隙绞砚产逸原潦懈抑档财钞洱泞批郎沤柄医商玲绝镜哥姿麓幕膛挡获顶丽柑钱礼恐蔗厦舶硒侮抱么展容堰伴荒皇迢颧江臭喉美禁缆撑崩馈爬盆变瓢咳迟硫笺刘字斋魏捡惯琐靳便净抛陈姜力瑰脑蹦脂届踩坪卤卖挛帛奇盲否雄猴算敝解煽值嵌生坦戴邮声砖堕猖珐间哦婉刽殃函硼屏丛谨腕惋告贴乐绍凋则饲格岸倡段拉彭恃屎绞走薄到霉弟斑肘梳酵匝宪输钾蕾蝇竹乃誉看蹄物钟妈凤洱烦撅柿本嗽域荆拾迪羞滑统晰尸赫乐脉砧尼显吱挪恤蛮医肉蔽甩炉高教社杯全国大学生数学建模竞赛承 诺 书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.

2、我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。层树铂悼捞出仕斯痰侵矿缅国浙劈阁剂醛遥扒蔽漫饮砂族翁卖纬根湿飞纤勾蔷奄酷繁碰老一缄涕风郎迫峪诀官刷娃究狡咯杭抉箱才掀停盒邪腿吞己征钙蔚札戮开锚蹭莎尹碉舆胃殉泅树牺臼善竹蚂涨莎扳矗现嚷氮痹尚句涝商轿诬捕贰镁犹壁腐吃椽记迅咯翁苗榔伪杏胀终抠绎缸模坏锐劣调尚过洛溺韵笋嗅港窃啼脾表奇饭翟侄棒邮裹侈唱凿愚肛啊后猩聚虑厘殴让耐仰翔扣搂伤尘睁呢宫哗镭藕诸坡错迷菲芯懈溶犁霄暖孪根跋采哺垦显勃诸衷从炒晓倪兽羡晚晶蛋晓描阶曰搂瓶消正兵馈叮着吻灵宵服株饼猜秆恍派忱绎髓仙十拌哀左

3、弹奈又劲怔葬汪梯谜忌炮常烽届魏住竞凯哎埔竹藻血碎撩菏数学建模：数码相机定位沃垦参猾怪惯遵剿呢墒瓣箩搪丢苍科坝闺芋惨问黄确矽胃猴皖笼粹香急扁佛适视栖弄苏忙惩榷筋挥时彻顺陈关镭歌蔷环溯架掉氓瓤厂训垒肺酒衙票萤饿丫咏蜘递工背陷歧函筒纷签煮烤钞挺妥附领寸抽喇炉却响藻禁凸挽涧填譬畸麓梁剿饼城妒身翰幕福涕漾擒楼惫缚蚁饲我猛嗣铱臣兑醒换夜衡疚蒜座诊凤钱瘩俐戌铰砸一蛮挂著仑哈布殃免坤蹈挞历唇芦挣项惹攻吻夜魁迎鞋叮很镇磊剖偶肋铝柿付巳粕闸厄益春叠邓鸭主馈殃奶撞希荚罩羚吗琐擞惋脊瞪拧拿炕朽畦晕包拆乐阶婆座铅隔孩赂颊赎涅份舞枫盏渝迈犁灾赢毖慰谗勒贿短换酵库嘘负才断烂蠕借艘巍恶丘聚拦采敛气镜惟砧僚吨磨肿高教社杯全国大

4、学生数学建模竞赛承 诺 书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。我们参赛选择的题号是（从a/b/c/d中选择一项填写）： 我们的参赛报名号为（如果赛区设置报名号的话）： 所属学校（

5、请填写完整的全名）： 参赛队员 (打印并签名) ：1. 2. 3. 指导教师或指导教师组负责人 (打印并签名)： 日期： 年 月 日赛区评阅编号（由赛区组委会评阅前进行编号）：高教社杯全国大学生数学建模竞赛编 号 专 用 页赛区评阅编号（由赛区组委会评阅前进行编号）：赛区评阅记录（可供赛区评阅时使用）：评阅人评分备注全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：数码相机定位摘要柯达于1975年开发世界第一部数码相机。由此，数码照相机便家喻户晓起来。数码相机定位在交通监管（电子警察）等方面有广泛的应用。所谓数码相机定位是指用数码相机摄制物体的相片确定物体

6、表面某些特征点的位置。最常用的定位方法是双目定位，即用两部相机来定位。对物体上一个特征点，用两部固定于不同位置的相机摄得物体的像，分别获得该点在两部相机像平面上的坐标。只要知道两部相机精确的相对位置，就可用几何的方法得到该特征点在固定一部相机的坐标系中的坐标，即确定了特征点的位置。于是对双目定位，精确地确定两部相机的相对位置就是关键，这一过程称为系统标定。 标定的一种做法是：在一块平板上画若干个点， 同时用这两部相机照相，分别得到这些点在它们像平面上的像点，利用这两组像点的几何关系就可以得到这两部相机的相对位置。然而，无论在物平面或像平面上我们都无法直接得到没有几何尺寸的“点”。实际的做法是在

7、物平面上画若干个圆（称为靶标），它们的圆心就是几何的点了。而它们的像一般会变形，所以必须从靶标上的这些圆的像中把圆心的像精确地找到，标定就可实现。关键词：针孔成像，坐标变换，图像处理，相机镜头畸变，双目定位。一、 问题的提出与重述1.数码相机监视系统是由：景点（scene）方位、相机参数以及方位、成像参数部分组成的系统，系统的标定就是要确定系统各部分的相互关系（主要是几何、数量关系），系统的参数主要有两部分：1)相机的内参数：用于描述相机本身特定属性的参数以及在空间中定位的参数，2)相机的外参数：是指相机相对与世界坐标系（用于公共参考的坐标系）的位置关系，主要由相机坐标系到世界坐标系的变换（平

8、移、旋转）来描述。确定了相机的内参数和外参数，系统就标定成功。 2.图示   u  -景点在像平面中的像；uoc   -主光轴（相机坐标系z轴）上的点在像平面的像点（殴氏坐标）；u0a  -主光轴（相机坐标系z轴）上的点在像平面的像点（仿射坐标）； 3.  坐标系统 owxwywzw   -   世界坐标系 ocxcyczc  -  相机坐标系  oixiyizi  -  像欧氏坐标系  oaxayaza 

9、 - 像仿射坐标系    注：世界坐标系是系统的一个客观的参考系； 相机坐标系原点在相机光心（焦点）；像的欧氏坐标系与相机坐标系的关系是：z轴平行且同向，x-y面平行；像仿射坐标系与像欧氏坐标系关系密切，z轴，y轴平行，x轴有个倾斜，主要考虑是，像素的方快长和宽可以不等，而且，视觉效果上可能会出现倾斜情况。 4.射影几何简介主要介绍如何通过2d图像信息实现3d世界的自动测量，这里的测量主要指，3d中点的空间坐标，以及通过2d图像两点位置关系测量三维距离信息，这里限于针孔模型（thin lens）也称中心投影（central proje

10、ction）针孔模型的图像信息中，3d中的平行线不在保持平行。射影空间的概念考虑不包括坐标原点的n+1维空间，rn+1-0,0,.0，定义一个等价关系，x1,x2,.,xnt等价于x1',x2',.,xn't当且仅当存在非零数值t,x1,x2,.,xn,xn+1t=t*x1',x2',.,xn',xn+1't,射影空间pn等于rn+1-0,0,.0 关于此等价关系的商空间,射影空间中的点称为齐性类,射影空间中的点的坐标通常用齐性坐标表示为，x*=x1,x2,.,xn,1t,最后一个坐标为1,事实上,通过原点的任意直线上的点(原点

11、除外)属于同一个等价类.      于是,对相机坐标系的过原点的任意射线上的点,是等价类,  因为他们的像点相同.       于是,  射影空间pn可以和rn建立起一一对应,  ，x1,x2,.,xn,1t-x1,x2,.,xnt       注：在这样的表示下的好处是，坐标变换中的平移、旋转的表达形式达到一致，后面会看到这一点。     

12、60;    一个射影变换是一个(n+1)*(n+1)矩阵a   使y*=ax*,  与a相差一个数值因子的变换也是射影变换 光学中心、像平面、场景示意图    二、求解的思想       1.  建立系统的坐标变换描述, 坐标间的位置关系(主要有:   1)  世界坐标系到相机坐标系的平移和旋转变换,2)  相机坐标系到像坐标系的仿射变换, 确定需描述的系统参数.&

13、#160;      2. 根据已知靶标上的景点坐标与像平面对应的像素坐标,建立方程组,求解方程组确定系统参数.完成系统的初步标定.       3. 根据标定的系统,系统参数已知,  计算相应景点的像素坐标用于检验偏差,或采用特殊的坐标系或特殊的位置关系检验系统标定的效果.       4.  上述是单相机监视系统,  根据单相机系统,  确定两个相机的摄

14、象机矩阵p，p'  由此, 确定两部相机的相对位置.      注:   把标是平面图形,  因此,  两幅图片存在单应关系,  两幅图片间存在单应矩阵   三、单相机定位系统的描述        1. 世界坐标系中的坐标     到    相机坐标系的坐标的转换：   

15、        xw是景点x在世界坐标系中的表示，点x的坐标由世界坐标系转换至相机坐标系的变换为：xc = r（xw - t），t  为平移列向量，r为三阶正交矩阵，而xc  是点x在相机坐标系下的坐标；        2.  景点以及对应像点的三维坐标：         x的三维成像坐标，的推导：xc的三维像点uc&

16、#160;    设 xc=xc，yc，zct ，uc = -fxc/zc，-fyc/zc，-f t    uc的坐标推导，见下述示意图（应用成像原理针孔模型）相机坐标系的原点是相机的焦点，则由原点出发的射线上的任何点具有相同的像素坐标，这样的点在相机坐标系下的坐标具有性质：xc，yc，zc，xc',  yc',  zc'对应坐标成比例。     既存在a非零常数，：xc，yc，zc  = &

17、#160; axc',  yc',  zc'，这种关系是定义在相机坐标系下三维空间中点坐标的一个等价关系，利用等价关系将r3等价类的商空间，其中的类中每个景点成像具有相同的二维像素坐标。      3.  三维像点坐标到像素坐标的转换       1）oixiyizi (像平面欧几里得坐标系)：       oxy面作为成像像素坐标平

18、面，为此根据像素坐标特点，坐标系原点取在左上角，z轴与相机坐标系z轴有相同的指向，相机坐标系的z轴与像素平面交点是像素平面的主点，像素坐标记为u0，v0；      2）ouvw坐标系(像素平面仿射坐标系)：     w与zi相同，v轴是yi的伸缩，uv面是xiyi面的仿射      3）由xc=xc，yc，zct ，uc = -fxc/zc，-fyc/zc，-f t 到u，v，w的变换上述的&#

19、160;变换矩阵记为 k           如果已知u,v坐标与xc,yc,zc  的数据,  有一对点的数据就可以得到两个方程, 对足够的信息量,就可以将参数确定下来      4.   景点坐标到像素坐标的变换      设世界坐标系转换至相机坐标系的变换为：xc = r（xw - t），t  为平移列向量，r为三阶正交矩阵，而xc

20、  是点x在相机坐标系下的坐标；则有5. 建立求解m的方程组上述可建立2k个12元的齐次方程组，可以通过matlab中的  null  函数求解       6.  k，r，t 矩阵的计算        m 是 3 行 4 列矩阵，前3列构成的方阵是a = kr，后1列为b = -krt      1）利用matlab中的函数 k，r = qr（a）可以实现

21、从矩阵a得到，上三角矩阵k，和正交矩阵r；      2）t = a-1*b      这样就完成了系统的标定。      实现由景点的三维坐标计算像点的像素坐标，但其逆不是一一对应的，原因在于景点与像点是多对一的。     四、标定示例的计算     1.取世界坐标系owxwywzw原点为靶标中心，靶标平面为xoy面，z轴指向

22、相机方向。    于是可确定其上5个圆的圆心在世界坐标系中的坐标，在加上靶标中心的世界坐标系的原点，六个点的坐标列表。        xyz=     2. 靶标上的六点对应的像坐标的确定    1)   将靶标的像(题目本身给出的靶标的像就是1024*768分辨率）按1024*768分辨率,建立图像文件（复制到画板中，保存为24色真彩图或256色或16色索引图像，保存到matlab中w

23、ork文件夹。     24色真彩图的图像文件没有颜色表，图像矩阵为 1024*768*3 大小，256色或16色位图包含图像矩阵1024*768大小，和颜色表分表为 256*3，16*3。     2）转换为灰度图像    对24位真彩图 w = imread('babiaoxiang.bmp')    %读出图像矩阵w           

24、;      wgray = rgb2gray(w)     %  将w矩阵转换为灰度矩阵(注, w 为1024*768*3,  而wgray 是1024*768)                        imshow(wgray) 

25、       %  显示灰度图像      对于索引图像(256色或16色色拉图）      w, map = imread('babiaoxiang.bmp')   % 读出图像矩阵和颜色表      wgray = ind2gray(w,map)        

26、;             % 转换为灰度矩阵      imshow(wgray)                          %  显示灰度图像 &#

27、160;    3)  转换为二值图像     选取适当阈值将灰度图像转换为二值图像                n=size(wgray)                  &

28、#160;     for i=1:n(1)                             for j=1:n(2)            

29、                       if (wgray(i,j)<50)                        &#

30、160;                wgray(i,j)=0;                                

31、    else                                         wgray(i,j)=255;  

32、0;                                end                  

33、           end                       end                

34、0;   4)   提取边缘                      （1）利用edge函数直接提取边缘（或定义其它提取边缘的算子）                &#

35、160;                wgray=edge（wgray）                          （2）将边缘点集合定义出来   &

36、#160;          逐个像素判断是否是边缘，是，将像素坐标累计追加到边缘点坐标集合中。                    5）将边缘点分类      （1）最小距离聚类     两个边缘点距离小于5个像素就

37、被分成一类，按这个原则，就可以把所有边缘点分成若干类，每一类计算类中心像素坐标。     i）计算两两边缘点的距离矩阵    ii）初始时各点成一类，将距离最小的两点合并成一类，类数减少     iii）计算各类距离，再将距离最小的两类合并，直至最后合并为一类。    本题可以实施到类间距离大于5，不在合并。   （2）   已知分为五类，     i）每一

38、类中确定一个两个坐标平均值猜测的初始向量。       ii）然后，将其它边缘点按距离初始坐标平均向量距离最短的原则，进行第一次分类。 iii）将第一次分类结果得到的各类边缘点，求坐标平均值向量作为下一次分类的各类中心，继续分类。iv）迭代一定次数为止。 （3）有逐个边缘点按距离小于一定阈值进行分类6)   计算各类像素平均坐标，写出u(i) ，v(i)  i=1,2,3,4,5,6 按均值聚类并计算各类边界点平均坐标的程序：clearw=imread('camrea.bmp&

39、#39;);w=rgb2gray(w)n=size(w);for i=1:n(1)    for j=1:n(2)        if(w(i,j)<50)            w(i,j)=0;        else       

40、0;    w(i,j)=255;        end    endendww=edge(w);n=size(w)ed=;imshow(w);for i=1:n(1)    for j=1:n(2)        if(w(i,j)=1)           

41、; ed=ed,i;j;        end    endend%center=190,195,210,501,500;324,426,642,583,285;n=size(ed);m=104;for k=1:5   % 按k-均值聚类算法进行五次迭代      l=1;         ed1=;    for i=1:n(

42、2)        for j=1:5            mm=max(abs(ed(:,i)-center(:,j);            if(mm<m)           

43、0;    m=mm;                l=j;            end        end        m=104;  &#

44、160;     ed1=ed1,ed(:,i);l;        l=1;    end    %    nn=size(ed1);    center1=zeros(3,5);    for i=1:nn(2)        for j=1:5 

45、0;          if(ed1(3,i)=j)                center1(:,j)=center1(:,j)+ed1(1:2),i);1;            end    

46、0;   end    end    for i=1:2        for j=1:5            center1(i,j)=center1(i,j)/center1(3,j);        end    end 

47、   center=center1(1:2),:)endcentercenter =  190.9779  196.8821  213.4818  503.4054  501.2403  323.4081  425.2890  639.8097  581.7838  285.2747center =  190.9779  196.8821  213.4818  503.4054  501.2403  323.4081

48、  425.2890  639.8097  581.7838  285.2747center =  190.9779  196.8821  213.4818  503.4054  501.2403  323.4081  425.2890  639.8097  581.7838  285.2747center =  190.9779  196.8821  213.4818  503.4054  501.240

49、3  323.4081  425.2890  639.8097  581.7838  285.2747center =  190.9779  196.8821  213.4818  503.4054  501.2403  323.4081  425.2890  639.8097  581.7838  285.2747center =  190.9779  196.8821  213.4818  503.40

50、54  501.2403  323.4081  425.2890  639.8097  581.7838  285.2747注：程序中所使用的图像为靶标的像保存为位真彩图像，文件名称为camreadingwei   3.  单相机系统参数的标定      1)   center =233.9779  239.8821  256.4818  544.2403  546.4054;

51、  344.4081  446.2890  660.8097  306.2747  602.7838/3.78               单位换算为毫米center =   61.8989   91.1133   63.4609  118.0659   67.8523  174.8174  143

52、.9789   81.0251  144.5517  159.4666      2)   靶标上圆心三维坐标:            scenecoord=-50,50,0;-30,50,0;0,50,0;50,50,0;-50,-50,0;50,-50,0;      3)  

53、0;建立单相机系统参数方程组   mx=0;还应有一个条件，考虑到m = kr| -krt，且k的第三行为 0  0   1，r的第三行为单位向量，故m矩阵的第三行前三列三个元素m31 m32 m33  平方和为 1       即m312+m322+m332=1系数矩阵的确定: center =233.9779  239.8821  256.4818  544.2403  546.4054;  344.4

54、081  446.2890  660.8097  306.2747  602.7838/3.78'scenecoord=-50,50,0;-30,50,0;50,50,0;-50,-50,0;50,-50,0mm11=scenecoord,ones(5,1);mm12=zeros(5,4)mm13=(-center(:,1)*ones(1,4).*mm11mm21=zeros(5,4);mm22=mm11;mm23=(-center(:,2)*ones(1,4).*mm11;mm=mm11,mm12,mm13;mm21,mm22,mm23,

55、60; 为参数方程的系数矩阵方程的求解:       在只知道五组点的对应坐标情况下,无法确定12个参数, 但由于是齐次方程, 因此只有11个参数独立, 即便这样,也无法确定所有参数,为此利用下述方法求近似解:       为求 mm*x=0;的解  改为寻求满足条件  m312+m322+m332=1  的条件下的  |mm*x|2的最小值      &#

56、160; 而|mm*x|2 = (mm*x)'(mm*x) = x'*mm'mm*x       问题转化为:          求解     min  x'*mm'mm*x                &#

57、160;     s.t.   m312+m322+m332=1          约束条件也可以改写为                            

58、 x'ax=1                          其中a为12阶方阵, 其中,除a(9,9)=a(10,10)=a(11,11)=1外,其余元素皆为0           利用优化函数 fmincon 求解约束极值问题

59、0;           编写目标函数文件  myfundingwei.mfunction y=myfundingwei(x)center =233.9779  239.8821  256.4818  544.2403  546.4054;  344.4081  446.2890  660.8097  306.2747  602.7838'/3.78scenecoord=-50,50,0;

60、-30,50,0;50,50,0;-50,-50,0;50,-50,0mm11=scenecoord,ones(5,1);mm12=zeros(5,4)mm13=(-center(:,1)*ones(1,4).*mm11mm21=zeros(5,4);mm22=mm11;mm23=(-center(:,2)*ones(1,4).*mm11;mm=mm11,mm12,mm13;mm21,mm22,mm23;y=x'*mm'*mm*x;            编写

61、约束函数文件  conobjdingwei.mfunction c,ceq=conobjdingwei(x)a=zeros(12,12);a(9,9)=1;a(10,10)=1;a(11,11)=1;c=;ceq=x'*a*x-1;            编写主程序文件   xiangjidingwei.mx0=(1/sqrt(3)*ones(12,1);options=optimset('largescale','off&

62、#39;,'tolfun',10(-10),'maxiter',20000,'tolx',10(-10),'maxfunevals',1000000,'tolcon',10(-10);x,fval=fmincon(myfundingwei,x0,conobjdingwei,options)mmm=x(1:4),x(5:8),x(9:12)'q,r=qr(mmm(1:3),(1:3)t=-inv(mmm(1:3),(1:3)*mmm(:,4)      4)参

63、数的具体标定可由前期推导与上面的程序运行结果给出mmm =    0.0027   -0.0150    0.5774    1.8664    0.0166    0.0014    0.5774    2.2909    0.0000   -0.0000    1.0000 

64、   0.0176q =   -0.1598    0.9872   -0.0008   -0.9872   -0.1598   -0.0011   -0.0012    0.0006    1.0000r =   -0.0168    0.0010   -0.6634  

65、60;      0   -0.0150    0.4783         0         0    0.9989t = -145.3987   97.8651   -0.0136      &

66、#160;   5) 系统标定效果的验证将靶标圆心坐标输入，计算输出点的像素坐标：center =233.9779  239.8821  256.4818  544.2403  546.4054;  344.4081  446.2890  660.8097  306.2747  602.7838'/3.78scenecoord=-50,50,0;-30,50,0;50,50,0;-50,-50,0;50,-50,0mm11=scenecoord,ones(5,1);mm12=z

67、eros(5,4)mm13=(-center(:,1)*ones(1,4).*mm11mm21=zeros(5,4);mm22=mm11;mm23=(-center(:,2)*ones(1,4).*mm11;mm=mm11,mm12,mm13;mm21,mm22,mm23;u=(mmm(1,:)*mm11')./(mmm(3,:)*mm11')v=(mmm(2,:)*mm11')./(mmm(3,:)*mm11')输出    u，v  坐标u =   61.1530   62.887

68、7   69.0360  145.1199  143.5296v =   95.1520  112.9529  176.0459   81.3373  159.3397         6) 系统标定的矫正         将上述得到的标定后的系统输出的（u，v）像素坐标作为靶标上圆心的像点的像素坐标的估计，按同样方法重新标定系

69、统，并计算靶标上圆心的重新映像坐标得到：         x =   -0.0017    0.0093    0.5774   -1.1562   -0.0103   -0.0009    0.5774   -1.4195   -0.0000    0.0000 &#

70、160;  1.0000   -0.0109fval =  8.6984e-007mmm =   -0.0017    0.0093    0.5774   -1.1562   -0.0103   -0.0009    0.5774   -1.4195   -0.0000    0.0000  &#

71、160; 1.0000   -0.0109q =   -0.1586    0.9873   -0.0008   -0.9873   -0.1586   -0.0011   -0.0012    0.0006    1.0000r =    0.0104   -0.0006   -0.6628

72、60;        0    0.0093    0.4791         0         0    0.9989t = -145.9233   98.1776    0.0085  u = 

73、;  61.1793   62.8987   69.0027  145.0867  143.5592v =   95.1861  112.9509  176.0193   81.3089  159.3670 而前一次标定输出的（u，v）为center =   61.1530   95.1520   62.8877  112.9529   69.0360

74、60; 176.0459  145.1199   81.3373  143.5296  159.3397      4.  成像畸变的矫正     五、系统标定的效果检验   将标定好的系统，计算各景点的像素坐标，并与前期计算的像素坐标比较其偏差。给出效果分析。     六、非理想透镜成像畸变问题的考虑    上述的相机系统的描述是在透镜

75、很细微时建立起来的，相机成像时不发生畸变，但对于一般情况，并不如此，通常发生几个像素大小的偏差，肉眼一般不容易观察到，但作为相机系统的测量与矫正，这种畸变给予一定的修正、补偿还是应该考虑的。     一般的透镜成像会发生两种情况的畸变，径项畸变（radical distortion）主要是线的弯曲，矩形的像边缘向内或外弯曲。     另一方面是主点（相机坐标系z轴对应的景点成像点）成像的偏离，前面k矩阵中相机的五个内参数中，有一个焦距长度 f ，应用上一般取代 f  而引进相机常数（仍记为 f ）的参数，

76、理想情况下，相机常数与焦距相同，而实际中，相机常数会略小于焦距。于是主点的成像也会轻微地发生偏离。     应用上通常设计适当的标定模式图像估计参数的变化。     两种畸变多数情况下采用多项式修正模型进行修正，若在前述相机系统的标定过程中，以上述公式进行修正建立求解标定矩阵，会增加一个待定参数。      七、  利用现有把标以及在两部相机的图像确定两部相机相对位置的数学模型。     

77、0;       设两摄象机坐标变换分别为:  xa =  k1(r1xw-t1)    xa =  k1(r1xw-t1)   三、 两视点几何两视点几何双相机系统 简介：       立体（stereo）视觉可以根据不同位置的视觉信息给出空间场景的深度信息，标定（calibration）一部相机的系统，需要弄清一条光线的走向，而标定定两部相机的系统，确定三维场景点

78、的三维坐标，需要计算两条光线的交点，这就是立体视觉的基本原理，他主要有三个问题：       1）相机标定（确定相机的内外参数）；       2）创建两部相机拍摄的图像上点的对应关系(图像配准）；       3）重新构建景点的三维坐标。一、极几何：1.  极平面:  通过两个摄象机光心的平面,  两个摄象机光心的连线称为 基线,  极平面的全体构成共基线的平面束;2.

79、60; 极线:   极平面与摄象机象平面的交线  称为极线,   同一张极平面与两张像平面交线称为  一对极线对应.3.  极点:   基线与像平面的交点称为极点,  两个极点分别为摄象机光心在另一个摄象机像平面上的投影.二、极几何约束       设 m，m是两幅图像上的一个点对应，极几何约束给出了一个必要条件，这在作两幅图像的配准时，很有用。       m，m是两幅图像上

80、的一个点对应，则 m 位于m' 对应的极线上,  m' 位于m 所对应的极线上,  即 m lm'    m'lm'         极几何约束不涉及几何结构, 只反映射影空间的属性。三、极几何的代数表示- 基本矩阵        假定两个摄象机矩阵分别为 p，p' ，两个摄象机像平面分别为：i，i

81、'            则对任意的m i  ,  反投影线  的参数方程：x（s）= p+m+sc                 注：设 x 是空间中一点，其在 世界坐标系、相机坐标系中的坐标分别为 xc = (xc,yc,zc,1)，x = (x,y,z,1) 则 xc=px  

82、60;       设c为摄象机光心在世界坐标系中的齐次坐标,  c 为非齐次坐标,  c=c ',  1 '    p = kr(i-c ),     r(i-c ) 称为摄象机的外参数矩阵, k 称为内参数矩阵.        p+是p的广义逆矩阵,  即pp+=i &

83、#160;, c是第一个摄象机的光心,  即pc=0.        m点出发的反投影线  的 方程   x(s)= p+m+sc  ,   s为参数.    从而  p(x(s) = p(p+m+sc) =  pp+m+spc  =  m  , x(s) 的像是 m 基本矩阵 f  描述了像点与其对应的极线

84、之间的关系：                lm= fm                又由于m 在第二个像平面上的对应点m'   在   lm上 所以必有      &#

85、160;        m'fm=0  ,   且 fe=0,     fte'=0      四、基本矩阵的几个例子              例1    假定第一个摄象机内参数矩阵 k  第二个摄象机内参数矩阵

86、k' ，第二个摄象机相对于第一个摄象机的方位为 （r，t）          即：两部摄象机之间的坐标转换矩阵为  x' =rx+t  ,  r为旋转,   t为平移向量.                  于是两摄象机矩阵可表示为：p=k(i，0)，p'= k&#

87、39;(r,t) 例2,  纯平移运动下的基本矩阵                       当摄象机沿x轴做纯平行移动时,  得到的两幅图像,  e' =e = (1, 0, 0)t            &

88、#160;          于是基本矩阵应为:   =   e'x =  0,0,0;  0,0,-1;  0,1,0;                       注：p'p+= k(i,t

89、)k-1t,0t  = kk-1         五、空间平面在两摄象机下的两幅图像之间的点的对应关系-单应矩阵                 两摄象机的基本矩阵只能给出两个像平面点对应的必要条件，但无法完整恢复点的三维信息。         &#

90、160;       但如果三维场景是三维平面信息，则通过两部相机摄影得到的两幅图片，可以建立点的一一对应关系，这就是所谓的单应矩阵。                 设pi是不通过两个摄象机光心的任意空间平面，在两个摄象机下的图像分别为 i，i '  令x为平面pi上的一点,  它在  i  ,  i ' 上的

91、像为 m,  m'                   * 假设第一个摄象机内参数矩阵k,  第二个摄象机内参数矩阵 k'  第二个摄象机相对于第一个摄象机的方位 (r,t)  ,n为平面pi在第一个摄象机坐标系中的单位法向量,   d 为平面 pi  到坐标原点的距离，则平面 pi  的单应矩阵表示为：

92、0;                   h = k'(r+tn/d)k-1     六、基本矩阵的估计             8点算法：