CN101770583A - 一种基于场景全局特征的模板匹配方法 - Google Patents
一种基于场景全局特征的模板匹配方法 Download PDFInfo
- Publication number
- CN101770583A CN101770583A CN201010028992A CN201010028992A CN101770583A CN 101770583 A CN101770583 A CN 101770583A CN 201010028992 A CN201010028992 A CN 201010028992A CN 201010028992 A CN201010028992 A CN 201010028992A CN 101770583 A CN101770583 A CN 101770583A
- Authority
- CN
- China
- Prior art keywords
- template
- pixel
- image
- scene
- subgraph
- 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
本发明公开了一种基于场景全局特征的模板匹配方法,用于根据目标模板检测图像中同类目标。该方法对场景图像逐像素点遍历提取与模板图像相同区域范围的搜索子图,计算每个搜索子图与模板图像的相关性,依据搜索子图与模板图像的相关性越大,该搜索子图对应的遍历像素点成为目标点可能性越大的准则选取候选目标像素点,在候选目标像素点中确定最终目标点。本发明关键之处在于利用全局特征来度量搜索子图与模板图像的相关性,全局特征是通过把模板和子图的相关性放到场景后计算整个场景的特征得到的,全局特征是整个场景的特征,带有一定的语义信息,使得目标检测的准确性提高。
Description
技术领域
本发明属于图像处理方法,具体涉及一种基于场景全局特征的模板匹配方法,运用于根据目标模板检测图像中同类目标。
背景技术
模板匹配是图像处理中经典的方法,其中模板就是一幅已知的小图像,在目标检测任务中模板是一个目标的实例;模板匹配就是在一幅大图像中搜索目标,确定其坐标位置,在目标检测任务中就是要找到所有的该类目标的实例,这个大图像称为场景图像。比如,我们在执行汽车检测的任务,模板就是一幅汽车的图像,模板匹配就是根据这个模板的信息在场景图像中找到所有的汽车,并确定其位置。
模板匹配的搜索过程就是将模板叠放在大图像上平移,模板覆盖被搜索图的那个区域叫做子图。传统的模板匹配方法是基于灰度值的,即通过度量模板与子图之间的灰度相关性来确定目标的位置,这种方法有很大的局限性,比如光照变化等情况对于图像灰度值的影响是很大的。因此有很多基于特征的模板匹配方法,利用模板与子图的特征之间的相关性来确定目标位置,因为特征的鲁棒性,在一定程度上提高了匹配的性能。基于特征的模板匹配方法已经在目标检测任务中得到广泛的应用。
模板匹配的方法非常简单,同时又非常的有效,但是模板匹配在目标检测时还有很多地方可以改进。模板匹配的性能很大程度上依赖于特征的性能,特征的分类性能越好模板匹配的效果也越好。以上所提及的方法都是基于局部特征的,而局部特征在目标检测时也有一定的局限性,比如目标被遮挡、目标姿态变化、同类目标的不同外观,目标视点变化和成像条件很差等情况。
发明内容
本发明的目的是提出一种基于全局特征的模板匹配方法,把模板与子图放到场景上下文中来度量模板与子图的相关性,提高目标检测的准确性。
一种基于全局特征的模板匹配方法,具体步骤为:
(1)计算场景图像的特征H(k),1≤k≤K0,K0为特征维度;
(2)对场景图像逐像素点遍历提取与模板图像相同区域范围的搜索子图,计算每个搜索子图与模板图像的相关性;
相关性计算方法具体为:在遍历像素点(i,j),将其对应搜索子图S(i,j)的内容替换为模板图像中的内容,并计算替换后对应场景图像的全局特征Gi,j(k),计算搜索子图S(i,j)与模板图像的相关性:
(3)在场景图像中,依据搜索子图与模板图像的相关性越大,该搜索子图对应的遍历像素点成为目标点可能性越大的准则选取候选目标像素点;
(4)对于每一个候选目标像素点,在场景图像中以其为中心选取局部感兴趣区,若局部感兴趣区内该候选目标像素点对应的相关性大于该区内其它像素点对应的相关性,则确认该候选目标点为最终目标点。
本发明采用的搜索过程中,利用全局特征来度量模板和子图的相关性,其中全局特征是通过把模板和子图的相关性放到场景后计算整个场景的特征得到的,基于全局特征的模板匹配方法比传统的使用局部特征(模板和子图范围内的特征)的模板匹配方法可以更准确的检测到目标,本发明中使用的全局特征是整个场景的特征,带有一定的语义信息,使得目标检测的准确性提高;本发明采用的确定位置的过程充分利用了周围的相关性的信息,得到了更精确的目标位置信息。
附图说明
图1是本发明流程图;
图2是模板图像;
图3是场景图像;
图4搜索过程中,位置为(20,30)时的示意图;
图5是位置为(20,30)时替换后的场景图像;
图6是将相关性结果显示为位图;
图7是对相关性结果阈值分割后结果;
图8是检测结果图像;
图9是梯度方向角量化示意图;
图10是图像分为两层的示意图。
具体实施方式
本发明步骤流程如图1所示,现举例说明:
给出实际的模板图像T(m,n)和场景图像F(x,y),分别如图2和图3所示,其中1≤m≤M,1≤n≤N,1≤x≤W,1≤y≤H,模板图像的宽度M=138和高度N=69,场景图像的宽度W=256和高度H=256。
(1)计算场景图像的特征
计算整个场景的特征,一般采用统计低级特征的方法,这里的的低级特征包括灰度值,边缘,梯度幅值,梯度方向,滤波器(Gabor滤波器,DoG滤波器等)等。本发明中使用梯度方向的统计,其统计结果表示为直方图H(k),其中1≤k≤K0,特征的维度K0=80;
(2)对场景图像逐像素点遍历提取搜索子图,计算各搜索子图与模板图像的相关性。
在每个搜索像素点位置(i,j),计算模板图像与搜索子图Si,j(m,n)的相关性D(i,j),其中(i,j)为搜索子图左上角像素点在场景图像F(x,y)上的坐标,搜索范围为1≤i≤W-M=118,1≤j≤H-N=187,搜索子图大小与模板大小相同,即为138×69,图4给出了在256×256的场景中搜索的示意图,此时搜索的像素点位置(i,j)=(20,30)。
(2.1)把子图Si,j(m,n)的内容替换为模板T(m,n),称为替换后场景图像,记为Ci,j(x,y),i=20,j=30时的替换后场景图像如图5所示;
(2.2)计算替换后场景图像Ci,j(x,y)的全局特征,这里我们同样计算场景图像的梯度方向直方图Gi,j(k);
(2.3)通过计算替换前场景图像和替换后场景图像的全局特征的相关性来度量子图Si,j(m,n)与模板T(m,n)的相关性,计算公式为:
搜索完成后得到的D(i,j)是一个118×187的矩阵。将D(i,j)归一化到0~255,显示为位图,如图6所示。
(3)选取候选目标点。
依据相关性结果对场景图像的像素点进行阈值分割,以突出目标存在可能性大的像素点,将其作为候选目标点。实例中具体操作方法为:若D(i,j)大于等于相关性阈值r,则场景图像中像素(i,j)的值置为1,反之(i,j)的值置为零。相关性阈值r的选择依据为:r=0.95×(Vmax-Vmin)+Vmin,其中Vmax和Vmin分别为D(i,j)的最大值和最小值。实例中Vmax=0.9579和Vmin=0.7962,所以r=0.95×(Vmax-Vmin)+Vmin=0.9498。将D(i,j)大于等于0.9498的位置对应场景图像中的位置的像素值置为1,小于0.9498的位置的像素值置为零,分割结果显示如图7所示,右下方的缺口白色区域内的像素点即为候选目标点。
(4)确认最终目标点。
为了精确的确定目标的位置,确认最终的目标点需要考虑到候选目标邻域的信息,该邻域称为感兴趣区。感兴趣区的宽度和高度的范围分别为[M/2,M]和[N/2,N],本实例中使用的感兴趣区的宽度和高度分别为M/2=69和N/2=34。对于每一个候选目标点逐像素搜索,以其为中心选取局部感兴趣区,如果该候选目标点对应的相关性是该感兴趣区内的最大值点,说明这个候选目标点有更大的可能成为最终目标点。
为了排除孤立的噪声点的干扰,在极大值判断的基础上进一步计算该极值点的邻域支持度,即该局部感兴趣内候选目标点的个数。实例中,在上述的69×34的窗口中只有一个最大值点,位置为(118,153),其邻域支持度为148;对于该实例中得到的唯一的最大值点,其邻域支持度大于支持度阈值q=(M/2)×(N/2)/100=23,因此该位置即是最后的目标位置,图8标出了最后的检测结果。
本发明所依据的全局特征梯度方向直方图的计算原理如下:图像采用灰度值F(x,y)表征,(x,y)表示图像内像素点坐标
1)计算图像F(x,y)各点的梯度
当(x,y)是图像F(x,y)最右边一列或者最下边一行时,F(x,y)在该位
置的梯度μ(x,y)=0;
否则,图像F(x,y)在(x,y)位置的梯度定义为下列向量:
2)计算F(x,y)各点的梯度方向角
用α(x,y)表示向量μ(x,y)在(x,y)处的梯度方向角:
3)量化梯度方向角
(4)计算梯度方向直方图
对于图像中除去最右一列和最下一行所有位置的梯度方向量化值统计每个值出现的个数,即得到一个16维向量。
本发明采用分层描述的特征,将图像F(x,y)分为两层,如图10所示,图10(a)表示第一层,是整个图像,表示为FL1,图10(b)表示第二层,平均分为四部分,分别为FL21,FL22,FL23和FL23。采用上述的梯度直方图的计算方法,对两层总共五个部分计算各自的梯度直方图,然后依次组合为一个向量,即为该图像的全局特征H(k),1≤k≤K0,K0=16×5=80表示特征的维度。
Claims (4)
1.一种基于全局特征的模板匹配方法,具体步骤为:
(1)计算场景图像的特征H(k),1≤k≤K0,K0为特征维度;
(2)对场景图像逐像素点遍历提取与模板图像相同区域范围的搜索子图,计算每个搜索子图与模板图像的相关性;
相关性计算方法具体为:在遍历像素点(i,j),将其对应搜索子图S(i,j)的内容替换为模板图像中的内容,并计算替换后对应场景图像的全局特征Gi,j(k),计算搜索子图S(i,j)与模板图像的相关性:
(3)在场景图像中,依据搜索子图与模板图像的相关性越大,该搜索子图对应的遍历像素点成为目标点可能性越大的准则选取候选目标像素点;
(4)对于每一个候选目标像素点,在场景图像中以其为中心选取局部感兴趣区,若局部感兴趣区内该候选目标像素点对应的相关性大于该区内其它像素点对应的相关性,则确认该候选目标点为最终目标点。
2.根据权利要求1所述的基于全局特征的模板匹配方法,其特征在于,所述步骤(3)具体为:若搜索子图与模板图像的相关性D(i,j)大于等于相关性阈值r,则像素点(i,j)为候选目标点,相关性阈值r=0.95×(Vmax-Vmin)+Vmin,Vmax和Vmin分别为各搜索子图与模板图像的相关性中的最大值和最小值。
3.根据权利要求1或2所述的基于全局特征的模板匹配方法,其特征在于,还包括步骤(5):计算满足步骤(4)的候选目标点对应的感兴趣区内包含的候选目标点的个数,若其大于支持度阈值,确认该候选目标点为最终目标点,其中支持度阈值q=(M/2)×(N/2)/100,M和N分别为模板图像的宽度和高度。
4.根据权利要求1或2所述的基于全局特征的模板匹配方法,其特征在于,所述局部感兴趣区的宽度和高度范围分别为[M/2,M]和[N/2,N],M和N分别为模板图像的宽度和高度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010100289920A CN101770583B (zh) | 2010-01-15 | 2010-01-15 | 一种基于场景全局特征的模板匹配方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010100289920A CN101770583B (zh) | 2010-01-15 | 2010-01-15 | 一种基于场景全局特征的模板匹配方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101770583A true CN101770583A (zh) | 2010-07-07 |
CN101770583B CN101770583B (zh) | 2012-06-13 |
Family
ID=42503435
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2010100289920A Expired - Fee Related CN101770583B (zh) | 2010-01-15 | 2010-01-15 | 一种基于场景全局特征的模板匹配方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101770583B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103020613A (zh) * | 2013-01-07 | 2013-04-03 | 信帧电子技术(北京)有限公司 | 一种基于视频的信号灯识别方法及装置 |
CN103263278A (zh) * | 2013-01-23 | 2013-08-28 | 郑末晶 | 一种从超声图像上自动测量胎儿颈背透明物厚度的图像处理方法 |
CN105469052A (zh) * | 2015-11-25 | 2016-04-06 | 东方网力科技股份有限公司 | 一种车辆检测跟踪方法和装置 |
WO2018068304A1 (zh) * | 2016-10-14 | 2018-04-19 | 深圳配天智能技术研究院有限公司 | 一种图像匹配的方法及装置 |
CN110309785A (zh) * | 2019-07-03 | 2019-10-08 | 孙启城 | 一种基于图像识别技术的导盲机器人控制方法 |
CN111951211A (zh) * | 2019-05-17 | 2020-11-17 | 株式会社理光 | 一种目标检测方法、装置及计算机可读存储介质 |
CN112085033A (zh) * | 2020-08-19 | 2020-12-15 | 浙江华睿科技有限公司 | 一种模板匹配方法、装置、电子设备及存储介质 |
CN112241975A (zh) * | 2019-07-17 | 2021-01-19 | 杭州海康机器人技术有限公司 | 特征模板的匹配定位方法和匹配定位装置 |
CN114022510A (zh) * | 2021-10-20 | 2022-02-08 | 浙江大立科技股份有限公司 | 一种基于内容检索的目标长时跟踪方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101488224B (zh) * | 2008-01-16 | 2011-01-19 | 中国科学院自动化研究所 | 基于相关性度量的特征点匹配方法 |
-
2010
- 2010-01-15 CN CN2010100289920A patent/CN101770583B/zh not_active Expired - Fee Related
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103020613A (zh) * | 2013-01-07 | 2013-04-03 | 信帧电子技术(北京)有限公司 | 一种基于视频的信号灯识别方法及装置 |
CN103263278A (zh) * | 2013-01-23 | 2013-08-28 | 郑末晶 | 一种从超声图像上自动测量胎儿颈背透明物厚度的图像处理方法 |
CN103263278B (zh) * | 2013-01-23 | 2015-05-13 | 珠海艾博罗生物技术有限公司 | 一种从超声图像上自动测量胎儿颈背透明物厚度的图像处理方法 |
CN105469052B (zh) * | 2015-11-25 | 2018-10-19 | 东方网力科技股份有限公司 | 一种车辆检测跟踪方法和装置 |
CN105469052A (zh) * | 2015-11-25 | 2016-04-06 | 东方网力科技股份有限公司 | 一种车辆检测跟踪方法和装置 |
CN109348731B (zh) * | 2016-10-14 | 2022-05-17 | 深圳配天智能技术研究院有限公司 | 一种图像匹配的方法及装置 |
CN109348731A (zh) * | 2016-10-14 | 2019-02-15 | 深圳配天智能技术研究院有限公司 | 一种图像匹配的方法及装置 |
WO2018068304A1 (zh) * | 2016-10-14 | 2018-04-19 | 深圳配天智能技术研究院有限公司 | 一种图像匹配的方法及装置 |
CN111951211A (zh) * | 2019-05-17 | 2020-11-17 | 株式会社理光 | 一种目标检测方法、装置及计算机可读存储介质 |
CN111951211B (zh) * | 2019-05-17 | 2024-05-14 | 株式会社理光 | 一种目标检测方法、装置及计算机可读存储介质 |
CN110309785A (zh) * | 2019-07-03 | 2019-10-08 | 孙启城 | 一种基于图像识别技术的导盲机器人控制方法 |
CN110309785B (zh) * | 2019-07-03 | 2023-10-20 | 孙启城 | 一种基于图像识别技术的导盲机器人控制方法 |
CN112241975A (zh) * | 2019-07-17 | 2021-01-19 | 杭州海康机器人技术有限公司 | 特征模板的匹配定位方法和匹配定位装置 |
CN112241975B (zh) * | 2019-07-17 | 2024-03-29 | 杭州海康机器人股份有限公司 | 特征模板的匹配定位方法和匹配定位装置 |
CN112085033A (zh) * | 2020-08-19 | 2020-12-15 | 浙江华睿科技有限公司 | 一种模板匹配方法、装置、电子设备及存储介质 |
CN112085033B (zh) * | 2020-08-19 | 2024-04-09 | 浙江华睿科技股份有限公司 | 一种模板匹配方法、装置、电子设备及存储介质 |
CN114022510A (zh) * | 2021-10-20 | 2022-02-08 | 浙江大立科技股份有限公司 | 一种基于内容检索的目标长时跟踪方法 |
Also Published As
Publication number | Publication date |
---|---|
CN101770583B (zh) | 2012-06-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101770583B (zh) | 一种基于场景全局特征的模板匹配方法 | |
Li et al. | Automatic pavement crack detection by multi-scale image fusion | |
CN108510467B (zh) | 基于深度可变形卷积神经网络的sar图像目标识别方法 | |
Yin et al. | Hot region selection based on selective search and modified fuzzy C-means in remote sensing images | |
CN110287826B (zh) | 一种基于注意力机制的视频目标检测方法 | |
CN103530600B (zh) | 复杂光照下的车牌识别方法及系统 | |
CN110598609A (zh) | 一种基于显著性引导的弱监督目标检测方法 | |
CN105279772B (zh) | 一种红外序列图像的可跟踪性判别方法 | |
CN109460764B (zh) | 一种结合亮度特征与改进帧间差分法的卫星视频船舶监测方法 | |
CN102176208B (zh) | 基于三维空时特征的鲁棒视频指纹方法 | |
CN107784663A (zh) | 基于深度信息的相关滤波跟踪方法及装置 | |
CN103871039B (zh) | 一种sar图像变化检测差异图生成方法 | |
CN104657980A (zh) | 一种改进的基于Meanshift的多通道图像分割算法 | |
CN105976376B (zh) | 一种基于部件模型的高分辨sar图像目标检测方法 | |
CN105225226A (zh) | 一种基于图像分割的级联可变形部件模型目标检测方法 | |
CN107092890A (zh) | 基于红外视频的舰船检测及跟踪方法 | |
CN110008900B (zh) | 一种由区域到目标的可见光遥感图像候选目标提取方法 | |
CN110866926B (zh) | 一种红外遥感图像快速精细海陆分割方法 | |
CN104182985A (zh) | 遥感图像变化检测方法 | |
CN111208479B (zh) | 一种降低深度网络检测中虚警概率的方法 | |
CN106815323A (zh) | 一种基于显著性检测的跨域视觉检索方法 | |
CN110633727A (zh) | 基于选择性搜索的深度神经网络舰船目标细粒度识别方法 | |
CN112308883A (zh) | 一种基于可见光和红外图像的多船舶融合跟踪方法 | |
CN111833353B (zh) | 一种基于图像分割的高光谱目标检测方法 | |
JP2009163682A (ja) | 画像識別装置及びプログラム |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20120613 Termination date: 20160115 |
|
EXPY | Termination of patent right or utility model |