CN103400393B - 一种图像匹配方法和系统 - Google Patents

Abstract

本申请提供一种图像匹配方法和系统。该图像匹配方法得到当前时刻一种尺度下可能包含指定物体的搜索区域图像，该图像的中心点为前一时刻指定物体的目标图像的中心点，然后以所述目标图像为划分模板划分所述搜索区域图像，得到与所述目标图像的横向分辨率和纵向分辨率相同的至少两个候选图像，从多个候选图像中选出相似度最大的候选图像，依据目标图像和相似度最大的候选图像的弹性参数调整相似度最大的候选图像，将调整后的候选图像作为当前时刻指定物体的目标图像。因此在得到当前时刻指定物体的目标图像过程中，可以将前一时刻指定物体的目标图像和当前时刻一种尺度下指定物体的图像进行匹配，降低匹配复杂度和匹配数据量，从而提高匹配实时性。

Description

一种图像匹配方法和系统

技术领域

本申请涉及图像识别技术领域，特别涉及一种图像匹配方法和系统。

背景技术

在视频分析过程中，快速准确的跟踪视频中指定物体，具有重要的意义。随着空间位置和形态的变化，指定物体的颜色、纹理、尺度等特征都会发生变化。其中尺度特征变化的计算最为复杂。

为计算指定物体的尺度变化，现有的算法依赖于多尺度特征提取。算法过程是：获取当前时刻不同尺度下指定物体的图像，将每个尺度下指定物体的图像与前一时刻指定物体的目标图像进行匹配，得到当前时刻指定物体的目标图像，其中当前时刻指定物体的目标图像是当前时刻含有指定物体的区域的图像。

虽然上述多尺度特征提取匹配算法可以得到当前时刻指定物体的目标图像，但是其算法复杂度高，匹配数据量大，降低匹配的实时性。

发明内容

本申请所要解决的技术问题是提供一种图像匹配方法和系统，用以解决现有多尺度特征提取匹配算法的复杂度高，匹配数据量大导致的匹配实时性降低的问题。技术方案如下：

一方面，本申请提供一种图像匹配方法，包括：

将前一时刻指定物体的目标图像的中心点作为当前时刻图像的中心点，构建当前时刻一种尺度下搜索区域图像，其中所述搜索区域图像可能包含指定物体，所述搜索区域图像的横向分辨率大于所述目标图像的横向分辨率和/或所述搜索区域图像的纵向分辨率大于所述目标图像的纵向分辨率；

以所述目标图像为划分模板划分所述搜索区域图像，得到与所述目标图像的横向分辨率和纵向分辨率相同的至少两个候选图像；

计算每个所述候选图像和所述目标图像的相似度，选取相似度最大的候选图像；

计算所述目标图像和所述相似度最大的候选图像的弹性参数；

依据所述弹性参数调整所述相似度最大的候选图像；

将调整后的候选图像作为当前时刻指定物体的目标图像。

优选地，所述计算每个所述候选图像和所述目标图像的相似度包括：以m*m邻域对所述目标图像进行划分，其中m是领域的分辨率；

计算每个邻域中各个像素点的梯度，并选取梯度最大值的像素点作为所述像素点所在邻域的特征点；

选取所述候选图像中与每个所述特征点各自对应的相似点，其中所述相似点是以对应的所述特征点为圆心，以R为半径得到的圆形邻域中，梯度数值与所述特征点的梯度数值最接近的像素点；

计算每个所述特征点和对应的所述相似点的相似度；

对同一个候选图像中的相似度进行求和，得到该候选图像和目标图像的相似度。

优选地，所述计算每个领域中各个像素点的梯度包括：

依据公式：

$\mathrm{GRAD}\left(x\right)={\mathrm{\Σ}}_{x=1}^{m/2}{\mathrm{\Σ}}_{y=-m/2}^{m/2}P(x,y)-{\mathrm{\Σ}}_{x=-m/2}^{-1}{\mathrm{\Σ}}_{y=-m/2}^{m/2}P(x,y)/m\×(\frac{m}{2}-1)$

$\mathrm{GRAD}\left(y\right)={\mathrm{\Σ}}_{y=1}^{m/2}{\mathrm{\Σ}}_{x=-m/2}^{m/2}P(x,y)-{\mathrm{\Σ}}_{y=-m/2}^{-1}{\mathrm{\Σ}}_{x=-m/2}^{m/2}P(x,y)/m\×(\frac{m}{2}-1)$

得到每个领域中各个像素点的梯度GRAD=(GRAD(x),GRAD(y))，其中P(x,y)表示像素点，x为横向分辨率，y为横向分辨率。

优选地，计算每个所述特征点和对应的所述相似点的相似度包括：

计算所述特征点和对应的所述相似点的均值相似度D_EX(P₁,P₂)=W_EX×|EX₁-EX₂|，其中P₁为所述特征点，P₂为所述特征点对应的所述相似点，W_EX为均值权值，EX₁为P₁的均值：EX₂为P₂的均值： ${\mathrm{EX}}_2=\frac{{\mathrm{\Σ}}_{i=1}^n{P}_i}{n};$

计算所述特征点和对应的所述相似点的方差相似度D_DX(P₁,P₂)=W_DX×|DX₁-DX₂|，其中，W_DX为方差权值，DX₁为P₁的方差： ${\mathrm{DX}}_1=\frac{{\mathrm{\Σ}}_{i=1}^n{P}_i*P_i}{n}-{\left({\mathrm{EX}}_1\right)}^2,$ DX₂为P₂的方差： ${\mathrm{DX}}_2=\frac{{\mathrm{\Σ}}_{i=1}^n{P}_i*P_i}{n}-{\left({\mathrm{EX}}_2\right)}^2;$

计算所述特征点和对应的所述相似点的梯度相似度D_GRAD(P₁,P₂)=W_GRAD×|GRAD₁-GRAD₂|，其中，W_GRAD为梯度权值，GRAD₁为P₁的梯度，GRAD₂为P₂的梯度；

依据公式S(P₁,P₂)=D_EX(P₁,P₂)+D_DX(P₁,P₂)+D_GRAD(P₁,P₂)得到所述特征点和对应的所述相似点的相似度。

优选地，所述计算所述目标图像和所述相似度最大的候选图像的弹性参数包括：

依据公式得到所述弹性参数的位置偏移，其中T_i-1是前一时刻，T_i是当前时刻，是前一时刻指定物体的目标图像中所有特征点的中心： 是当前时刻指定物体的候选图像中所有相似点的中心：n是特征点和相似点的总个数，P_k是特征点，P′_k是相似点；

依据公式得到所述弹性参数的尺度变化，其中D_k为所述特征点到中心的距离：D″_k为所述相似点到中心的距离： $D_k^{\′}=\left|\right|P_k-\stackrel{\‾}{C^{\′}}\left|\right|.$

优选地，依据所述弹性参数调整所述相似度最大的候选图像，得到当前时刻指定物体的目标图像包括：

调整所述相似度最大的候选图像的中心点为其中CT_i-1(x_i,y_i)为前一时刻指定物体的目标图像的中心点；

调整所述相似度最大的候选图像的横向分辨率为M×Z(T_i-1,T_i)，纵向分辨率为N×Z(T_i-1,T_i)，其中M为前一时刻指定物体的目标图像的横向分辨率，N为前一时刻指定物体的目标图像的纵向分辨率；

调整后的所述相似度最大的候选图像作为当前时刻指定物体的目标图像。

另一方面，本申请还提供一种图像匹配系统，包括：

构建模块，用于将前一时刻指定物体的目标图像的中心点作为当前时刻图像的中心点，构建当前时刻一种尺度下搜索区域图像，其中所述搜索区域图像可能包含指定物体，所述搜索区域图像的横向分辨率大于所述目标图像的横向分辨率和/或所述搜索区域图像的纵向分辨率大于所述目标图像的纵向分辨率；

划分模块，用于以所述目标图像为划分模板划分所述图像，得到与所述目标图像的横向分辨率和纵向分辨率相同的至少两个候选图像；

第一计算模块，用于计算每个所述候选图像和所述目标图像的相似度，选取相似度最大的候选图像；

第二计算模块，用于计算所述目标图像和所述相似度最大的候选图像的弹性参数；

调整模块，用于依据所述弹性参数调整所述相似度最大的候选图像；

图像获取模块，用于将调整后的候选图像作为当前时刻指定物体的目标图像。

优选地，所述第一计算模块包括：

划分单元，用于以m*m邻域对所述目标图像进行划分，其中m是领域的分辨率；

第一计算单元，用于计算每个邻域中各个像素点的梯度，并选取梯度最大值的像素点作为所述像素点所在邻域的特征点；

选取单元，用于选取所述候选图像中与每个所述特征点各自对应的相似点，其中所述相似点是以对应的所述特征点为圆心，以R为半径得到的圆形邻域中，梯度数值与所述特征点的梯度数值最接近的像素点；

第二计算单元，用于计算每个所述特征点和对应的所述相似点的相似度；

求和单元，用于对同一个候选图像中的相似度进行求和，得到该候选图像和目标图像的相似度。

优选地，所述第一计算单元具体用于：

依据公式：

$\mathrm{GRAD}\left(x\right)={\mathrm{\Σ}}_{x=1}^{m/2}{\mathrm{\Σ}}_{y=-m/2}^{m/2}P(x,y)-{\mathrm{\Σ}}_{x=-m/2}^{-1}{\mathrm{\Σ}}_{y=-m/2}^{m/2}P(x,y)/m\×(\frac{m}{2}-1)$

$\mathrm{GRAD}\left(y\right)={\mathrm{\Σ}}_{y=1}^{m/2}{\mathrm{\Σ}}_{x=-m/2}^{m/2}P(x,y)-{\mathrm{\Σ}}_{y=-m/2}^{-1}{\mathrm{\Σ}}_{x=-m/2}^{m/2}P(x,y)/m\×(\frac{m}{2}-1)$

得到每个领域中各个像素点的梯度GRAD=(GRAD(x),GRAD(y))，其中P(x,y)表示像素点，x为横向分辨率，y为横向分辨率。

优选地，所述第二计算单元包括：

第一计算子单元，用于计算所述特征点和对应的所述相似点的均值相似度D_EX(P₁,P₂)=W_EX×|EX₁-EX₂|，其中P₁为所述特征点，P₂为所述特征点对应的所述相似点，W_EX为均值权值，EX₁为P₁的均值：EX₂为P₂的均值： ${\mathrm{EX}}_2=\frac{{\mathrm{\Σ}}_{i=1}^n{P}_i}{n};$

第二计算子单元，用于计算所述特征点和对应的所述相似点的方差相似度D_DX(P₁,P₂)=W_DX×|DX₁-DX₂|，其中，W_DX为方差权值，DX₁为P₁的方差： ${\mathrm{DX}}_1=\frac{{\mathrm{\Σ}}_{i=1}^n{P}_i*P_i}{n}-{\left({\mathrm{EX}}_1\right)}^2,$ DX₂为P₂的方差： ${\mathrm{DX}}_2=\frac{{\mathrm{\Σ}}_{i=1}^n{P}_i*P_i}{n}-{\left({\mathrm{EX}}_2\right)}^2;$

第三计算子单元，用于计算所述特征点和对应的所述相似点的梯度相似度D_GRAD(P₁,P₂)=W_GRAD×|GRAD₁-GRAD₂|，其中，W_GRAD为梯度权值，GRAD₁为P₁的梯度，GRAD₂为P₂的梯度；

第四计算子单元，用于依据公式S(P₁,P₂)=D_EX(P₁,P₂)+D_DX(P₁,P₂)+D_GRAD(P₁,P₂)得到所述特征点和对应的所述相似点的相似度。

优选地，所述第二计算模块包括：

位置偏移计算单元，用于依据公式得到所述弹性参数的位置偏移，其中T_i-1是前一时刻，T_i是当前时刻，是前一时刻指定物体的目标图像中所有特征点的中心： 是当前时刻指定物体的候选图像中所有相似点的中心：n是特征点和相似点的总个数，P_k是特征点，P′_k是相似点；

尺度变化计算单元，用于依据公式得到所述弹性参数的尺度变化，其中D_k为所述特征点到中心的距离：D″_k为所述相似点到中心的距离：

优选地，所述调整模块包括：

中心点调整单元，用于调整所述相似度最大的候选图像的中心点为其中CT_i-1(x_i,y_i)为前一时刻指定物体的目标图像的中心点；

分辨率调整单元，用于调整所述相似度最大的候选图像的横向分辨率为M×Z(T_i-1,T_i)，纵向分辨率为N×Z(T_i-1,T_i)，其中M为前一时刻指定物体的目标图像的横向分辨率，N为前一时刻指定物体的目标图像的纵向分辨率。

与现有技术相比，本申请包括以下优点：

在本申请中，首先得到当前时刻一种尺度下可能包含指定物体的搜索区域图像，该图像的中心点为前一时刻指定物体的目标图像的中心点，然后以所述目标图像为划分模板划分所述搜索区域图像，得到与所述目标图像的横向分辨率和纵向分辨率相同的至少两个候选图像，再从多个候选图像中选取出相似度最大的候选图像，依据目标图像和相似度最大的候选图像的弹性参数调整相似度最大的候选图像，将调整后的候选图像作为当前时刻指定物体的目标图像。因此在得到当前时刻指定物体的目标图像过程中，可以将前一时刻指定物体的目标图像和当前时刻一种尺度下指定物体的图像进行匹配，相对于现有多尺度特征提取匹配算法来说，降低匹配复杂度和匹配数据量，从而提高匹配实时性。

附图说明

为了更清楚地说明本申请实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。

图1是本申请提供的一种图像匹配方法的流程图；

图2是图像划分的一种示意图；

图3是图像划分的另一种示意图；

图4是本申请提供的一种图像匹配方法的子流程图；

图5是图像划分的再一种示意图；

图6是相似点选取的示意图；

图7是本申请提供的一种图像匹配系统的结构示意图；

图8是本申请提供的一种图像匹配系统的子结构示意图。

具体实施方式

首先对本申请提供的一种图像匹配方法进行说明，该方法可以包括以下步骤：

将前一时刻指定物体的目标图像的中心点作为当前时刻图像的中心点，得到当前时刻一种尺度下可能搜索区域图像，其中所述搜索区域图像包含指定物体，所述搜索区域图像的横向分辨率大于所述目标图像的横向分辨率，所述搜索区域图像的纵向分辨率大于所述目标图像的纵向分辨率；

以所述目标图像为划分模板划分所述搜索区域图像，得到与所述目标图像的横向分辨率和纵向分辨率相同的至少两个候选图像；

计算每个所述候选图像和所述目标图像的相似度，选取相似度最大的候选图像；

计算所述目标图像和所述相似度最大的候选图像的弹性参数；

依据所述弹性参数调整所述相似度最大的候选图像；

将调整后的候选图像作为当前时刻指定物体的目标图像。

该方法提出了匹配当前时刻单一尺度下指定物体的图像与前一时刻的目标图像，以得到当前时刻的目标图像的图像匹配方法，与现有多尺度特征提取匹配算法来说，匹配复杂度和匹配数据量降低，从而提高匹配的实时性。

进一步在匹配过程中，选取目标图像的每个邻域中梯度最大值的像素点作为特征点进行匹配，再次降低了匹配数据量，提高匹配的实时性。

下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本申请保护的范围。

请参阅图1，其示出了本申请实施例提供的一种图像匹配方法的一种流程图，可以包括以下步骤：

步骤101：将前一时刻指定物体的目标图像的中心点作为当前时刻图像的中心点，得到当前时刻一种尺度下搜素区域图像。

在本实施例中，指定物体是视频分析过程中所跟踪的物体，跟踪视频中该指定物体需要匹配到当前时刻该指定物体所在区域，并获取到该区域的图像作为指定物体的目标图像。

在获取当前时刻指定物体的目标图像时，需要选取一个匹配基准。本申请的一种可行方式是：以前一时刻指定物体的目标图像的中心点作为当前时刻图像的中心点，构建当前时刻一种尺度下搜索区域图像。其中该搜索区域图像可能包含指定物体。

在本申请的另一种可行方式中：还可以以其他时刻指定物体的目标图像的中心点作为当前时刻图像的中心点，构建当前时刻一种尺度下可能包含指定物体的搜索区域图像。由于当前时刻指定物体的轨迹和状态更加接近于前一时刻目标图像中指定物体的轨迹和状态，所以本实施例优选以前一时刻指定物体的目标图像的中心点作为当前时刻图像的中心点，构建当前时刻一种尺度下可能包含指定物体的搜索区域图像。

需要注意的一点是：实际视频中，指定物体会不断移动，所以为了跟踪该指定物体，构建的可能包含指定物体的搜索区域图像要比前一时刻的目标图像大。如搜索区域图像的横向分辨率大于目标图像的横向分辨率，或者搜索区域图像的纵向分辨率大于目标图像的纵向分辨率。当然搜索区域图像的横向分辨率和纵向分辨率可以分别大于目标图像的横向分辨率和纵向分辨率。

步骤102：以目标图像为划分模板划分图像，得到与目标图像的横向分辨率和纵向分辨率相同的至少两个候选图像。

为了便于以目标图像为划分模板划分图像，当前时刻搜索区域图像的横向分辨率优选是目标图像的横向分辨率的整数倍，搜索区域图像的纵向分辨率优选是目标图像的纵向分辨率的整数倍，两者的关系是Mc*Nc=（F*M）*（F*N）。其中Mc是搜索区域图像的横向分辨率，Nc是搜索区域图像的纵向分辨率，M是目标图像的横向分辨率，N是目标图像的纵向分辨率。

在目标跟踪中，相邻帧间指定物体移动的范围有限，所以在本申请中F的取值优选2，如目标图像的尺寸是M*N，则指定物体的图像的尺寸是2M*2N。使用目标图像对该种尺寸的图像进行划分时，可以将该种尺寸的图像进行均分，得到4个与目标图像的横向分辨率和纵向分辨率相同的候选图像。

在划分图像时，目标图像需要全部位于搜索区域图像内，而不能使目标图像部分位于搜索区域图像外，如图2所示示意图，其中填充有白色的图像是目标图像，未填充颜色的图像是搜索区域图像。当目标图像对搜索区域图像划分后，如果剩余图像区域不能全部包括目标图像，则将目标图像向指定物体的图像内部移动直至目标图像全部包括在搜索区域图像中，如图3所示。

步骤103：计算每个候选图像和目标图像的相似度，选取相似度最大的候选图像。

其中相似度是指候选图像和目标图像的相似程度。相似度越大，说明候选图像和目标图像含有更多相同的内容，候选图像中包括指定物体的可能性也越大。

在本实施例中，计算候选图像和目标图像的相似度的过程请参阅图4所示，包括：

步骤1031：以m*m邻域对目标图像进行划分，其中m是领域的分辨率。

在对目标图像进行邻域划分时，需要将目标图像划分成多个m*m的邻域，如图5所示，12*9的目标图像被均分成12个3*3的邻域。

需要说明的是：目标图像若不能被均分成多个m*m的邻域，则在目标图像划分过程中，选择使未能划分至邻域的边缘部分越少的邻域划分。如7*7的目标图像，可以被均分成4个3*3的邻域，相对于被分成1个4*4邻域来说，被划分为3*3的邻域的目标图像中未划分至邻域的边缘部分小于采用4*4方式划分。

进一步舍弃目标图像的边缘部分。由于舍弃的边缘部分在目标图像内所占比例较小，不会影响匹配过程。

步骤1032：计算每个邻域中各个像素点的梯度，并选取梯度最大值的像素点作为像素点所在邻域的特征点。

其中像素点采用P(x,y)表示，梯度的表示方式是GRAD=(GRAD(x),GRAD(y))，x为横向分辨率，y为横向分辨率，GRAD(x)表示x方向的梯度，GRAD(y)表示y方向的梯度。计算公式如下：

$\mathrm{GRAD}\left(x\right)={\mathrm{\Σ}}_{x=1}^{m/2}{\mathrm{\Σ}}_{y=-m/2}^{m/2}P(x,y)-{\mathrm{\Σ}}_{x=-m/2}^{-1}{\mathrm{\Σ}}_{y=-m/2}^{m/2}P(x,y)/m\×(\frac{m}{2}-1)$

$\mathrm{GRAD}\left(y\right)={\mathrm{\Σ}}_{y=1}^{m/2}{\mathrm{\Σ}}_{x=-m/2}^{m/2}P(x,y)-{\mathrm{\Σ}}_{y=-m/2}^{-1}{\mathrm{\Σ}}_{x=-m/2}^{m/2}P(x,y)/m\×(\frac{m}{2}-1).$

步骤1033：选取候选图像中与每个特征点各自对应的相似点，其中相似点是以对应的特征点为圆心，以R为半径得到的圆形邻域中，梯度数值与特征点的梯度数值最接近的像素点。

下面结合附图6对选取过程进行说明，左边为前一时刻的目标图像，右边为当前时刻的候选图像，对于前一时刻目标图像的特征点P_i，在对应当前时刻的候选图像中，以P_i为中心、R为半径得到圆形邻域。其中R的取值根据不同应用场景可以取不同数值，对此本实施例不加以限制。

计算该圆形邻域中各个像素点的梯度，然后匹配出梯度数值与特征点P_i的梯度数值最接近的像素点，如图6中的像素点p′_i为特征点P_i的相似点。圆形邻域中各个像素点的梯度计算公式可以参阅特征点的梯度计算公式。

步骤1034：计算每个特征点和对应的相似点的相似度。具体过程如下：

计算特征点和对应的相似点的均值相似度D_EX(P₁,P₂)=W_EX×|EX₁-EX₂|，其中P₁为特征点，P₂为特征点对应的相似点，W_EX为均值权值，EX₁为P₁的均值：EX₂为P₂的均值： ${\mathrm{EX}}_2=\frac{{\mathrm{\Σ}}_{i=1}^n{P}_i}{n};$

计算特征点和对应的相似点的方差相似度D_DX(P₁,P₂)=W_DX×|DX₁-DX₂|，其中，W_DX为方差权值，DX₁为P₁的方差： ${\mathrm{DX}}_1=\frac{{\mathrm{\Σ}}_{i=1}^n{P}_i*P_i}{n}-{\left({\mathrm{EX}}_1\right)}^2,$ DX₂为P₂的方差： ${\mathrm{DX}}_2=\frac{{\mathrm{\Σ}}_{i=1}^n{P}_i*P_i}{n}-{\left({\mathrm{EX}}_2\right)}^2.$

计算特征点和对应的相似点的梯度相似度D_GRAD(P₁,P₂)=W_GRAD×|GRAD₁-GRAD₂|，其中，W_GRAD为梯度权值，GRAD₁为P₁的梯度，GRAD₂为P₂的梯度。

依据公式S(P₁,P₂)=D_EX(P₁,P₂)+D_DX(P₁,P₂)+D_GRAD(P₁,P₂)得到特征点和对应的相似点的相似度。

步骤1035：对同一个候选图像中的相似度进行求和，得到该候选图像和目标图像的相似度。

步骤104：计算目标图像和相似度最大的候选图像的弹性参数。

其中弹性参数包括位置偏移和尺度变化。位置偏移的计算公式是：其中T_i-1是前一时刻，T_i是当前时刻，是前一时刻指定物体的目标图像中所有特征点的中心： 是当前时刻指定物体的候选图像中所有相似点的中心：n是特征点和相似点的总个数，P_k是特征点，P′_k是相似点。

尺度变化的计算公式是：得到所述弹性参数的尺度变化，其中D_k为所述特征点到中心的距离：D″_k为所述相似点到中心的距离：

步骤105：依据弹性参数调整相似度最大的候选图像。

在本实施例中，调整候选图像主要包括调整候选图像的中心点、横向分辨率和纵向分辨率，具体过程可以是：

调整所述相似度最大的候选图像的中心点为其中CT_i-1(x_i,y_i)为前一时刻指定物体的目标图像的中心点；调整所述相似度最大的候选图像的横向分辨率为M×Z(T_i-1,T_i)，纵向分辨率为N×Z(T_i-1,T_i)，其中M为前一时刻指定物体的目标图像的横向分辨率，N为前一时刻指定物体的目标图像的纵向分辨率。

步骤106：将调整后的候选图像作为当前时刻指定物体的目标图像。

应用本申请实施例提供的图像匹配方法，首先得到当前时刻一种尺度下搜索区域图像，该图像的中心点为前一时刻指定物体的目标图像的中心点，然后以所述目标图像为划分模板划分所述图像，得到与所述目标图像的横向分辨率和纵向分辨率相同的至少两个候选图像，再从多个候选图像中选取出相似度最大的候选图像，依据目标图像和相似度最大的候选图像的弹性参数调整相似度最大的候选图像，调整后的候选图像作为当前时刻指定物体的目标图像。因此在得到当前时刻指定物体的目标图像过程中，可以将前一时刻指定物体的目标图像和当前时刻一种尺度下指定物体的图像进行匹配，相对于现有多尺度特征提取匹配算法来说，降低匹配复杂度和匹配数据量，从而提高匹配实时性。

进一步，在具体匹配过程中，通过单一特征点和相似点的相似度计算得到最有可能包括指定物体的候选图像，再次降低匹配数据量和匹配复杂度，提高匹配实时性。

与上述方法实施例相对应，本申请实施例还提供一种图像匹配系统，其结构示意图请参阅7所示，可以包括：构建模块11、划分模块12、第一计算模块13、第二计算模块14、调整模块15和图像获取模块16。其中，

构建模块11，用于将前一时刻指定物体的目标图像的中心点作为当前时刻图像的中心点，构建当前时刻一种尺度下搜索区域图像，其中所述搜索区域图像可能包含指定物体，所述搜索区域图像的横向分辨率大于所述目标图像的横向分辨率和/或所述图像的纵向分辨率大于所述目标图像的纵向分辨率。

划分模块12，用于以所述目标图像为划分模板划分所述图像，得到与所述目标图像的横向分辨率和纵向分辨率相同的至少两个候选图像。

第一计算模块13，用于计算每个所述候选图像和所述目标图像的相似度，选取相似度最大的候选图像。其中相似度是指候选图像和目标图像的相似程度。相似度越大，说明候选图像和目标图像含有更多相同的内容，候选图像中包括指定物体的可能性也越大。

在本实施例中，第一计算模块13的结构示意图请参阅图8所示，可以包括：划分单元131、第一计算单元132、选取单元133、第二计算单元134和求和单元135。其中，

划分单元131，用于以m*m邻域对所述目标图像进行划分，其中m是领域的分辨率。

第一计算单元132，用于计算每个邻域中各个像素点的梯度，并选取梯度最大值的像素点作为所述像素点所在邻域的特征点。具体第一计算单元132具体用于：

依据公式：

$\mathrm{GRAD}\left(x\right)={\mathrm{\Σ}}_{x=1}^{m/2}{\mathrm{\Σ}}_{y=-m/2}^{m/2}P(x,y)-{\mathrm{\Σ}}_{x=-m/2}^{-1}{\mathrm{\Σ}}_{y=-m/2}^{m/2}P(x,y)/m\×(\frac{m}{2}-1)$

$\mathrm{GRAD}\left(y\right)={\mathrm{\Σ}}_{y=1}^{m/2}{\mathrm{\Σ}}_{x=-m/2}^{m/2}P(x,y)-{\mathrm{\Σ}}_{y=-m/2}^{-1}{\mathrm{\Σ}}_{x=-m/2}^{m/2}P(x,y)/m\×(\frac{m}{2}-1)$

得到每个领域中各个像素点的梯度GRAD=(GRAD(x),GRAD(y))，其中P(x,y)表示像素点，x为横向分辨率，y为横向分辨率。

选取单元133，用于选取所述候选图像中与每个所述特征点各自对应的相似点，其中所述相似点是以对应的所述特征点为圆心，以R为半径得到的圆形邻域中，梯度数值与所述特征点的梯度数值最接近的像素点。

第二计算单元134，用于计算每个所述特征点和对应的所述相似点的相似度。求和单元135，用于对同一个候选图像中的相似度进行求和，得到该候选图像和目标图像的相似度。

在本实施例中，第二计算单元134可以包括：第一计算子单元、第二计算子单元、第三计算子单元和第四计算子单元。其中，

第一计算子单元，用于计算所述特征点和对应的所述相似点的均值相似度D_EX(P₁,P₂)=W_EX×|EX₁-EX₂|，其中P₁为所述特征点，P₂为所述特征点对应的所述相似点，W_EX为均值权值，EX₁为P₁的均值：EX₂为P₂的均值： ${\mathrm{EX}}_2=\frac{{\mathrm{\Σ}}_{i=1}^n{P}_i}{n}.$

第二计算子单元，用于计算所述特征点和对应的所述相似点的方差相似度D_DX(P₁,P₂)=W_DX×|DX₁-DX₂|，其中，W_DX为方差权值，DX₁为P₁的方差： ${\mathrm{DX}}_1=\frac{{\mathrm{\Σ}}_{i=1}^n{P}_i*P_i}{n}-{\left({\mathrm{EX}}_1\right)}^2,$ DX₂为P₂的方差： ${\mathrm{DX}}_2=\frac{{\mathrm{\Σ}}_{i=1}^n{P}_i*P_i}{n}-{\left({\mathrm{EX}}_2\right)}^2.$

第三计算子单元，用于计算所述特征点和对应的所述相似点的梯度相似度D_GRAD(P₁,P₂)=W_GRAD×|GRAD₁-GRAD₂|，其中，W_GRAD为梯度权值，GRAD₁为P₁的梯度，GRAD₂为P₂的梯度。

第四计算子单元，用于依据公式S(P₁,P₂)=D_EX(P₁,P₂)+D_DX(P₁,P₂)+D_GRAD(P₁,P₂)得到所述特征点和对应的所述相似点的相似度。

第二计算模块14，用于计算所述目标图像和所述相似度最大的候选图像的弹性参数。其中第二计算模块14可以包括：位置偏移计算单元和尺度变化计算单元。其中，

位置偏移计算单元，用于依据公式得到所述弹性参数的位置偏移，其中T_i-1是前一时刻，T_i是当前时刻，是前一时刻指定物体的目标图像中所有特征点的中心： 是当前时刻指定物体的候选图像中所有相似点的中心：n是特征点和相似点的总个数，P_k是特征点，P′_k是相似点。

尺度变化计算单元，用于依据公式得到所述弹性参数的尺度变化，其中D_k为所述特征点到中心的距离：D″_k为所述相似点到中心的距离：

调整模块15，用于依据所述弹性参数调整所述相似度最大的候选图像。在本实施例中，所述调整模块15包括：中心点调整单元和分辨率调整单元。

中心点调整单元，用于调整所述相似度最大的候选图像的中心点为其中CT_i-1(x_i,y_i)为前一时刻指定物体的目标图像的中心点。

分辨率调整单元，用于调整所述相似度最大的候选图像的横向分辨率为M×Z(T_i-1,T_i)，纵向分辨率为N×Z(T_i-1,T_i)，其中M为前一时刻指定物体的目标图像的横向分辨率，N为前一时刻指定物体的目标图像的纵向分辨率。

图像获取模块16，用于将调整后的候选图像作为当前时刻指定物体的目标图像。

需要说明的是，本说明书中的各个实施例均采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似的部分互相参见即可。对于系统类实施例而言，由于其与方法实施例基本相似，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。

最后，还需要说明的是，在本文中，诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。

为了描述的方便，描述以上装置时以功能分为各种模块/单元分别描述。当然，在实施本申请时可以把各单元的功能在同一个或多个软件和/或硬件中实现。

以上对本申请所提供的一种图像匹配方法和系统进行了详细介绍，本文中应用了具体个例对本申请的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本申请的方法及其核心思想；同时，对于本领域的一般技术人员，依据本申请的思想，在具体实施方式及应用范围上均会有改变之处，综上所述，本说明书内容不应理解为对本申请的限制。

Claims (8)

1.一种图像匹配方法，其特征在于，包括：

将前一时刻指定物体的目标图像的中心点作为当前时刻图像的中心点，构建当前时刻一种尺度下搜索区域图像，其中所述搜索区域图像可能包含指定物体，所述搜索区域图像的横向分辨率大于所述目标图像的横向分辨率和/或所述搜索区域图像的纵向分辨率大于所述目标图像的纵向分辨率；

以所述目标图像为划分模板划分所述搜索区域图像，得到与所述目标图像的横向分辨率和纵向分辨率相同的至少两个候选图像；

计算每个所述候选图像和所述目标图像的相似度，选取相似度最大的候选图像；

计算所述目标图像和所述相似度最大的候选图像的弹性参数；

依据所述弹性参数调整所述相似度最大的候选图像；

将调整后的候选图像作为当前时刻指定物体的目标图像；

其中，所述计算每个所述候选图像和所述目标图像的相似度包括：以m*m邻域对所述目标图像进行划分，其中m是领域的分辨率；

计算每个邻域中各个像素点的梯度，并选取梯度最大值的像素点作为所述像素点所在邻域的特征点；

选取所述候选图像中与每个所述特征点各自对应的相似点，其中所述相似点是以对应的所述特征点为圆心，以R为半径得到的圆形邻域中，梯度数值与所述特征点的梯度数值最接近的像素点；

计算每个所述特征点和对应的所述相似点的相似度；

对同一个候选图像中的相似度进行求和，得到该候选图像和目标图像的相似度；

其中，计算每个所述特征点和对应的所述相似点的相似度包括：

计算所述特征点和对应的所述相似点的均值相似度D_EX(P₁,P₂)＝W_EX×|EX₁-EX₂|，其中P₁为所述特征点，P₂为所述特征点对应的所述相似点，W_EX为均值权值，EX₁为P₁的均值：EX₂为P₂的均值： ${\mathrm{EX}}_2=\frac{{\mathrm{\Σ}}_{i=1}^n{P}_i}{n};$

计算所述特征点和对应的所述相似点的方差相似度D_DX(P₁,P₂)＝W_DX×|DX₁-DX₂|，其中，W_DX为方差权值，DX₁为P₁的方差： ${\mathrm{DX}}_1=\frac{{\mathrm{\Σ}}_{i=1}^n{P}_i*P_i}{n}-{\left({\mathrm{EX}}_1\right)}^2,$ DX₂为P₂的方差： ${\mathrm{DX}}_2=\frac{{\mathrm{\Σ}}_{i=1}^n{P}_i*P_i}{n}-{\left({\mathrm{EX}}_2\right)}^2;$

计算所述特征点和对应的所述相似点的梯度相似度：

D_GRAD(P₁,P₂)＝W_GRAD×|GRAD₁-GRAD₂|

其中，W_GRAD为梯度权值，GRAD₁为P₁的梯度，GRAD₂为P₂的梯度；

依据公式S(P₁,P₂)＝D_EX(P₁,P₂)+D_DX(P₁,P₂)+D_GRAD(P₁,P₂)得到所述特征点和对应的所述相似点的相似度。

2.根据权利要求1所述的图像匹配方法，其特征在于，所述计算每个领域中各个像素点的梯度包括：

依据公式：

$GRAD\left(x\right)={\mathrm{\Σ}}_{x=1}^{m/2}{\mathrm{\Σ}}_{y=-m/2}^{m/2}P(x,y)-{\mathrm{\Σ}}_{x=-m/2}^{-1}{\mathrm{\Σ}}_{y=-m/2}^{m/2}P(x,y)/m\×(\frac{m}{2}-1)$

$GRAD\left(y\right)={\mathrm{\Σ}}_{y=1}^{m/2}{\mathrm{\Σ}}_{x=-m/2}^{m/2}P(x,y)-{\mathrm{\Σ}}_{y=-m/2}^{-1}{\mathrm{\Σ}}_{x=-m/2}^{m/2}P(x,y)/m\×(\frac{m}{2}-1)$

得到每个领域中各个像素点的梯度GRAD＝(GRAD(x),GRAD(y))，其中P(x,y)表示像素点，x为横向分辨率，y为横向分辨率。

3.根据权利要求1所述的图像匹配方法，其特征在于，所述计算所述目标图像和所述相似度最大的候选图像的弹性参数包括：

依据公式得到所述弹性参数的位置偏移，其中T_i-1是前一时刻，T_i是当前时刻，是前一时刻指定物体的目标图像中所有特征点的中心： 是当前时刻指定物体的候选图像中所有相似点的中心：n是特征点和相似点的总个数，P_k是特征点，P′_k是相似点；

依据公式得到所述弹性参数的尺度变化，其中D_k为所述特征点到中心的距离：D′_k为所述相似点到中心的距离： $D_k^{\′}=\left|\right|P_k-\stackrel{\‾}{C^{\′}}\left|\right|.$

4.根据权利要求3所述的图像匹配方法，其特征在于，依据所述弹性参数调整所述相似度最大的候选图像，得到当前时刻指定物体的目标图像包括：

调整所述相似度最大的候选图像的中心点为其中CT_i-1(x_i,y_i)为前一时刻指定物体的目标图像的中心点；

调整所述相似度最大的候选图像的横向分辨率为M×Z(T_i-1,T_i)，纵向分辨率为N×Z(T_i-1,T_i)，其中M为前一时刻指定物体的目标图像的横向分辨率，N为前一时刻指定物体的目标图像的纵向分辨率；

调整后的所述相似度最大的候选图像作为当前时刻指定物体的目标图像。

5.一种图像匹配系统，其特征在于，包括：

构建模块，用于将前一时刻指定物体的目标图像的中心点作为当前时刻图像的中心点，构建当前时刻一种尺度下搜索区域图像，其中所述搜索区域图像可能包含指定物体，所述搜索区域图像的横向分辨率大于所述目标图像的横向分辨率和/或所述搜索区域图像的纵向分辨率大于所述目标图像的纵向分辨率；

划分模块，用于以所述目标图像为划分模板划分所述图像，得到与所述目标图像的横向分辨率和纵向分辨率相同的至少两个候选图像；

第一计算模块，用于计算每个所述候选图像和所述目标图像的相似度，选取相似度最大的候选图像；

第二计算模块，用于计算所述目标图像和所述相似度最大的候选图像的弹性参数；

调整模块，用于依据所述弹性参数调整所述相似度最大的候选图像；

图像获取模块，用于将调整后的候选图像作为当前时刻指定物体的目标图像；

其中，所述第一计算模块包括：

划分单元，用于以m*m邻域对所述目标图像进行划分，其中m是领域的分辨率；

第一计算单元，用于计算每个邻域中各个像素点的梯度，并选取梯度最大值的像素点作为所述像素点所在邻域的特征点；

选取单元，用于选取所述候选图像中与每个所述特征点各自对应的相似点，其中所述相似点是以对应的所述特征点为圆心，以R为半径得到的圆形邻域中，梯度数值与所述特征点的梯度数值最接近的像素点；

第二计算单元，用于计算每个所述特征点和对应的所述相似点的相似度；

求和单元，用于对同一个候选图像中的相似度进行求和，得到该候选图像和目标图像的相似度；

其中，所述第二计算单元包括：

第一计算子单元，用于计算所述特征点和对应的所述相似点的均值相似度D_EX(P₁,P₂)＝W_EX×|EX₁-EX₂|，其中P₁为所述特征点，P₂为所述特征点对应的所述相似点，W_EX为均值权值，EX₁为P₁的均值：EX₂为P₂的均值： ${\mathrm{EX}}_2=\frac{{\mathrm{\Σ}}_{i=1}^n{P}_i}{n};$

第二计算子单元，用于计算所述特征点和对应的所述相似点的方差相似度D_DX(P₁,P₂)＝W_DX×|DX₁-DX₂|，其中，W_DX为方差权值，DX₁为P₁的方差： ${\mathrm{DX}}_1=\frac{{\mathrm{\Σ}}_{i=1}^n{P}_i*P_i}{n}-{\left({\mathrm{EX}}_1\right)}^2,$ DX₂为P₂的方差： ${\mathrm{DX}}_2=\frac{{\mathrm{\Σ}}_{i=1}^n{P}_i*P_i}{n}-{\left({\mathrm{EX}}_2\right)}^2;$

第三计算子单元，用于计算所述特征点和对应的所述相似点的梯度相似度D_GRAD(P₁,P₂)＝W_GRAD×|GRAD₁-GRAD₂|，其中，W_GRAD为梯度权值，GRAD₁为P₁的梯度，GRAD₂为P₂的梯度；

第四计算子单元，用于依据如下公式得到所述特征点和对应的所述相似点的相似度：

S(P₁,P₂)＝D_EX(P₁,P₂)+D_DX(P₁,P₂)+D_GRAD(P₁,P₂)。

6.根据权利要求5所述的图像匹配系统，其特征在于，所述第一计算单元具体用于：

依据公式：

$GRAD\left(x\right)={\mathrm{\Σ}}_{x=1}^{m/2}{\mathrm{\Σ}}_{y=-m/2}^{m/2}P(x,y)-{\mathrm{\Σ}}_{x=-m/2}^{-1}{\mathrm{\Σ}}_{y=-m/2}^{m/2}P(x,y)/m\×(\frac{m}{2}-1)$

$GRAD\left(y\right)={\mathrm{\Σ}}_{y=1}^{m/2}{\mathrm{\Σ}}_{x=-m/2}^{m/2}P(x,y)-{\mathrm{\Σ}}_{y=-m/2}^{-1}{\mathrm{\Σ}}_{x=-m/2}^{m/2}P(x,y)/m\×(\frac{m}{2}-1)$

得到每个领域中各个像素点的梯度GRAD＝(GRAD(x),GRAD(y))，其中P(x,y)表示像素点，x为横向分辨率，y为横向分辨率。

7.根据权利要求5所述的图像匹配系统，其特征在于，所述第二计算模块包括：

位置偏移计算单元，用于依据公式得到所述弹性参数的位置偏移，其中T_i-1是前一时刻，T_i是当前时刻，是前一时刻指定物体的目标图像中所有特征点的中心： 是当前时刻指定物体的候选图像中所有相似点的中心：n是特征点和相似点的总个数，P_k是特征点，P′_k是相似点；

尺度变化计算单元，用于依据公式得到所述弹性参数的尺度变化，其中D_k为所述特征点到中心的距离：D′_k为所述相似点到中心的距离：

8.根据权利要求7所述的图像匹配系统，其特征在于，所述调整模块包括：

中心点调整单元，用于调整所述相似度最大的候选图像的中心点为其中CT_i-1(x_i,y_i)为前一时刻指定物体的目标图像的中心点；

分辨率调整单元，用于调整所述相似度最大的候选图像的横向分辨率为M×Z(T_i-1,T_i)，纵向分辨率为N×Z(T_i-1,T_i)，其中M为前一时刻指定物体的目标图像的横向分辨率，N为前一时刻指定物体的目标图像的纵向分辨率。