数码相机定位探讨

1、河北经贸大学毕业论文数码相机定位探讨 专 业：信息与计算科学 摘 要数码相机的应用十分广泛，它不仅应用于人们的日常生活中，还大量应用于科学研究中。本文所要讨论的数码相机定位问题，就是指用数码相机摄制物体的相片以确定物体表面某些特征点的位置。在数码相机定位中最常用的定位方法是双目定位法，即用两部相机来定位，只要知道两部相机精确的相对位置，就可用几何的方法得到该特征点在固定一部相机的坐标系中的坐标，即确定了特征点的位置。本文对一个数学建模问题利用小孔成像理论，来求得给定相片中一点在像平面中的成像点。首先，建立了三个坐标系：相机坐标系（即世界坐标系）、像平面坐标系和物体坐标系，其中相机坐标系原点在相

2、机光学中心，像平面坐标系和相机坐标系相互平行；其次，在物体坐标系中任取一点，经过坐标旋转和坐标平移得到它在相机坐标系中的坐标；然后，利用点和相机坐标系原点建立空间直线方程，直线交像平面坐标于一点，即为所要求点在像平面中的成像点。而对于双目标定问题，首先通过坐标系的旋转和平移得到两个像点在世界坐标系中的坐标，然后将这两个点分别与其所对应的光学中心建立两条直线，这两条直线交于一点，即为所求物体的中心。本文在图像处理中大量应用MATLAB软件，在MATLAB中直接读取图像，调用imread函数将图像转化为二进制矩阵，其中分别用0、1表示成像点和非像点；同时引入像素坐标系，即图像每个像素点位置组成的坐

3、标系，这大大简化了图像信息向数据信息转化的难度。通过对二进制矩阵编程寻找各个成像椭圆的中心坐标，利用所求坐标代入成像模型中，求得旋转矩阵和平移矩阵中所涉及的参数，方便后续问题的分析。在模型检测中同时考察误差分析和灵敏度检验，确保了模型的精确性。关键词： 旋转矩阵；小孔成像；双目定位；像素坐标系ABSTRACTThe application of digital camera pretty much extensive. It is not only applied in people's daily life, but also applied in science research.

4、This text the digital camera fixed position problem for discussing wanted is the photograph that points to shoot to make object with the digital camera with the position that makes sure that the object surface some characteristics order.Positions in the digital camera is a pair of eyes positioning i

5、n the inest common use fixed position method, then position with two cameras, as long as you know the opposite position of two camera precisions, can use several what of method's getting the characteristic's point be mediumly fixing sitting of a camera mark to fasten to sit a mark, then make

6、 sure the position that the characteristic orders.This text learns a model problem to make use of eyelet to become to be like theories to a piece, beg to settle photograph in at 1:00 in the middle of being like flat surface of become be like a point.First, built up 3 to fasten by mark:The camera fas

7、ten(namely the world fastens by mark), is like flat surface by mark by mark to fasten to fasten with object by mark, among them, the camera fasten by mark a little bit original in the camera optical center, be like flat surface to fasten by mark to fasten mutual parallelism with camera by mark;Faste

8、n by mark in the object secondly in allow to take 1:00, has been revolved and sat to mark by mark even move to get it to fasten to mediumly sit a mark by mark in the camera;Then, the exploitation orders to fasten with camera by mark original the point build up a space equation of a straight line and

9、 straightly hand over to be like flat surface to sit a mark at 1:00, for request to order in the middle of being like flat surface of become be like a point.But settle a problem to the double target, pass to sit to mark to fasten first of revolve peace move to get 2 to be like a point to fasten to m

10、ediumly sit a mark by mark in the world, then these two point differences with it to should of the optical center build up two straight lines, these 2 straightly hand over at 1:00, is the center of object begged.This text handles in the in great quantities applied MATLAB software in the picture and

11、directly read in MATLAB picture, adjust to convert the picture as binary system matrix with the imread function, use respectively among them 0, 1 mean be like a point and don't be like a point;Lead into a pixel to fasten by mark in the meantime, namely picture each pixel order position to consti

12、tute of fasten by mark, this difficulty that consumedly simplified a picture information to convert toward the data information.Pass to weave a distance to look for to the binary system matrix each become be like oval center to sit a mark, make use of beg to go into by mark generation be like model

13、in, beg to revolve matrix to peaceably move matrix in the parameter involved, the analysis of convenient follow-up problem.Examines in the model in investigate error margin analysis in the meantime the precision that examine, ensured model is with intelligent degree.Keywords: Revolve matrix;The eyel

14、et becomes to be like;The double eyes position;The pixel fastens by mark目 录一、提出问题1二、基本假设与符号约定3（一）基本假设3（二）符号约定3三、问题分析与模型建立3（一）根据小孔成像原理求空间中任意一点的像平面坐标3（二）利用图像范围，计算靶标上椭圆的圆心在像平面内的坐标7（三）双目标定模型和方法9四、模型求解与检验11（一）模型参数求解11（二）模型检验和误差分析11（三）稳定性检验12五、模型评估13附录14参考文献18 数码相机定位探讨一、提出问题数码相机定位在交通监管（电子警察）等方面有广泛的应用。所谓数码

15、相机定位是指用数码相机摄制物体的相片确定物体表面某些特征点的位置。最常用的定位方法是双目定位9，即用两部相机来定位。对物体上一个特征点，用两部固定于不同位置的相机摄得物体的像，分别获得该点在两部相机像平面上的坐标。只要知道两部相机精确的相对位置，就可用几何的方法得到该特征点在固定一部相机的坐标系中的坐标，即确定了特征点的位置。于是对双目定位，精确地确定两部相机的相对位置就是关键，这一过程称为系统标定。系统标定的一种方法是在物平面上画若干个圆（称为靶标），它们的圆心就是几何的点。同时用两部相机照相，分别得到这些圆的圆心在它们像平面上的像点，利用这两组像点的几何关系就可以得到这两部相机的相对位置。

16、但是这些圆的像一般会变形，所以必须从靶标上的这些圆的像中把圆心的像精确地找到，标定就可实现。设计靶标如下，取1个边长为100mm的正方形，分别以四个顶点（对应为A、C、D、E）为圆心，12mm为半径作圆。以AC边上距离A点30mm处的B为圆心，12mm为半径作圆，如图1所示。用一位置固定的数码相机摄得其像，如图2所示。问题要求：（1） 建立数学模型和算法以确定靶标上圆的圆心在该相机像平面的像坐标（这里坐标系原点取在该相机的光学中心，x-y平面平行于像平面）。（2） 对由图1、图2分别给出的靶标及其像，计算靶标上圆的圆心在像平面上的像坐标, 假设该相机的像距（即光学中心到像平面的距离）是1577

17、个像素单位(1毫米约为3.78个像素单位)，相机分辨率为1024×768；（3） 建立用此靶标给出两部固定相机相对位置的数学模型和方法。图 1 靶标示意图图2 靶标的像二、基本假设与符号约定（一）基本假设 1. 假设由相机所拍摄靶标圆的像平面上圆为椭圆；2. 假设相机为单凸透镜；3. 由于物距（即物体到相机平面的距离）远大于像距，假设物体通过相机成像只考虑小孔成像模型，不考虑相机透镜发生畸变情况。（二）符号约定现对本文中用到的符号进行如下说明：1. ：表示像距，即光学中心到像平面中心的距离。2.三、问题分析与模型建立1对于靶标平面上任意一点，通过相机光学中心映射到像平面上的一点，借助

18、旋转矩阵和空间直线方程找到这两点之间的坐标关系。（一）根据小孔成像原理求空间中任意一点的像平面坐标8、13相机成像原理与小孔成像原理相同，即光线是沿直线传播的。由题目中的靶标和相机所拍图像所知，相机的拍摄位置与靶标所在平面的法向量存在一定的角度，而不是垂直于靶标拍摄的。取靶标上任意点P，把P的坐标经过坐标旋转2和平移转化为镜头坐标系（定义为世界坐标系）中的坐标，过P点和镜头坐标原点建立一条直线L，直线L交像平面于一点Q，即P点在像平面上所对应的像点。再由Q点在Z轴上的分坐标Qz(像距),代入直线方程L，得到Q点的横、纵坐标即所要求的像坐标。建立如下坐标系YZXAqoEDCBycxcDcEc图3

19、 小孔成像原理图如图3所示，就是依题目要求所做的相机坐标系， 为P点的自身坐标系，在坐标系下A、B、C、D、E点的坐标依次为，设P点坐标为（）。在下的坐标为在下的方向余弦,在下的方向余弦,在下的方向余弦。则新旧坐标转换公式2为：推得 （1）把代入上式求得P点在世界坐标系下的坐标为：。在世界坐标系下，过P点与点建立直线方程是：， （2）其中t为参数。设像平面方程为：， （3）联立（2）、（3）得到： （4），此点即为P点所对应像点在相机坐标系下的坐标。由（1）式可以求得A、B、C、D、E在相机坐标下的坐标为：A：， B：, C ：，D：， E：成像平面方程为 则的方程为，A的像点坐标为的方程为，

20、B的像点坐标为的方程为，C的像点坐标为的坐标为，D的像点为的坐标为，E的像点为（二）利用图像范围，计算靶标上椭圆的圆心在像平面内的坐标3首先利用imread函数把像图导入MATLAB软件3、4中，生成一个三维矩阵K，调用MATLAB图像处理函数im2bw生成一个由0和1组成的二维矩阵，其中，（=1,2,768； =1,2,1024），T中的行列为该点的像素坐标。其次对于二进制矩阵T，我们先找出5个成像椭圆圆环的坐标，可以采取以下方法方法12、14：遍历矩阵T，寻找T中相邻的两个元素值之和为1的点（若两个元素值之和为1，即表示这两个元素在0、1分界的边缘上，也就是处在成像椭圆的圆环上），记录下值

21、为0的点的行列坐标,而由向量即得到原矩阵T中5个椭圆的圆环坐标。 图4 靶标圆所对应的成像图片而由向量组成的矩阵把5个椭圆的坐标杂合在一起，现我们分开这5个椭圆，利用5个椭圆的近似位置（如图4），我们发现图像上的椭圆（a，b），c，d，e分别处在像素平面的四个区域中，我们可以先把图像按照行列分别二等分，可以分别把椭圆（a、b），椭圆c，椭圆d，椭圆e分开。而对于椭圆a、b可以先找出椭圆a和椭圆b的坐标中距离最长的两个点，然后再次根据中点把椭圆a、椭圆b分开。接着我们寻找各个椭圆的几何中心。以椭圆a为例，我们先找出椭圆a中距离最长的两个点的坐标，即为椭圆的长轴，再由椭圆的中心对称性质求出椭圆长轴

22、两个点的中点坐标（即为中心坐标）。求得这5个椭圆圆心的像素坐标分别是：a(190.5,322.5),b(197.5,422.5),c(214,640),d(502.5,284),e(503.5,582).最后我们由上述方法所求出的各个椭圆圆心的像素坐标（r,c），根据像素坐标和像平面坐标转换公式7， （5）其中分别表示在行列轴上单位像素的距离；为像平面坐标原点在像素平面下的坐标如图5所示：URVC 图5 像素坐标和像平面坐标转化关系 可以求出各个椭圆圆心的像平面坐标(u,v):圆圆心：（-51.19，-50.132）圆：（-49.339，-23.677）圆（-44.974，33.862）圆（3

23、1.349，-60.317）圆（31.614，18.519） （三）双目标定模型和方法建立坐标系如图6所示，靶标与平面平行，其中一部相机置于原点另一部相机置于,相机坐标系分别是和，每个相机中的像坐标系分别为和。设和在坐标系下的坐标分别是和,P的坐标为,P点在和的像坐标分别是和。YZXOOP 图6 双目标定模型和。A是从到的旋转矩阵，B是从到的旋转矩阵， （6） （7）其中是P1点在中的坐标，是P2在中的坐标，过P1点和O点做一条直线，直线方程是： （8） 过P2和点做一条直线，其方程是： （9）联立方程（8）和（9），得到 (10)即为物体P在中的坐标，其中满足（6）式，满足（7）式。四、模型

24、求解与检验（一）模型参数求解1、5由（1）（4）（5）可知所建立的数学模型中共有12个未知参数，而我们从成像平面上所得的椭圆的坐标共有5个，即含有15个等式信息。因此从其中取出4个椭圆中心，即列12个等式来求解参数。我们取（b,c,d,e）这4个椭圆的圆心来确定参数。然而我们观察发现A、B、C、D、E点这5个点经过相机投影之后，所得像点信息是和无关的变量，因此我们在本题中只找到其他9个参数。通过MATLAB求解所得参数为：， （11）（二）模型检验和误差分析6、15对于模型的检验我们可以考虑改变C点圆心的坐标，即C点移动到点。为了研究方便，我们讨论相对误差。A、C所对应的像平面点分别为。不妨设

25、坐标为（10，100，0），则由上述建立的模型以及MATLAB软件可以求得坐标分别为：为考虑方便，由于像点z轴坐标相同，故只考虑x、y轴坐标表1 物点相对误差表物点AC1C相对误差|CC1|/|AC|坐标（0，0，0）（10，100，0）（0，100，0）0.1表2 像点相对误差表像点相对误差坐标（-0.278,-0.163）（6.958,-48.804）（2.068,-48.119）0.1004由上表可见，上述建立的模型误差较小，比较理想。（三）稳定性检验尝试改变公式（1）中所求得的9个参数中的一个数值，来分析物点所对应像点的变化，原来参数为：式（11）所表值，调整=0.6，得到修正后的像变

26、化表格，如表3所式：表3 物点对就像点变化表XYXYB投影坐标0.4435-14.914C投影坐标2.0682-48.1197修改后B坐标0.9894-15.0059修改后C投影坐标3.7708-48.29840.54950.09191.70260.1787通过对比发现当发生微小变化时，原物点所对应像点就会发生比较大的变化，因此该模型能很好的反映微小变化所引起的变化，灵敏度比较高。五、模型评估1. 本模型巧妙地借助坐标旋转和坐标平移将不同坐标下的空间点转化到同一个坐标系中（即世界坐标系），大大简化了问题的求解，同时借助小孔成像原理建立模型，模型简单实用。2. 本模型在图像处理中，将大量的图像信

27、息通过MATLAB软件处理为数字化信息，同时借助编程寻找圆心，提高了模型的数度。3. 本模型灵敏度良好，较小的变化也能在模型中显示出来。4. 模型在探究成像原理时，只考虑了小孔成像原理，没有考虑非线性畸变情况，模型的应用性不广。附录（本文中用到的程序源代码）：1导入像图，转化为二进制矩阵：z,mat1=imread('G:tupianb.bmp');q=im2bw(z);2找环形坐标：m,n=size(y);s=1;for i=1:m for j=1:n-1 if (y(i,j)+y(i,j+1)=1 if y(i,j)=0 a(s)=i; b(s)=j; s=s+1; els

28、e a(s)=i; b(s)=j-1; s=s+1; end end endenda=a'b=b'c=a b;3分割坐标（5个椭圆）：分割行：function A,B=fengehang(C,z) m,n=size(C); s=1;t=1;for i=1:m if C(i,1)<=z/2 for j=1:n A(s,j)=C(i,j); end s=s+1; else for j=1:n B(t,j)=C(i,j); end t=t+1; endend 分割列：function A,B=fengelie(C,z) m,n=size(C); s=1;t=1;for i=1:

29、m if C(i,2)<=z/2 for j=1:n A(s,j)=C(i,j); end s=s+1; else for j=1:n B(t,j)=C(i,j); end t=t+1; endend 分割（a,b）：function A11,A12=fengeA1(q)o=zuizhi(q);%见下A11,A12=fengelie(q,2*o(2); %plot(A11(:,1),A11(:,2),'*',A12(:,1),A12(:,2),'o');grid on;text(o(1),o(2),'o'); %hold on; %plot

30、(A11(:,1),A11(:,2),'bo')最值函数：function o=zuizhi(q)m,n=size(q);d=sqrt(q(1,1)-q(2,1)2+(q(1,2)-q(2,2)2);h=1,2;for i=1:m for j=1:m v=sqrt(q(i,1)-q(j,1)2+(q(i,2)-q(j,2)2); if v>=d h=i,j; d=v; end endendo=(q(h(1),1)+q(h(2),1)/2,(q(h(1),2)+q(h(2),2)/2;%plot(q(:,1),q(:,2),'*');hold on;%tex

31、t(o(1),o(2),'o');%grid on;4求中点：function o=zhongdian(q)m,n=size(q);d=sqrt(q(1,1)-q(2,1)2+(q(1,2)-q(2,2)2);h=1,2;for i=1:m for j=1:m v=sqrt(q(i,1)-q(j,1)2+(q(i,2)-q(j,2)2); if v>=d h=i,j; d=v; end endendo=(q(h(1),1)+q(h(2),1)/2,(q(h(1),2)+q(h(2),2)/2;plot(q(:,1),q(:,2),'*');hold on;

32、text(o(1),o(2),'o');grid on;参考文献1 姜启源 谢金星 叶俊，数学模型（第三版），高等教育出版社出版，2006年5月；2 陈志杰，高等代数与解析几何（下），高等教育出版社和施普林格出版社出版，2006年4月； 3 王家文 李仰军，MATLAB7.0 图形图像处理，国防工业出版社，2006年6月；4 曾建军 李世航 王永国 叶仁玉 夏慧异，MATLAB语言与数学建模，安徽大学出版社，2005年10月；5 马颂德 张正友计算机视觉计算理论与算法基础北京：科学出版社，1998；6 哈尔滨工业大学机器人研究所，机器人双目视觉系统的标定与定位算法，哈尔滨工业大

33、学学报，2007年11期；7 孔令富等，并联机器人双目主动视觉目标定位的研究，计算机集成制造系统，2007年11期；8 毛剑飞 邹细勇 诸静，改进的平面模板两步法标定摄像机，中国图象图形学报，2004,9(7)；9 罗翔 席文明 颜景平，一种双目主动立体视觉系统的目标定位算法，东南大学学报(自然科学版)， 2002年01期；10 何绪堂等，用于交通事故处理的图像测距研究与实现，现代电子技术，2007 年第24 期；11 熊锋等，一种简单的基于共面的摄像机参数标定方法，计算机工程与应用，2008年第44期；12 张洪涛等，基于网格靶标快速高精度摄像机标定方法，光电工程，第33卷第11期；13 R

34、ENAUD P，ANDREFF N,MARTINET G,et al， Kinemat-ic calibration of parallel mechanism:a novel approach usinglegs observation，IEEE Transactions on Robotics， 2005年；14 ANDREFF N，MARTINET P，Vision servoing of a Gough-Stewart parallel robot without proprioceptive sensors，Proceedings of the 5th International W

35、orkshop on Robot Mo-tion and Control， 2005年；15 Points l-J，IEEE Transactions on Pattern Analysis and Machine Intelligence，2000年。河北经贸大学本科毕业论文（设计）开题报告申请书姓名：刘士龙 专业：信息与计算科学 所属学院：数学与统计学学院指导老师：蔡红艳 开题时间：2009年4月11日1.拟选课题国内、外的研究动态、水平、存在问题，与本人实习、社会实践、调研的关系，并附主要的参考文献：数码相机定位在交通监管（电子警察）等方面有广泛的应用。所谓数码相机定位是指用数码相机摄制物体