CN102654902A - 基于轮廓向量特征的嵌入式实时图像匹配方法 - Google Patents
基于轮廓向量特征的嵌入式实时图像匹配方法 Download PDFInfo
- Publication number
- CN102654902A CN102654902A CN2012100128010A CN201210012801A CN102654902A CN 102654902 A CN102654902 A CN 102654902A CN 2012100128010 A CN2012100128010 A CN 2012100128010A CN 201210012801 A CN201210012801 A CN 201210012801A CN 102654902 A CN102654902 A CN 102654902A
- Authority
- CN
- China
- Prior art keywords
- template
- point
- mapping
- image
- layer
- 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
本发明提出一种基于轮廓向量特征的嵌入式实时图像匹配方法。采用基于X、Y方向向量的线特征,抗图像畸变、噪声、遮挡、光照变化、极性反转等能力强。采用图像金字塔搜索策略,在高层小分辨率待测图中快速匹配模板,然后逐级向下精确找出目标位置,大大减少匹配时间。根据模板图像具体信息,自动计算出最佳金字塔分层数及每层金字塔模板匹配最佳旋转角度步长。图像金字塔最高层三级筛选匹配策略根据待测图具体内容进行处理,第一、二级筛选,仅利用加减运算和判断语句筛除非目标位置,在嵌入式***中比乘除运算更高效,第三级则仅对满足上两级要求的较少位置进行处理,极大提高匹配速度。整体方法可实现任意角度、坐标下的目标匹配定位工作。
Description
技术领域
本发明涉及一种图像匹配方法,具体是指一种基于轮廓向量特征的嵌入式实时图像匹配方法。
背景技术
图像匹配,是指根据给定目标图像所包含的颜色、纹理、形状或它们的组合特征等信息建立索引,并且通过图像特征向量间的相似性度量,在一幅(或一批)图像中寻找与给定目标相同或相似的图像区域(子图像)的过程。通常将给定目标图像称为模板图像,而将后续待搜索的序列图像称为待测图像。
图像匹配技术在诸多领域内有着广泛的应用,其中包括如下方面。医学图像处理中的CT、MRI、PET、SPECT图像结构信息融合。遥感图像领域中的多电磁波段图像(如微波、红外、可见光或多光谱图像)信息融合,定位和识别已定义或已知特征的如飞机场、高速路、车站等场景,矿物寻找、特定区域搜索、农作物估产、森林和水源保护、台风追踪、天气预报等方面。在机器视觉领域中的字符识别、标志确认、波形分析、运动跟踪、序列图像分析等研究方向。在工业检测领域中的目标产品定位、表面缺陷检测、工业测量、产品质量评估、构建三维信息等。
根据图像匹配的基本原理可以将图像匹配算法分为二类:基于灰度相关的匹配和基于特征的匹配。其中,基于灰度相关的匹配方法简单,匹配精度较高。但是需要处理的数据量多,计算复杂,很难达到实时匹配的要求,对噪声较为敏感,当外界条件的变化如光照条件的影响或是目标图像缺损、遮挡都会很大地影响匹配精度。基于特征的匹配对图像畸变、噪声、遮挡、光照变化等具有一定的鲁棒性,计算量小,匹配速度高,目前先进的图像匹配都采用了这一类的方法。
Luigi Di stefano等提出了一种基于零均值归一化积相关的边界局部相关算法。Wang等首先利用灰度相关法进行出匹配,再用MAPSAC算法对变换矩阵进行修正,以消除误匹配。Lowe等利用尺度空间性质,提出了一种尺度不变关键点的检测方法,以同时在尺度域和空间域取得极值的点作为关键点,并以关键点作为特征区域的中心。同时,关键点所在的尺度还用以确定特征区域的大小。这种方法较好地同时解决了特征区域定位和大小选择的问题。在每个特征区域内,Lowe以梯度方向的直方图为基础构造了一种称为SIFT的局部不变特征提取性能的衡量标准。YKe等提出了PCA-SIFT算法,通过采用PCA方法把SIFT的128维特征向量降到36维,从而加快匹配速度。Herbert Bay等提出了一种新的特征提取算法SURF,它通过快速Hessian矩阵求极值点,并通过非极大值抑制进行插值方法得到特征点,计算特征描述向量是通过计算区域的Harr响应得到。Keller等提出的基于Fourier-Mellin变换的算法从原理上很好地解决了具有平移、旋转和尺度差异的图像匹配问题。随着极坐标Fourier变换计算方法的改进,Averbuch等利用伪极坐标Fourier变换技术改进了对图像平移、旋转和尺度的检测稳定性。
在商业领域,很多公司开发出了成熟的图像匹配定位算法包,如美国Cognex公司的PatMax根据模板定位原理进行几何特征模板匹配来取代基于像素灰度的匹配方式。但是几何模式匹配的基础是图像的几何特征提取,因而此法的定位精度依赖于图像几何特征的边缘检测精度,在特征提取过程中此方法比NCC法对随机噪声更敏感。比利时Euresys公司开发的基于点特征的图像分析软件EasyFind,采用创新的特征点技术,EasyFind可迅速在图像中找到一个或多个参考模型的例证。它不是将样本图像的像素层次与参考模型相对比,而是精心选择模型中的重要特征。Adept公司开发的HexSight机器视觉软件***,其定位工具是根据几何特征,采用最先进的轮廓检测技术来识别对象和模式。Matrox公司的Inspector软件则先提取模板和待匹配对象的边缘点,再将边缘点进行匹配,从而计算出模板和对象之间的几何位姿关系。日本Keyence公司的定位软件ShapeTrax II则同样通过物体轮廓抽取的形状为特征基础,进行目标定位。德国MVTec公司的HALCON软件基于形状的匹配模块则使用边缘特征定位物体。此外,加拿大的Dalsa公司、丹麦的JAI/Pulnix公司、加拿大的Coreco Imaging公司、美国的Navitar公司、瑞士的Photonfocus公司、日本的CCS公司等世界一流的图像和机器视觉技术领先企业对图像匹配定位技术都有深入的研究和应用。
由于拍摄时间、拍摄角度、拍摄位置、自然环境的变化、多种传感器的使用、传感器本身的缺陷及噪声、相同规格不同个体之间的略微差别等因素的影响,对同一或相同规格的不同目标拍摄的图像之间会存在灰度失真和几何畸变等差别。同时,图像预处理过程又会引入误差,这都会导致模板图像与待测图像之间通常存在着一定程度上的差异。其差异主要表现为:不同的分辨率、不同的灰度属性、不同的位置(平移和旋转)、不同的大小(缩小和放大)、以及其他的非线性变化(部分物体被遮挡、光照不均)等。在这种情况下,如何使匹配算法精度高、正确匹配率高、速度快、鲁棒性和抗干扰性强成为人们关心的问题。
发明内容
本发明目的在于针对以往图像匹配方法计算量大、耗时严重、抗图像畸变、噪声、遮挡和光照变化等外界环境影响力小、模板极性反转时不易匹配的等多方面缺点,提出了一种在嵌入式***上利用轮廓点方向向量为特征,考虑目标平移和旋转,根据待测图像具体图像内容,实现图像目标快速、准确匹配定位。
按照本发明提供的技术方案,所述基于轮廓向量特征的嵌入式实时图像匹配方法包括以下步骤:
第一步,模板图像特征提取模块在输入的模板灰度图中,利用2*2均值法构建图像金字塔、高斯滤波、Sobel边缘提取、大津法阈值分割、非极大值抑制、高低阈值处理算法,提取用于图像匹配的模板图像金字塔各层X、Y方向向量特征图和二值化轮廓图;
第二步,待测图像特征提取模块在输入的待测灰度图中,利用均值法构建图像金字塔、高斯滤波、Sobel边缘提取、大津法阈值分割算法,提取与模板特征进行相似度吻合比较的待测图图像金字塔各层X、Y方向向量特征图和二值化轮廓图;
第三步,匹配参数自动设置模块通过获取模板二值化轮廓图中距质心最远边缘点和统计边缘点总数,计算最佳旋转角度步长、金字塔分层数;
第四步,图像金字塔最高层匹配模块分析不同待测图像的具体内容,利用三级筛选匹配策略,剔除待测图像中非目标位置,在图像金字塔最高层实现快速准确的模板粗匹配;
第五步,图像金字塔非最高层匹配模块根据上一层图像金字塔传递来的匹配信息,通过图像金字塔构建逆向工程原理,定位出本层目标位置,然后将本层定位信息传递到下一层图像金字塔进行更精确匹配,直至图像金字塔第一层。
进一步的,所述第一步具体为:
(1.1)利用2*2均值法构建多层模板图像金字塔,模板原始灰度图作为第一层,高一层金字塔图像则根据其下一层灰度图相邻四个像素的灰度值的均值求得;
(1.2)对图像金字塔各层模板灰度图利用高斯滤波去噪,高斯卷积模板如下:
继而应用Sobel边缘检测算子提取图像中目标轮廓点的X、Y方向向量及对应的梯度值,X、Y方向检测子如下:
(1.3)对图像金字塔各层模板图像对应的梯度图进行大津法全局阈值分割,剔除部分梯度较小的边缘点,以删除背景信息,保留前景目标信息;
(1.4)对经阈值分割的模板图像金字塔各层梯度图利用非极大值抑制的方法,比较边缘点梯度值与该边缘点梯度正、逆方向上相邻两个点的大小,若大于这两点梯度值,则保留该边缘点,否则去除该边缘点;
(1.5)对经非极大值抑制的模板梯度图利用高低阈值处理方法,提取模板前景目标的关键特征,去除细小、琐碎边缘,将此时模板梯度取直方图,当累加梯度值由小到大排列所对应的边缘点数达到0.4~0.6倍总边缘点数时,所对应的梯度值取为高阈值,所述高阈值的0.3~0.5倍取为低阈值;
判断经非极大值抑制的模板梯度图中梯度值与高低阈值的关系,若一点大于高阈值则保留,小于低阈值则剔除,介于高低阈值之间且与大于高阈值的点8方向连通则保留,否则剔除;
(1.6)根据上一步处理所得的模板梯度图与步骤1.2所得的X、Y方向向量图对应,获得图像金字塔各层模板X、Y方向向量图和二值化轮廓图。
所述第二步具体为:
(2.1)利用2*2均值法构建多层待测图图像金字塔,待测图原始灰度图作为第一层,高一层金字塔图像则根据其下一层灰度图相邻四个像素的灰度值的均值求得;
(2.2)对图像金字塔各层待测灰度图利用高斯滤波去噪,高斯卷积模板如下:
继而应用Sobel边缘检测算子提取图像中目标轮廓点的X、Y方向向量及对应的梯度值,X、Y方向检测子如下:
(2.3)对图像金字塔各层待测图像对应的梯度图进行大津法全局阈值分割,剔除部分梯度较小的边缘点,以删除背景信息,保留前景目标信息;
(2.4)根据上一步处理所得的待测图的梯度图与步骤2.2所得的X、Y方向向量图对应,获得图像金字塔各层对应的待测图的X、Y方向向量图和二值化轮廓图。
所述第三步具体为:
(3.1)根据处理过的图像金字塔各层模板二值化轮廓图计算其质心坐标及轮廓总点数,并求解距质心点欧式距离最远的模板轮廓点,以下为质心点求解:
模板某层金字塔质心坐标为(XC,XC),边缘总点数为K,边缘点X、Y方向坐标分别为XC、XC;
(3.2)在模板图像金字塔各层上,假设距离质心最远的轮廓点绕质心旋转1~2个像素距离后距离不变,则根据余弦定理求得的该点此时所转过的角度认为是该模板旋转角度步长;设点L为模板质心点,点p为模板最远轮廓点,点w、q为点p旋转过1、2个像素距离后对应的点,点L、p、w构成三角形Lpw,三边边长分别为Lp=Lw=a,pw=1,则由余弦定理求∠wLp得:
此时,∠wLp为保证模板旋转后至少有一个轮廓点旋转出原来像素点坐标所对应的最小角度,模板的旋转角度步长需保证模板旋转前后能够区分出来;
(3.3)当模板图像金字塔某一层轮廓点数大于20个,而其更高一层轮廓点数小于20时,则认为该层为图像金字塔最高层,能够进行整体有效的图像匹配工作。
所述第四步具体为:
(4.1)构建搜索框:根据图像金字塔最高层的模板二值化轮廓图,以其质心点N为旋转基准点,构建包围模板在旋转角度范围内任意角度位置的最小外接矩形,设置旋转角度范围为0~360°,并且使该矩形四条边与原模板图像坐标轴平行,该矩形即为搜索框,纵长M个像素,横长H个像素,旋转基准点即为搜索框基准点N;
(4.2)待测图搜索框范围内边缘总点数统计:假设该层模板轮廓总点数为K,相似度阈值为S,以搜索框基准点N在待测图中逐像素遍历移动搜索框,判断搜索框覆盖区域中待测图边缘点数,若统计点数大于K*S,则进入第4.3步;否则,则认为待测图该像素上不存在以它为质心的任意角度目标,而后搜索框成“S”形在待测图中移动再次执行本步操作;若搜索框沿X轴移动,只需在上一步统计的边缘总点数基础上,加上沿移动方向向前新进入搜索框范围的一列M个像素点中的待测图边缘点数,并减去相应移出搜索框范围的一列M个像素点中的待测图边缘点数;若搜索框沿Y轴移动,只需在上一步统计的边缘总点数基础上,加上沿移动方向向前新进入搜索框范围的一行H个像素点中的待测图边缘点数,并减去相应移出搜索框范围的一行H个像素点中的待测图边缘点数;
(4.3)模板旋转与待测图对应边缘点数统计:以搜索框基准点N与模板质心对齐并以搜索框基准点N为旋转中心,按该层旋转角度步长旋转模板二值化轮廓图,统计每旋转一个角度步长的模板边缘点对应到待测图的二值化轮廓图相应坐标下同时也存在边缘点的个数,当模板旋转至某角度该统计数大于K*S时,则进行第4.4步,否则,继续旋转至下一个角度;
其中轮廓点旋转采用如下公式:
(x0 x0)T为边缘点原坐标,(x1 x1)T为边缘点旋转后坐标,θ为旋转角度步长;
(4.4)基于方向向量特征进行匹配判断:利用模板和待测图边缘点的方向向量通过如下公式,对旋转至某角度的模板与待测图进行匹配判断:
其中,di=(ai,bi)由模板边缘点上的X、Y方向向量组成,ei=(wi,ri)由待测图边缘点上的X、Y方向向量组成, 分别为方向向量对应的梯度强度值;为归一化的模板边缘点上的X、Y方向向量,为归一化的待测图边缘点上的X、Y方向向量;s表示模板与待测图的相似度,取值范围为0~1,取s=0时不匹配,随着s增加相似度提高,s=1时完全匹配;因使用轮廓线特征,只要能够提取出轮廓,即可进行图像匹配操作;
s′T<(K-1)S+T
则这T个点相似度未达到最低阈值,即使剩余的K-T个完全相同也不可能达到模板相似度阈值S,因此结束本次相似度计算;
若在某位置匹配出了满足相似度阈值的图形,则记录此时模板质心点坐标和旋转角度作为位置信息,当遍历完整个待测图后,将所述位置信息传递到下一层金字塔中。
所述第五步具体为:
(5.1)根据图像金字塔上一层在待测图中定位的目标位置信息,即:
(质心点坐标(x0 y0)T,旋转角度θ0),
通过图像金字塔构建逆向工程原理,计算本层目标位置信息(质心点坐标(x y)T,旋转角度θ)如下:
在以(x y)T为中心、周围2~4个像素范围内,以θ为中心角度、左右1~2个该层对应的角度步长范围内,直接采用基于方向向量特征进行匹配判断的方法,即利用模板和待测图边缘点的方向向量通过如下公式,对旋转至某角度的模板与待测图进行匹配判断:
s′T<(K-1)S+T ,
则这T个点相似度未达到最低阈值,结束本次相似度计算;
若在某位置匹配出了满足相似度阈值的图形,则记录此时模板质心点坐标和旋转角度作为位置信息,当遍历完整个待测图后,将所述位置信息传递到下一层金字塔中;
(5.2)根据以下两步筛选策略进行无效位置信息剔除:
当目标位置信息中在相同质心点坐标上检测有不同角度目标位置时,找出相似度最高所对应的角度,然后计算剩余角度与其相差的度数,当差值的绝对值大于10°时,则认为该角度对应的目标属于该质心点坐标下的第二个目标,保留该角度值,否则不保留该位置信息;继而再比较剩余角度与上述两个角度的差值,同样若差值的绝对值大于10°时保留为第三个认可角度,其它角度依次类推;处理完成后,即能够保证每个质心点坐标上,只会对应真实的目标角度;
再考虑临近质心点坐标对应同一真实目标的情况,首先,找出相似度最高的位置对应的质心点坐标,然后比较该质心点坐标与剩余位置信息的质心点坐标的关系,当二者X、Y方向坐标差值的绝对值大于设定阈值,则保留该质心点坐标对应的目标位置信息作为第二目标,否则不保留该目标位置信息;剩余目标位置信息同样与记录的两个认可位置信息进行质心坐标比较,依此类推直至处理完所有位置;
此时所得所有目标位置信息即认为是该层图像金字塔待测图中所有真实目标对应的信息,将其传递到下一层图像金字塔。
本发明具有以下优点:采用基于图像线特征之一的轮廓作为特征,并用轮廓线上每个点的方向向量对其描述,相对于使用图像的点特征、面特征来说,对图像畸变、噪声、遮挡、光照变化、极性反转等具有一定的鲁棒性,并且计算量小,匹配速度高。采用图像金字塔搜索策略,在图像高层小分辨率待测图中可快速匹配模板,然后逐级向下精确匹配出模板在原始待测图中的位置,大大减少匹配时间。图像金字塔最高层三级筛选匹配策略,根据待测图具体内容,逐渐提高每级达标要求,大量剔除待测图中非模板目标位置,如第一级筛选,仅利用加法运算统计搜索框大小范围内待测图对应区域中边缘点数,第二级筛选,则在满足第一级筛选的位置上进行模板二值化轮廓图与待测图逐角度相应边缘点有无判断,在嵌入式***中进行加法和判断比乘除运算更高效,第三级则为基于轮廓点方向向量的匹配,此时虽然进行乘除运算,但整幅待测图中能达到该级的位置并不多,相对于穷举遍历法的匹配策略,提高匹配速度。整体算法可实现目标任意角度旋转、坐标情况下的图像匹配定位。
附图说明
图1是图像金字塔最高层搜索匹配示意图。
图2是图像金字塔构建示意图。
图3是轮廓点旋转角度步长计算示意图。
图4是本发明流程图。
具体实施方式
本发明是可用于嵌入式***对给定目标图案进行轮廓向量特征提取并据此在后续输入图像序列中进行快速相同或相似目标图案匹配定位的方法。目的在于针对以往图像匹配方法计算量大、耗时严重、抗图像畸变、噪声、遮挡和光照变化等外界环境影响力小的多方面缺点,提出在嵌入式***上实现图像快速、精确匹配的方法。
下面结合附图和实施例对本发明作进一步说明。
如图1所示,本发明所述待测图像区域为矩形OXZY,并以O点为原点,OX方向为X轴,OY方向为Y轴,建立坐标系,图像像素坐标以该坐标系表示。并且规定,后述图像及点坐标均以此坐标系为准。矩形abcd为根据模板efg及旋转角度范围(在此设全角度旋转即360°)构建的搜索框,点N为搜索框基准点,该搜索框保证模板efg以N点为基准点360°范围内任意旋转,距点N距离最远点f始终落在其内部或边上,并且该框四边各自平行于X、Y轴,该图状态为搜索框游走到待测图中某一位置时的情况。
如图4所示,本发明的工作过程具体说明如下:
首先,模板图像特征提取模块在输入的模板灰度图中,利用2*2均值法构建图像金字塔、高斯滤波、Sobel边缘提取、大津法阈值分割、非极大值抑制、高低阈值处理算法,提取用于图像匹配的模板图像金字塔各层X、Y方向向量特征图和二值化轮廓图;
第二步,待测图像特征提取模块在输入的待测灰度图中,利用均值法构建图像金字塔、高斯滤波、Sobel边缘提取、大津法阈值分割算法,提取与模板特征进行相似度吻合比较的待测图图像金字塔各层X、Y方向向量特征图和二值化轮廓图;
第三步,匹配参数自动设置模块通过获取模板二值化轮廓图中距质心最远边缘点和统计边缘点总数,计算最佳旋转角度步长、金字塔分层数等参数;
第四步,图像金字塔最高层匹配模块分析不同待测图像的具体内容,利用三级筛选匹配策略,剔除待测图像中非目标位置,在图像金字塔最高层实现快速准确的模板粗匹配;
第五步,图像金字塔非最高层匹配模块根据上一层图像金字塔传递来的匹配信息,通过图像金字塔构建逆向工程原理,快速准确定位出本层目标位置,然后将本层定位信息传递到下一层图像金字塔进行更精确匹配,直至图像金字塔第一层。
具体的,第一步包括以下步骤:
(1.1)如图2所示,利用2*2均值法构建多层模板图像金字塔,模板原始灰度图作为第一层,高一层金字塔图像则根据其下一层灰度图相邻四个像素的灰度值的均值求得,这样每上升一层,分辨率变为其下一层的1/4倍;
(1.2)对图像金字塔各层模板灰度图利用高斯滤波去噪,高斯卷积模板如下所示:
继而应用Sobel边缘检测算子提取图像中目标轮廓点的X、Y方向向量及对应的梯度值,X、Y方向检测子如下所示:
(1.3)对图像金字塔各层模板图像对应的梯度图进行大津法全局阈值分割,剔除部分梯度较小的边缘点,以删除背景信息,保留前景目标信息。
(1.4)对经阈值分割的模板图像金字塔各层梯度图利用非极大值抑制的方法,比较某边缘点梯度值与其梯度正、逆方向上相邻两个点的大小,若大于这两点梯度值,则保留该边缘点,否则去除该点。至此实现模板图大部分轮廓细化到一个像素,保留局部梯度最大值,提高特征提取的精度。
(1.5)对经非极大值抑制的模板梯度图利用高低阈值处理方法,提取模板前景目标的关键特征,去除细小、琐碎边缘,将此时模板梯度取直方图,当累加梯度值由小到大排列所对应的边缘点数达到0.4~0.6倍总边缘点数时,所对应的梯度值取为高阈值,高阈值的0.3~0.5倍取为低阈值。
判断经非极大值抑制的模板梯度图中梯度值与高低阈值的关系,若某点大于高阈值则保留,小于低阈值则剔除,介于高低阈值之间且与大于高阈值的点8方向连通则保留,否则剔除。
(1.6)又经过上一步的处理后的模板梯度图基本已满足图像匹配的需要,根据处理所得的模板梯度图与步骤(1.2)所得的X、Y方向向量图对应,获得图像金字塔各层模板X、Y方向向量图和二值化轮廓图。
第二步包括以下步骤:
(2.1)同样利用2*2均值法构建多层待测图图像金字塔,待测图原始灰度图作为第一层,高一层金字塔图像则根据其下一层灰度图相邻四个像素的灰度值的均值求得。
(2.2)对图像金字塔各层待测灰度图利用高斯滤波去噪,高斯卷积模板如下所示:
继而应用Sobel边缘检测算子提取图像中目标轮廓点的X、Y方向向量及对应的梯度值,X、Y方向检测子如下所示:
(2.3)对图像金字塔各层待测图像对应的梯度图进行大津法全局阈值分割,剔除部分梯度较小的边缘点,以删除背景信息,保留前景目标信息。
(2.4)经过以上处理后的待测图的梯度图基本已满足图像匹配的需要,根根据处理所得的待测图的梯度图与步骤(2.2)所得的X、Y方向向量图对应,获得图像金字塔各层对应的待测图的X、Y方向向量图和二值化轮廓图。
第三步包括以下步骤:
(3.1)根据处理过的图像金字塔各层模板二值化轮廓图计算其质心坐标及轮廓总点数,并求解距质心点欧式距离最远的模板轮廓点,如下为质心求解:
模板某层金字塔质心坐标为(XC,YC),边缘总点数为K,边缘点X、Y方向坐标分别为XC、YC。
(3.2)在图像金字塔各层上,假设距离质心最远的轮廓点绕质心旋转1~2个像素距离后距离不变,则根据余弦定理求得的该点此时所转过的角度可认为是该模板旋转角度步长。如图3为点旋转示意图,点L为质心点,点p为最远轮廓点,点w、q为点p旋转过1、2个像素距离后对应的点,假设三角形两边长度Lp=Lw=a,pw=1,则由余弦定理求∠wLp得:
此时,∠wLp为保证模板旋转后至少有一个不轮廓点旋转出原来像素点坐标所对应的最小角度,模板的旋转角度步长θ必须满足θ≥∠wLp。但由于取θ=∠wLp时将线段Lp绕点L旋转θ,可发现在与点L距离在0.6*Lp长度范围内的部分点,未旋转出原来像素点坐标,若在图像匹配时,模板上有大量点存在这种情况,特别是当这种情况的点数占总点数比例达到相似度阈值时,前后两个角度的图像可认为未分辨开,匹配结果是冗余的。因此,旋转角度步长需保证模板旋转前后能够区分出来。
(3.3)图像金字塔构建层数设定的原则是保证每一层图像未发生较大畸变,形状上有一定相似性,通过分析得到,当模板图像金字塔某一层轮廓点数大于20个,而其更高一层轮廓点数小于20时,则认为该层为图像金字塔最高层,可进行整体有效的图像匹配工作。
第四步包括以下步骤:
(4.1)构建搜索框:如图1中矩形abcd,根据图像金字塔最高层的模板二值化轮廓图,以其质心点N为旋转基准点,构建包围模板在旋转角度范围内任意角度位置的最小外接矩形(此处设置旋转角度范围为0~360°),并且使该矩形四条边与原模板图像坐标轴平行,该矩形即为搜索框,纵长M个像素,横长H个像素,旋转基准点即为搜索框基准点N。
(4.2)待测图搜索框范围内边缘总点数统计:假设该层模板轮廓总点数为K,相似度阈值为S,以搜索框基准点N在待测图中逐像素遍历移动搜索框,判断搜索框覆盖区域中待测图边缘点数,若统计点数大于K*S,则进入第(4.3)步。否则,则认为待测图该像素上不存在以它为质心的任意角度目标,而后搜索框成“S”形在待测图中移动再次执行本步操作。若沿X轴移动,只需在上一步统计的边缘总点数基础上,加上沿移动方向向前新进入搜索框范围的一列M个像素点中的待测图边缘点数,并减去相应移出搜索框范围的一列M个像素点中的待测图边缘点数;若沿Y轴移动,只需在上一步统计的边缘总点数基础上,加上沿移动方向向前新进入搜索框范围的一行H个像素点中的待测图边缘点数,并减去相应移出搜索框范围的一行H个像素点中的待测图边缘点数。
(4.3)模板旋转与待测图对应边缘点数统计:以搜索框基准点N与模板质心对齐并以其为旋转中心,按该层旋转角度步长旋转模板二值化轮廓图,统计每旋转一个角度步长的模板边缘点对应到待测图的二值化轮廓图相应坐标下同时也存在边缘点的个数,当模板旋转至某角度该统计数大于K*S时,则进行第(4.4)步,否则,继续旋转至下一个角度。
其中轮廓点旋转采用如下公式:
(x0 x0)T为边缘点原坐标,(x1 x1)T为其旋转后坐标,θ为旋转角度步长。
(4.4)基于方向向量特征进行匹配判断:利用模板和待测图边缘点的方向向量通过如下公式,对旋转至某角度的模板与待测图进行匹配判断:
其中,di=(ai,bi)由模板边缘点上的X、Y方向向量组成,ei=(wi,ri)由待测图边缘点上的X、Y方向向量组成, 分别为方向向量对应的梯度强度值。为归一化的模板边缘点上的X、Y方向向量,为归一化的待测图边缘点上的X、Y方向向量,s表示模板与待测图的相似度,取值范围为0~1,取s=0时不匹配,随着s增加相似度提高,s=1时完全匹配。因使用轮廓线特征,只要能够提取出轮廓,即可进行图像匹配操作,因上式中两端取绝对值,即使待测图中目标极性反转,也可进行匹配。
s′T<(K-1)S+T
则这T个点相似度未达到最低阈值,即使剩余的K-T个完全相同也不可能达到模板相似度阈值S,因此结束本次相似度计算。
若在某位置匹配出了满足相似度阈值的图形,则记录此时模板质心点坐标和旋转角度作为位置信息,当遍历完整个待测图后,将这些位置信息传递到下一层金字塔中。
第五步包括以下步骤:
(5.1)根据图像金字塔上一层在待测图中定位的目标位置信息(质心点坐标(x0 y0)T,旋转角度θ0),通过图像金字塔构建逆向工程原理,计算本层目标位置信息(质心点坐标(x y)T,旋转角度θ)如下:
考虑在上述图像金字塔逆向推理过程中存在坐标和角度信息丢失计算的情况,因此需在以(x y)T为中心、周围2~4个像素范围内,以θ为中心角度、左右1~2个该层对应的角度步长范围内,进行全匹配,匹配方法直接采用步骤(4.4),以能完全匹配所有可能位置。
(5.2)由于相似度的大小、模板和待测图特征提取的边缘信息不完全一致性等原因,可能在待测图极小范围内匹配出多个目标位置,它们具有相近的质心点坐标或旋转角度,而实际上这些位置只对应一个真实目标。需根据以下两步筛选策略进行无效位置信息剔除。
当目标位置信息中在相同质心点坐标上检测有不同角度目标位置时,找出相似度最高所对应的角度,然后计算剩余角度与其相差的度数,当差值的绝对值大于10°时,则认为该角度对应的目标属于该质心点坐标下的第二个目标,保留该角度值,否则不保留该位置信息。继而再比较剩余角度与上述两个角度的差值,同样若差值的绝对值大于10°时保留为第三个认可角度,其它角度依次类推。处理完成后,即可保证每个质心点坐标上,只会对应真实的目标角度。
另外,经以上处理还需考虑临近质心点坐标对应同一真实目标的情况,首先,找出相似度最高的位置对应的质心点坐标,然后比较其与剩余位置信息的质心点坐标的关系,当二者X、Y方向坐标差值的绝对值大于一定阈值(分辨率为640*480的待测图设为12~15,其余分辨率图相应比例增减),则保留该质心点坐标对应的目标位置信息作为第二目标,否则不保留该目标位置信息。剩余目标位置信息同样与记录的两个认可位置信息进行质心坐标比较,依此类推直至处理完所有位置。
此时所得所有目标位置信息即可认为是该层图像金字塔待测图中所有真实目标对应的信息,将其传递到下一层图像金字塔。
Claims (6)
1.基于轮廓向量特征的嵌入式实时图像匹配方法,其特征是,包括以下步骤:
第一步,模板图像特征提取模块在输入的模板灰度图中,利用2*2均值法构建图像金字塔、高斯滤波、Sobel边缘提取、大津法阈值分割、非极大值抑制、高低阈值处理算法,提取用于图像匹配的模板图像金字塔各层X、Y方向向量特征图和二值化轮廓图;
第二步,待测图像特征提取模块在输入的待测灰度图中,利用均值法构建图像金字塔、高斯滤波、Sobel边缘提取、大津法阈值分割算法,提取与模板特征进行相似度吻合比较的待测图图像金字塔各层X、Y方向向量特征图和二值化轮廓图;
第三步,匹配参数自动设置模块通过获取模板二值化轮廓图中距质心最远边缘点和统计边缘点总数,计算最佳旋转角度步长、金字塔分层数;
第四步,图像金字塔最高层匹配模块分析不同待测图像的具体内容,利用三级筛选匹配策略,剔除待测图像中非目标位置,在图像金字塔最高层实现快速准确的模板粗匹配;
第五步,图像金字塔非最高层匹配模块根据上一层图像金字塔传递来的匹配信息,通过图像金字塔构建逆向工程原理,定位出本层目标位置,然后将本层定位信息传递到下一层图像金字塔进行更精确匹配,直至图像金字塔第一层。
2.如权利要求1所述基于轮廓向量特征的嵌入式实时图像匹配方法,其特征是,所述第一步具体为:
(1.1)利用2*2均值法构建多层模板图像金字塔,模板原始灰度图作为第一层,高一层金字塔图像则根据其下一层灰度图相邻四个像素的灰度值的均值求得;
(1.2)对图像金字塔各层模板灰度图利用高斯滤波去噪,高斯卷积模板如下:
继而应用Sobel边缘检测算子提取图像中目标轮廓点的X、Y方向向量及对应的梯度值,X、Y方向检测子如下:
(1.3)对图像金字塔各层模板图像对应的梯度图进行大津法全局阈值分割,剔除部分梯度较小的边缘点,以删除背景信息,保留前景目标信息;
(1.4)对经阈值分割的模板图像金字塔各层梯度图利用非极大值抑制的方法,比较边缘点梯度值与该边缘点梯度正、逆方向上相邻两个点的大小,若大于这两点梯度值,则保留该边缘点,否则去除该边缘点;
(1.5)对经非极大值抑制的模板梯度图利用高低阈值处理方法,提取模板前景目标的关键特征,去除细小、琐碎边缘,将此时模板梯度取直方图,当累加梯度值由小到大排列所对应的边缘点数达到0.4~0.6倍总边缘点数时,所对应的梯度值取为高阈值,所述高阈值的0.3~0.5倍取为低阈值;
判断经非极大值抑制的模板梯度图中梯度值与高低阈值的关系,若一点大于高阈值则保留,小于低阈值则剔除,介于高低阈值之间且与大于高阈值的点8方向连通则保留,否则剔除;
(1.6)根据上一步处理所得的模板梯度图与步骤1.2所得的X、Y方向向量图对应,获得图像金字塔各层模板X、Y方向向量图和二值化轮廓图。
3.如权利要求2所述基于轮廓向量特征的嵌入式实时图像匹配方法,其特征是,所述第二步具体为:
(2.1)利用2*2均值法构建多层待测图图像金字塔,待测图原始灰度图作为第一层,高一层金字塔图像则根据其下一层灰度图相邻四个像素的灰度值的均值求得;
(2.2)对图像金字塔各层待测灰度图利用高斯滤波去噪,高斯卷积模板如下:
继而应用Sobel边缘检测算子提取图像中目标轮廓点的X、Y方向向量及对应的梯度值,X、Y方向检测子如下:
(2.3)对图像金字塔各层待测图像对应的梯度图进行大津法全局阈值分割,剔除部分梯度较小的边缘点,以删除背景信息,保留前景目标信息;
(2.4)根据上一步处理所得的待测图的梯度图与步骤2.2所得的X、Y方向向量图对应,获得图像金字塔各层对应的待测图的X、Y方向向量图和二值化轮廓图。
4.如权利要求3所述基于轮廓向量特征的嵌入式实时图像匹配方法,其特征是,所述第三步具体为:
(3.1)根据处理过的图像金字塔各层模板二值化轮廓图计算其质心坐标及轮廓总点数,并求解距质心点欧式距离最远的模板轮廓点,以下为质心点求解:
模板某层金字塔质心坐标为(XC,XC),边缘总点数为K,边缘点X、Y方向坐标分别为XC、XC;
(3.2)在模板图像金字塔各层上,假设距离质心最远的轮廓点绕质心旋转1~2个像素距离后距离不变,则根据余弦定理求得的该点此时所转过的角度认为是该模板旋转角度步长;设点L为模板质心点,点p为模板最远轮廓点,点w、q为点p旋转过1、2个像素距离后对应的点,点L、p、w构成三角形Lpw,三边边长分别为Lp=Lw=a,pw=1,则由余弦定理求∠wLp得:
此时,∠wLp为保证模板旋转后至少有一个轮廓点旋转出原来像素点坐标所对应的最小角度,模板的旋转角度步长需保证模板旋转前后能够区分出来;
(3.3)当模板图像金字塔某一层轮廓点数大于20个,而其更高一层轮廓点数小于20时,则认为该层为图像金字塔最高层,能够进行整体有效的图像匹配工作。
5.如权利要求4所述基于轮廓向量特征的嵌入式实时图像匹配方法,其特征是,所述第四步具体为:
(4.1)构建搜索框:根据图像金字塔最高层的模板二值化轮廓图,以其质心点N为旋转基准点,构建包围模板在旋转角度范围内任意角度位置的最小外接矩形,设置旋转角度范围为0~360°,并且使该矩形四条边与原模板图像坐标轴平行,该矩形即为搜索框,纵长M个像素,横长H个像素,旋转基准点即为搜索框基准点N;
(4.2)待测图搜索框范围内边缘总点数统计:假设该层模板轮廓总点数为K,相似度阈值为S,以搜索框基准点N在待测图中逐像素遍历移动搜索框,判断搜索框覆盖区域中待测图边缘点数,若统计点数大于K*S,则进入第4.3步;否则,则认为待测图该像素上不存在以它为质心的任意角度目标,而后搜索框成“S”形在待测图中移动再次执行本步操作;若搜索框沿X轴移动,只需在上一步统计的边缘总点数基础上,加上沿移动方向向前新进入搜索框范围的一列M个像素点中的待测图边缘点数,并减去相应移出搜索框范围的一列M个像素点中的待测图边缘点数;若搜索框沿Y轴移动,只需在上一步统计的边缘总点数基础上,加上沿移动方向向前新进入搜索框范围的一行H个像素点中的待测图边缘点数,并减去相应移出搜索框范围的一行H个像素点中的待测图边缘点数;
(4.3)模板旋转与待测图对应边缘点数统计:以搜索框基准点N与模板质心对齐并以搜索框基准点N为旋转中心,按该层旋转角度步长旋转模板二值化轮廓图,统计每旋转一个角度步长的模板边缘点对应到待测图的二值化轮廓图相应坐标下同时也存在边缘点的个数,当模板旋转至某角度该统计数大于K*S时,则进行第4.4步,否则,继续旋转至下一个角度;
其中轮廓点旋转采用如下公式:
(x0 x0)T为边缘点原坐标,(x1 x1)T为边缘点旋转后坐标,θ为旋转角度步长;
(4.4)基于方向向量特征进行匹配判断:利用模板和待测图边缘点的方向向量通过如下公式,对旋转至某角度的模板与待测图进行匹配判断:
其中,di=(ai,bi)由模板边缘点上的X、Y方向向量组成,ei=(wi,ri)由待测图边缘点上的X、Y方向向量组成, 分别为方向向量对应的梯度强度值;为归一化的模板边缘点上的X、Y方向向量,为归一化的待测图边缘点上的X、Y方向向量;s表示模板与待测图的相似度,取值范围为0~1,取s=0时不匹配,随着s增加相似度提高,s=1时完全匹配;因使用轮廓线特征,只要能够提取出轮廓,即可进行图像匹配操作;
s′T<(K-1)S+T
则这T个点相似度未达到最低阈值,即使剩余的K-T个完全相同也不可能达到模板相似度阈值S,因此结束本次相似度计算;
若在某位置匹配出了满足相似度阈值的图形,则记录此时模板质心点坐标和旋转角度作为位置信息,当遍历完整个待测图后,将所述位置信息传递到下一层金字塔中。
6.如权利要求5所述基于轮廓向量特征的嵌入式实时图像匹配方法,其特征是,所述第五步具体为:
(5.1)根据图像金字塔上一层在待测图中定位的目标位置信息,即:
(质心点坐标(x0 y0)T,旋转角度θ0),
通过图像金字塔构建逆向工程原理,计算本层目标位置信息(质心点坐标(x y)T,旋转角度θ)如下:
在以(x y)T为中心、周围2~4个像素范围内,以θ为中心角度、左右1~2个该层对应的角度步长范围内,直接采用基于方向向量特征进行匹配判断的方法,即利用模板和待测图边缘点的方向向量通过如下公式,对旋转至某角度的模板与待测图进行匹配判断:
s′T<(K-1)S+T ,
则这T个点相似度未达到最低阈值,结束本次相似度计算;
若在某位置匹配出了满足相似度阈值的图形,则记录此时模板质心点坐标和旋转角度作为位置信息,当遍历完整个待测图后,将所述位置信息传递到下一层金字塔中;
(5.2)根据以下两步筛选策略进行无效位置信息剔除:
当目标位置信息中在相同质心点坐标上检测有不同角度目标位置时,找出相似度最高所对应的角度,然后计算剩余角度与其相差的度数,当差值的绝对值大于10°时,则认为该角度对应的目标属于该质心点坐标下的第二个目标,保留该角度值,否则不保留该位置信息;继而再比较剩余角度与上述两个角度的差值,同样若差值的绝对值大于10°时保留为第三个认可角度,其它角度依次类推;处理完成后,即能够保证每个质心点坐标上,只会对应真实的目标角度;
再考虑临近质心点坐标对应同一真实目标的情况,首先,找出相似度最高的位置对应的质心点坐标,然后比较该质心点坐标与剩余位置信息的质心点坐标的关系,当二者X、Y方向坐标差值的绝对值大于设定阈值,则保留该质心点坐标对应的目标位置信息作为第二目标,否则不保留该目标位置信息;剩余目标位置信息同样与记录的两个认可位置信息进行质心坐标比较,依此类推直至处理完所有位置;
此时所得所有目标位置信息即认为是该层图像金字塔待测图中所有真实目标对应的信息,将其传递到下一层图像金字塔。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2012100128010A CN102654902B (zh) | 2012-01-16 | 2012-01-16 | 基于轮廓向量特征的嵌入式实时图像匹配方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2012100128010A CN102654902B (zh) | 2012-01-16 | 2012-01-16 | 基于轮廓向量特征的嵌入式实时图像匹配方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102654902A true CN102654902A (zh) | 2012-09-05 |
CN102654902B CN102654902B (zh) | 2013-11-20 |
Family
ID=46730528
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2012100128010A Active CN102654902B (zh) | 2012-01-16 | 2012-01-16 | 基于轮廓向量特征的嵌入式实时图像匹配方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102654902B (zh) |
Cited By (63)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103278127A (zh) * | 2013-06-20 | 2013-09-04 | 爱司帝光电科技(苏州)有限公司 | 计算待测元件的轮廓数据的方法 |
CN103679702A (zh) * | 2013-11-20 | 2014-03-26 | 华中科技大学 | 一种基于图像边缘矢量的匹配方法 |
CN105261012A (zh) * | 2015-09-25 | 2016-01-20 | 上海瑞伯德智能***科技有限公司 | 基于Sobel向量的模板匹配方法 |
CN106355559A (zh) * | 2016-08-29 | 2017-01-25 | 厦门美图之家科技有限公司 | 一种图像序列的去噪方法及装置 |
CN106355573A (zh) * | 2016-08-24 | 2017-01-25 | 北京小米移动软件有限公司 | 图片中目标物的定位方法及装置 |
JP2017111638A (ja) * | 2015-12-17 | 2017-06-22 | キヤノン株式会社 | 画像処理方法、画像処理装置、画像処理システム、生産装置、プログラム及び記録媒体 |
CN107222726A (zh) * | 2014-05-12 | 2017-09-29 | 江苏省电力公司常州供电公司 | 电力设施防外破预警方案 |
CN107316315A (zh) * | 2017-05-04 | 2017-11-03 | 佛山市南海区广工大数控装备协同创新研究院 | 一种基于模板匹配的目标识别定位方法 |
CN107356610A (zh) * | 2017-09-08 | 2017-11-17 | 李尚东 | 一种识别led支架正方向的视觉检测装置和判断方法 |
CN107423770A (zh) * | 2017-08-03 | 2017-12-01 | 广东顺德西安交通大学研究院 | 一种面向高速生产线的机器人视觉定位方法 |
WO2017206099A1 (zh) * | 2016-06-01 | 2017-12-07 | 深圳配天智能技术研究院有限公司 | 一种图像模式匹配的方法及装置 |
CN107545572A (zh) * | 2016-06-23 | 2018-01-05 | 爱威科技股份有限公司 | 图像中目标的处理方法和装置 |
CN107665350A (zh) * | 2016-07-29 | 2018-02-06 | 广州康昕瑞基因健康科技有限公司 | 图像识别方法及***和自动聚焦控制方法及*** |
CN108010082A (zh) * | 2017-12-28 | 2018-05-08 | 上海觉感视觉科技有限公司 | 一种几何匹配的方法 |
CN108038495A (zh) * | 2017-12-04 | 2018-05-15 | 昆明理工大学 | 一种残缺汉字识别方法 |
CN108230327A (zh) * | 2016-12-14 | 2018-06-29 | 南京文采科技有限责任公司 | 一种基于mvp平台的包装定位与分类研究通用方法 |
CN108319979A (zh) * | 2018-02-05 | 2018-07-24 | 哈工大机器人(合肥)国际创新研究院 | 一种基于缩放和旋转匹配的图像定位识别检测方法 |
CN108399620A (zh) * | 2018-03-19 | 2018-08-14 | 嘉兴学院 | 一种基于低秩稀疏矩阵分解的图像质量评价方法 |
CN108416735A (zh) * | 2018-03-19 | 2018-08-17 | 深圳市深图医学影像设备有限公司 | 基于几何特征的数字化x线图像的拼接方法及装置 |
CN108537732A (zh) * | 2018-04-10 | 2018-09-14 | 福州大学 | 基于pca-sift的快速图像拼接方法 |
CN108550168A (zh) * | 2018-04-18 | 2018-09-18 | 湖南牛顺科技有限公司 | 智能化机械臂基于视觉的多目标快速高精度识别与定位方法 |
CN108573273A (zh) * | 2017-03-10 | 2018-09-25 | 欧姆龙株式会社 | 图像处理装置、图像处理方法、模板作成装置、物体识别处理装置以及记录介质 |
CN108701161A (zh) * | 2016-03-04 | 2018-10-23 | 微软技术许可有限责任公司 | 为搜索查询提供图像 |
CN108734706A (zh) * | 2018-05-21 | 2018-11-02 | 东南大学 | 一种融合区域分布特性与边缘尺度角度信息的转子绕线图像检测方法 |
CN108734250A (zh) * | 2018-05-29 | 2018-11-02 | 西安理工大学 | 基于Sobel算子的视觉二维码生成方法 |
CN108805220A (zh) * | 2018-07-25 | 2018-11-13 | 广东奥普特科技股份有限公司 | 一种基于梯度积分的快速模板匹配算法 |
CN108955901A (zh) * | 2017-07-06 | 2018-12-07 | 深圳市朗驰欣创科技股份有限公司 | 一种红外测温方法、***及终端设备 |
CN109034258A (zh) * | 2018-08-03 | 2018-12-18 | 厦门大学 | 基于特定物体像素梯度图的弱监督目标检测方法 |
CN109241807A (zh) * | 2018-08-17 | 2019-01-18 | 湖南大学 | 一种远距离二维码定位方法 |
CN109447083A (zh) * | 2018-08-31 | 2019-03-08 | 广州市派客朴食信息科技有限责任公司 | 一种基于图像边缘特征的餐具识别方法 |
CN109583368A (zh) * | 2018-11-28 | 2019-04-05 | 北京京东金融科技控股有限公司 | 特征比对方法及其***、计算机***及计算机可读介质 |
CN109596625A (zh) * | 2019-02-01 | 2019-04-09 | 东莞中科蓝海智能视觉科技有限公司 | 料盘中的工件缺陷检测识别方法 |
CN109801318A (zh) * | 2019-01-30 | 2019-05-24 | 东北大学 | 一种快速物体匹配算法 |
CN109883433A (zh) * | 2019-03-21 | 2019-06-14 | 中国科学技术大学 | 基于360度全景视图的结构化环境中车辆定位方法 |
CN110059101A (zh) * | 2019-04-16 | 2019-07-26 | 北京科基中意软件开发有限公司 | 一种基于图像识别的车辆数据查找***和查找方法 |
CN110136160A (zh) * | 2019-05-13 | 2019-08-16 | 南京大学 | 一种基于圆投影的快速图像匹配方法 |
CN110197232A (zh) * | 2019-06-05 | 2019-09-03 | 中科新松有限公司 | 基于边缘方向和梯度特征的图像匹配方法 |
CN110717497A (zh) * | 2019-09-06 | 2020-01-21 | 中国平安财产保险股份有限公司 | 图像相似度匹配方法、装置及计算机可读存储介质 |
CN110717546A (zh) * | 2019-10-15 | 2020-01-21 | 上海眼控科技股份有限公司 | 图片分布相似度计算方法、装置、计算机设备和存储介质 |
CN110728326A (zh) * | 2019-10-15 | 2020-01-24 | 易思维(杭州)科技有限公司 | 一种带旋转的边缘模板匹配方法 |
CN111311673A (zh) * | 2018-12-12 | 2020-06-19 | 北京京东尚科信息技术有限公司 | 一种定位方法及装置、存储介质 |
CN111415378A (zh) * | 2020-02-27 | 2020-07-14 | 湖南大学 | 用于汽车玻璃检测的图像配准方法及汽车玻璃检测方法 |
CN111553425A (zh) * | 2020-04-29 | 2020-08-18 | 广州大学 | 一种用于视觉定位的模板匹配lsp算法、介质和设备 |
CN112037134A (zh) * | 2020-09-10 | 2020-12-04 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种背景同质处理的图像拼接方法、存储介质及终端 |
CN112085033A (zh) * | 2020-08-19 | 2020-12-15 | 浙江华睿科技有限公司 | 一种模板匹配方法、装置、电子设备及存储介质 |
CN112132783A (zh) * | 2020-08-21 | 2020-12-25 | 成都飞机工业(集团)有限责任公司 | 一种基于数字图像处理技术的零件识别方法 |
CN112184785A (zh) * | 2020-09-30 | 2021-01-05 | 西安电子科技大学 | 基于mcd度量和vtm的多模态遥感图像配准方法 |
CN112418300A (zh) * | 2020-11-19 | 2021-02-26 | 重庆连芯光电技术研究院有限公司 | 一种加速图像匹配的方法 |
CN112801141A (zh) * | 2021-01-08 | 2021-05-14 | 吉林大学 | 基于模板匹配和孪生神经网络优化的异源图像匹配方法 |
CN112861983A (zh) * | 2021-02-24 | 2021-05-28 | 广东拓斯达科技股份有限公司 | 图像匹配方法、装置、电子设备以及存储介质 |
CN113111929A (zh) * | 2021-04-01 | 2021-07-13 | 广东拓斯达科技股份有限公司 | 一种模板匹配方法、装置、计算机设备及存储介质 |
CN113159103A (zh) * | 2021-02-24 | 2021-07-23 | 广东拓斯达科技股份有限公司 | 图像匹配方法、装置、电子设备以及存储介质 |
CN113158928A (zh) * | 2021-04-27 | 2021-07-23 | 浙江云奕科技有限公司 | 一种基于图像识别的混凝土试块防造假方法 |
CN113283442A (zh) * | 2020-02-19 | 2021-08-20 | 北京四维图新科技股份有限公司 | 特征点的提取方法和装置 |
CN113313766A (zh) * | 2021-03-09 | 2021-08-27 | 常州龙源智能机器人科技有限公司 | 机器人及其位姿计算方法、位姿计算装置 |
CN113658253A (zh) * | 2021-07-27 | 2021-11-16 | 广东省大湾区集成电路与***应用研究院 | 模板匹配定位方法、装置、图像采集设备及存储介质 |
CN114842212A (zh) * | 2022-04-19 | 2022-08-02 | 湖南大学 | 一种快速的多目标多角度模板匹配方法 |
CN115131587A (zh) * | 2022-08-30 | 2022-09-30 | 常州铭赛机器人科技股份有限公司 | 一种基于边缘轮廓的梯度向量特征的模板匹配方法 |
CN116188585A (zh) * | 2023-04-24 | 2023-05-30 | 成都垣景科技有限公司 | 一种基于无人机摄影测量的山区光伏目标定位方法 |
CN116386074A (zh) * | 2023-06-07 | 2023-07-04 | 青岛雅筑景观设计有限公司 | 一种园林工程设计数据智能处理及管理*** |
CN116433672A (zh) * | 2023-06-15 | 2023-07-14 | 山东九思新材料科技有限责任公司 | 基于图像处理的硅片表面质量检测方法 |
CN117274112A (zh) * | 2023-11-20 | 2023-12-22 | 张家港乐达汽车电器有限公司 | 一种智能化全景监视方法及*** |
CN113689378B (zh) * | 2021-07-07 | 2024-04-05 | 杭州未名信科科技有限公司 | 试纸带精准定位的判定方法、装置、存储介质及终端 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1691065A (zh) * | 2004-04-23 | 2005-11-02 | 中国科学院计算技术研究所 | 一种视频运动对象分割方法 |
CN101131728A (zh) * | 2007-09-29 | 2008-02-27 | 东华大学 | 一种基于Shape Context的人脸形状匹配方法 |
-
2012
- 2012-01-16 CN CN2012100128010A patent/CN102654902B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1691065A (zh) * | 2004-04-23 | 2005-11-02 | 中国科学院计算技术研究所 | 一种视频运动对象分割方法 |
CN101131728A (zh) * | 2007-09-29 | 2008-02-27 | 东华大学 | 一种基于Shape Context的人脸形状匹配方法 |
Non-Patent Citations (1)
Title |
---|
贾佳等: "嵌入式机器视觉实用彩色图像二值化方法", 《江南大学学报(自然科学版)》 * |
Cited By (97)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103278127A (zh) * | 2013-06-20 | 2013-09-04 | 爱司帝光电科技(苏州)有限公司 | 计算待测元件的轮廓数据的方法 |
CN103679702A (zh) * | 2013-11-20 | 2014-03-26 | 华中科技大学 | 一种基于图像边缘矢量的匹配方法 |
CN103679702B (zh) * | 2013-11-20 | 2016-08-31 | 华中科技大学 | 一种基于图像边缘矢量的匹配方法 |
CN107222726B (zh) * | 2014-05-12 | 2020-01-10 | 江苏省电力公司常州供电公司 | 电力设施防外破预警方法 |
CN107222726A (zh) * | 2014-05-12 | 2017-09-29 | 江苏省电力公司常州供电公司 | 电力设施防外破预警方案 |
CN105261012A (zh) * | 2015-09-25 | 2016-01-20 | 上海瑞伯德智能***科技有限公司 | 基于Sobel向量的模板匹配方法 |
JP2017111638A (ja) * | 2015-12-17 | 2017-06-22 | キヤノン株式会社 | 画像処理方法、画像処理装置、画像処理システム、生産装置、プログラム及び記録媒体 |
CN108701161A (zh) * | 2016-03-04 | 2018-10-23 | 微软技术许可有限责任公司 | 为搜索查询提供图像 |
CN108701161B (zh) * | 2016-03-04 | 2022-04-01 | 微软技术许可有限责任公司 | 为搜索查询提供图像 |
CN107851196B (zh) * | 2016-06-01 | 2020-02-14 | 深圳配天智能技术研究院有限公司 | 一种图像模式匹配的方法及装置 |
CN107851196A (zh) * | 2016-06-01 | 2018-03-27 | 深圳配天智能技术研究院有限公司 | 一种图像模式匹配的方法及装置 |
WO2017206099A1 (zh) * | 2016-06-01 | 2017-12-07 | 深圳配天智能技术研究院有限公司 | 一种图像模式匹配的方法及装置 |
CN107545572A (zh) * | 2016-06-23 | 2018-01-05 | 爱威科技股份有限公司 | 图像中目标的处理方法和装置 |
CN107665350A (zh) * | 2016-07-29 | 2018-02-06 | 广州康昕瑞基因健康科技有限公司 | 图像识别方法及***和自动聚焦控制方法及*** |
CN106355573A (zh) * | 2016-08-24 | 2017-01-25 | 北京小米移动软件有限公司 | 图片中目标物的定位方法及装置 |
CN106355559A (zh) * | 2016-08-29 | 2017-01-25 | 厦门美图之家科技有限公司 | 一种图像序列的去噪方法及装置 |
CN108230327A (zh) * | 2016-12-14 | 2018-06-29 | 南京文采科技有限责任公司 | 一种基于mvp平台的包装定位与分类研究通用方法 |
CN108573273B (zh) * | 2017-03-10 | 2022-02-18 | 欧姆龙株式会社 | 图像处理装置、图像处理方法、模板作成装置、物体识别处理装置以及记录介质 |
CN108573273A (zh) * | 2017-03-10 | 2018-09-25 | 欧姆龙株式会社 | 图像处理装置、图像处理方法、模板作成装置、物体识别处理装置以及记录介质 |
CN107316315A (zh) * | 2017-05-04 | 2017-11-03 | 佛山市南海区广工大数控装备协同创新研究院 | 一种基于模板匹配的目标识别定位方法 |
CN108955901A (zh) * | 2017-07-06 | 2018-12-07 | 深圳市朗驰欣创科技股份有限公司 | 一种红外测温方法、***及终端设备 |
CN108955901B (zh) * | 2017-07-06 | 2020-03-31 | 深圳市朗驰欣创科技股份有限公司 | 一种红外测温方法、***及终端设备 |
CN107423770A (zh) * | 2017-08-03 | 2017-12-01 | 广东顺德西安交通大学研究院 | 一种面向高速生产线的机器人视觉定位方法 |
CN107356610A (zh) * | 2017-09-08 | 2017-11-17 | 李尚东 | 一种识别led支架正方向的视觉检测装置和判断方法 |
CN107356610B (zh) * | 2017-09-08 | 2023-08-22 | 深圳市祥云达光电设备有限公司 | 一种识别led支架正方向的视觉检测装置和判断方法 |
CN108038495A (zh) * | 2017-12-04 | 2018-05-15 | 昆明理工大学 | 一种残缺汉字识别方法 |
CN108038495B (zh) * | 2017-12-04 | 2021-08-20 | 昆明理工大学 | 一种残缺汉字识别方法 |
CN108010082B (zh) * | 2017-12-28 | 2021-07-23 | 上海觉感视觉科技有限公司 | 一种几何匹配的方法 |
CN108010082A (zh) * | 2017-12-28 | 2018-05-08 | 上海觉感视觉科技有限公司 | 一种几何匹配的方法 |
CN108319979A (zh) * | 2018-02-05 | 2018-07-24 | 哈工大机器人(合肥)国际创新研究院 | 一种基于缩放和旋转匹配的图像定位识别检测方法 |
CN108399620B (zh) * | 2018-03-19 | 2022-04-19 | 嘉兴学院 | 一种基于低秩稀疏矩阵分解的图像质量评价方法 |
CN108416735A (zh) * | 2018-03-19 | 2018-08-17 | 深圳市深图医学影像设备有限公司 | 基于几何特征的数字化x线图像的拼接方法及装置 |
CN108399620A (zh) * | 2018-03-19 | 2018-08-14 | 嘉兴学院 | 一种基于低秩稀疏矩阵分解的图像质量评价方法 |
CN108537732B (zh) * | 2018-04-10 | 2021-11-02 | 福州大学 | 基于pca-sift的快速图像拼接方法 |
CN108537732A (zh) * | 2018-04-10 | 2018-09-14 | 福州大学 | 基于pca-sift的快速图像拼接方法 |
CN108550168A (zh) * | 2018-04-18 | 2018-09-18 | 湖南牛顺科技有限公司 | 智能化机械臂基于视觉的多目标快速高精度识别与定位方法 |
CN108734706B (zh) * | 2018-05-21 | 2022-07-19 | 东南大学 | 一种融合区域分布特性与边缘尺度角度信息的转子绕线图像检测方法 |
CN108734706A (zh) * | 2018-05-21 | 2018-11-02 | 东南大学 | 一种融合区域分布特性与边缘尺度角度信息的转子绕线图像检测方法 |
CN108734250A (zh) * | 2018-05-29 | 2018-11-02 | 西安理工大学 | 基于Sobel算子的视觉二维码生成方法 |
CN108734250B (zh) * | 2018-05-29 | 2021-06-15 | 西安理工大学 | 基于Sobel算子的视觉二维码生成方法 |
CN108805220A (zh) * | 2018-07-25 | 2018-11-13 | 广东奥普特科技股份有限公司 | 一种基于梯度积分的快速模板匹配算法 |
CN109034258A (zh) * | 2018-08-03 | 2018-12-18 | 厦门大学 | 基于特定物体像素梯度图的弱监督目标检测方法 |
CN109241807A (zh) * | 2018-08-17 | 2019-01-18 | 湖南大学 | 一种远距离二维码定位方法 |
CN109241807B (zh) * | 2018-08-17 | 2021-04-20 | 湖南大学 | 一种远距离二维码定位方法 |
CN109447083A (zh) * | 2018-08-31 | 2019-03-08 | 广州市派客朴食信息科技有限责任公司 | 一种基于图像边缘特征的餐具识别方法 |
CN109583368A (zh) * | 2018-11-28 | 2019-04-05 | 北京京东金融科技控股有限公司 | 特征比对方法及其***、计算机***及计算机可读介质 |
CN111311673A (zh) * | 2018-12-12 | 2020-06-19 | 北京京东尚科信息技术有限公司 | 一种定位方法及装置、存储介质 |
CN111311673B (zh) * | 2018-12-12 | 2023-11-03 | 北京京东乾石科技有限公司 | 一种定位方法及装置、存储介质 |
CN109801318A (zh) * | 2019-01-30 | 2019-05-24 | 东北大学 | 一种快速物体匹配算法 |
CN109596625A (zh) * | 2019-02-01 | 2019-04-09 | 东莞中科蓝海智能视觉科技有限公司 | 料盘中的工件缺陷检测识别方法 |
CN109883433A (zh) * | 2019-03-21 | 2019-06-14 | 中国科学技术大学 | 基于360度全景视图的结构化环境中车辆定位方法 |
CN109883433B (zh) * | 2019-03-21 | 2023-07-18 | 中国科学技术大学 | 基于360度全景视图的结构化环境中车辆定位方法 |
CN110059101B (zh) * | 2019-04-16 | 2021-08-13 | 北京科基中意软件开发有限公司 | 一种基于图像识别的车辆数据查找***和查找方法 |
CN110059101A (zh) * | 2019-04-16 | 2019-07-26 | 北京科基中意软件开发有限公司 | 一种基于图像识别的车辆数据查找***和查找方法 |
CN110136160A (zh) * | 2019-05-13 | 2019-08-16 | 南京大学 | 一种基于圆投影的快速图像匹配方法 |
CN110197232A (zh) * | 2019-06-05 | 2019-09-03 | 中科新松有限公司 | 基于边缘方向和梯度特征的图像匹配方法 |
CN110717497B (zh) * | 2019-09-06 | 2023-11-07 | 中国平安财产保险股份有限公司 | 图像相似度匹配方法、装置及计算机可读存储介质 |
CN110717497A (zh) * | 2019-09-06 | 2020-01-21 | 中国平安财产保险股份有限公司 | 图像相似度匹配方法、装置及计算机可读存储介质 |
CN110717546A (zh) * | 2019-10-15 | 2020-01-21 | 上海眼控科技股份有限公司 | 图片分布相似度计算方法、装置、计算机设备和存储介质 |
CN110728326A (zh) * | 2019-10-15 | 2020-01-24 | 易思维(杭州)科技有限公司 | 一种带旋转的边缘模板匹配方法 |
CN113283442A (zh) * | 2020-02-19 | 2021-08-20 | 北京四维图新科技股份有限公司 | 特征点的提取方法和装置 |
CN113283442B (zh) * | 2020-02-19 | 2024-05-10 | 北京四维图新科技股份有限公司 | 特征点的提取方法和装置 |
CN111415378A (zh) * | 2020-02-27 | 2020-07-14 | 湖南大学 | 用于汽车玻璃检测的图像配准方法及汽车玻璃检测方法 |
CN111415378B (zh) * | 2020-02-27 | 2023-11-24 | 湖南大学 | 用于汽车玻璃检测的图像配准方法及汽车玻璃检测方法 |
CN111553425B (zh) * | 2020-04-29 | 2023-04-25 | 广州大学 | 一种用于视觉定位的模板匹配lsp算法、介质和设备 |
CN111553425A (zh) * | 2020-04-29 | 2020-08-18 | 广州大学 | 一种用于视觉定位的模板匹配lsp算法、介质和设备 |
CN112085033A (zh) * | 2020-08-19 | 2020-12-15 | 浙江华睿科技有限公司 | 一种模板匹配方法、装置、电子设备及存储介质 |
CN112085033B (zh) * | 2020-08-19 | 2024-04-09 | 浙江华睿科技股份有限公司 | 一种模板匹配方法、装置、电子设备及存储介质 |
CN112132783A (zh) * | 2020-08-21 | 2020-12-25 | 成都飞机工业(集团)有限责任公司 | 一种基于数字图像处理技术的零件识别方法 |
CN112132783B (zh) * | 2020-08-21 | 2022-09-20 | 成都飞机工业(集团)有限责任公司 | 一种基于数字图像处理技术的零件识别方法 |
CN112037134B (zh) * | 2020-09-10 | 2023-04-21 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种背景同质处理的图像拼接方法、存储介质及终端 |
CN112037134A (zh) * | 2020-09-10 | 2020-12-04 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种背景同质处理的图像拼接方法、存储介质及终端 |
CN112184785A (zh) * | 2020-09-30 | 2021-01-05 | 西安电子科技大学 | 基于mcd度量和vtm的多模态遥感图像配准方法 |
CN112184785B (zh) * | 2020-09-30 | 2023-03-24 | 西安电子科技大学 | 基于mcd度量和vtm的多模态遥感图像配准方法 |
CN112418300A (zh) * | 2020-11-19 | 2021-02-26 | 重庆连芯光电技术研究院有限公司 | 一种加速图像匹配的方法 |
CN112801141A (zh) * | 2021-01-08 | 2021-05-14 | 吉林大学 | 基于模板匹配和孪生神经网络优化的异源图像匹配方法 |
CN112861983A (zh) * | 2021-02-24 | 2021-05-28 | 广东拓斯达科技股份有限公司 | 图像匹配方法、装置、电子设备以及存储介质 |
CN113159103B (zh) * | 2021-02-24 | 2023-12-05 | 广东拓斯达科技股份有限公司 | 图像匹配方法、装置、电子设备以及存储介质 |
CN113159103A (zh) * | 2021-02-24 | 2021-07-23 | 广东拓斯达科技股份有限公司 | 图像匹配方法、装置、电子设备以及存储介质 |
WO2022179002A1 (zh) * | 2021-02-24 | 2022-09-01 | 广东拓斯达科技股份有限公司 | 图像匹配方法、装置、电子设备以及存储介质 |
CN113313766A (zh) * | 2021-03-09 | 2021-08-27 | 常州龙源智能机器人科技有限公司 | 机器人及其位姿计算方法、位姿计算装置 |
CN113111929A (zh) * | 2021-04-01 | 2021-07-13 | 广东拓斯达科技股份有限公司 | 一种模板匹配方法、装置、计算机设备及存储介质 |
CN113111929B (zh) * | 2021-04-01 | 2024-04-12 | 广东拓斯达科技股份有限公司 | 一种模板匹配方法、装置、计算机设备及存储介质 |
CN113158928A (zh) * | 2021-04-27 | 2021-07-23 | 浙江云奕科技有限公司 | 一种基于图像识别的混凝土试块防造假方法 |
CN113158928B (zh) * | 2021-04-27 | 2023-09-19 | 浙江云奕科技有限公司 | 一种基于图像识别的混凝土试块防造假方法 |
CN113689378B (zh) * | 2021-07-07 | 2024-04-05 | 杭州未名信科科技有限公司 | 试纸带精准定位的判定方法、装置、存储介质及终端 |
CN113658253A (zh) * | 2021-07-27 | 2021-11-16 | 广东省大湾区集成电路与***应用研究院 | 模板匹配定位方法、装置、图像采集设备及存储介质 |
CN114842212A (zh) * | 2022-04-19 | 2022-08-02 | 湖南大学 | 一种快速的多目标多角度模板匹配方法 |
CN114842212B (zh) * | 2022-04-19 | 2024-06-25 | 湖南大学 | 一种快速的多目标多角度模板匹配方法 |
CN115131587A (zh) * | 2022-08-30 | 2022-09-30 | 常州铭赛机器人科技股份有限公司 | 一种基于边缘轮廓的梯度向量特征的模板匹配方法 |
CN116188585A (zh) * | 2023-04-24 | 2023-05-30 | 成都垣景科技有限公司 | 一种基于无人机摄影测量的山区光伏目标定位方法 |
CN116386074B (zh) * | 2023-06-07 | 2023-08-15 | 青岛雅筑景观设计有限公司 | 一种园林工程设计数据智能处理及管理*** |
CN116386074A (zh) * | 2023-06-07 | 2023-07-04 | 青岛雅筑景观设计有限公司 | 一种园林工程设计数据智能处理及管理*** |
CN116433672B (zh) * | 2023-06-15 | 2023-08-25 | 山东九思新材料科技有限责任公司 | 基于图像处理的硅片表面质量检测方法 |
CN116433672A (zh) * | 2023-06-15 | 2023-07-14 | 山东九思新材料科技有限责任公司 | 基于图像处理的硅片表面质量检测方法 |
CN117274112A (zh) * | 2023-11-20 | 2023-12-22 | 张家港乐达汽车电器有限公司 | 一种智能化全景监视方法及*** |
CN117274112B (zh) * | 2023-11-20 | 2024-02-06 | 张家港乐达汽车电器有限公司 | 一种智能化全景监视方法及*** |
Also Published As
Publication number | Publication date |
---|---|
CN102654902B (zh) | 2013-11-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102654902B (zh) | 基于轮廓向量特征的嵌入式实时图像匹配方法 | |
CN104637073B (zh) | 一种基于太阳照射阴影补偿的带状地下结构探测方法 | |
CN107067415B (zh) | 一种基于图像匹配的目标定位方法 | |
CN105335973B (zh) | 运用于带钢加工生产线的视觉处理方法 | |
CN103714541B (zh) | 一种利用山体轮廓区域约束识别定位建筑物的方法 | |
CN107092871B (zh) | 基于多尺度多特征融合的遥感影像建筑物检测方法 | |
CN107330376A (zh) | 一种车道线识别方法及*** | |
CN105865344A (zh) | 一种基于机器视觉的工件尺寸测量方法和装置 | |
CN103839265A (zh) | 基于sift和归一化互信息的sar图像配准方法 | |
CN103218787B (zh) | 多源异构遥感影像控制点自动采集方法 | |
CN105761271A (zh) | 核安全壳表面缺陷自动检测方法及*** | |
CN102865859B (zh) | 一种基于surf特征的航空序列图像位置估计方法 | |
CN107230203A (zh) | 基于人眼视觉注意机制的铸件缺陷识别方法 | |
CN103136525B (zh) | 一种利用广义Hough变换的异型扩展目标高精度定位方法 | |
CN104536009A (zh) | 一种激光红外复合的地面建筑物识别及导航方法 | |
CN102903109B (zh) | 一种光学影像和sar影像一体化分割配准方法 | |
CN103400365A (zh) | 肺区ct序列自动分割方法 | |
CN103325126A (zh) | 一种在尺度变化和遮挡情况下的视频目标跟踪方法 | |
CN108596165A (zh) | 基于无人机低空航拍图像的道路交通标线检测方法及*** | |
CN107451982A (zh) | 一种基于无人机影像的高郁闭度林分树冠面积获取方法 | |
CN103065135A (zh) | 基于数字图像处理的车牌号码匹配算法 | |
CN109523516A (zh) | 一种基于双重约束条件的对象级土地覆盖变化检测方法 | |
CN101777181A (zh) | 基于脊波双框架***的sar图像机场跑道提取方法 | |
CN107727658A (zh) | 基于图像处理的受电弓裂纹在线检测装置及方法 | |
CN103198319A (zh) | 用于矿山井筒环境下的模糊图像角点提取方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C53 | Correction of patent of invention or patent application | ||
CB02 | Change of applicant information |
Address after: 1800 No. 214122 Jiangsu city of Wuxi Province Li Lake Avenue Applicant after: Jiangnan University Applicant after: Xinje Electronic Co., Ltd. Address before: 1800 No. 214122 Jiangsu city of Wuxi Province Li Lake Avenue Applicant before: Jiangnan University Applicant before: Wuxi Czech Automation Co., Ltd. |
|
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |