CN101556694A - 一种旋转图像的匹配方法 - Google Patents
一种旋转图像的匹配方法 Download PDFInfo
- Publication number
- CN101556694A CN101556694A CNA2009100394888A CN200910039488A CN101556694A CN 101556694 A CN101556694 A CN 101556694A CN A2009100394888 A CNA2009100394888 A CN A2009100394888A CN 200910039488 A CN200910039488 A CN 200910039488A CN 101556694 A CN101556694 A CN 101556694A
- Authority
- CN
- China
- Prior art keywords
- matching
- search window
- image
- angles
- theta
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Image Analysis (AREA)
Abstract
本发明涉及一种旋转图像的匹配方法,其包括:(1)确定输入角度搜索范围,以中间角度值及搜索范围的首末两个角度为搜索角度,旋转搜索窗口分别进行灰度归一化互相关匹配,获得各自匹配值;(2)将角度间距缩小一倍定为新角度搜索范围,以步骤(1)中取得最大匹配值角度为中心,左右再取两个角度,旋转搜索窗口进行灰度归一化互相关匹配;(3)重新比较步骤(2)中三个角度进行匹配时的匹配值,取最大者为中心,角度间距继续缩小一半,左右取两个角度,旋转搜索窗口进行灰度归一化互相关匹配;(4)重复步骤(3),直到角度间距达到指定的精度,进行匹配,获得最后的匹配角度、匹配位置及匹配值。本发明匹配时间更短、计算量更小。
Description
技术领域
本发明属于图形处理方法,尤其涉及一种旋转图像的匹配方法。
背景技术
随着科学技术的发展,图像匹配已成为图像信息处理领域中一项非常重要的技术。现在约40%的机器视觉应用中需要用到图像匹配技术,所涉及的应用领域从工业检测到导弹的地形匹配、光学和雷达的图像跟踪、工业流水线的自动监控、工业仪表的自动监控、资源分析、气象预报、医疗诊断、交通管理、文字识别以及图像检索等。图像匹配研究涉及到了图像采集、图像预处理、图像分割、特征提取等,并且与计算机视觉、多维信号处理和数值计算方法等紧密结合。它也是其它一些图像分析技术,如立体视觉、运动分析、数据融合等的基础:正因为其应用的广泛性,新的应用和新的要求逐步产生,使得匹配算法的研究逐步走向深入,出现了快速、稳定、鲁棒性好的匹配算法。因此,对现有匹配算法展开分析对于实际工程提高图像处理质量和识别精度具有非常重要的意义。
归一化互相关匹配算法通过计算模板和待匹配图像的灰度归一化互相关系数,即互相关值来确定匹配的程度。互相关法可以克服由于图像采集和系统误差带来的灰度畸变和几何畸变(平移变形和旋转)。
如PCB元件可能存在旋转,为了检测出旋转的角度,需要进行旋转的匹配。最简单的方法就是将在可能的旋转角度范围的所有角度依次代入到层叠二分算法中,求得最大匹配值的角度即实际元件的旋转角度。例如,假设旋转的范围在±15度之间,精确到1度,则将搜索窗口分别旋转-15、-14,直到15的31个角度应用匹配算法,取得最大匹配值的角度即元件的旋转角度。但这种方法随着匹配角度范围的增大,匹配时间也成倍增大,匹配时间太长,不适合应用于实际。
发明内容
针对现有技术存在的问题,本发明提供一种匹配时间更短、计算量更小的图像匹配方法。
为实现上述目的,本发明的技术方案为:一种旋转图像的匹配方法,包括以下步骤:
(1)确定输入角度搜索范围,以中间角度值及搜索范围的首末两个角度为搜索角度,旋转搜索窗口分别进行灰度归一化互相关匹配,获得各自匹配值;
(2)将角度间距缩小一倍定为新的角度搜索范围,以步骤(1)中取得的最大匹配值的角度为中心,左右再取两个角度,旋转搜索窗口进行灰度归一化互相关匹配;
(3)重新比较步骤(2)中三个角度进行匹配时的匹配值,取最大者为中心,角度间距继续缩小一半,左右取两个角度,旋转搜索窗口进行灰度归一化互相关匹配;
(4)重复步骤(3),直到角度间距达到指定的精度,进行匹配,获得的最后的匹配角度、匹配位置及匹配值。
灰度归一化互相关匹配运用层叠二分互相关匹配方法,其包括以下步骤:
(1)将待匹配图像和模板图像在宽度和高度方向均缩小一倍;
(2)重复步骤(1),直到匹配图像和模板图像达到指定的尺寸,在缩小的图像上使用灰度归一化互相关匹配获得最佳匹配位置(x,y);
(3)返回到上一级缩小一倍的图像上在(2x-1~2x+1,2y-1~2y+1)范围内的9个搜索位置再进行进一步的匹配,获得新的最佳匹配位置(x,y);
(4)重复步骤(3),直到返回到第一级图像,即原图像,进行最后的匹配,获得的最佳匹配位置及匹配值。
上述旋转搜索窗口的方法进行灰度归一化互相关值的计算,首先计算第一个搜索位置的搜索窗口内所有像素的坐标,其它搜索位置的坐标通过第一个搜索位置的坐标平移得到。
旋转搜索窗口的相关值计算公式为:
使用的旋转搜索窗口图像坐标系为原点在左上角,向右为x轴正方向,向下为y轴正方向,假设搜索窗口的尺寸为w*h像素大小,搜索窗口顺时针旋转θ时得到第一个搜索位置搜索窗口坐标,搜索窗口顺时针旋转θ时的坐标变换公式为:
搜索窗口逆时针旋转θ时的第一个搜索位置搜索窗口坐标,为使坐标值非负,再向上平移(w-1)*sinθ得到,搜索窗口顺时针旋转θ时的坐标变换公式为:
匹配过程中互相关值的计算采用增量方法进行计算,先计算出第一个搜索位置的灰度平均值,搜索窗口向右移动一个位置,新位置的灰度平均值的计算只需要在上一个平均值的基础上参考刚移入的一列和移出的一列,同理搜索窗口向下移动时,只需要参考刚移入的一行和刚移出的一行。
匹配过程中互相关值的计算采用增量方法进行计算,先计算出第一个搜索位置的灰度平均值,对于搜索窗口旋转的情况,当旋转角度为顺时针旋转时,向右移动一个像素时,最上一行和最右一列都新移入的,同时移出了最下一行及最左一列,新位置的灰度平均值的计算只需要在上一个平均值的基础上参考新移入的最上一行和最右一列;当旋转角度为逆时针旋转,向右移动一个像素时,只需参考新移入了最下一行和最右一列。
本发明相对于现有技术的优势:
本发明通过旋转搜索窗口的方法进行灰度归一化互相关值的计算,并结合运用层叠二分互相关匹配方法,大大减小了匹配的时间。
附图说明
图1是搜索角度和匹配值的关系曲线;
图2是层叠二分互相关匹配算法示意图;
图3是搜索窗口顺时针旋转时的坐标系变换示意图;
图4增量法计算搜索窗口平均灰度值示意图。
具体实施方式
本发明提供了一种图像匹配方法,如图3所示,其包括以下步骤:
(1)确定输入角度搜索范围,以中间角度值及搜索范围的首末两个角度为搜索角度,旋转搜索窗口分别进行灰度归一化互相关匹配,获得各自匹配值;
(2)将角度间距缩小一倍定为新的角度搜索范围,以步骤(1)中取得的最大匹配值的角度为中心,左右再取两个角度,旋转搜索窗口进行灰度归一化互相关匹配;
(3)重新比较步骤(2)中三个角度进行匹配时的匹配值,取最大者为中心,角度间距继续缩小一半,左右取两个角度,旋转搜索窗口进行灰度归一化互相关匹配;
(4)重复步骤(3),直到角度间距达到指定的精度,进行匹配,获得的最后的匹配角度、匹配位置及匹配值。
灰度归一化互相关匹配运用层叠二分互相关匹配方法,如图2所示,其包括以下步骤:
(1)将待匹配图像和模板图像在宽度和高度方向均缩小一倍;
(2)重复步骤(1),直到匹配图像和模板图像达到指定的尺寸,在缩小的图像上使用灰度归一化互相关匹配获得最佳匹配位置(x,y);
(3)返回到上一级缩小一倍的图像上在(2x-1~2x+1,2y-1~2y+1)范围内的9个搜索位置再进行进一步的匹配,获得新的最佳匹配位置(x,y);
(4)重复步骤(3),直到返回到第一级图像,即原图像,进行最后的匹配,获得的最佳匹配位置及匹配值。
上述旋转搜索窗口的方法进行灰度归一化互相关值的计算,首先计算第一个搜索位置的搜索窗口内所有像素的坐标,其它搜索位置的坐标通过第一个搜索位置的坐标平移得到。
旋转搜索窗口的相关值计算公式为:
使用的旋转搜索窗口图像坐标系为原点在左上角,向右为x轴正方向,向下为y轴正方向,假设搜索窗口的尺寸为w*h像素大小,搜索窗口顺时针旋转θ时得到第一个搜索位置搜索窗口坐标,搜索窗口顺时针旋转θ时的坐标变换公式为:
搜索窗口逆时针旋转θ时的第一个搜索位置搜索窗口坐标,为使坐标值非负,再向上平移(w-1)*sinθ得到,搜索窗口顺时针旋转θ时的坐标变换公式为:
由于模板图像是不变的,因此模板图像的平均灰度值P只需要计算一次,也只需要计算一次;其余部分每移动一个搜索位置都需要重新计算,但其中的搜索框窗口的图像的灰度平均值C(i,j)可以用增量法进行计算,减少计算时间。由于搜索过程是从左到右或从上到下一个一个像素移动搜索窗口的,因此可以这样简化计算搜索窗口图像的灰度平均值:先计算出第一个搜索位置的灰度平均值,搜索窗口向右移动一个位置,新位置的灰度平均值的计算只需要在上一个平均值的基础上考虑刚移入的一列和移出的一列,同理搜索窗口向下移动时,只需要考虑刚移入的一行和刚移出的一行。对于搜索窗口旋转的情况,当旋转角度为正时(顺时针旋转),向右移动一个像素时,最上一行和最右一列都新移入的,同时移出了最下一行及最左一列;当旋转角度为负时(逆时针旋转),向右移动一个像素时,新移入了最下一行和最右一列,同时移出了最上一行及最左一列,如图4所示。
例如,假设上一个搜索位置的平均值为Cpre,那么计算新移入的像素灰度之和为Gin,刚移出的像素灰度之和为Gout,搜索窗口的大小为m*n像素,则新位置搜索窗口的灰度平均值为:
假设匹配对象实际旋转了6度,在待匹配图像中的位置为(30,31),角度精度为1度,角度搜索范围为±12,则使用本发明进行匹配的实施方法如下:
(1)图1是各个搜索角度对应的最大互相关值,实际图像旋转了6度,搜索角度越接近6度互相关值越大。根据这一特性,本发明进行角度匹配的过程为:首先分别以-12、0和12度三个角度进行匹配,结果在0度时的匹配值最大;再以0为中心间距缩小一半,得到-6度和6度两个角度,分别以这两个角度进行匹配,并将得到的匹配值和0度时的匹配值比较,结果在6度时的匹配值最大;再以6度为中心,间距继续缩小一半,得到3度和9度,并分别以这两个角度进行匹配,将匹配值和6度时的匹配值比较,结果仍然在6度时的匹配值最大;如此不断进行下去最后是以6度为中心,间距为0.75,得到5.25度和6.75,比较三个角度的匹配值,最后在6度时的匹配值最大,此时三个角度的间距小于要求的1度精度,停止继续下去,匹配过程结束。角度匹配计算过程为0,-12,12,-6,6,3,9,4.5,7.5,5.25,6.75。最后的匹配结果为6度。
(2)在角度匹配过程中,每个角度的匹配值使用层叠二分互相关匹配求得,并同时返回最佳匹配位置,以0度的匹配过程为例进行说明,如图2所示。将待匹配图像和模板图像称为第一级图像,对一级图像进行连续两次缩小一半的操作,获得第二级和第三级图像,第三级图像的尺寸已经很小,对其使用传统灰度互相关匹配算法在整个搜索范围查找最佳匹配位置,假设是在位置(7,8)。然后返回到上一级图像,即第二级图像,则只需要在包括位置(2*7,2*8)在内及其周围的九个搜索位置进行进一步的精匹配,这九个位置为(13,15),(14,15),(15,15),(13,16),(14,16),(15,16),(13,17),(14,17),(15,17)。假设在位置(15,15)找到最大互相关值,则继续返回到上一级图像即第一级图像进行最后的精确匹配查找,同理也只需要在包括(2*15,2*15)在内及其周围的九个搜索位置进行查找,这九个位置为(29,29),(30,29),(31,29),(29,30),(30,30),(31,30),(29,31),(30,31),(31,31)。最后在(30,31)求得最大互相关值,这个位置即匹配位置,对应的互相关值即0度时的匹配值。
(3)第一个搜索位置的搜索窗口内的图像像素坐标和互相关值的计算根据式(5)、式(7)和式(8)计算,即通过旋转搜索窗口的方式进行计算,如图3所示。
(4)匹配过程下一个搜索位置搜索窗口内的图像像素坐标根据第一个搜索窗口的坐标平移得到,搜索窗口内的灰度平均值根据式(9)使用增量法计算,以最大限度减少计算互相关值的计算时间,如图4所示。
Claims (6)
1、一种旋转图像的匹配方法,其特征在于包括以下步骤:
(1)确定输入角度搜索范围,以中间角度值及搜索范围的首末两个角度为搜索角度,旋转搜索窗口分别进行灰度归一化互相关匹配,获得各自匹配值;
(2)将角度间距缩小一倍定为新的角度搜索范围,以步骤(1)中取得的最大匹配值的角度为中心,左右再取两个角度,旋转搜索窗口进行灰度归一化互相关匹配;
(3)重新比较步骤(2)中三个角度进行匹配时的匹配值,取最大者为中心,角度间距继续缩小一半,左右取两个角度,旋转搜索窗口进行灰度归一化互相关匹配;
(4)重复步骤(3),直到角度间距达到指定的精度,进行匹配,获得最后的匹配角度、匹配位置及匹配值。
2、根据权利要求1所述的旋转图像的匹配方法,其特征在于:灰度归一化互相关匹配运用层叠二分互相关匹配方法,其包括以下步骤:
(1)将待匹配图像和模板图像在宽度和高度方向均缩小一倍;
(2)重复步骤(1),直到匹配图像和模板图像达到指定的尺寸,在缩小的图像上使用灰度归一化互相关匹配获得最佳匹配位置(x,y);
(3)返回到上一级缩小一倍的图像上在(2x-1~2x+1,2y-1~2y+1)范围内的9个搜索位置再进行进一步的匹配,获得新的最佳匹配位置(x,y);
(4)重复步骤(3),直到返回到第一级图像,即原图像,进行最后的匹配,获得的最佳匹配位置及匹配值。
3、根据权利要求2所述的旋转图像的匹配方法,其特征在于:上述旋转搜索窗口的方法进行灰度归一化互相关值的计算,首先计算第一个搜索位置的搜索窗口内所有像素的坐标,其它搜索位置的坐标通过第一个搜索位置的坐标平移得到。
4、根据权得要求3所述的旋转图像的匹配方法,其特征在于:旋转搜索窗口的相关值计算公式为:
使用的旋转搜索窗口图像坐标系为原点在左上角,向右为x轴正方向,向下为y轴正方向,假设搜索窗口的尺寸为w*h像素大小,搜索窗口顺时针旋转θ时得到第一个搜索位置搜索窗口坐标,搜索窗口顺时针旋转θ时的坐标变换公式为:
搜索窗口逆时针旋转θ时的第一个搜索位置搜索窗口坐标,为使坐标值非负,再向上平移(w-1)*sinθ得到,搜索窗口顺时针旋转θ时的坐标变换公式为:
5、根据权得要求4所述的旋转图像的匹配方法,其特征在于:匹配过程中互相关值的计算采用增量方法进行计算,先计算出第一个搜索位置的灰度平均值,搜索窗口向右移动一个位置,新位置的灰度平均值的计算只需要在上一个平均值的基础上参考刚移入的一列和移出的一列,同理搜索窗口向下移动时,只需要参考刚移入的一行和刚移出的一行。
6、根据权得要求4所述的旋转图像的匹配方法,其特征在于:匹配过程中互相关值的计算采用增量方法进行计算,先计算出第一个搜索位置的灰度平均值,对于搜索窗口旋转的情况,当旋转角度为顺时针旋转时,向右移动一个像素时,最上一行和最右一列都新移入的,同时移出了最下一行及最左一列,新位置的灰度平均值的计算只需要在上一个平均值的基础上参考新移入的最上一行和最右一列;当旋转角度为逆时针旋转,向右移动一个像素时,只需参考新移入了最下一行和最右一列。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009100394888A CN101556694B (zh) | 2009-05-15 | 2009-05-15 | 一种旋转图像的匹配方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009100394888A CN101556694B (zh) | 2009-05-15 | 2009-05-15 | 一种旋转图像的匹配方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101556694A true CN101556694A (zh) | 2009-10-14 |
CN101556694B CN101556694B (zh) | 2011-07-06 |
Family
ID=41174799
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2009100394888A Active CN101556694B (zh) | 2009-05-15 | 2009-05-15 | 一种旋转图像的匹配方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101556694B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101840570A (zh) * | 2010-04-16 | 2010-09-22 | 广东工业大学 | 一种快速图像拼接方法 |
CN102915115A (zh) * | 2012-09-25 | 2013-02-06 | 上海华勤通讯技术有限公司 | 根据眼部调整显示画面的方法 |
CN107004139A (zh) * | 2014-12-01 | 2017-08-01 | 欧司朗股份有限公司 | 借助于互相关的图像处理 |
CN107305681A (zh) * | 2016-04-22 | 2017-10-31 | 北京仿真中心 | 一种用于光学场景图像的快速扩充方法 |
CN107851196A (zh) * | 2016-06-01 | 2018-03-27 | 深圳配天智能技术研究院有限公司 | 一种图像模式匹配的方法及装置 |
CN108010074A (zh) * | 2017-10-19 | 2018-05-08 | 宁波蓝圣智能科技有限公司 | 一种基于机器视觉的工件检测方法及系统 |
CN108992089A (zh) * | 2013-03-15 | 2018-12-14 | 敦提大学校董事会 | 医疗设备及其可视化 |
CN112669411A (zh) * | 2020-12-31 | 2021-04-16 | 一诺仪器(中国)有限公司 | 一种矢量网络分析仪迹线平滑方法 |
-
2009
- 2009-05-15 CN CN2009100394888A patent/CN101556694B/zh active Active
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101840570A (zh) * | 2010-04-16 | 2010-09-22 | 广东工业大学 | 一种快速图像拼接方法 |
CN102915115A (zh) * | 2012-09-25 | 2013-02-06 | 上海华勤通讯技术有限公司 | 根据眼部调整显示画面的方法 |
CN108992089A (zh) * | 2013-03-15 | 2018-12-14 | 敦提大学校董事会 | 医疗设备及其可视化 |
CN107004139A (zh) * | 2014-12-01 | 2017-08-01 | 欧司朗股份有限公司 | 借助于互相关的图像处理 |
CN107305681A (zh) * | 2016-04-22 | 2017-10-31 | 北京仿真中心 | 一种用于光学场景图像的快速扩充方法 |
CN107305681B (zh) * | 2016-04-22 | 2020-09-18 | 北京仿真中心 | 一种用于光学场景图像的快速扩充方法 |
CN107851196A (zh) * | 2016-06-01 | 2018-03-27 | 深圳配天智能技术研究院有限公司 | 一种图像模式匹配的方法及装置 |
CN107851196B (zh) * | 2016-06-01 | 2020-02-14 | 深圳配天智能技术研究院有限公司 | 一种图像模式匹配的方法及装置 |
CN108010074A (zh) * | 2017-10-19 | 2018-05-08 | 宁波蓝圣智能科技有限公司 | 一种基于机器视觉的工件检测方法及系统 |
CN112669411A (zh) * | 2020-12-31 | 2021-04-16 | 一诺仪器(中国)有限公司 | 一种矢量网络分析仪迹线平滑方法 |
Also Published As
Publication number | Publication date |
---|---|
CN101556694B (zh) | 2011-07-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101556694B (zh) | 一种旋转图像的匹配方法 | |
US10867189B2 (en) | Systems and methods for lane-marker detection | |
CN101216941B (zh) | 剧烈光照变化下基于角点匹配与光流法的运动估计方法 | |
CN102519481B (zh) | 一种双目视觉里程计实现方法 | |
CN103456022B (zh) | 一种高分辨率遥感图像特征匹配方法 | |
CN101556695A (zh) | 一种图像匹配方法 | |
CN102541954B (zh) | 一种商标检索方法及系统 | |
CN103136525B (zh) | 一种利用广义Hough变换的异型扩展目标高精度定位方法 | |
CN102426019A (zh) | 一种无人机景象匹配辅助导航方法及系统 | |
CN101984463A (zh) | 全景图合成方法及装置 | |
CN103679702A (zh) | 一种基于图像边缘矢量的匹配方法 | |
CN102006425A (zh) | 一种基于多摄像机的视频实时拼接方法 | |
CN109712071B (zh) | 基于航迹约束的无人机图像拼接与定位方法 | |
CN105389774A (zh) | 对齐图像的方法和装置 | |
CN102800097A (zh) | 多特征多级别的可见光与红外图像高精度配准方法 | |
CN102903109B (zh) | 一种光学影像和sar影像一体化分割配准方法 | |
CN111028292A (zh) | 一种亚像素级图像匹配导航定位方法 | |
CN104820997A (zh) | 一种基于分块稀疏表达与hsv特征融合的目标跟踪方法 | |
CN111429344B (zh) | 基于感知哈希的激光slam闭环检测方法及系统 | |
CN103761739A (zh) | 一种基于半全局能量优化的影像配准方法 | |
CN115267724B (zh) | 一种基于激光雷达可估位姿的移动机器人位置重识别方法 | |
CN104457758A (zh) | 基于视频采集的Visual Map数据库建立方法及利用该数据库的室内视觉定位方法 | |
CN103700082A (zh) | 基于对偶四元数相对定向的图像拼接方法 | |
CN103733225A (zh) | 特征点对等系统、特征点对等方法、以及记录介质 | |
CN102324043B (zh) | 基于dct的特征描述算子及优化空间量化的图像匹配方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |