Sift图像快速匹配的尺寸预处理方法
【专利摘要】一种SIFT图像快速匹配的尺寸预处理方法，其特征是当有M对缩放比例相近的大尺寸图像和小尺寸图像进行匹配时，以一对图像中景物相同部分分辨率尽可能相等为原则，选择一对原始图像做预匹配，确定大尺寸图像的缩小比例。之后，或按照该缩小比例将后续M-1幅大尺寸图像在匹配前缩小，或利用实时更新的缩小比例将下一幅大尺寸图像在匹配前缩小。除预匹配之外，其余图像匹配的速度随着大尺寸图像的缩小而提升，而匹配的性能并未下降。本发明便于与其他SIFT改进算法相结合，以得到更快的匹配速度。
【专利说明】SIFT图像快速匹配的尺寸预处理方法
【技术领域】
[0001]本发明涉及计算机视觉领域，特别是涉及一种图像匹配的快速方法。
【背景技术】
[0002]图像匹配在视觉导航、目标定位、识别与跟踪、遥感图像处理、图像检索、立体视觉测距和三维重建等应用上不可或缺，故已成为计算机视觉领域的研究热点之一。SIFT (scale invariant feature transform)是一种著名的图像匹配算法，由 D.G.Lowe 于1999年提出，2004年完善。SIFT算法的突出优点是具有尺度、旋转、光照等不变性。换句话说，SIFT能够实现具有不同尺度、不同旋转参数、不同光照强度的异源图像匹配。
[0003]然而，运行速度慢是SIFT算法的一个公认问题。为此，许多学者给出了各具特色的改进方案，但鲜有涉及预处理环节。众所周知，在能够匹配的条件下，两幅待匹配图像的尺寸越小，图像匹配的速度越快。所以，在匹配前对图像做尺寸缩小预处理是提高匹配速度的一种途径。然而，在实际应用中，如果同时缩小两幅待匹配图像，往往存在着降低匹配性能的风险，使得原本能够匹配的两幅图像失去匹配关系。尤其当两幅待匹配图像的尺寸不同时，小尺寸图像的缩小意味着特征点的减少，进而导致匹配点的减少。当匹配点数少于一定阈值时，图像匹配失败。
【发明内容】
[0004]本发明的目的在于，针对存在对应景物但对应部分分辨率不同(即尺寸不同)的两幅待匹配图像，给出一种大尺寸图像的预缩小方法，一方面使图像匹配时间得以减少，另一方面使图像匹配性能得以保持。
[0005]本发明的技术方案是，当有M对具有相近缩放比例的大尺寸图像和小尺寸图像进行匹配时，在大批量进行匹配之前，以两幅图像中景物相同部分分辨率尽可能相等为原则，通过选择一对原始图像做预匹配，确定大尺寸图像的缩小比例。之后，按照该缩小比例将后续M-1幅大尺寸图像在匹配前做缩小预处理。具体步骤如下:
[0006]第一步，基于SIFT算法进行一组原始尺寸的大尺寸图像和小尺寸图像匹配，即图像预匹配。如果M幅大尺寸图像是序列图像，则选取第一对原始图像进行预匹配。
[0007]第二步，从所有匹配的特征点对中随机选取两对，如图1中的A-A'和B-B'，计算大尺寸图像中特征点A和B之间的距离I,以及小尺寸图像中对应特征点A'和B'之间的距离Γ，得到第一个缩小比例:
[0008]rl = I' /I
[0009]第三步，重复第二步N-1次，得到N-1个缩小比例r2,…，rN。N可取3?10。
[0010]第四步，求取N个缩小比例的平均值，得到最终的大尺寸图像缩小比例:
[0011]r = (rl+r2+…+rN)/N
[0012]第五步，在后续的M-1次图像匹配任务中，先将每幅大尺寸图像缩小到原来的r倍(图像的宽和高均缩小到原来的r倍)，然后再和对应的小尺寸图像进行匹配。
[0013]在上述尺寸预处理中，第一次预匹配所获得的缩小比例r与后续某些图像对的实际缩小比例之间可能存在较大偏差。如果M对缩放比例相近的大尺寸图像和小尺寸图像是序列图像，可采用缩小比例跟踪校正策略消除该偏差，即在第k+Ι次匹配时，采用第k次的匹配结果更新缩小比例r, k = 1，...，Μ-1,具体如下:
[0014]步骤1:令k= l,r = I,基于SIFT算法进行第一对原始尺寸的大尺寸图像和小尺寸图像预匹配。
[0015]步骤2:从第k对图像的所有匹配特征点对中随机选取两对计算缩放比例，k = I时，缩放比例=小尺寸图像两特征点间距离/大尺寸图像两特征点间距离；k>l时，缩放比例=小尺寸图像两特征点间距离/缩小的大尺寸图像两特征点间距离，重复该做法N次(N=3?10),共得至Ij N个缩放比例rl,…，rN,求取其平均值r' = (rl+r2+…+rN)/N，进而得至Li更新的大尺寸图像缩小比例r = rXr'。
[0016]步骤3:令 k = k+l。
[0017]步骤4:将第k幅大尺寸图像缩小到原来的r倍，然后再和第k幅小尺寸图像进行匹配。
[0018]步骤5:若k = M，停止匹配过程。否则，转到第二步。
[0019]因为缩放比例的计算时间较之图像匹配时间可以忽略不计，所以缩小比例的跟踪校正并不影响匹配速度的提升。
[0020]本发明所达到的效果和益处是，除了第一次图像预匹配的速度保持不变之外，其余图像匹配的速度都将随着大尺寸图像的缩小而提升。同时，图像匹配的性能并未下降。其原因有二，一是缩小后的大尺寸图像与小尺寸图像在匹配内容上具有近似相等的分辨率，在尺度空间上的匹配对应关系也变得更加明确；二是本专利中大尺寸图像的高斯模糊程度较小，图像的细节特征保留得更好。以一幅大尺寸航拍图(1024X1024)和一幅小尺寸Google Earth图(640X640)为例，见图2。航拍图缩小前(r = 1)，匹配率=1.31% (总特征点数=7883，匹配点数=103)，匹配时间为11.8s(计算机CPU为E6550，主频2.34GHz,内存3.25GB)。航拍图缩小后(r = 0.23)，匹配率=15.99% (总特征点数=613，匹配点数=98)，匹配时间为4.6s。可见，航拍图缩小后，匹配速度约提升为原来的3倍，但匹配点数基本不变，匹配率上升了 12倍。需要说明的是，本发明处于图像预处理环节，便于与其他的SIFT改进算法相结合，以得到更加快速的匹配速度。
【专利附图】
【附图说明】
[0021]图1是基于图像预匹配确定大尺寸图像缩小比例的示意图。
[0022]图2是一幅大尺寸航拍图和一幅小尺寸Google Earth图的预匹配示意图。
【具体实施方式】
[0023]下面结合技术方案和图2，详细叙述本发明的一个具体实施例。
[0024]现有100幅大尺寸航拍图(1024 X 1024)和100幅小尺寸Google Earth图(640X640)，每对图像中相同景物具有相近的缩放比例，其中第一组图像对如图2所示。
[0025]航拍图和Google Earth图的来源不同，但存在共同景物,如图2中的广场部分。以下是本发明的图像匹配步骤:
[0026]第一步，基于SIFT算法进行第一组图像对(即图2中两幅图像)的预匹配，得到103对匹配的特征点。
[0027]第二步，从103对匹配的特征点中随机选取两对，如图2中的A-a与B_b，根据这些特征点的坐标(见表1)，计算航拍图中特征点A和B之间的距离1AB,以及Google Earth图中特征点a和b之间的距离lab，得到第一个缩小比例:
[0028]rl = lab/IAB ^ 0.23147
[0029]其中，Ipq表示点P与点q之间的距离。设p点坐标为(xp，yp)，q点坐标为(xq，yq)，
【权利要求】
1.一种SIFT图像快速匹配的尺寸预处理方法，其特征是，当有M对具有相近缩放比例的大尺寸图像和小尺寸图像进行匹配时，在大批量进行匹配之前，以两幅图像中景物相同部分分辨率尽可能相等为原则，通过选择一对原始图像做预匹配，确定大尺寸图像的缩小比例；之后，按照该缩小比例将后续M-1幅大尺寸图像在匹配前做缩小预处理；具体步骤如下: 第一步，基于SIFT算法进行一组原始尺寸的大尺寸图像和小尺寸图像预匹配；如果M幅大尺寸图像是序列图像，选取第一对原始图像进行预匹配； 第二步，从所有匹配的特征点对中随机选取两对，计算得到第一个缩小比例rl =小尺寸图像两特征点间距离/大尺寸图像两特征点间距离； 第三步，重复第二步N-1次，得到N-1个缩小比例r2，…，rN ； 第四步，求取N个缩小比例的平均值，得到最终的大尺寸图像缩小比例: r = (rl+r2+...+rN) /N 第五步，在后续的M-1次图像匹配任务中，先将每幅大尺寸图像缩小到原来的r倍，再和对应的小尺寸图像进行匹配。
2.根据权利要求1所述的SIFT图像快速匹配的尺寸预处理方法，其特征是，如果M对缩放比例相近的大尺寸图像和小尺寸图像是序列图像，采用缩小比例跟踪校正策略，在第k+Ι次匹配时，采用第k次的匹配结果更新缩小比例r，k= 1，...，Μ-1，消除第一次预匹配所获得的缩小比例与后续某些图像对的实际缩小比例之间可能存在的偏差；具体如下: 步骤1:令k= l，r= 1，基于SIFT算法进行第一对原始尺寸的大尺寸图像和小尺寸图像预匹配； 步骤2:从第k对图像的所有匹配特征点对中随机选取两对计算缩放比例，k=I时，缩放比例=小尺寸图像两特征点间距离/大尺寸图像两特征点间距离；k>l时，缩放比例=小尺寸图像两特征点间距离/缩小的大尺寸图像两特征点间距离，重复该做法N次，共得到N个缩放比例rl,…，rN,求取其平均值r' = (rl+r2+…+rN)/N，得到更新的大尺寸图像缩小比例 T = TXr'； 步骤3:令k = k+Ι ; 步骤4:将第k幅大尺寸图像缩小到原来的r倍，再和第k幅小尺寸图像进行匹配； 步骤5:若k = M，停止匹配过程；否则，转到第二步。
3.根据权利要求1或2所述的SIFT图像快速匹配的尺寸预处理方法，其特征是，缩放比例的计算次数N取3?10。
【文档编号】G06T3/40GK104182930SQ201410407096
【公开日】2014年12月3日 申请日期:2014年8月16日 优先权日:2014年8月16日
【发明者】林秋华, 张存柱, 曹建超 申请人:大连理工大学