一种动态背景目标跟踪方法与流程

1.本发明涉及图像目标跟踪技术领域，特别是一种动态背景目标跟踪方法。


背景技术：

2.运动目标的跟踪作为军事制导、安全监控等系统中的核心技术，是计算机视觉领域的一个研究热点。国内对视频运动目标检测和跟踪技术的研究相对较晚。近年来，国内许多研究机构在该领域也开展了大量研究工作，但是目前对视频目标跟踪理论和应用的研究还处于在一定的限制条件下设计特定方法的阶段，许多目标跟踪算法对复杂的背景或者变化的目标的跟踪常常会遇到很大的困难，甚至丢失目标导致跟踪失败，通常只适用于一些背景不太复杂、目标变换不大或者其他一些特定的工作场合。因此，如何有效、精准地表征目标，如何在运动目标发生变化和背景复杂地情况下实时精确地跟踪目标，是目前视频目标跟踪研究中的主要难题。
3.目前比较经典的跟踪算法有meanshift、camshift算法、粒子滤波、卡尔曼滤波算法等。camshift比较简单，计算量较少，但是容易陷入局部最优，在背景比较复杂的视频中容易受到干扰，跟踪效果较差；粒子滤波运算复杂度较高，但是对于遮挡后再跟踪有比较好的效果。在复杂动态背景下，单一的传统跟踪算法并不能实现稳定的跟踪效果，因此多种跟踪算法融合的方式被提出并应用于实现稳定跟踪。
4.由于单一的传统目标跟踪算法在背景环境复杂时比较容易受到外界干扰，并且不可避免的存在目标被遮挡后难以继续跟踪的情况。此外，一些较为先进的基于深度学习的目标跟踪算法，需要事先训练大量样本集，而机场围界经过的运动目标类型较多，不利于统计样本，这些算法对平台性能要求较高，增加了应用成本。


技术实现要素：

5.本发明的目的在于提供一种系统稳健性强、跟踪效果好、抗干扰能力强、计算量小、实时性强的基于粒子滤波与camshift相结合的动态背景目标跟踪方法。
6.实现本发明目的的技术解决方案为：一种动态背景目标跟踪方法，包括以下步骤：
7.步骤1、手动选择跟踪目标作为初始的目标状态模板，并在跟踪目标所在区域放置粒子集；
8.步骤2、通过状态转移模型得到k时刻的粒子集；
9.步骤3、计算初始化后的目标状态模板的h
‑
s联合颜色直方图和梯度方向直方图，并利用camshift算法的迭代对k时刻粒子集中的每个粒子进行漂移，使粒子集中的粒子移动到目标当前所在的区域，得到漂移后的k时刻粒子集；
10.步骤4、根据目标的h
‑
s联合颜色直方图和梯度方向直方图建立观测模型，通过观测模型判断漂移后的k时刻粒子集中每个粒子状态与实际状态的相似度，然后计算粒子权值并归一化；
11.步骤5、对权值小于阈值e的粒子，利用随机采样法复制从粒子集中随机选取的粒
子取代原粒子，完成粒子重采样，其中e＝1/n；
12.步骤6、对重采样后的粒子加权求和，得到当前目标的估计状态，确定目标位置。
13.本发明与现有技术相比，其显著优点在于：(1)解决了遮挡物与目标颜色相似时对跟踪效果的干扰问题，系统有较强的稳健性；(2)在网格围栏与遮挡物较多的复杂背景中，跟踪过程基本不受网格干扰，跟踪效果好；(3)能够实时反映目标框位于遮挡物体处时目标所处的位置；(4)利用camshift迭代每个粒子使其漂移到目标区域，计算量小，实时性强。
附图说明
14.图1为本发明动态背景目标跟踪方法的流程示意图。
具体实施方式
15.下面结合附图和具体实施例对本发明做进一步的详细说明。
16.结合图1，本发明一种动态背景目标跟踪方法，包括如下步骤：
17.步骤1、手动选择跟踪目标作为初始的目标状态模板，并在跟踪目标所在区域放置粒子集，具体如下：
18.在视频中手动选择需要跟踪的目标，作为初始的目标状态模板，在此目标附近随机放置粒子集取共有n个粒子样本。该粒子集为在目标附近通过加高斯噪声方式产生的n个粒子其中为各个粒子的初始权值，表示初始状态下编号为i的粒子。
19.步骤2、通过状态转移模型得到k时刻的粒子集，具体如下：
20.通过状态转移模型得到的k时刻的粒子集记为其中代表粒子的随机传播半径，表示k时刻第i个粒子的可能状态。
21.状态转移模型为：
22.x
k
＝2x
k
‑1‑
x
k
‑2+ru
k
‑123.其中，x
k
为k时刻的目标状态，u
k
‑1为高斯噪声且u
k
‑1∈[
‑
1,1]，代表粒子的随机传播半径。
[0024]
步骤3、计算初始化后的目标状态模板的h
‑
s联合颜色直方图和梯度方向直方图，并利用camshift算法的迭代对k时刻粒子集中的每个粒子进行漂移，使粒子集中的粒子移动到目标当前所在的区域，得到漂移后的k时刻粒子集，具体如下：
[0025]
步骤3.1、计算初始目标模型的h
‑
s联合颜色直方图，具体如下：
[0026]
将rgb图像空间转换成hsv图像空间，再对其中的h
‑
s色调分量计算直方图，作为目标的颜色特征参数。由于在框选目标时，选择框的边界部分有可能是背景图像，因此利用核函数k(r)对初始目标模型的颜色分布直方图进行加权处理，使得随着像素远离目标模型中心，该像素所占权重逐渐减小。当像素在目标框边缘部分视，像素权重为0。增强了颜色特征模型的可靠性；
[0027]
核函数k(r)为：
[0028][0029]
式中r为像素点到目标框中心的距离；
[0030]
设定目标框中的像素点为x
i
，中心点为y，则目标框内的颜色分布直方图为：
[0031][0032][0033][0034]
其中m为目标的颜色分布直方图的维数，n为目标区域内像素总数；h(x
i
)为像素点x
i
点对应颜色分布直方图中的值，δ[h(x
i
)
‑
u]为单位冲击函数，根据x
i
点的颜色值是否位于第u个颜色分布直方图的区间内进行取值，若在则为1，否则为0；a为目标区域自适应变量，为目标框的宽和高；f为归一化因子；
[0035]
分别构造跟踪目标的h分量直方图和s分量直方图，然后将上述两种直方图加权融合得到跟踪目标的h
‑
s联合颜色直方图：
[0036][0037]
其中α∈(0.5,1)为加权系数；分别为h分量直方图、s分量直方图和h
‑
s联合颜色直方图，α∈(0.5,1)为加权系数；
[0038]
步骤3.2、利用sobel算子计算框选目标区域的梯度方向直方图，具体如下：
[0039]
sobel算子的横向和纵向梯度模板分别为g
x
、g
y
：
[0040][0041]
利用sobel算子的梯度模板计算图中目标区域所有像素点在x、y方向的梯度值，得到该框选目标在边缘区域的梯度幅值g和边缘方向θ为：
[0042][0043]
式中
[0044]
然后对各个方向上的梯度幅值进行统计，构造目标区域的梯度方向直方图，公式如下：
[0045][0046]
其中u为目标的颜色分布直方图的维数，n为目标区域内像素总数；θ(x
i
)为像素点x
i
点对应颜色分布直方图中的值，δ[θ(x
i
)
‑
u]为单位冲击函数；
[0047]
步骤3.3、利用camshift算法的迭代对k时刻粒子集中的每个粒子进行漂移，使粒子移动到目标当前所在的区域，得到漂移后的k时刻粒子集，具体如下：
[0048]
步骤3.3.1、计算方向投影，根据步骤3.1得到的目标区域的h
‑
s联合颜色直方图，计算输入图像的颜色概率密度分布图，即用输入图像的一位置上像素值对应在直方图上的值来代替该像素值，并将该数值归一化到0
‑
255范围；通过对输入图像中的每个像素均采用以上步骤，得到当前图像的颜色概率密度分布图，在颜色概率密度分布图中，图像的像素越大，代表该像素值出现的概率越大；
[0049]
步骤3.3.2、在当前帧图像的联合概率密度分布图中设置粒子i的初始搜索窗口，如果是第一帧图像则通过手动选择的方式确定跟踪目标，如果是其他帧图像则利用上一帧图像最终得到的粒子i窗口的质心位置和大小来设置粒子i在当前帧图像中搜索窗口的位置和尺寸；
[0050]
步骤3.3.3、计算粒子i搜索区域的零阶矩m
00
和一阶矩m
10
、m
01
，公式为：
[0051][0052][0053][0054]
其中(x,y)为粒子i的在搜索区域中的像素位置，i(x,y)为粒子i所在像素位置在图像的颜色概率密度分布图中的值；
[0055]
步骤3.3.4、计算粒子i的搜索区域的质心位置，并将粒子i新的搜索窗口的中心移动到质心(x
c
,y
c
)处，其中，
[0056]
步骤3.3.5、调整粒子i新的搜索窗口大小，具体如下：
[0057]
计算(x,y)的二阶矩m
20
、m
02
、m
11
，公式为：
[0058][0059][0060][0061]
根据(x,y)的二阶矩m
20
、m
02
、m
11
调整粒子i新搜索窗口，公式为：
[0062]
[0063][0064][0065]
其中l和w分别为粒子i新的搜索窗口的高和宽，θ为方向角；
[0066]
步骤3.3.6、判断此时点(x0,y0)与质心(x
c
,y
c
)的距离，若小于ε，或者n＞n
ε
，则停止迭代，并将当前搜索框的大小与质心位置作为下一帧图像的粒子i位置输入，然后进入步骤3.3.7；反之，则返回步骤3.3.3；其中ε和n
ε
分别为事先设定好的距离阈值和最大迭代次数；
[0067]
步骤3.3.7、判断是否每个粒子是否都进行漂移，若是，则输出漂移后的k时刻粒子集其中表示k时刻下编号为i的粒子；反之，则返回步骤3.3.2。
[0068]
步骤4、根据目标的h
‑
s联合颜色直方图和梯度方向直方图建立观测模型，通过观测模型判断漂移后的k时刻粒子集中每个粒子状态与实际状态的相似度，然后计算粒子权值并归一化，具体如下：
[0069]
步骤4.1、基于目标的h
‑
s联合颜色直方图和梯度方向直方图建立观测模型，具体如下：
[0070]
计算每个粒子的颜色直方图与初始目标框区域的颜色直方图q
(u)
之间的bhattacharyya距离d：
[0071][0072]
其中，u＝1,
…
,m，m为目标区域的颜色直方图的维数，ρ为bhattacharyya系数：
[0073][0074]
由bhattacharyya距离d来衡量直方图和q
(u)
之间的相似程度，d越小表示两种颜色直方图越相似，则赋予该粒子比较大的权值，d越大则相反；
[0075]
由此得粒子滤波的观测函数为：
[0076][0077]
其中z
k
为通过观测来判断每个粒子状态，粒子i实际状态，σ为距离d对应的方差。
[0078]
根据h
‑
s联合颜色直方图和目标区域的梯度方向直方图分别求得其对应的观测函数和通过自适应融合目标的h
‑
s联合颜色直方图的观测函数和目标区域的梯度方向直方图的得到最终改进后的观测函数：
[0079][0080]
其中β为加权系数；
[0081]
步骤4.2、根据观测函数更新漂移后的k时刻粒子集中粒子的权值为：
[0082][0083]
步骤4.3、归一化漂移后的k时刻粒子集中粒子的权值为：
[0084][0085]
其中，归一化漂移后粒子k时刻粒子i的权值，为步骤4.1所述最终改进后的观测函数。
[0086]
步骤5、对权值小于阈值e的粒子，利用随机采样法复制从粒子集中随机选取的粒子取代原粒子，完成粒子重采样，其中e＝1/n；
[0087]
对步骤4中权值小的粒子，利用随机采样法复制权值较大的粒子取代原粒子，完成粒子重采样，避免退化；
[0088]
步骤6、对重采样后的粒子加权求和，得到当前目标的估计状态，确定目标位置；
[0089]
当前时刻目标的估计状态的公式为：
[0090][0091]
其中，归一化漂移后粒子k时刻粒子i的权值，其中表示k时刻下编号为i的粒子，e[x
k
]为k时刻跟踪目标的估计位置。
[0092]
本发明解决了遮挡物与目标颜色相似时对跟踪效果的干扰问题，系统有较强的稳健性；在网格围栏与遮挡物较多的复杂背景中，跟踪过程基本不受网格干扰，跟踪效果好；能够实时反映目标框位于遮挡物体处时目标所处的位置；利用camshift迭代每个粒子使其漂移到目标区域，计算量小，实时性强。