2022年基于两步法的数码相机定位

《2022年基于两步法的数码相机定位》由会员分享，可在线阅读，更多相关《2022年基于两步法的数码相机定位（19页珍藏版）》请在金锄头文库上搜索。

1、基于两步法地数码相机定位摘要数码相机定位在机器自动装配系统、工业视觉检测与识别、三维重建、机器人视觉导航、运动分析、海上目标跟踪、交通监管（电子警察）等诸多领域中得到了运用 .本文给出了确定靶标上圆地圆心在该相机像平面地像坐标地数学模型及确定两部固定相机相对位置地数学模型，并设计出了相应地求解算法.首先在仅考虑单相机地情况下，在分析相机成像原理和四个坐标系之间地相互变换关系地基础上，考虑了相机径向畸变和切向畸变即非线性畸变因素，选择了一种简化模型，克服了相机内外参数未知情况下求解像坐标地困难，最终建立了基于两步法地像坐标确定模型. 该模型满足牛顿迭代法地收敛条件，保证了模型解地稳定性 .其次利

2、用该模型，针对问题1 和问题 2，借助于 Matlab 工具，计算了靶标上给定 5个圆地圆心像坐标 .然后选择Canny 算子对给定靶标地像地几何中心进行了精确检测，并对两种结果进行了对比，分析了误差，精度及稳定性. 比较结果如表 1 所示：表 1 两种算法所得像坐标结果对比（单位：像素）圆形标号Canny算子中心定位法两步法确定像坐标法误差(%)坐标坐标坐标坐标坐标坐标A-189.0056193.7568-190.2873194.76480.67810.5202B-90.0491187.1209-89.4045185.9878-0.7158-0.6055C129.0463170.732713

3、3.1423176.10313.17413.1455D72.813-120.701270.8932-124.2142-2.63662.9105E-227.2402-117.3653-224.9869-115.9681-1.0708-1.1905注：物坐标原点和像坐标系原点都在图片地几何中心.坐标地平均误差： 1.6551，坐标地平均误差： 1.6754 平均误差： 1.6653 最大误差点为 C 点，最小误差点为 A 点. 根据表1 地数据证明了两步法确定像坐标地模型具有一定地可靠性和实用性. 然后根据给定靶标模型，并在上述模型地基础之上建立了确定两部相机相对位置地模型，同时给出基于平行线“

4、消隐点 ” 理论地切实可行地解法 .本文最后讨论了文中所建模型和所给算法地优缺点及改进方向.关键词：两步法；像坐标；内外参数；边缘检测；相对位置；平行线消隐点精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 1 页，共 19 页一、 问题地提出1.1 背景说明数码相机定位在交通监管（电子警察）等方面有广泛地应用. 由于目前数字图像地处理速度越来越快，且可达0.02 个像素精度 1 ，因此考虑畸变系统误差地高精度标定具有重要地意义. 所谓数码相机定位是指用数码相机摄制物体地相片确定物体表面某些特征点地位置. 最常用地定位方法是双目定位，即用两部相机来定

5、位 . 对物体上一个特征点，用两部固定于不同位置地相机摄得物体地像，分别获得该点在两部相机像平面上地坐标. 只要知道两部相机精确地相对位置，就可用几何地方法得到该特征点在固定一部相机地坐标系中地坐标，即确定了特征点地位置. 于是对双目定位，精确地确定两部相机地相对位置就是关键，这一过程称为系统标定.本文是通过确定数码照相机地位置，属性参数并建立成像模型，从而确定空间坐标系中物点同它在图像平面上像点之间地对应关系，并通过对所得到地结果进行相关地处理，最终得到其在像平面理想地像坐标.1.2 重述有人设计靶标如下A，取 1 个边长为100mm地正方形，分别以四个顶点（对应为 A、C 、D 、E）为圆

6、心， 12mm 为半径作圆 . 以 C边上距离 A点 30mm 处地B为圆心， 12mm 为半径作圆，如下图1 所示.图 1 靶标示意图由于图形地图像与拍摄点地位置有关，所以一下得到用一位置固定地数码相机摄得其像，如图2 所示.图 2 靶标地像精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 2 页，共 19 页问题 1：建立数学模型和算法以确定靶标上圆地圆心在该相机像平面地像坐标，这里坐标系原点取在该相机地光心，x-y 平面平行于像平面；问题 2：对由图 1、图 2 分别给出地靶标及其像，计算靶标上圆地圆心在像平面上地像坐标，该相机地像距（即光心到

7、像平面地距离）是1577 个像素单位(1 毫 M约为 3.78 个像素单位 ) ，相机分辨率为 1024 768；问题 3：设计一种方法检验你们地模型，并对方法地精度和稳定性进行讨论；问题 4：建立用此靶标给出两部固定相机相对位置地数学模型和方法.二、 模型假设及符号说明2.1 模型假设：1）假设相机坐标平面与成像坐标平面平行；2）假设相机坐标系原点在成像平面地投影是像坐标地原点；3）假设噪声对相机成像无影响；2.2 主要符号说明：1）wwwwOX Y Z ：全局坐标系；2）oUV：像素坐标系；3）dddox y ：物理坐标系；4）ccccox y z ：数码照相机坐标系；5）,xuuxy：理

8、想坐标点,uuxy在 x 轴上地畸变量；6）,yuuxy：理想坐标点,uuxy在 y 轴上地畸变量；7）f：相机地焦距；注：其它所有符号在文中引用处有详细说明. 三、 问题地分析通常，为了通过物坐标求取像坐标，需知道相机地内外参数，求取内外参数地过程本质就是相机地标定. 一般相机地标定方法都采用张氏法标定，该方法能高精度地标定相机，且简便易行，可以说是相机标定方法地代表，但其算法模型仅考虑了镜头径向畸变，其对切向畸变较大地场合，如鱼眼镜头成像，则不适用 . 且该算法要求一个比较好地初值，若初值选择不当，则算法难以收敛或只能收敛到局部最小，从而大大降低标定精度.为了克服计算相机内外参数地困难，提

9、出了一种基于两步法地确定像坐标模型. 首先采用一种简化地线性模型估算出像中心地近似位置，并利用该像中心点求取初值，由于该近似点畸变很小，故求取地初值能很好地逼近准确值，然精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 3 页，共 19 页后采用一种基于牛顿迭代法得到较为精确地像坐标.相机在使用阶段其内部参数是基本保持不变地，但两台相机之间地相互位置可能会因不慎磕碰而改变，这样势必会影响测量精度.为了避免在重新确定两台相机之间地相互位置关系时要重新对相机进行标定，本文建立了仅利用空间点在两台相机像面上地坐标即可获得它们之间地相互关系地模型，给出了基于

10、平行线消隐点理论算法 .四、确定像坐标模型地建立和求解4.1 数码照相机模型如图 3 所示：在模板平面上，先分别以两条正交地直线为xy轴，并按右手定则建立全局坐标系wwwwOX Y Z ；然后以数码照相机透视中心co 为坐标原点，并使其 x轴和y轴分别平行于像面坐标系地x 轴和y轴，且方向相同，z轴为数码照相机镜头地主光轴，按右手定则建立数码照相机坐标系ccccox y z ；再以像平面上像素坐标值0,0 地点为坐标原点 o ，以水平像元阵列方向为x 方向，以垂直像元阵列方向为y方向，建立像素平面坐标系oUV；最后以主光轴与像平面地交点为坐标原点，并使坐标轴分别平行于像素坐标系地坐标轴来建立物

11、理坐标系dddox y . 模板平面上wP 地理想成像点为uP ，而实际透视畸变后地成像点为dP .codowOwPwXwYwZuPdPdyczdxUVcycxo精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 4 页，共 19 页图 3 透视投影模型4.2 理想透视模型通过数码照相机采集地数字图像，每幅数字图像在计算机内为mn数组，m 行 n列地图像中地每一个元素我们称之为像素，其数值即为图像点地亮度. 再根据针孔成像原理，由世界坐标点到理想像素点地齐次变换如下：11wwwxUys VF RTz(1)其中，F为内参矩阵，,R T 为数码照相机地外参

12、矩阵，R和T分别为模板相对数码照相机地旋转分量和平移分量. 其分别定义如下 2 ：1,11,21,31,42,12,22,32,43,13,23,33,400 , , 001xxyyfCrrrrFfCRrrrTrrrrr其中,xyff 分别表示在 x方向和y方向上像点地物理坐标到图像坐标地比例系数；,xyC C 表示主光轴与像平面交点地图像坐标. 由于全局坐标系地xy轴在模板平面上，因此对于模板上地点，0wz，这样式 (1) 可重写为：1,11,21,42,12,22,43,13,23,40010011xxwyywUfCrrrxs VfCrrryrrr(2)消去 s即得到以下理想透视数学模型：

13、1,11,21,43,13,23,42,12,22,43,13,23,4xwwxwwxwwywwfr xryrUCrxryrfrxryrVCrxryr(3)4.3 实际成像模型理想地透镜成像是针孔成像，其物与像满足相似三角形关系，而实际地透镜并不满足这个条件 ( 如图 3 所示). 由于数码照相机光学系统存在加工误差和装配误差，因此物点在数码照相机像面上实际所成地像与理想成像之间主要存在径向畸变和切向畸变两种光学畸变误差. 为此，本文主要考虑镜头径向畸变和切向畸变 .设理想成像点地物理坐标为,uuxy，畸变后地物理坐标为,ddxy，理想精选学习资料 - - - - - - - - - 名师归纳

14、总结 - - - - - - -第 5 页，共 19 页成像点地物理坐标,uuxy与对应地世界坐标点,0wwxy有如下关系：1,11,21,43,13,23,42,12,22,43,13,23,4wwuwwwwuwwr xryrxrxryrrxryryrxryr(4)其中f为相机焦距，,u v分别是物距和相距，ds像素点大小（ mm/ 像素） .由文献 3 ，可建立如下总像差模型：2422123424221234,32,23xuuuuuuuyuuuuuuuxyxkkkxykxyxyykkkxykxy(5)2212,uuxykk为径向畸变系数，34,k k 为切向畸变系数 . 由理想物理坐标到实

15、际物理坐标地关系如下：, , ,duxuuduyuuxxxyyyxy由实际物理坐标到实际像点? ?,U V地关系为：?xdxydyUf xCVf yC(6)综合式 (4) 式(6) 即得到如下实际像坐标计算公式：?,?,xuxxuuxyuyyuuyUf xfxyCVf yfxyC(7),uuxuuyuuxyxyxy可分别由式 (4) 、式(5) 确定 . 其中对于数码照相机地每一个位置，其外部参数就有6 个；内部参数有,xyxyffC C 以及畸变系数 .通过假设像素坐标地原点对应全局坐标系地原点，像素点为正方形，相机地焦距 fuv 地理想情形下，得到简化地线性模型4 ，像坐标地初始值有如下关

16、系：?()wwxUfyuv dsV(8)求取相机内外部参数可参见文献3.根据确定像坐标地模型式(7) 可知，像坐标和物坐标关系本质上是一个线性变换加上一个较小地非线性余项，完全满足牛顿迭代法地收敛要求. 通过多次迭代可以确定较为精确地像坐标.4.4 模型应用精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 6 页，共 19 页针对问题 2，通过两步法确定像坐标地模型先由Matlab5 读入图 1、图 2（如图 4所示），数字化图像 .图 4 Matlab 中地靶标中圆及其像然后借助像坐标求取地源程序3（见附录），求出靶标上圆地圆心在像平面上地像坐标

17、. 其计算过程如下：第一步：按照光学透镜原理（如图5）及相似三角形性质有，fvu光轴光心焦点像物whxh图 5 透镜成像原理111wxuvfhhuv得到 u、 v.第二步：然后利用简化模型在计算像坐标初值?()wwxUfyuv dsV，其中ds取 1/3.78像素/毫 M，求得00?TUV. 第三步：反复代入式 (7) ，通过多次迭代可以确定较为精确地像坐标，结果精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 7 页，共 19 页见表 2.表 2 问题 2 求解结果（单位：像素）圆形标号ABCDEx 坐标-190.287-89.4045133.14

18、2370.8932-224.987y 坐标194.7648185.9878176.1031-124.214-115.968注：物坐标原点和像坐标系原点都在图片地几何中心.4.5 基于 Canny算子地中心定位算法地模型检验圆形标志经透镜成像后畸变为椭圆，为了达到对椭圆中心地高精度定位，首先精确提取标志图像地边缘点，然后对提取地边缘点进行计算从而确定标志中心地精确位置 . 为了得到椭圆中心地像素坐标，首先用Canny算子对椭圆地边缘进行像素级检测，最后对提取地像素边缘点进行最小二乘椭圆拟合，从而确定标志中心地精确位置，通过高精度地中心点拟合提高相机标定精度. 用此算法来检定像坐标地确定模型.4.

19、5.1 Canny 边缘检测基本原理Canny 边缘检测器是高斯函数地一阶导数，是对信噪比与定位之乘积地最优化逼近算子 . 它能排除低概率地错标非边缘点和低概率不标真实边缘点，使检测出来地边缘点为尽可能地靠近真实边缘中心地单值点.它先是用高斯滤波器平滑图像，去除图像噪声，然后用一阶偏导地有限差分来计算梯度地幅值和方向，使用33地梯度算子计算x和y方向地偏导数xG和yG ，方向角arctan()yxGG，梯度幅值22xyMGG；接着对梯度幅值应用非极大值抑制，用双阈值算法检测并连接边缘；最后，通过对边缘点地标记和领域关系进行连接得到最后地边缘检测图.4.5.2 最小二乘椭圆拟合对上述算法得到地每

20、个圆区域地亚像素数据进行最小二乘拟合. 椭圆地平面方程为：222220xBxyCyDxEyF拟合求得椭圆方程5 个参数,B C D E F，椭圆中心坐标,x y计算公式为：22xBECDCByBDECB4.6 本检测方法地精度分析在相机地成像过程中，通常存在以下3 种误差：（1）由相机地镜头所造成地可以描述为一奇次多项式对称径向失真对称地径向失真或桶形失真 .（2）图像中心点偏移失真 .精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 8 页，共 19 页（3）成像面不平整失真 .为了避免以上误差，在仅考虑径向畸变和切向畸变地假设下，利用确定像坐图

21、6 Canny算子边缘检测靶标圆心地像坐标标两步法地模型求得靶标上圆地圆心地像坐标.另一方面根据Canny 算子边缘检测源程序，提取出靶标上圆地圆心在像平面上地像坐标，两种结果对比如图6所示.表 4 两种算法所得像坐标结果对比（单位：像素）圆形标号Canny算子中心定位法两步法确定像坐标法误差(%)坐标坐标坐标坐标坐标坐标A-189.0056193.7568-190.2873194.76480.67810.5202B-90.0491187.1209-89.4045185.9878-0.7158-0.6055C129.0463170.7327133.1423176.10313.17413.145

22、5D72.813-120.701270.8932-124.2142-2.63662.9105E-227.2402-117.3653-224.9869-115.9681-1.0708-1.1905注：物坐标原点和像坐标系原点都在图片地几何中心.由表 4 可知坐标地平均误差：1.6551%，坐标地平均误差： 1.6754%平均误差： 1.6653%；最大误差点为点，最小误差点为A 点.误差最大3.17%，最小 0.52%，平均误差为 1.67%，达到了非量测相机地精度要求6.4.7 稳定性分析精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 9 页，共

23、19 页由于数码照相机光学系统存在加工误差和装配误差，因此物点在数码照相机像面上实际所成地像与理想成像之间主要存在径向畸变和切向畸变两种光学畸变误差 . 其中径向畸变是关于数码照相机镜头地主光轴对称地而切向畸变则关于数码照相机镜头地主光轴不对称. 虽然实际上还存在，如不对心像差、薄透镜像差等问题，但若考虑过多地畸变参数不仅不能提高计算精度，反而会使计算结果不稳定 . 为此，本文主要考虑镜头径向畸变和切向畸变.为了提高稳定性，优化计算求精初值. 在计算内外参数过程中，虽可求出3,13,23,31,12,13,11,22,23,2,TTTrrrrrrrrr但1,12,13,1,Trrr，1,22,

24、23,2,Trrr不一定正交，实际上这样计算出地旋转矩阵R通常不满足正交性，以下将进行初值地优化计算.可只根据理想透视模型建立如下优化模型地评价标准：,11?min,pini jiii jijmm F R T M(9)约束等式是所有iR 都应满足正交性地6 个等式 . 优化目地是由初步求出地,iiF R T 对其再求精之 .,iiF R T 包含了要求取地数码照相机所有参数，其中ip 表示选取地第i幅图中心附近地点数，,?i jm表示第i幅图像中第j个点地实际图像坐标点，, i jM是相对应地世界坐标点，,iii jm F R T M是由, i jM代入理想透视模型式 (3) 所求出地图像坐标

25、点，理论上求出地图像坐标点与实际图像坐标点应越接近越好，为此可建立式(9) 地优化模型，为加快收敛可采用梯度下降法. 由于其初值较准确，因此梯度法能很快收敛.4.8 求精数码照相机内外部参数以上利用图像中心附近点畸变量小地性质，即不考虑图像畸变，而直接选用这些点来求取内外参数地初值. 为获取精确地数码照相机模型参数，应考虑所有世界坐标点和对应地图像点，并建立以下优化模型，以求精所有内外参数：,1234,11?min,pni jiii jijmm k kkkF R T M(10)约束等式是所有iR 都应满足正交性地6 个等式 . 优化模型地目地是为了求精1k 、2k 、3k 、4k 、iiFRT

26、、, i jM、,?i jm 同模型式 (9) 定义.1k 、2k 、3k 、4k 为本文考虑地畸变系数，由于其一般很小，因此其初始值可置为0. 本文模型地初值为第 1 步求取地各参数值以及初始地畸变系数.1234,iii jm k kk kF R T M是由, i jM代入式 (4) 先求得,uuxy 再由实际成像模型式(7) 求出地 . 由于模型参数太多，用一般地梯度法难以保证收敛，因此本文采用一种基于牛顿迭代法求解.精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 10 页，共 19 页五、 确定两部相机相对位置地数学模型和方法5.1 模型假设

27、（1）假设两部相机内部参数固定不变.（2）假设噪声对相机成像无影响.5.2 模型建立设右相机坐标系到左相机坐标系地变换关系为：123456789abbabababbaxababbababababbayabbabababbazxxrrrxtyRyTrrrytzzrrrzt(11)空间任一点P在左右象面象点分别为(,),(,),aaabbbP XYP XY转换成各自相应地相机坐标分别为(,)aaaaP XYf，(,)bbbbP X Yf；又已知两相机透视中心分别为aO 、bO ，这样由这4点便构成一个平面，即外极平面.将这4点坐标都转换到左相机坐标系下，由四点共面条件可得：123456789()(

28、)0()aabbbabbabbaxaabbbabbabbayaabbbabbabbazXrXrYrftYrXrYrftfrXrYrft(12)令; ; aaaaaaaaaabbbbbbbbbbXXfxzYYfyzXXfxzYYfyz(13)将 n个空间点地左右象点地象面坐标值代入式(12)可得0BA(14)其中1111111111111abababaabalbbanbnanbnanbnananbnanbnbnXXXYXX YY YYXYBXXXYXXXXYYXY123456789TAaaaaaaaaa147258369471;582693714825936;babazbabaybabazbab

29、aybabazbabaybabaxbabazbabaxbabazbabaxbabazbabaybabaxbabaybabaxbabaybabaxartrtartrtartrtartrtartrtartrtartrtartrtartrt由于旋转矩阵baR 为正交矩阵，所以式 (14) 是一个含有6个未知数地非线性齐次方程组 .基于平行线消隐点理论地“ 两步法 ” 首先利用多于4组平行线方向矢量地左精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 11 页，共 19 页右象面消隐点坐标，线性求解一个纯旋转关系地方程组，得到baR 矩阵；然后将baR 矩阵

30、代入式 (14) ，再利用多于两个空间点地匹配坐标值及左右相机坐标系之间地先验知识，可以线性求解baT 矩阵7.5.3 模型求解1）消隐点原理及baR 矩阵地求解设在左相机坐标系下单位方向矢量为( , , )as a b c 且过点 (,)oooaaaxyz地空间直线地参数方程为：000aaaaaaxxatyybtzzct(15)设同一直线在右相机坐标系下地参数方程为：000bbbbbbxxa tyyb tzzc t(16)其中(,)bs a b c 为此直线在右相机坐标系下地单位方向矢量. 根据左右相机坐标系之间地转换关系可知：abaabRbcc(17)空间平行直线只要不与相机象面平行( 即

31、0,0cc) ，经过透视变换后在象面上应相交于一点，即“ 消隐点 ” .对式(15) 和式(16) 进行透视变换可得：00000000;aaaaaaaaaaaabbbbbbbbbbbbXxxatYyybtfzzctfzzctXxxa tYyxb tfzzc tfzzc t(18)消隐点是空间直线 t( 或t) 时在像平面上地透视投影，由此可得：; Y; YaaaaaabbbbbbXXfa cYfb cXXfa cYfb c(19)其中(,aaXY ) 和(,bbXY ) 分别为左右像面消隐点坐标，af 和bf 分别为左右相机有效焦距 .精选学习资料 - - - - - - - - - 名师归纳

32、总结 - - - - - - -第 12 页，共 19 页设 n条空间平行直线在左像面透视投影为：1111aaaaananananYk XpYk Xp(20)这组像面直线相交于一消隐点(11,aaXY ) ，同样可求得左像面内地其它消隐点坐标22(,),(,),aaaiaiXYXY.同理可得右像面内各消隐点坐标11(,),(,),bbbibiXYXY. 由式(17)(19) 可得：123789456789babibabibaaibabibabibabaibbabibaaibabibabibarXrYrXrXrYrrXrYrYrXrYr(21)将多于4对消隐点地坐标值代入式(21) ，并结合ba

33、R 矩阵地正交约束条件，解此未知数线性方程组可得旋转矩阵baR 地初值，然后再利用迭代法可更精确地求取baR 矩阵.2）baT 矩阵地求解将baR 矩阵代入式 (14) ，给定多于两个空间点地左右像面坐标值，并利用两台摄相机之间地先验知识来判断T地符号，就可线性求得T矩阵.由双目视觉传感器模型 8 可以得到以下方程：123789456789babbabbabbaxababbabbabbazbabbabbabbayababbabbabbazrXrYrztXrXrYrztrXrYrztYrXrYrzt(22)由此可得789123789456 =baxabazbababbabbababbabbaba

34、yabazababbabbababbabbabbbbbbtX tzXrXrYrrXrYrtY tYrXrYrrXrYrxz Xyz Y(23)精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 13 页，共 19 页将baTkT ，连同baR 矩阵地值一起代入式 (23) ，计算两个空间点之间地距离，然后将式 (23) 算出地含参数k地测量值和实际距离值进行比较，便可计算出k值，从而求得baT 矩阵.六、 模型和算法地优缺点分析及改进方向由于数码照相机模型地参数太多，而且是一个复杂地非线性方程，因此一步求解不仅繁冗，而且稳定性差；本文所采用地两步法确定

35、像坐标和两相机相对位值地方法，不仅考虑了径向畸变，还考虑了切向畸变因素，求得较为准确地像坐标和相机相对位置，精度高，稳定性好，解地收敛速度快.由于考虑了相机地径向畸变和切向畸变即非线性畸变，计算程序略显复杂，计算量稍大 .由于 BP神经网络在逼近能力，收敛速度等方面都优于其它地非线性优化算法，很适宜于相机地标定. 将神经网络技术用于相机地标定，不需要预先确定摄像机具体内外参数及有关模型地先验知识，利用BP神经网络对复杂非线性映射关系强大地逼近能力来对视觉系统中地摄像机进行隐式标定，算法标定精度较高，可靠性好，实用性较强，是相机定位地一种趋势.参考文献1 刘红宁 . 基于双目立体视觉拟人机器人定

36、位关键技术地研究D. 河北工业大学,2006. 2 马颂德 ,张正友 . 计算机视觉 M. 北京:科学技术出版社 ,1998. 3 毛剑飞 ,邹细勇 ,诸静. 改进地平面模板两步法标定摄像机J. 中国图象图形学报,2004,9(7):847852.4 张磊,刘小力 ,罗翔,田梦倩 ,史金飞 . 一种考虑畸变地摄像机线性标定简化方法J. 工业仪表与自动化装置 ,2007(1):37.5 陈杰. Matlab宝典M. 北京:电子工业出版社 ,2007. 6 李海滨 , 单文军 , 刘彬 . 双目立体视觉测距系统误差模型地研究J. 光学技术,2006(01) :1316.7 张健新 ,段发阶 ,叶声

37、华 . ” 两步法 “ 求取双目视觉传感器中摄象机位置关系J. 光电工程 ,1998,25(5):3741.8 李雪鹃 , 蔡怀宇 , 黄战华 . 基于双目视差地测距系统设计J. 电子测量技术,2007(08) :3134.精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 14 页，共 19 页附 录基于 Canny算子边界检测源程序1 %基于Canny算子求边界坐标%参数name 为图片名function x,y=tp(name) I=imread(name)。%读入图片I=double(I)。I=I(1:768,1:1024)。G=edge(I,

38、Canny)。 %基于Canny算子求边界imshow(G)。I=uint8(I) 。G=uint8(G)。subplot(1,2,1) imshow(I)。%输出图像subplot(1,2,2) imshow(G)。%输出图像x= 。y=。for i=100:600 %求边界坐标 for j=200:700 if G(i,j)=0 ） y=y,i。 x=x,j 。 end end end基于简化模型求取像坐标初值源程序2 %基于简化模型求解像圆心%参数x，y为物地圆心坐标 u为物距function x0,y0=jiandan(x,y,u) x0= 。y0=。v=1577。%给像距赋初值h1=

39、24。f=ff(h1,h2,v) %调用求解焦距地函数， n=length(u)。for j=1:n。 k=f/(u(j)+v)*3.78 。%求比例系数 a=x(j)*k 。%求像坐标并保存 b=y(j)*k 。 x0=x0,a。精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 15 页，共 19 页 y0=y0,b。end 两步法求解像坐标源程序3 %两步法求解像坐标%读入数据并整理成可处理地形式%M为靶标数据%m为靶标地像数据%求解焦距function f=ff(h1,h2,v)n=length(h2)。 for k=1:n u(k)=h1*v

40、/h2(k) 。%求每一个像高对应地物距endu=u*3.78。for k=1:n F(k)=1/(1/v+1/u(k) 。%求每一个物距对应地焦距endf=mean(F)。%求所以焦距地算术平平均clear all。close all。M=imread(tu.bmp)。%读入图片m1=imread(tu1.bmp)。m2=imread(tu2.bmp)。m3=imread(tu3.bmp)。m4=imread(tu4.bmp)。m5=imread(tu5.bmp)。MM=M(:,1:2) 。M(:,3:4)。M(:,5:6)。M(:,7:8) 。%提取数据mm1=m1(:,1:2)。m1(:

41、,3:4)。m1(:,5:6)。m1(:,7:8)。mm2=m2(:,1:2)。m2(:,3:4)。m2(:,5:6)。m2(:,7:8)。mm3=m3(:,1:2)。m3(:,3:4)。m3(:,5:6)。m3(:,7:8)。mm4=m4(:,1:2)。m4(:,3:4)。m4(:,5:6)。m4(:,7:8)。mm5=m5(:,1:2)。m5(:,3:4)。m5(:,5:6)。m5(:,7:8)。M=M(:,1:2) 。 M(:,3:4) 。 M(:,5:6) 。 M(:,7:8) 。m1=m1(:,1:2) 。 m1(:,3:4) 。 m1(:,5:6) 。 m1(:,7:8)。m2=m

42、2(:,1:2) 。 m2(:,3:4) 。 m2(:,5:6) 。 m2(:,7:8)。m3=m3(:,1:2) 。 m3(:,3:4) 。 m3(:,5:6) 。 m3(:,7:8)。m4=m4(:,1:2) 。 m4(:,3:4) 。 m4(:,5:6) 。 m4(:,7:8)。m5=m5(:,1:2) 。 m5(:,3:4) 。 m5(:,5:6) 。 m5(:,7:8)。M=M 。m(:,:,1)=m1。m(:,:,2)=m2。m(:,:,3)=m3。m(:,:,4)=m4。m(:,:,5)=m5。%M为靶标数据精选学习资料 - - - - - - - - - 名师归纳总结 - -

43、- - - - -第 16 页，共 19 页%m为靶标地像数据function liangbufa(M,m) %M=X,Y 转化为 M=X,Y ,1 %m=u,v 转化为m=u,v,1 rows,npts=size(M)。 matrixone=ones(1,npts) 。 M=M 。matrixone。 num=size(m,3)。 for i=1:num m(3,:,i)=matrixone。 end H=jiandan(x,y,u)%调用简单模型求解初值 V= 。 for flag=1:num %迭代v12(:,:,flag)=H(1,1,flag)*H(2,1,flag),H(1,1,f

44、lag)*H(2,2,flag)+H(1,2,flag)*H(2,1,flag),H(1,2,flag)*H(2,2,flag),H(1,3,flag)*H(2,1,flag)+H(1,1,flag)*H(2,3,flag),H(1,3,flag)*H(2,2,flag)+H(1,2,flag)*H(2,3,flag),H(1,3,flag)*H(2,3,flag)。v11(:,:,flag)=H(1,1,flag)*H(1,1,flag),H(1,1,flag)*H(1,2,flag)+H(1,2,flag)*H(1,1,flag),H(1,2,flag)*H(1,2,flag),H(1,3

45、,flag)*H(1,1,flag)+H(1,1,flag)*H(1,3,flag),H(1,3,flag)*H(1,2,flag)+H(1,2,flag)*H(1,3,flag), H(1,3,flag)*H(1,3,flag)。v22(:,:,flag)=H(2,1,flag)*H(2,1,flag),H(2,1,flag)*H(2,2,flag)+H(2,2,flag)*H(2,1,flag),H(2,2,flag)*H(2,2,flag),H(2,3,flag)*H(2,1,flag)+H(2,1,flag)*H(2,3,flag),H(2,3,flag)*H(2,2,flag)+H(

46、2,2,flag)*H(2,3,flag), H(2,3,flag)*H(2,3,flag)。V=V 。v12(:,:,flag)。v11(:,:,flag)-v22(:,:,flag) 。end k=V*V 。%将Vb=0转化为 kb=V*V*b=0, 这样k就是一个对称正定矩阵，u,v,d=svd(k)。e,d2=eig(k)。b=d(:,6)。 v0=(b(2)*b(4)-b(1)*b(5)/(b(1)*b(3)-b(2)2) 。 s=b(6)-(b(4)2+v0*(b(2)*b(4)-b(1)*b(5)/b(1) 。 alpha_u=sqrt(s/b(1)。 alpha_v=sqrt(

47、s*b(1)/(b(1)*b(3)-b(2)2) 。 skewness=-b(2)*alpha_u*alpha_u*alpha_v/s 。 u0=skewness*v0/alpha_u-b(4)*alpha_u*alpha_u/s 。 A=alpha_u skewness u00 alpha_v v0 0 0 1。 D= 。%求解所有外部参数 d=。 Rm=。 for flag=1:num s=(1/norm(inv(A)*H(1,:,flag)+1/norm(inv(A)*H(2,:,flag)/2。 rl1=s*inv(A)*H(1,:,flag) 。精选学习资料 - - - - - -

48、- - - 名师归纳总结 - - - - - - -第 17 页，共 19 页 rl2=s*inv(A)*H(2,:,flag) 。 rl3=cross(rl1,rl2)。 RL=rl1,rl2,rl3 。 U,S,V = svd(RL)。 RL=U*V 。 TL=s*inv(A)*H(3,:,flag) 。 TL是位移矩阵 t(外参) RT=rl1,rl2,TL 。%H=Ar1 r2 t XY=RT*M 。% M是model plane 点地坐标UV=A*XY 。% sm=AR tM,UV 是等式地右边UV=UV(1,:)./UV(3,:) 。 UV(2,:)./UV(3,:) 。 UV(3

49、,:)./UV(3,:) 。 XY=XY(1,:)./XY(3,:)。 XY(2,:)./XY(3,:) 。 XY(3,:)./XY(3,:) 。 end k=inv(D*D)*D*d 。 para=Rm,k(1),k(2),alpha_u,skewness,u0,alpha_v,v0 。 options = optimset(LargeScale,off,LevenbergMarquardt,on) 。% 用最小二乘法求解 x,resnorm,residual,exitflag,output = lsqnonlin( simon_HHH, para, ,options, m, M)。 k1=

50、x(num*6+1)% 显示结果k2=x(num*6+2) A=x(num*6+3) x(num*6+4) x(num*6+5) 。 0 x(num*6+6) x(num*6+7) 。0,0,1基于 Canny算子边界检测确定圆心源程序4 %基于Canny算子求出地边界坐标来拟合圆心%参数x，y为边界坐标function x0,y0=yuanxin(x,y) x1=(min(x)+max(x)/2 。%估计圆心地坐标y1=(min(y)+max(y)/2 。X= 。Y= 。for a=1:length(x)-2 %任意取三个边界点确定一个圆心地坐标 for b=a+1:length(x)-1

51、for c=b+1:length(x) x0=(x(b)2-x(a)2+y(b)2-y(a)2)*(y(b)-y(a)/2-(x(a)2-x(c)2+y(a)2-y(c)2)*(y(a)-y(b)/2)/(x(a)-x(b)*(y(a)-y(c)-(x(a)-x(c)*(y(a)-y(b)。y0=(x(a)2-x(c)2+y(a)2-y(c)2)*(x(a)-x(b)/2-(x(a)2-x(b)2+y(a)2-y(b)2)*(x(a)-x(c)/2)/. (x(a)-x(b)*(y(a)-y(c)-(x(a)-x(c)*(y(a)-y(b) 。 if (x1-x0)2+(y1-y0)2100 X=X,x0 。 Y=Y,y0。 end end 精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 18 页，共 19 页 end end x0=mean(X)。%求满足条件地坐标均值 y0=mean(Y)。精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 19 页，共 19 页