数学建模-数码相机(高一)

数码相机定位的几何方法摘要本文在数码相机定位问题中，首先建立起4个空间坐标系，并给出它们之间的影射关系，把空间中的三维世界坐标通过平面投影成像过程转换为二维坐标。

对于第一问，根据观测到的靶标的像C B A ''',,三点的中心在同一直线上，利用旋转变换和平移变换以及仿射变换，建立由世界坐标系到像平面坐标系的影射关系。

利用坐标变换，根据已知参考坐标系坐标和图像坐标系坐标的控制点来确定相机的光学和几何参数以及相机相对坐标系的位置和方向变换，得出由物平面到像平面的函数关系式。

对于第二问，由给出的靶标及其像，通过刻度尺量出靶标的像的坐标，利用分辨率和测量值计算出像平面大小的比例系数k ，利用第一问的表达式算出靶标上圆的圆心在像平面上的像坐标。

第三问中为了纠正相机拍摄的几何畸变，本文采用仿射变换来检验模型的准确度。

建立了平面直角坐标变换下的仿射变换，即两直角坐标系w w Y X -和u -v 坐标系，利用不在同一直线上的E C A E C A ''',,,,,六点坐标代入建立的仿射变换方程组中，就可得到各点在u -v 坐标系中的预测值，通过与观测值的对比计算中，即可对所建立的模型进行精度和稳定性的讨论。

第四问中，给出了两种求解方法。

方法一是假定两部相机是平行于物平面放置，经过空间坐标转化，根据几何比例关系得出结论。

方法二则建立空间三维直角坐标系，给出“正像平面”的法向量和“反像平面”的平面向量，且不考虑镜头点的位置改变，从对称及非对称两不同方法分别考虑两相机的相对位置，对任取靶标的像的两点及其正像平面的中心点，分别求得正像平面与直线的交点，利用相似三角形性质、交点坐标、像距从而来确定出两部固定相机的相对位置。

两方法都利用了数形结合，比较清晰的给出表达式。

关键词：相机标定 世界坐标系 小孔成像 仿射变换一、问题的提出1、背景随着计算机硬件、软件、图像采集、处理技术的迅速发展，计算机视觉的理论和技术已被广泛地应用于医学图像处理、机器人技术、文字识别、工业检测、军事侦察、地理勘察和现场测量等。

计算机视觉在多种测量中的应用是一种定量分析系统，有确定的精度要求。

一般用于现场不可到达、现场复杂不便直接测量、保留现场状态以便复测、补测和事后重构现场等场合。

通常人们把摄像机和胶片相机作为获取原始图像的主要设备，但随着数码相机的出现和发展，在许多场合人们已逐渐采用数码相机作为主要的图像采集设备。

由于数码相机具有很多优点，比如可以根据测量的需要选择不同的像素大小、可直接与计算机进行通讯、可选定多种焦距进行定焦距拍摄等，这些特点给图像采集和图像处理带来了许多方便之处。

因此，照相测量也就成为近年来发展起来的一种可用于交通事故勘察与信息采集等方面的智能化技术。

在照相测量技术里，相机标定是最基本的任务。

因此利用物体在相机像平面的位置求出物体在现实世界中的具体相对位置，其实质就是建立像平面与物平面之间的对应关系，这是我们研究与讨论的重点。

2、问题的重述数码相机定位在交通监管(电子警察)等方面有广泛的应用。

所谓数码相机定位是指用数码相机拍摄物体的相片来确定物体表面某些特征点的位置。

最常用的定位方法是双目定位，即用两部相机来定位。

本文将对相机拍摄时的相对位置进行确定性的解决及讨论。

对物体上一个特征点，用两部固定于不同位置的相机摄得物体的像，分别获得该点在两部相机像平面上的坐标。

只要知道两部相机精确的相对位置，就可用几何的方法得到该特征点在固定一部相机的坐标系中的坐标，即确定了特征点的位置。

对双目定位，精确地确定两部相机的相对位置就是关键，这一过程称为系统标定。

标定的一种做法是：在一块平板上画若干个点，同时用这两部相机照相，分别得到这些点在它们像平面上的像点，利用这两组像点的几何关系就可以得到这两部相机的相对位置。

然而，无论在物平面或像平面上我们都无法直接得到没有几何尺寸的“点”。

实际的做法是在物平面上画若干个圆(称为靶标)，它们的圆心就是几何的点了。

而它们的像一般会变形，所以必须从靶标上的这些圆的像中把圆心的像精确地找到，标定就可实现。

有人设计靶标如下，取1个边长为100mm的正方形，分别以四个顶点(对应为A，C，D，E)为圆心，12mm为半径作圆。

以AC边上距离A点30mm处的B为圆心，12mm为半径作圆，就是靶标示意图。

用一位置固定的数码相机摄得其像，就是靶标的像。

就上述所述：(1)建立数学模型以确定靶标上圆的圆心在该相机像平面的像坐标，这里坐标系原点取在该相机的光学中心，x-y平面平行于像平面；(2)对由图1、图2分别给出的靶标及其像，计算靶标上圆的圆心在像平面上的像坐标，该相机的像距(即光学中心到像平面的距离)是1577个像素单位(1毫米约为3。

78个像素单位)，相机分辨率为1024×768；(3)设计一种方法检验你们的模型，并对方法的精度和稳定性进行讨论；图1 靶标示意图 图2 靶标的像二、模型的假设1.假设相机的内部结构正常，尤其是镜片不会出现故障。

2.拍摄时由于物体与相机的距离很远，因此由于倾斜产生的小角度所产生的误差可忽略不计。

3.由于测量误差是难免的，因此假设我们在利用图2进行测量时误差是允许的。

4.假设采用的是定焦数码相机，即相机的焦距是一定的。

三、模型的分析及建立[问题一]相机定位的目的是确定相机的图像坐标系与物体空间中的三维参考坐标系之间的对应关系。

为此需要知道相机的内部参数(即光学和几何参数),及相机的相对外部参数(即外部参考坐标系的位置和方向)。

数码相机图像拍摄实际上是一个光学成像过程，它涉及四个坐标系，这四个坐标系[1]分别是：(1)世界坐标系——根据自然环境所选的坐标系，原点为W ,坐标值用),,(w w w Z Y X 表示。

(2)图像坐标系——坐标原点在像平面的中心，X 轴，Y 轴分别为平行于图像平面的两条垂直边，原点为O ,坐标值用),(y x 表示，是一个平面坐标系。

(3)像素坐标系——坐标原点在图像平面的左上角，X 轴，Y 轴分别平行于图像坐标系的X 轴和Y 轴，原点为'O ,坐标值用()v u ,来表示，是一个平面坐标系。

(4)光心坐标系——以相机的光心为坐标原点，X 轴，Y 轴分别平行于图像坐标系的X 轴和Y 轴，相机的光轴为Z 轴，原点为C ,坐标值用()c c c Z Y X ,,表示。

这四个坐标系的空间位置关系如图3所示由这四个空间坐标系分别将世界坐标系中的信息转换到图像坐标系，最后由图像坐标系转换到像素坐标系。

光学成像的理论模型是针孔模型，根据这个模型由光心坐标系向图像坐标系转换过程符合中心影射或透视影射，可用齐次坐标与矩阵表示⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡100010000001c c c C Z Y X h h y x Z (1) 其中，()c c c Z Y X ,,是光心坐标系中一空间点P 的坐标，()y x ,是P 点在对应图像坐标系中的坐标，h 是拍摄相机的像距。

由世界坐标系到光心坐标系的转换关系为⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⎥⎦⎤⎢⎣⎡=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡111W W W T C C C Z Y X O T R Z Y X (2) 其中，R 为旋转矩阵，T 为位移向量，T O 是元素为0的列向量。

由图像坐标系到像素坐标系的转换关系是⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡∆∆=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡11001001100y x v y u x v u (3) 其中，()00,v u 是图像坐标系原点在像素坐标系中的坐标，y x ∆∆,分别是像素坐标系在X 方向和Y 方向相邻像素间的距离。

由以上三式，整理化简可得⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⎥⎦⎤⎢⎣⎡⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡∆∆=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡11010000001100w w w T C Z Y X O T R v y h u x h Z v u (4) 以上是针对针孔模型下的成像系统的基本标定原理来进行求解的。

其中C Z 及y x ∆∆,可以由以下方法给出[2]。

(ⅰ)确定CZ 1 C Z 可由物平面上的点与光心的连线在光轴方向的投影得到。

(ⅱ)确定()00,v u图像坐标系原点在像素坐标系中的坐标()00,v u 是相机光轴与像平面的交点。

(ⅲ)相邻像素点的距离y x ∆∆,的确定对于给定的靶标，它经过针孔成像，其影射在像平面的图像形状与大小会发生变化，我们称这种变化为图像的失真。

现在，就是要确定这种图像失真的程度，利用所提供的图1与图2中靶标的四条边所在的空间位置与尺寸变化，建立起一一对应的关系，求出较精确的值。

以下就是具体的求解方法。

由题给的正方形ACDE ，边长为a 。

设s 为物平面到像平面的距离，在忽略小误差的情况下，可视为h l s +=。

则当h 一定时，有h s ha x -=1 hs ha y -=1 (5) 又由靶标图，可测得正方形成像后的四个圆A, C, D, E 的坐标分别是()()()()),(,,,,,,,,5544332211v u v u v u v u v u ，则()()()()()()()()⎪⎪⎪⎩⎪⎪⎪⎨⎧-+-=-+-=-+-=-+-=2322324241241324324322212211v v u u s v v u u s v v u u s v v u u s (6) 其中4321,,,s s s s 是正方形的像的四条边的边长，取其中边长最大的两个，不妨令21,s s 是四条边中最大的两条，那么就有2112s s x x +=∆ 4312s s y y +=∆ (7) 由(5)、(7)式，可得()()()⎪⎪⎩⎪⎪⎨⎧-+=∆-+=∆h s s s ha y h s s s ha x 43212)(2 (8) (ⅳ)旋转矩阵R 及位移向量T 的确定假设目标图形围绕原点顺时针旋转，分别围绕X 轴,Y 轴,Z 轴旋转角度ψϕθ,,，因此总的转换矩阵可以表示成绕单个坐标轴旋转时的转换矩阵的乘积⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=ψψψψϕϕϕϕθθθθcos 0sin 010sin 0cos 1000cos sin 0sin cos cos sin 0sin cos 0001R 将以上数据及等式代入⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⎥⎦⎤⎢⎣⎡⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡∆∆=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡11010000001100w w w T C Z Y X O T R v y h u x h Z v u 中，即可得出靶标上圆的圆心在该相机像平面的像坐标。

[问题二]因为相机分辨率为1024×768(1毫米约为3.78个像素单位)，即可得到像平面的大小为()()mm mm 78.376878.31024⨯，即()()mm mm 17.20390.270⨯，利用厘米刻度尺可量得靶标大小为()()cm cm 39.789.9⨯，因为3333.17681024=，3383.139.789.9=，在忽略误差的情况下可近似两者相等，因此可得像平面大小的比例系数4926.2739.717.203==k 。