利用仿射变换建立双目定位模型-Read

1、利用仿射变换建立双目定位模型摘要：本文介绍了利用仿射变换理论来建立双目定位的数学模型，以仿射变换为理论基础，综合应用仿射矩阵、坐标变换和 MATLAB 编程作图建立两个坐标系之间的关系的数学模型。可以得到物坐标系与像坐标系之间的点坐标的相互转换，利用两台相机拍得同一实物的相片即可实现精确定位。相机照相时发生形状变形，但直线还是直线，符合仿射变换。这里我们采取两种方法来建立模型：方法一建立三个空间坐标系：物坐标系，物原像坐标系和像坐标系，由仿射变换得到物坐标系与物原像坐标系相对应的 3*3的矩阵，矩阵可由最小二乘法再进行求导得到表达式，而物原像坐标系与像坐标系可由透镜成像公式得到，即可建立一个确

2、定像平面坐标的数学模型的算法。方法二直接得到物坐标系与像坐标的仿射矩阵变换关系。解决了问题一。再用 MATALAB 得像坐标系中作为特征点的圆心坐标，可以求得模型的 3*3 的仿射矩阵，解决了问题二中求解物圆心坐标对应的像平面上的坐标。精度是由一个像素点对应的实物坐标系下的大小和算法中通过矩阵范数求得的两部分组成，进行分析比较得到：方法一的误差主要出现在仿射变换中，误差小；而方法二的误差在经过放大的情况下还有仿射变换带来的误差，误差大。所以可以得到方法一的精度与稳定性明显高于方法二。这就解决了问题三中的精度与稳定性的讨论。在这里双目定位是由两部相机对同一实物拍照，得到两个相应的坐标系。将像平面

3、的坐标代入建立的模型可以得到物坐标系中的坐标，这里可以得到两部相机的相对距离，由此问题四得到了解决。同时我们还能根据这两部相机实现定位，达到了双目定位的目的。本文中的模型具有简单，直观的将物坐标定位到像坐标，而且精确度高，误差小的特点。关键词：仿射变换，物坐标系，物原像坐标系，像坐标系，圆心坐标，特征点1 1.问题的重述当今社会，利用数码相机进行双目定位在交通监管等方面有着广泛的应用。双目定位是指用两台固定在不同位置的相机摄得物体的像，从而获得物体特征点在像平面上的坐标。这里运用仿射变换可以建立对实物坐标与像平面坐标的对应关系，只要知道了两部相机精确的相对位置，就可由几何方法得到特征点在固定的

4、一部相机坐标中的坐标。由此确定两部相机的相对位置就成了双目定位的关键，而这一过程即为系统标定。系统标定的实际做法：取 1 个边长为 100mm 的正方形，分别以四个顶点（对应为 A、C、D、E）为圆心,12mm 为半径作圆。以 AC 边上距离 A 点 30mm 处的 B 为圆心，12mm 为半径作圆。（1）建立数学模型和算法以确定靶标上圆的圆心在该相机像平面的像坐标，这里坐标系原点取在该相机的焦点，x-y 平面平行于像平面；（2）计算靶标上圆的圆心在像平面上的像坐标，该相机的像距（即透镜中心到像平面的距离）是 1577 个像素单位（1 毫米约为 3.78 个像素单位），相机分辨率为 1024X

5、768;（3）设计一种方法检验模型，并对方法的精度和稳定性进行讨论；（4）建立用此靶标给出两部固定相机相对位置的数学模型和方法。2 2 . .问题分析这里主要研究的是双目定位的问题，即利用两台固定位置的相机拍得实物的像，建立了三个坐标系：物坐标系 X,物原像坐标系 Y 和像平面坐标系乙实际的圆在物原像上会发生变形成为椭圆（非标准椭圆）。直线还是直线，这符合仿射变换特征。由仿射变换原理，可以假设存在一个 3*3 的仿射变换矩阵，将物坐标变换成物原像坐标（或像坐标）。由标定的特征点可以求得一个从物坐标系到物原像（或像）坐标系的仿射变换的 3*3 的矩阵。再由像距可以得到物原像与像的相似关系，即可得

6、到一个确定像平面坐标的数学模型，可以求得物平面上圆的圆心在像平面上的坐标。而在相机分辨率一定的情况下，可由算法求得像上一个像素对应的物的实际距离就是精度，在本文建立了两种求解方法中进行比较，并对稳定性和误差进行分析。对于求解两部相机的相对位置。可以由另一部相机用同样的方法建立一个数学模型，最后由两个模型的像之间对应的实际距离可以求得两相机的相对位置。3 3.假设与符号3.13.1假设照相中实物发生变形，但直线变直线，符合仿射变换的条件.所以模型中可以运用仿射变换原理。可以将各空间 XYZ 坐标系投影到平面 XOY标系中，Z 坐标统一取为 00假设模型中相机的参数为：实物到相机的距离 D（即拍摄

7、距离）是 5m.镜头尺寸取为 1/2 英寸。3.23.2符号说明X 表示靶标平面上的坐标,X=（i=1的圆心；Y 表示物原像上的坐标 Y=（i=15）为圆心；Z 表示靶标像平面上的坐标 Z=（i=15）为圆心；A 表示从靶标平面到物镜平面的仿射矩阵；Q 表示矩阵 A 的范数；m 表示靶标像的中心到焦点的距离；n 表示焦点到透镜中心的距离；K 表示靶标平面到靶标像平面的仿射矩阵；F 为镜头的焦距长度；V 为拍摄对象的纵向尺寸；H 为拍摄对象的横向尺寸；D 为镜头至拍摄对象之间的距离；V 为镜头纵向尺寸；H 为镜头横向尺寸；d 为两相机的相对距离；h,h,h 为圆心连线的中点。4 4.模型的建立与

8、求解4.14.1问题（1 1）的模型与算法这里可以由两种方法来建立模型。图 1 相机成像原理图图 1 是反映物坐标系，物原像坐标系，像坐标系之间的关系方法一在原实物图即靶标上建立 XiX2平面坐标系，并且取 25 个特征点Xi=(Xii,Xi2,0)i=125普杼而仿射变火韧带由蛙的Y5 5 而透镜外相和流像聃 7 7#而由于我们先是在平面上面讨论变换问题，所以设 ai3=0823=0a33=1,故Yi=AXi(1IY25=AX252TQi=|Yi-AXi|=Yi-AXiY-AXi(1.4)(1.4)式是用精度误差较小的最下二乘法来解出矩阵 A25我们先来求 Q 的最小值 minQ,其中 Q=

9、ZQiEQi100首先用 Q 对 aj求导得并且令=0二 am工-2Y；TXi+XTAXiXTAT-XiT 二 aijy二 aijiJ二 aij首先对一个分量求导即二=0。(即各圆中心和上下左右距点为特征点)在物镜所成的像上建立 Y1Y2平面坐标系，取与 X 轴上的特征点相对应的 25 个特征点Yi=(Yii,Yi2,0)(即各圆中心和上下左右距点为特征点)由仿射变换建立如下对应关系i=12SY=AX(1.1)其中矩阵 A 就是我们建立的仿射变换矩阵。设A=aiia213iai2a22a32ai3a23a33所以aii=a211a31ai2a22a320101FQ：5：:Qi=乙:aijy(1

10、.5)二 aj-2yiixii2aiixii2xiixi2ai2=07.7.2C2c8.8.2yiixi22xiixi2aii2xi2ai2=0入c2cc9.9.2yi2xii2a2ixii-2为函2a22-0a22a3i021可以看到经过求导后展开式中有三个变量在计算过程中消去原因是我们建立模型的时候是在平面坐标系中讨论变换，所以第三维的坐标为 0.由于要保证矩阵 A 满秩可逆，所以我们令a7=0,a23=0,a33=i。所有特征点对应的式子相加便是能够解出 A 中未知变量的线性方程组2525ZPiL=ZBiiiii-2yi2Xi2-2a2iXiiXi22xi2a22=0(1.6)22a3i

11、xii2xiixi2a322a3ixiixi22xi22a32建立方程组的系数矩阵一 2 冷22xiixi200002xiixi22xi220000002xii22xiixi2001002xiixi22x220000002X：2xiixi2:.00002xiixi22xi22常数项矩阵Bi变量矩阵L=a2iP=2yiixii|2yiixi22yi2xii2yi2xi202525运用 matlab 中的矩阵命令计算出矩阵R,Bj的值。i1i1即可解得L=/BJ(1.8)前面说过，由于第三维的坐标假设是零而且要保证矩阵满秩可逆，所以 A 矩阵可以写为allai20IA=a2ia220(1J01一由

12、于现在题目没有给出 Y 的坐标值，所以现在在相机所成的像平面上建立一个坐标系 Z_,-Z2来确定 Y 和 Z 的关系。同样在其上确立 25 个点，Zj=(Zii,Zi2,0),i=125。我们来找这个坐标系和物镜所称坐标系的对应比例关系。照相机的成像原理可得靶标的像与物原像的比例为 m:n即就是m、，z1:y1=m:n,z2:y2=m:n,且2=Y。n可得yi1=n*zi1/myi2=n*zi2/m。(1.10)现在可以得到 Z 的坐标值Zi=(zi1,zi2,0)i=125,可以求得Yi=Wi1,yi2,0)=(n*zi1/m,n*zi2/m.0),将计算结果代入(1.8)式一 2x212x

13、i1xi200002n*zi1xi1/m12 乂1的22x2200002n*zi1xi2/mL=.25Z002x；2xi1xi200 x2n*zi2xi1/mi1=002xi1xi22xil002n*zi2xi2/m00002x：2xixi2000002xixi22x221,10-解出 4 个坐标值，得到矩阵a11a1201A=a21a2200001J方法二在方法一中保持坐标系不变，改为Z=KX(1.12)同理因为我们是讨论平面上的问题所以第三维坐标为零，且为了保证矩阵可逆，所以置k33=1k11k120K=k21k220k31k3211直接求得物坐标系与像坐标系的关系令Wi=|Zi-AXi|

14、25W 八 Wii1同样运用最小二乘法求 MIN(W)解出 Ko 对 W，求导同样的方法求得 K 矩阵4.24.2问题(2)(2)的解法将靶标的像导入软件 MATALAB 中， 由程序 qiuyuan.m(见附录 4.2)得到五个圆的圆心像素坐标， 如图2 所示。/0得出的 25 个方程组相加把系数矩阵求逆的到25J=JY=BX2X2=BYX1是第一部相机在实物坐标系中的位置，X2是第二部相机在实物坐标系中的位置。两相机的相对距离为：d=J|X1_X2|2=J|A,YBAH2=d(A，YB，YT(A,YB丫5.模型的检验根据仿射变换具有定比的特性，我们可以在图 2 和图 3 选取适当的比例点来

15、检验我们模型的正确性。如图 4 所示图 4 对于模型精确度检验，各两个圆心连线的中点分得线段比例不变符合仿射变换，可以代入检验。我们取点 A 和点 E 的中点，点 E 和点 D 的中点以及点 B 和点 C 的中点并计算出它们的坐标值，同样，在像平面上一样的处同样取三个对应的点，计算出其坐标值。AEm=(0,50,0)线段 AE 的中点；(4.4)(4.3)图 4 取两圆心中点进行检验EDm=(50,0,0)线段 ED 的中点；BCm=(65,100,0)线段 BC 的中点。对应像平面上的点h1=(0,41.68135。h2=(34.9885,0,0)h3=(55.55485,83.3627,0

16、)0.71520.12480G=mA=0.12960.78900n-0011将实物图中的三个中点坐标代入变换 Z=GX 中，求得相对应的三个点的坐标hn=(6.4800,39.4500,0)h21=(35.7600,6.2400,0)h31=(59.4480,87.0120,0)=|h-hn|=6.8534.由2=|h2-h21|=6.2875h=|h-h31|=5.3361求平均值.：h=.：h1：h2m3/3=6.1590计算出误差率PI=七/|=0.16p2fA.h2/|h2七 0.17p3fA-/|巾3|=0.053p=0.127这个就是我们的模型在像平面上产生误差的平均值，由此我们可

17、以计算出模型的精确率是 12.7%,由于是 6mm 误差左右，而像平面是 1024*786 分辨率的图，换算后是 271*203(mm)的图， 这样误差差了两个数量级， 所以很小， 说明我们的模型基本可以解决这个问题。6.模型的优化(1)对于变换中我们采用的 Z 坐标为 0,可以说是比较优化了， 但更好的办法是利用三部相机在不同的方位照相，采用在每一个候选位置，详细地模型拟合以确定位置和尺度。根据稳定性选择关键点通过拟合三维二次函数以精确确定关键点的位置和尺度，同时消除低对比度的关键点和不稳定的边缘响应点以增强匹配稳定性、提高抗噪声能力。(2)对于稳定性， 每个关键点周围在选定的尺度下测量局部

18、图像梯度。 为了增强匹配的稳健性，对每个关键点使用 4W 共 16 个种子点来描述，这样对于一个关键点就可以产生 128 个数据，最终形成 128 维的 SIFT(尺度不变特征变换)特征向量(注：SIFT 特征是图像的局部特征，其对旋转、尺度缩放、亮度变化保持不变性，对视角变化、仿射变换、噪声也保持一定程度的稳定)。(3)对于精度分析，要求大测量范围和较高测量精度的场合，采用基于双摄像机的双目立体视觉系统比较合适；对测量范围要求比较小，对视觉系统体积和质量要求严格，需要高速度实时测量对象，基于光学成像的单摄像机双目立体视觉系统便成为最佳选择。另外基于双摄像机的双目立体视觉系统必须安装在一个稳定的平台上，在进行双目视觉系统标定以及应用该系统进行测量时，要确保摄像机的内参(比如焦距)和两个摄像机相对位置关系不能够发生变化，如果任何一项发生变化，则需要重新对双目立体视觉系统进行标定。为了得到更高的精度，应该使摄像机的焦距以及基线长度增大，同时应该使被测物体尽可能的靠近立体视觉系统。另外这个精度和视差的精度有直接的关系。7,模型的优缺点1 1 模型的优点(1)本题采用的两种方法，法一在利用三个坐标系的变换时，物坐标和物原像坐标之间的仿射变换其线形宽度是不变的，只是在角度方面存在变换，并且在物坐标