您的当前位置:首页正文

改进的RANSAC匹配点提纯算法

2021-12-13 来源:好走旅游网
第45卷第6期 西 安 建 筑 科技 大 学 学 报(自然科学版) J.Xi ai1 Univ.of Arch.&Tech.(Natural Science Edition) Vo1.45 No.6 Dec.2O13 2013年12月 改进的RANSAC匹配点提纯算法 介 军 ,李智杰 ,姚 鹏 (1.西安建筑科技大学信息与控制工程学院,陕西西安710055; 2.长庆油田分公司机械制造总厂,陕西西安710201) 摘 要:针对在图像匹配中,随机抽样一致性(RANSAC)算法对匹配点提纯存在计算量大、效率低的问题,采 用将基本矩阵作为模型参数估计对象的方法,对RANSAC匹配点提纯算法进行了改进.在改进的算法中,运 用Bucket分割技术抽取粗匹配点对,进行两幅图像的检测角点和粗匹配,利用视差梯度对匹配点样本预检 验.实验结果表明,此方法在保证较高精度和鲁棒性的情况下,运算量大幅度减少,提高了图像匹配的速度. 关键词:角点检测;图像匹配;视差梯度;RANSAC 中图分类号:TP 391.4l3 文献标志码:A 文章编号:1006—7930(2013)06—0896—06 图像匹配是当前研究的热点问题之一,也是计算机视觉领域中相关理论的基础.在图像匹配的众多 方法中,目前应用最广泛、研究最多的方法是基于特征的匹配[1 方法,它根据原图像中特征趋于稳定、数 量不是太多等具有明显标记的的一些点、线、边缘等进行图像的匹配,使得在处理过程中能尽可能的减 小运算数据量,具有较高的运算速度和效率,且对图像灰度变化具有很强的适应性.美中不足的是,基于 特征的匹配仍需要在不同图像中进行大量的遍历性匹配运算.这样会使得匹配过程中存在计算量很大, 且精度不高的问题. 本文采用Harris算法提取不同角度的左右两幅图像的角点,利用灰度归一化互相关对角点进行粗 匹配,然后基于改进的随机抽样一致性(RANdom SAmple Consensus,RANSAC)算法进行精匹配.在 改进的算法中,使用反映两幅图像像素坐标关系的基本矩阵作为模型参数估计对象,采用Bucket平面 分割技术分块抽取粗匹配点对避免分布过于集中,利用视差梯度对参与模型计算的粗匹配点样本预检 验以减少不必要的模型验证,在确保精度的前提下,大大减少了计算量,提高了运算速度. 1 角点提取和粗匹配 角点是指在至少两个方向上图像灰度变化均较大的点.在实际的图像中,角点可以是两条直线的交 汇处或者图像轮廓的拐角点、甚至是线段的末端等都可以称之为角点.本文采用Harris算法检测角点, 采用使用最多的基于窗口的归一化互相关函数对角点粗匹配. 1.1角点提取 Harris算法l_3 定义的基础是如下定义的图像灰度强度的二阶导数矩阵: a I 8x a I 8xOv M:= a0I OxOy a 工 a、1 其中:J为图像灰度强度;5C与 为图像像素点的横坐标和纵坐标.Harris角点的响应函数为: H=== 1 2——k( l+ 2) (2) 其中 , 为矩阵为式(1)中M的两个特征值,根据Harris建议,是的取值为0.04,特征点就是此 收稿日期:2013-01—14 修改稿日期:2013—1卜25 基金项目:国家自然科学基金资助项目(50878176) 作者简介:介军(1971一),男,陕西彬县人,工程师,硕士,主要从事计算机应用和教学研究 第6期 介军等:改进的RANSAC匹配点提纯算法 897 函数的局部最大值. 1.2粗匹配 角点粗匹配是利用角点附近的灰度信息,采用归一化互相关方法,对于左图像中每一个角点,建立 一个局部匹配的准则.角点粗匹配计算其所有角点的归一化互相关系数. 归一化互相关函数定义如下: k Z ∑∑EI(u+ , + )一 ]×[ ( + , + )一了7 ] (3) 公式中,k、f是两个互相关窗口的宽度,一般取k—l, 表示的是灰度平均值(以点(“, )为中 心). 计算结果为相关系数C,C取值越大,则证明角点 表1 角点互相关匹配矩阵 相关程度越大,这里C的取值为:一1<C<1. Tab.1 Matching matrix of 本文通过计算角点互相关匹配矩阵,来获取匹配 点.查找原则为:遍历矩阵中的每一行数据,找出最大 值,然后再找出该大值所对应的一列,如果此最大值也 是这一列的最大值,那么就可以确定这是一对匹配点. 互相关矩阵如表1所示,其中左侧为左图像中角点 编号,右侧为右图像中角点编号. 2 图像的精匹配 精匹配的目的就是消除粗匹配中的伪匹配点,目前比较有效而且通用的方法就是由Fischler和 Bolles提出的RANSAC算法嘲,它对超过5O%的伪匹配点仍然有效,是一种模型参数估计算法.但是 该算法计算量大并且模型参数不稳定等因素导致了较低的匹配率,基于此问题,本文采取基于视差梯度 和Bucket平面分割技术的RANSAC算法来估计基本矩阵,通过基本矩阵剔除伪匹配点,达到图像的 精匹配. 2.1基本矩阵及其估计 极几何的代数形式通常用基本矩阵表示,基本矩阵描述了一幅图像上的像素点与另一幅图像上对 应极线之间关系,以及摄像机的内参数信息[6],对极关系描述形式如式(4)所示: ,rF f一0 (4) 其中,F是的3×3的基本矩阵,m 和m 是左右两幅图像的一对匹配点,m 一(z ,Y ,1),m 一 (z ,Y ,1)是像素的齐次坐标形式来展示. 基本矩阵的估计方法有多种,典型的也是最常用的有7点法、8点法和基于几何误差的非线性优化 算法等等 引,一般的方法就是以下的8点法,对于 和m .,按照式(4)可以展开为 73 z -厂1l+z fYff12+z 厂13+y'ixlf21+ {f22+3,t 厂23+z f31+Yi厂32+f33—0 (5) 用矢量,表示由F的元素组成并按行优先顺序排列的9维矢量,式(5)可以表示为一个矢量的内 积 (z f , Y ,z , z , t Y ,Y f, ,Y ,1)f=0 (6) 如果有n组点对,可以得到如下的线性方程组: Pl1 1;711- z 1Y z 1 Yl 1  Y t1z1z1  Y Ya1Y1  Y 11 z1 Yl  1] ] Af—l i i i ; ; i i ! ;I f一0 (7) f z t nz"ttz  nX z"z  Y nX n Y Y Y n z Y 1z” nY  .jf 根据矩阵描述的数据可以看出,假如矩阵A的秩等于8,那么就可以进行线性求解该方程组,并且 存在惟一的解.然而在实际情况中,点的坐标含有一定的噪声,造成矩阵的秩在有些情况下不等于8.对 898 西安建筑科技大学学报(自然科学版) 第45卷 矩阵A进行最小二乘解,然后再进行分解,即A—UDV ,V的最后一列矢量,在lJ f『l一1且{l AF l l的值最小的条件下,即为最小二乘解. 2.2 RANSAC算法 RANSAC算法是一种简单而巧妙的算法.RANSAC算法的计算量由两部分组成口]:(1)抽样选择 和模型估计需要的时间;(2)模型参数检验需要的时间. RANSAC算法可总结如下: (1)求出满足要求的最小抽样数D:根据计算模型参数需要的最小数据量d、数据错误率£和置信 概率P,利用式(8), 1一(1一(1一£) )。一P (8) (2)计算参数化模型:从原始数据S中随机选取一定量的数据点样本,一般采用数据量为d作为标准. (3)计算当前模型的支撑点集:使用原始数据作为参考,对模型数据进行检验计算,得到模型的支撑 点的数量. (4)返回第(2)、(3)步,重复进行计算D次,找出最优模型. (5)根据第四步得到的最优模型,获得其所对应得支撑点集合,得到最终模型参数. 2.3 改进的RANSAC算法 从上文RANSAC算法描述来看,算法还有可以改进优化之处.可以从抽样选择以及模型参数验证 这两个方面进行优化改进,本文利用平面分割技术和视差梯度分别解决匹配点在空间上分布不均的问 题和减少模型参数验证次数,这样获得的模型参数更具有代表性. 从相关参考文献中视差梯度的定义_8 得出结论:在一幅图像中,如果P,q和P ,q 分别是两幅图像 中两个相邻角点,若角点P,q分别与角点P ,q 匹配,并且是相容的,那么通过计算公式得出的视差梯 度G 应不大于2;否则,两个角点就被认为是不匹配的.视差梯度的公式为: G —z牦每寿 其中,II P lI表示向量P的模,( ,p)和(q ,q)分别是图像对应角点的坐标向量. 法,从不同方面进行描述可以有效地弥补灰度相关的不足之处,并做出改善. ㈤ 视差梯度区别于灰度相关的一般性质,此种方法从另外一个角度来描述角点之间的匹配模式与方 在使用RANSAC方法进行模型参数估计的过程中,如果在空间中使用的随机算法随机抽取点集 有较大的聚合性,经过计算得出的基础矩阵F将会有很差的稳定性,最后得出的结果也会有较大的误 差.根据经验,可以将图像平面进行网格划分,这里本文采用Buckets平面分割技术嗍|1o 形成若干个 Buckets来减少误差,用归于每个Bucket中的所有点的数目占总体点的数目比率的大小作为权值,采 用轮盘赌的方式,不重复的选择点集,能有效的解决点集分布过分集中的问题,并且增强基础矩阵的准 确性和稳定性. 改进的RANSAC算法描述: (1)根据图像所有角点性质,利用归一化互相关方法建立N个粗匹配角点对. (2)根据角点在图像上的位置将匹配点对等大小的划分成m个不同的Bucket,计算权值,统计并计 算所有Bucket中的角点数N ,i一1,2,3,…,m,这样每个Bucket的权值为P :N /N. (3)根据式(8)计算最小随机抽样数D. (4)根据每个Bucket的权值P ,采用轮盘赌的选择方式抽取一个没有重复匹配点对的大小为8的 样本. (5)设定阈值maxTD,如果计算所得样本中视差梯度总和,若大于maxTD,转回第(4)步骤进行重 新抽样;否则,继续执行下一步算法. (6)计算抽样基础矩阵 ,利用8点算法. (7)将所有的初始匹配点对代入公式(5),使得Y =m'iTF ITt ,如果y 的值在区间[口, 中,记录统 计点的数目N ,以及所对应的角点,重复第(4)~(7)步骤,一直到完成D次随机抽样为止. 900 西安建筑科技大学学报(自然科学版) 第45卷 与原始匹配点集的抽样数量大大减少,匹配的时间减少到1.812s.上文中图2显示了使用改进后的算 法所得到的结果,图像中存在的粗匹配点用“o”标出,可以看到在“o”中标“口”符号,表示的是算法剔 除的伪匹配点,将会被处理掉. 表2改进的算法与原算法时间的比较 Tab.2 Comparison of time consamed with different algorithms 为了检验本文算法的通用性,选用3组具有代表性的图像, 将本文的算法与原RANSAC算法进行了对比,由表3可以看 出,改进的算法在保持了精度的同时,提高了计算速度.同时,在 最小抽样数一定的前提下,通过不同分辨率,不同场景的图片进 行匹配,得出了在两种算法中粗匹配点数与匹配时间的关系,由 图3可知,随着粗匹配点数的增加,两种算法的匹配时间都在增 加,但是匹配时间差则越来越大,在粗匹配点较多时,改进算法 的优势更为明显. 4结论 三 针对RANSAC方法进行图像角点精匹配时速度慢、实时points and the match time 性差的问题,本文改进了RANSAC算法.在改进的算法中,以 基本矩阵作为模型参数估计对象,采用Bucket图像分割技术分块抽取粗匹配点对,利用视差梯度对抽 取的匹配点对预检验.实验证明,此方法在保证较高精度和鲁棒性的情况下,运算量大幅度减少,提高图 像匹配的速度. 参考文献 References [1]ZHU Q,WU B,XU Z.Seed point selection method for triangle con strained image matching propagation[J]. IEEE Geoscience and Remote Sensing Letters,2006,3(2):207—211. [2]曲天伟,安波,陈桂兰.改进的RANSAC算法在图像配准中的应用EJ].计算机应用,2010,30(7):1849—1851. QU Tian-wei,AN Bo,CHEN Gui—lan.Application of improved RANSAC algorithm to image registration[J]. Journal of Computer Applications,2010,30(7):1849—1851. [3]HARRIS C,STEPHENS M.A combined corner and edge detection[J].Image Vision Computing,1998,15(6): 121-127. [4]徐玮.一种基于角点匹配的视图合成方法[J-I.系统仿真学报,2007,19(14):3263—3265. XU Wei.Corner Matching—based Approach of View Synthesis[J].Journal of System Simulation,2007,1 9(14): 3263—3265. [5]FISCHLER M A,BOLLES R C.Random sample consensus:A paradigm for model fitting with applications to im— age analysis and automated cartography[J].CACM,1981,24(6):381—395. [63 HARTLEY R,ZISSERMAN A.Multiple View Geometry in Computer Vision[M].2nd ed.Cambridge Universi— 第6期 介 军等:改进的RANSAC匹配点提纯算法 901 ty,2003. [7]陈付幸,王润生.基于预检验的快速随机抽样一致性算法[J].软件学报,2005,16(8):1431—1437. CHEN Fu—xing,WANG Run-sheng.Fast RANSAC with Preview Model Parameters Evaluation[J].Journal of Software,2005,16(8):143卜1437. [8]马颂德,张正友.计算机视觉EM].北京:科学出版社,1998:82—83. MA Song—de,ZHANG Zheng-you.Computer vision[M].Beijing:Science Press,1998:82—83. [9]CHOUKROUN A,CHARVILLAT V.Bucketing techniques in robust regression for computer vision[C]∥In Pro— ceedings of SCIA 2003 Lecture Notes in Computer Science,Halmsted Sweden,2003,2749:609—616. [1O]黄以君,刘伟军.基于LQS的基本矩阵计算方法[J].中国图象图形学报,2009,14(10):2069—2073. HUANG Yi-jun,LIU Wei-jun.A Method for Fundamental Matrix Estimation Using LQS[J].Journal of Image and Graphics,,2009,14(10):2069—2073. Improved RANSAC algorithm of matched points purifying JIEJun ,LI Zhi-jie ,YAO Peng (1.School of Information and Control Engineering,Xi an Univ.of Arch.&Tech.,Xi an 710055,China; 2.Changqing Oilfield Company,Machine Manufacture Plant,Xi an 710201,China) Abstract:Improvement has been made for RANSAC alogorithm of matched points purifying in the image matching process by using fundamental matrix as object of model parameter in finding a solution to the problem with the large amount of calculation and the low efficienty that eonerge.By extracting rough matched points in sub—block generate from Bucketing techniques,algorithm test two images with corner detection and rough matching has also been improved.Moreover,it pretests the sample of matching points using disparity constraint.The experiment shows that this algorithm reduces the a— mount of computation largely,improves the speed of image matching and keeps high precision and robust. Key words:corner detection;image matching;disparity constraint;RANSAC Biography:JIE Jun,Engineer,Master,Xi an 710055,P.R.China,Tel:0086—29—82202537,E—mail:mouzi一1997@163.corn (上接第884页) Planning and design research of village emplacement based on‘‘original states’’residential conception ZHAO Hong—bin ’ (1.School of Arch.,Xi an Univ.of Arch.&Tech.,Xi an 710055,China; 2.State Key Laboratory of Architecture and Techology in West China(XAUAT),Xi an 710055,China) Abstract:The paper summarized a series of problem which arise in land—lost farmer village emplacement resettlement area. The research starts with the view of farmer whose 1and was commandeered for establishing fl village emplacement planning idea of“Original States”,and emphasize the protection of maintanance the traditional lifestyle of the people.Investigation is made for the basic settlement spatial form of a suitable emplacement community.The paper will explore a new planning concept and new design strategy of emplacement community. Key words:newly urbanization;land-lost farmer;emplacement planning Biography:ZHAO Hong—bin,Lecturer,Candidate for Ph.D.,Xi an 710055,P.R.China,Tel:0086—13991139113,E-mail:550998658 @q q.eom 

因篇幅问题不能全部显示,请点此查看更多更全内容