# Online Hierarchical Sparse Representation of Multifeature for Robust Object Tracking

Authors: Honghong Yang; Shiru Qu
Journal: Computational Intelligence and Neuroscience
 (2016)

---

## Abstract

Object tracking based on sparse representation has given promising tracking results in recent years. However, the trackers under the framework of sparse representation always overemphasize the sparse representation and ignore the correlation of visual information. In addition, the sparse coding methods only encode the local region independently and ignore the spatial neighborhood information of the image. In this paper, we propose a robust tracking algorithm. Firstly, multiple complementary features are used to describe the object appearance; the appearance model of the tracked target is modeled by instantaneous and stable appearance features simultaneously. A two-stage sparse-coded method which takes the spatial neighborhood information of the image patch and the computation burden into consideration is used to compute the reconstructed object appearance. Then, the reliability of each tracker is measured by the tracking likelihood function of transient and reconstructed appearance models. Finally, the most reliable tracker is obtained by a well established particle filter framework; the training set and the template library are incrementally updated based on the current tracking results. Experiment results on different challenging video sequences show that the proposed algorithm performs well with superior tracking accuracy and robustness.


---

## Body

## 1. Introduction

The task of visual tracking is to find the interested object and track it. It is an important research in computer vision due to its widespread applications in traffic monitoring, vehicle navigation, and visual surveillance. Robust object tracking in dynamic environment is still a challenging problem. This is mainly because the factors such as occlusion, pose variation, illumination change, and clutter background cause large appearance change [1, 2].A robust appearance model is important for dealing with occlusions or other interferences in the tracking process. A target object is represented by its visual information like color, edge, or texture features extracted from the target region. However, there are numerous trackers only that rely on single feature to build target appearance, ignore the complementary representation of different features, usually lack of robustness, and are sensitive to interferences in dynamic environment [3]. For example, Ross et al. [4] use the intensity feature to represent the appearance model of the target object and integrate incremental learning to obtain a low-dimensional subspace representation. Babenko et al. [5] propose Multiple Instance Learning (MIL), which employs the Haar-like feature to build the discriminative appearance model for robust tracking. Mei and Ling [6] introduce a l
1 minimization robust visual tracking method, it uses the intensity feature to represent the target appearance, and the target appearance is represented by sparse linear combination of the appearance template and trivial template in template space. However, the single feature ignores the complementary characteristic of different visual information; it is insufficient to describe the drastic changes of target appearance in complicated environment. Therefore, the representation ability will decline when there are occlusions or other interferences in complex background [3]. As a result, numerous trackers are proposed to represent the object by fusing the multiple features to describe the target object and build the object appearance model, which can better describe the appearance changes and is beneficial to improve the robustness of trackers in dynamic environments [7–11]. However, how to effectively use and integrate multiple features for robust tracking should be tackled urgently.Numerous trackers based on sparse representation have been proposed in recent years [12–17]. Mei et al. [6, 17] proposed a l
1 minimization robust tracking method that regards the tracking as a sparse approximation problem. Zhang et al. [14] proposed a low-rank sparse representation tracking method. Liu et al. [15] developed a robust tracking algorithm using a local sparse appearance model to balance the requirements of stability and flexibility in the process of tracking. These trackers are all solved as the sparse approximation problem by l
1 regularized least squares method and show promising results against many existing trackers. However, this sparse coding based on l
1 minimization provides very sparse representation but ignores the importance of collaboration representation and the correlation of visual information; it is vulnerable to interferences and the l
1 minimization is very time-consuming. In [16], Zhang et al. emphasized the role of collaborative representation with l
2 regularized least squares, which shows that l
2
_
RLS is beneficial for reducing the computation burden. Shi et al. [18] also demonstrated that l
2
_
RLS is more accurate, robust, and faster. In [19], Yu et al. demonstrated that traditional sparse coding methods ignore the spatial neighborhood structure of the image because they only encode the local patches independently; then they proposed an efficient discriminative image representation method by using a two-layer sparse coding scheme at the pixel level.Inspired by the challenges mentioned above, this paper proposes an object tracking algorithm that combines the multiple visual features with hierarchical sparse coding to realize the tracking. As shown in Figure1, a multiple complementary feature representation [20] is used to robustly represent the object; the target appearance is modeled by exploiting a two-stage sparse-coded method, which is based on l
2 regularized least squares to solve the sparse approximation problem. Then, each tracker is based on the different features to estimate the object state and build the multiple observation models. The corresponding reliability of each tracker is computed by the tracking likelihood function of instantaneous and reconstructed appearance models that take the transient and stable appearance changes into consideration. Finally, the most reliable tracker is obtained by a well established particle filter framework; the training set and template library are incrementally updated based on the current tracking result.Figure 1
The framework of the proposed tracking algorithm.The main contributions of the proposed tracking algorithm are as follows:(
1
) we construct the target appearance by taking account of instantaneous and stable appearance features; then the transient appearance model and reconstructed object appearance model are built independently. (
2
) A two-stage sparse-coded method is employed to obtain the reconstructed coefficient vector used to construct the reconstructed appearance model. The two-stage sparse-coded method takes the temporal correlation between target templates and spatial neighborhood structure of the image patches into consideration and solves the sparse approximation problem by l
2
_
RLS. This is beneficial for reducing the computational burden and improving the tracking performance. (
3
) To better describe the object appearance changes, the reliability of each tracker is measured by the tracking likelihood function of instantaneous and reconstructed appearance model that take transient and stable appearance changes into consideration. Experimental results on challenging sequences show that the proposed method performs well compared to state-of-the-art methods.


## 2. The Proposed Tracking Algorithm

### 2.1. Particle Filter Tracking Framework

Given the object observations in previoustth frame z
1
:
t
=
z
1
,
z
2
,
…
,
z
t, the object state in tth frame is defined as x
t
=
[
x
t
,
y
t
,
s
t
,
θ
t
,
ε
t
,
ϕ
t
]
T, where x
t
,
y
t denote the coordinates, s
t
,
θ
t are the scale and aspect, and ε
t
,
ϕ
t are the rotation angle and skew. In order to robustly represent the object, we use multifeatures to build observation models of multiple trackers, let k
t
∈
1
,
…
,
K denote the index of K trackers from K features, then the ith tracker’s index is k
t
i.Theith tracker’s posterior distribution of the state x
t is(1)
p
x
t
∣
z
1
:
t
,
k
t
i
=
p
z
t
∣
x
t
,
k
t
i
p
x
t
∣
z
1
:
t
-
1
,
k
t
i
,where p
z
t
∣
x
t
,
k
t
i is the observation model and p
x
t
∣
z
1
:
t
-
1
,
k
t
i is the predicted distribution of ith tracker:(2)
p
x
t
∣
z
1
:
t
-
1
,
k
t
i
=
∫
p
x
t
∣
x
t
-
1
,
k
t
i
p
x
t
-
1
∣
z
1
:
t
-
1
,
k
t
i
d
x
t
-
1
,
(3)
p
x
t
-
1
∣
z
1
:
t
-
1
,
k
t
i
=
∑
j
=
1
K
p
x
t
-
1
∣
z
1
:
t
-
1
,
k
t
-
1
j
P
k
t
-
1
j
∣
k
t
i
,
z
1
:
t
-
1
,where p
x
t
∣
x
t
-
1
,
k
t
i is the motion model of ith tracker between the tth and (
t
-
1
)th frame, which is restricted to Gaussian distribution N
x
t
∣
x
t
-
1
σ. p
x
t
-
1
∣
z
1
:
t
-
1
,
k
t
i denotes the prior distribution up to frame t
-
1 and P
k
t
i
∣
z
1
:
t is the probability of the ith tracker.The crossover probability ofith tracker for multiple features is (4)
P
k
t
-
1
j
∣
k
t
i
,
z
1
:
t
-
1
=
P
k
t
i
∣
k
t
-
1
j
,
z
1
:
t
-
1
P
k
t
-
1
j
∣
z
1
:
t
-
1
∑
l
=
1
K
P
k
t
i
∣
k
t
-
1
l
,
z
1
:
t
-
1
P
k
t
-
1
l
∣
z
1
:
t
-
1
.In addition, theith tracker probability P
· satisfies(5)
∑
i
P
k
t
i
∣
z
1
:
t
=
1
,
∑
j
P
k
t
-
1
j
∣
k
t
i
,
z
1
:
t
-
1
=
1
.Then, we sparsely represent the candidate samplez
t with state x
t
i from the template library f
t
k; the likelihood of the observation model is(6)
p
z
t
∣
x
t
,
k
t
i
=
exp
⁡
-
ε
i
,where ε
i
=
min
⁡
f
α
-
z
t is the sparse reconstruction error of candidate sample z
t and α is the sparse coefficients.Therefore, the tracking resultx
^
t at the tth frame is the most reliable tracker with the highest tracker probability:(7)
x
^
t
=
arg
⁡
max
x
t
⁡
p
x
t
∣
z
1
:
t
,
k
^
t
,
k
^
t
=
arg
⁡
max
k
t
i
⁡
P
k
t
i
∣
z
1
:
t
,
i
=
1
,
…
,
K
.


### 2.2. Multiple Features Representation for Object Appearance

The different features have complementary characteristics to cope with appearance changes, such that the HOG features are robust to pose variations [21], Haar-like features can effectively deal with occlusions [22] as the single appearance model is insufficient to represent the target in a complicated environment. Therefore, we exploit different types of the features to build the multiple appearance models to represent the object robustly. The multiple features with complementary characteristics are used to handle various appearance changes, which is beneficial for tracking the target object robustly.In the proposed method, we use three trackers based on HOG, Haar-like feature, and intensity feature to represent the object appearance, which can effectively deal with occlusions, illumination changes, and pose variations. For thetth frame, we extract the multiple features to form feature sets as f
t
k
∈
R
m
k, where k is the index of the feature and m
k is the dimension of the kth feature. Normalize the feature sets f
t
k
∈
R
m
k
×
n to form the target template f
t
k
∈
R
d
k and d
k denotes the dimension of kth multiple features.


### 2.3. Object Representation by Hierarchical Sparse Coding

In the proposed method, we use the transient and stable features to describe the abrupt and stable object appearance changes. The stable features are sparsely represented by the current template with hierarchical sparse coding. Then, the reliability of each tracker is measured by the tracking likelihood function of instantaneous and reconstructed appearance models.The transient features up totth frame is f
I
,
t
k
=
[
f
I
,
t
-
l
k
,
…
,
f
I
,
t
-
1
k
]. Then the transient appearance model f
-
I
,
t
k is achieved by averaging the recent L appearance features as (8)
f
-
I
,
t
k
=
1
L
∑
l
=
1
L
f
I
,
t
-
l
k
.The stable object appearancez
t
i
,
k is represented by sparse coding the stable features f
R
,
t
k as (9)
z
t
i
,
k
≈
f
R
,
t
k
α
t
i
,
k
+
ε
i
,
k
=
f
1
,
t
k
α
1
,
t
i
,
k
+
f
2
,
t
k
α
2
,
t
i
,
k
+
⋯
+
f
r
,
t
k
α
r
,
t
i
,
k
.Because the tracking algorithm based on sparse representation is to find samples with minimal reconstruction errors from the templates library, a target can be reconstructed from several templates [23]. Therefore, there are only some features having the discriminative capability to separate the target from its background. In order to achieve the goals that discriminatively separate the target from its background and minimal reconstruction errors from its template library, we utilize the hierarchical sparse coding to minimize reconstruction errors and maximize the discriminative capability of features. In addition, we use l
2
_
RLS to solve the sparse approximation problem, which is beneficial for reducing the computation burden.Letf
1
=
f
R
,
t
k
I
k, α
1
=
α
t
i
,
k
β
t
i
,
k, where α
t
i
,
k
=
[
α
1
,
t
i
,
k
,
…
,
α
r
,
t
i
,
k
]
T
∈
R
r is the sparse coefficient vector, β
t
i
,
k
=
[
β
1
,
t
i
,
k
,
…
,
β
d
k
,
t
i
,
k
]
T
∈
R
d
k is the noise coefficient vector, and I
k
∈
R
d
k
×
d
k is an identity matrix. The candidate sample z
t
i
,
k is sparsely represented by linear combination of the features f
t
i
,
k as(10)
z
t
i
,
k
=
f
R
,
t
k
α
t
i
,
k
+
ε
i
,
k
=
f
R
,
t
k
I
α
t
i
,
k
β
t
i
,
k
.For a new arriving frame, we can achieveK tracking results x
^
t
i
∣
i
=
1
,
…
,
K. For ith tracker, z
t
i
,
k denotes the candidate image patch represented by kth features and f
R
,
t
k
α
t
i
,
kis the reconstructed appearance for z
t
i
,
k.Then, a two-stage sparse-coded method byl
2
_
RLS is used to obtain the coefficient vectors α
t
i
,
k and β
t
i
,
k as follows:(11)
α
1
=
a
r
g
m
i
n
α
t
i
,
k
,
β
t
i
,
k
f
1
α
1
-
z
t
i
,
k
2
,
s.t.
α
t
i
,
k
2
≤
K
1
,
β
t
i
,
k
2
≤
K
2
,where K
1 and K
2 are nonzero components.To effectively tackle the high-dimensional data in feature space, we use the diagonal matrixW to decrease the dimension of the feature space. For a set of samples X
=
x
t
i
∈
R
1
×
p
∣
i
=
1
,
…
,
K, the joint sparse solution is shown as follows:(12)
α
1
,
W
=
arg
⁡
min
α
1
,
W
⁡
λ
W
f
1
α
1
-
W
z
t
i
,
k
2
2
+
γ
F
W
,
X
+
τ
1
α
1
2
2
+
τ
2
diag
⁡
W
2
2
,where F
(
W
,
X
) is the loss function and τ
1
,
τ
2 are the sparse parameters. If W
i
i
≠
0, the ith feature is activated.The loss function is computed as(13)
F
W
,
X
=
e
-
∑
i
=
1
K
x
t
i
w
t
i
,where w
t
i
∈
R
p
×
1
∣
i
=
1
,
…
,
K is the sparse vector. If w
i
≠
0, the ith feature is selected.Then, the solution to the minimum loss functionF
(
W
,
X
) is achieved by solving the sparse problem as(14)
w
t
i
∗
=
a
r
g
m
i
n
w
t
i
X
w
t
i
2
,
s.t.
w
t
i
2
≤
K
0
,where K
0 denotes the maximum number of features that can be selected.Considering the spatial neighborhood information of the image patch, letN
w
t
i
(
i
,
j
) denote the jth neighbor of ith feature; then the vector set is (15)
z
t
i
=
w
t
i
2
+
∑
j
=
1
τ
θ
j
2
N
w
t
i
2
i
,
j
,
i
=
1
,
…
,
p
,where θ is the weight of the neighbors.The diagonal matrixW is formed as(16)
W
t
i
j
,
j
=
1
,
w
t
i
j
∗
≠
0
0
,
otherwise
.From the above first-stage sparse representation coding, we take account of the spatial relationship of neighborhood features, which is beneficial for selecting a set of discriminative features to separate target from its background and reducing the computational burden byl
2
_
RLS to solve the sparse approximation problem, as the target templates always contain some features from background, which is not the same as its neighbors. By doing discriminative feature selection as above, we can efficiently eliminate the features from background in the target templates. Therefore, we can construct a more efficient and robust target template library.In the second sparse reconstruction stage,α
t
i
,
k and β
t
i
,
k in (12) can be computed as follows: (17)
α
t
i
,
k
,
β
t
i
,
k
=
a
r
g
m
i
n
α
t
i
,
k
,
β
t
i
,
k
W
f
1
α
1
-
W
z
t
i
,
k
2
,
s.t.
α
t
i
,
k
2
≤
K
1
,
β
t
i
,
k
2
≤
K
2
.The nonzero row of matrixW forms the matrix W
′
∈
R
K
0
×
p; let f
1
′
=
W
f
1, z
t
′
=
W
′
z
t, and β
′
=
W
′
β.Then,(18)
α
t
i
,
k
∗
,
β
t
i
,
k
∗
=
a
r
g
m
i
n
α
t
i
,
k
,
β
t
i
,
k
f
1
′
,
W
′
α
t
i
,
k
β
t
i
,
k
-
z
t
′
2
,
s.t.
α
2
≤
K
1
,
β
2
≤
K
2
,where K
1
,
K
2 are the sparsity parameters that control the sparse representation of the target template and the tolerance of interference in complicated environment, respectively.Therefore, the reconstructed object appearancef
-
R
,
t
i
,
k for z
t
i
,
k is represented as(19)
f
-
R
,
t
i
,
k
=
f
R
,
t
i
,
k
α
t
i
,
k
.After above sparse reconstruction, the feature dimension reduced fromp
×
L to K
0
×
L and L is the number of templates in the target library.The predicted reliable object state forith tracker at frame t is(20)
x
^
t
i
=
arg
⁡
max
x
t
⁡
p
x
t
∣
z
1
:
t
,
k
t
i
.Then the corresponding tracking likelihood function of theith tracker at frame t is (21)
p
z
t
∣
k
t
i
,
z
1
:
t
-
1
=
p
z
t
∣
x
^
t
i
.In the proposed method, we use instantaneous and reconstructed features to describe the transient and stable appearance changes. The reliability of each tracker is(22)
p
z
t
∣
x
^
t
i
=
p
I
z
t
∣
x
^
t
i
p
R
z
t
∣
x
^
t
i
=
∏
k
=
1
K
p
z
t
∣
x
^
t
i
,
f
I
,
t
k
p
z
t
∣
x
^
t
i
,
f
R
,
t
k
,where p
I
z
t
∣
x
^
t
i is the instantaneous appearance likelihood based on the transient object appearance f
-
I
,
t
k, which is formed by a set of recent frames features f
I
,
t
k. p
R
z
t
∣
x
^
t
i is the reconstructed object appearance likelihood based on the reconstructed object appearance f
-
R
,
t
i
,
k and f
-
R
,
t
i
,
k comes from the stable object appearance, which is formed by kth feature and the tracking result z
t
i
,
k from the ith tracker:(23)
p
I
z
t
∣
x
^
t
i
,
f
R
,
t
i
,
k
=
exp
⁡
-
ρ
f
-
I
,
t
k
-
z
t
i
,
k
2
,
p
R
z
t
∣
x
^
t
i
,
f
R
,
t
i
,
k
=
exp
⁡
-
ρ
f
-
R
,
t
i
,
k
-
z
t
i
,
k
2
,where ρ is the control parameter.


### 2.4. Predication and Update

To robustly track the target object, we update the tracker probability of the multiple trackers and the reliability of each tracker.The tracker probability is updated as follows:(24)
P
k
t
i
∣
z
1
:
t
=
p
z
t
∣
k
t
i
,
z
1
:
t
-
1
p
z
t
∣
z
1
:
t
-
1
P
k
t
i
∣
z
1
:
t
-
1
,where(25)
P
k
t
i
∣
z
1
:
t
-
1
=
∑
j
=
1
K
P
k
t
i
∣
k
t
-
1
j
,
z
1
:
t
-
1
P
k
t
-
1
j
∣
z
1
:
t
-
1
p
z
t
∣
z
1
:
t
-
1
=
∑
i
=
1
K
p
z
t
∣
k
t
i
,
z
1
:
t
-
1
∑
j
=
1
K
P
k
t
i
∣
k
t
-
1
j
,
z
1
:
t
-
1
P
k
t
-
1
j
∣
z
1
:
t
-
1
.The corresponding observation modelp
z
t
∣
x
t
,
k
t
i for each tracker is updated based on the incremental subspace model in [4]. Then, the particle filter is used to approximate the state posterior distribution p
x
t
∣
z
1
:
t
,
k
t
i by a set of N particles [24], the particles size N
=
600: (26)
p
x
t
∣
z
1
:
t
,
k
t
i
≈
∑
q
=
1
N
w
q
,
t
i
δ
x
q
,
t
i
-
x
t
,where δ
(
·
) is a delta function and w
q
,
t
i
q
=
1
N is the sample weight associated with x
q
,
t
i
q
=
1
,
…
,
N.The particlesx
q
,
t
i are obtained from state prediction p
x
t
∣
x
t
-
1
,
k
t
i, which is simplified by first-order Markov model(27)
x
q
,
t
i
~
p
x
t
∣
x
t
-
1
,
k
t
i
.The weights are updated as(28)
w
q
,
t
i
=
w
q
,
t
-
1
i
p
z
t
∣
x
t
,
k
t
i
.Then, we achieve a set ofK reliable states by maximizing the posterior estimates:(29)
x
^
t
i
=
x
q
^
,
t
i
,
q
^
=
arg
⁡
max
q
⁡
w
q
,
t
i
∣
i
=
1
,
…
,
K
,
q
=
1
,
…
,
N
.In the proposed method, the target appearance is constructed by multiple features that take account of transient and stable appearance changes to cope with occlusion and other interferences in complicated environments. For example, in a dynamic environment with drastic occlusion or illumination changes, the stable features are rarely updated, but the transient features can effectively describe the frequent appearance changes, while in a static background, if a background sample is added into the template, it usually has a good reconstruction with high likelihood because background is static at most of the time. Because the incorrect template is nonlinear, which is not the same as its neighbors, the two-stage sparse coding method taking account the spatial relationship of neighborhood features can prevent it from being selected. Therefore, we can construct a more efficient and robust target template library.In addition, we update the template library based on the current tracking result as done in IVT method [4]; the samples with high likelihood and near the target are added to the template library. We repeat this procedure for each frame in the entire sequences. The tracking based on joint multiple feature representation and hierarchical sparse coding can provide a robust and accurate tracking result.


## 2.1. Particle Filter Tracking Framework

Given the object observations in previoustth frame z
1
:
t
=
z
1
,
z
2
,
…
,
z
t, the object state in tth frame is defined as x
t
=
[
x
t
,
y
t
,
s
t
,
θ
t
,
ε
t
,
ϕ
t
]
T, where x
t
,
y
t denote the coordinates, s
t
,
θ
t are the scale and aspect, and ε
t
,
ϕ
t are the rotation angle and skew. In order to robustly represent the object, we use multifeatures to build observation models of multiple trackers, let k
t
∈
1
,
…
,
K denote the index of K trackers from K features, then the ith tracker’s index is k
t
i.Theith tracker’s posterior distribution of the state x
t is(1)
p
x
t
∣
z
1
:
t
,
k
t
i
=
p
z
t
∣
x
t
,
k
t
i
p
x
t
∣
z
1
:
t
-
1
,
k
t
i
,where p
z
t
∣
x
t
,
k
t
i is the observation model and p
x
t
∣
z
1
:
t
-
1
,
k
t
i is the predicted distribution of ith tracker:(2)
p
x
t
∣
z
1
:
t
-
1
,
k
t
i
=
∫
p
x
t
∣
x
t
-
1
,
k
t
i
p
x
t
-
1
∣
z
1
:
t
-
1
,
k
t
i
d
x
t
-
1
,
(3)
p
x
t
-
1
∣
z
1
:
t
-
1
,
k
t
i
=
∑
j
=
1
K
p
x
t
-
1
∣
z
1
:
t
-
1
,
k
t
-
1
j
P
k
t
-
1
j
∣
k
t
i
,
z
1
:
t
-
1
,where p
x
t
∣
x
t
-
1
,
k
t
i is the motion model of ith tracker between the tth and (
t
-
1
)th frame, which is restricted to Gaussian distribution N
x
t
∣
x
t
-
1
σ. p
x
t
-
1
∣
z
1
:
t
-
1
,
k
t
i denotes the prior distribution up to frame t
-
1 and P
k
t
i
∣
z
1
:
t is the probability of the ith tracker.The crossover probability ofith tracker for multiple features is (4)
P
k
t
-
1
j
∣
k
t
i
,
z
1
:
t
-
1
=
P
k
t
i
∣
k
t
-
1
j
,
z
1
:
t
-
1
P
k
t
-
1
j
∣
z
1
:
t
-
1
∑
l
=
1
K
P
k
t
i
∣
k
t
-
1
l
,
z
1
:
t
-
1
P
k
t
-
1
l
∣
z
1
:
t
-
1
.In addition, theith tracker probability P
· satisfies(5)
∑
i
P
k
t
i
∣
z
1
:
t
=
1
,
∑
j
P
k
t
-
1
j
∣
k
t
i
,
z
1
:
t
-
1
=
1
.Then, we sparsely represent the candidate samplez
t with state x
t
i from the template library f
t
k; the likelihood of the observation model is(6)
p
z
t
∣
x
t
,
k
t
i
=
exp
⁡
-
ε
i
,where ε
i
=
min
⁡
f
α
-
z
t is the sparse reconstruction error of candidate sample z
t and α is the sparse coefficients.Therefore, the tracking resultx
^
t at the tth frame is the most reliable tracker with the highest tracker probability:(7)
x
^
t
=
arg
⁡
max
x
t
⁡
p
x
t
∣
z
1
:
t
,
k
^
t
,
k
^
t
=
arg
⁡
max
k
t
i
⁡
P
k
t
i
∣
z
1
:
t
,
i
=
1
,
…
,
K
.


## 2.2. Multiple Features Representation for Object Appearance

The different features have complementary characteristics to cope with appearance changes, such that the HOG features are robust to pose variations [21], Haar-like features can effectively deal with occlusions [22] as the single appearance model is insufficient to represent the target in a complicated environment. Therefore, we exploit different types of the features to build the multiple appearance models to represent the object robustly. The multiple features with complementary characteristics are used to handle various appearance changes, which is beneficial for tracking the target object robustly.In the proposed method, we use three trackers based on HOG, Haar-like feature, and intensity feature to represent the object appearance, which can effectively deal with occlusions, illumination changes, and pose variations. For thetth frame, we extract the multiple features to form feature sets as f
t
k
∈
R
m
k, where k is the index of the feature and m
k is the dimension of the kth feature. Normalize the feature sets f
t
k
∈
R
m
k
×
n to form the target template f
t
k
∈
R
d
k and d
k denotes the dimension of kth multiple features.


## 2.3. Object Representation by Hierarchical Sparse Coding

In the proposed method, we use the transient and stable features to describe the abrupt and stable object appearance changes. The stable features are sparsely represented by the current template with hierarchical sparse coding. Then, the reliability of each tracker is measured by the tracking likelihood function of instantaneous and reconstructed appearance models.The transient features up totth frame is f
I
,
t
k
=
[
f
I
,
t
-
l
k
,
…
,
f
I
,
t
-
1
k
]. Then the transient appearance model f
-
I
,
t
k is achieved by averaging the recent L appearance features as (8)
f
-
I
,
t
k
=
1
L
∑
l
=
1
L
f
I
,
t
-
l
k
.The stable object appearancez
t
i
,
k is represented by sparse coding the stable features f
R
,
t
k as (9)
z
t
i
,
k
≈
f
R
,
t
k
α
t
i
,
k
+
ε
i
,
k
=
f
1
,
t
k
α
1
,
t
i
,
k
+
f
2
,
t
k
α
2
,
t
i
,
k
+
⋯
+
f
r
,
t
k
α
r
,
t
i
,
k
.Because the tracking algorithm based on sparse representation is to find samples with minimal reconstruction errors from the templates library, a target can be reconstructed from several templates [23]. Therefore, there are only some features having the discriminative capability to separate the target from its background. In order to achieve the goals that discriminatively separate the target from its background and minimal reconstruction errors from its template library, we utilize the hierarchical sparse coding to minimize reconstruction errors and maximize the discriminative capability of features. In addition, we use l
2
_
RLS to solve the sparse approximation problem, which is beneficial for reducing the computation burden.Letf
1
=
f
R
,
t
k
I
k, α
1
=
α
t
i
,
k
β
t
i
,
k, where α
t
i
,
k
=
[
α
1
,
t
i
,
k
,
…
,
α
r
,
t
i
,
k
]
T
∈
R
r is the sparse coefficient vector, β
t
i
,
k
=
[
β
1
,
t
i
,
k
,
…
,
β
d
k
,
t
i
,
k
]
T
∈
R
d
k is the noise coefficient vector, and I
k
∈
R
d
k
×
d
k is an identity matrix. The candidate sample z
t
i
,
k is sparsely represented by linear combination of the features f
t
i
,
k as(10)
z
t
i
,
k
=
f
R
,
t
k
α
t
i
,
k
+
ε
i
,
k
=
f
R
,
t
k
I
α
t
i
,
k
β
t
i
,
k
.For a new arriving frame, we can achieveK tracking results x
^
t
i
∣
i
=
1
,
…
,
K. For ith tracker, z
t
i
,
k denotes the candidate image patch represented by kth features and f
R
,
t
k
α
t
i
,
kis the reconstructed appearance for z
t
i
,
k.Then, a two-stage sparse-coded method byl
2
_
RLS is used to obtain the coefficient vectors α
t
i
,
k and β
t
i
,
k as follows:(11)
α
1
=
a
r
g
m
i
n
α
t
i
,
k
,
β
t
i
,
k
f
1
α
1
-
z
t
i
,
k
2
,
s.t.
α
t
i
,
k
2
≤
K
1
,
β
t
i
,
k
2
≤
K
2
,where K
1 and K
2 are nonzero components.To effectively tackle the high-dimensional data in feature space, we use the diagonal matrixW to decrease the dimension of the feature space. For a set of samples X
=
x
t
i
∈
R
1
×
p
∣
i
=
1
,
…
,
K, the joint sparse solution is shown as follows:(12)
α
1
,
W
=
arg
⁡
min
α
1
,
W
⁡
λ
W
f
1
α
1
-
W
z
t
i
,
k
2
2
+
γ
F
W
,
X
+
τ
1
α
1
2
2
+
τ
2
diag
⁡
W
2
2
,where F
(
W
,
X
) is the loss function and τ
1
,
τ
2 are the sparse parameters. If W
i
i
≠
0, the ith feature is activated.The loss function is computed as(13)
F
W
,
X
=
e
-
∑
i
=
1
K
x
t
i
w
t
i
,where w
t
i
∈
R
p
×
1
∣
i
=
1
,
…
,
K is the sparse vector. If w
i
≠
0, the ith feature is selected.Then, the solution to the minimum loss functionF
(
W
,
X
) is achieved by solving the sparse problem as(14)
w
t
i
∗
=
a
r
g
m
i
n
w
t
i
X
w
t
i
2
,
s.t.
w
t
i
2
≤
K
0
,where K
0 denotes the maximum number of features that can be selected.Considering the spatial neighborhood information of the image patch, letN
w
t
i
(
i
,
j
) denote the jth neighbor of ith feature; then the vector set is (15)
z
t
i
=
w
t
i
2
+
∑
j
=
1
τ
θ
j
2
N
w
t
i
2
i
,
j
,
i
=
1
,
…
,
p
,where θ is the weight of the neighbors.The diagonal matrixW is formed as(16)
W
t
i
j
,
j
=
1
,
w
t
i
j
∗
≠
0
0
,
otherwise
.From the above first-stage sparse representation coding, we take account of the spatial relationship of neighborhood features, which is beneficial for selecting a set of discriminative features to separate target from its background and reducing the computational burden byl
2
_
RLS to solve the sparse approximation problem, as the target templates always contain some features from background, which is not the same as its neighbors. By doing discriminative feature selection as above, we can efficiently eliminate the features from background in the target templates. Therefore, we can construct a more efficient and robust target template library.In the second sparse reconstruction stage,α
t
i
,
k and β
t
i
,
k in (12) can be computed as follows: (17)
α
t
i
,
k
,
β
t
i
,
k
=
a
r
g
m
i
n
α
t
i
,
k
,
β
t
i
,
k
W
f
1
α
1
-
W
z
t
i
,
k
2
,
s.t.
α
t
i
,
k
2
≤
K
1
,
β
t
i
,
k
2
≤
K
2
.The nonzero row of matrixW forms the matrix W
′
∈
R
K
0
×
p; let f
1
′
=
W
f
1, z
t
′
=
W
′
z
t, and β
′
=
W
′
β.Then,(18)
α
t
i
,
k
∗
,
β
t
i
,
k
∗
=
a
r
g
m
i
n
α
t
i
,
k
,
β
t
i
,
k
f
1
′
,
W
′
α
t
i
,
k
β
t
i
,
k
-
z
t
′
2
,
s.t.
α
2
≤
K
1
,
β
2
≤
K
2
,where K
1
,
K
2 are the sparsity parameters that control the sparse representation of the target template and the tolerance of interference in complicated environment, respectively.Therefore, the reconstructed object appearancef
-
R
,
t
i
,
k for z
t
i
,
k is represented as(19)
f
-
R
,
t
i
,
k
=
f
R
,
t
i
,
k
α
t
i
,
k
.After above sparse reconstruction, the feature dimension reduced fromp
×
L to K
0
×
L and L is the number of templates in the target library.The predicted reliable object state forith tracker at frame t is(20)
x
^
t
i
=
arg
⁡
max
x
t
⁡
p
x
t
∣
z
1
:
t
,
k
t
i
.Then the corresponding tracking likelihood function of theith tracker at frame t is (21)
p
z
t
∣
k
t
i
,
z
1
:
t
-
1
=
p
z
t
∣
x
^
t
i
.In the proposed method, we use instantaneous and reconstructed features to describe the transient and stable appearance changes. The reliability of each tracker is(22)
p
z
t
∣
x
^
t
i
=
p
I
z
t
∣
x
^
t
i
p
R
z
t
∣
x
^
t
i
=
∏
k
=
1
K
p
z
t
∣
x
^
t
i
,
f
I
,
t
k
p
z
t
∣
x
^
t
i
,
f
R
,
t
k
,where p
I
z
t
∣
x
^
t
i is the instantaneous appearance likelihood based on the transient object appearance f
-
I
,
t
k, which is formed by a set of recent frames features f
I
,
t
k. p
R
z
t
∣
x
^
t
i is the reconstructed object appearance likelihood based on the reconstructed object appearance f
-
R
,
t
i
,
k and f
-
R
,
t
i
,
k comes from the stable object appearance, which is formed by kth feature and the tracking result z
t
i
,
k from the ith tracker:(23)
p
I
z
t
∣
x
^
t
i
,
f
R
,
t
i
,
k
=
exp
⁡
-
ρ
f
-
I
,
t
k
-
z
t
i
,
k
2
,
p
R
z
t
∣
x
^
t
i
,
f
R
,
t
i
,
k
=
exp
⁡
-
ρ
f
-
R
,
t
i
,
k
-
z
t
i
,
k
2
,where ρ is the control parameter.


## 2.4. Predication and Update

To robustly track the target object, we update the tracker probability of the multiple trackers and the reliability of each tracker.The tracker probability is updated as follows:(24)
P
k
t
i
∣
z
1
:
t
=
p
z
t
∣
k
t
i
,
z
1
:
t
-
1
p
z
t
∣
z
1
:
t
-
1
P
k
t
i
∣
z
1
:
t
-
1
,where(25)
P
k
t
i
∣
z
1
:
t
-
1
=
∑
j
=
1
K
P
k
t
i
∣
k
t
-
1
j
,
z
1
:
t
-
1
P
k
t
-
1
j
∣
z
1
:
t
-
1
p
z
t
∣
z
1
:
t
-
1
=
∑
i
=
1
K
p
z
t
∣
k
t
i
,
z
1
:
t
-
1
∑
j
=
1
K
P
k
t
i
∣
k
t
-
1
j
,
z
1
:
t
-
1
P
k
t
-
1
j
∣
z
1
:
t
-
1
.The corresponding observation modelp
z
t
∣
x
t
,
k
t
i for each tracker is updated based on the incremental subspace model in [4]. Then, the particle filter is used to approximate the state posterior distribution p
x
t
∣
z
1
:
t
,
k
t
i by a set of N particles [24], the particles size N
=
600: (26)
p
x
t
∣
z
1
:
t
,
k
t
i
≈
∑
q
=
1
N
w
q
,
t
i
δ
x
q
,
t
i
-
x
t
,where δ
(
·
) is a delta function and w
q
,
t
i
q
=
1
N is the sample weight associated with x
q
,
t
i
q
=
1
,
…
,
N.The particlesx
q
,
t
i are obtained from state prediction p
x
t
∣
x
t
-
1
,
k
t
i, which is simplified by first-order Markov model(27)
x
q
,
t
i
~
p
x
t
∣
x
t
-
1
,
k
t
i
.The weights are updated as(28)
w
q
,
t
i
=
w
q
,
t
-
1
i
p
z
t
∣
x
t
,
k
t
i
.Then, we achieve a set ofK reliable states by maximizing the posterior estimates:(29)
x
^
t
i
=
x
q
^
,
t
i
,
q
^
=
arg
⁡
max
q
⁡
w
q
,
t
i
∣
i
=
1
,
…
,
K
,
q
=
1
,
…
,
N
.In the proposed method, the target appearance is constructed by multiple features that take account of transient and stable appearance changes to cope with occlusion and other interferences in complicated environments. For example, in a dynamic environment with drastic occlusion or illumination changes, the stable features are rarely updated, but the transient features can effectively describe the frequent appearance changes, while in a static background, if a background sample is added into the template, it usually has a good reconstruction with high likelihood because background is static at most of the time. Because the incorrect template is nonlinear, which is not the same as its neighbors, the two-stage sparse coding method taking account the spatial relationship of neighborhood features can prevent it from being selected. Therefore, we can construct a more efficient and robust target template library.In addition, we update the template library based on the current tracking result as done in IVT method [4]; the samples with high likelihood and near the target are added to the template library. We repeat this procedure for each frame in the entire sequences. The tracking based on joint multiple feature representation and hierarchical sparse coding can provide a robust and accurate tracking result.


## 3. Tracking Based on Online Hierarchical Sparse Representation of Multifeature

As described above, the main step of the proposed tracking algorithm is shown in Algorithm1.Algorithm 1 (tracking based on hierarchical sparse representation of multifeature).
Input. There are the initial states of target x
0
i
=
x
0
∣
i
=
1
,
…
,
K
.
Initializing. Construct L training samples X
∈
R
L
×
p, the set of samples for particle filter is x
q
,
0
i
,
w
q
,
0
i
=
1
/
N
∣
q
=
1
,
…
,
N, and the tracker probability is P
0
i
·
=
1
/
K
∣
i
=
1
,
…
,
K. Fort
=
1 to the end of video sequence, consider the following:
Fori
=
1
:
K,
(
1
)
achieve the solution to minimize loss functionF
(
W
,
X
) by (14),
(
2
)
construct diagonal matrixW by (16),
(
3
)
for candidate samplez
i in state x
t
i, achieve the sparse representation coefficients α
t
i
,
k and β
t
i
,
k by performing (18),
(
4
)
predict the reliable object state of each tracker at framet by (20),
(
5
)
compute the state posterior distributionp
x
t
∣
z
1
:
t
,
k
t
i by a set of N particles as (26),
(
6
)
predict state samplesx
q
,
t
i
,
w
q
,
t
i
q
=
1
N by (27),
(
7
)
update the sample weightsw
q
,
t
i as (28),
(
8
)
achieve a set ofK reliable states x
^
t
i from the ith tracking by (29),
End
(
9
)
compute the tracking likelihood functionp
z
t
∣
k
t
i
,
z
1
:
t
-
1 of the ith tracker at frame t by (21),
(
10
)
update the tracker probabilityP
k
t
i
∣
z
1
:
t using (24),
(
11
)
the tracking result at thetth frame is achieved by (7),
(
12
)
update the training set and template library with the tracking results.
End.


## 4. Experiments

To analyze the performance of the proposed tracking method, we compared our method with other five state-of-the-art trackers [25] such as IVT [4], L1 [6], MIL [5], OAB [26], and VTD [27] on several challenging video sequences. The target objects in the test videos are either nonrigid or rigid objects that suffered significant pose variation, heavy occlusion, in-plane and out of plane rotation, or motion blur. The video sequences are available in https://sites.google.com/site/trackerbenchmark/benchmarks/v10. The proposed tracker algorithm is implemented in MATLAB, which is run on a PC with 2CPU, 2.5 GHz, and 3.1 GB RAM, at around 20 frames per second.


### 4.1. Parameters Setting

For all test video sequences, we manually select the initial target location. Each image patch is normalized to32
∗
32 pixels and sparsity parameters τ
1
=
τ
2
=
0.001 and γ
=
0.1 and the dimensions of intensity features, HOG features, and Haar-like feature is 1024, 1296, and 1760, respectively. The number of particles is N
=
600, and the number of template samples is L
=
16. Table 1 lists the characteristics of the evaluated sequences used in the experiments of this paper.Table 1
Tracking sequence in the experiments.
Sequence
Frame
Main challenge
(1) Car4
659
Occlusion and illumination change
(2) CarScale
252
Occlusion and scale change
(3) CarDark
393
Illumination change, occlusion, and motion blur
(4) FaceOcc2
812
Occlusion and in-plane rotation
(5) Freeman1
326
Scale change and view change
(6) Girl
500
Out of plane rotation, occlusion, and scale and pose change
(7) Shaking
365
Illumination change and pose change
(8) Woman
569
Scale change, view variation, and occlusion
(9) Jogging
307
Occlusion, abrupt motion, and scale change


### 4.2. Qualitative Comparison

Experiment 1 (illumination variation, occlusion, scale change, and fast motion of rigid object).
The sequence of Car4 is to track a car in an open road with illumination variation and partial occlusion as shown in Figure2(a). At frame 86, the OAB tracker appears to slightly drift due to the trees and bridge occlusion and fails to track the car at frame 233. The L1, MIL, and VTD trackers start to drift away from the target when drastic illumination changes occur at frame 195 and fail to track the target at frame 255. The IVT and the proposed method can successfully track the target because they dynamically updated the template, which is beneficial for coping with the occlusion and illumination changes. However, the result of IVT is less satisfied because the tracking box is larger than the target object from frame 195 to the end sequences.Figure 2
The tracking results of the rigid object undergoing severe occlusion, illumination, and scale change.
(a)
Car4
(b)
CarScale
(c)
CarDarkIn the CarScale sequence, the tracking target is a fast motion car in an open road. Compared with the Car4 sequence, this sequence is more challenging because the tracked car undergoes large scale changes and fast motion on the entire sequence. Due to the fast motion accompanied with the tree’s occlusion, IVT, L1, MIL, and VTD trackers drift with different degree at frame 164 and finally lost the target at frame 171. The proposed method gives the best results followed by the OAB tracker.The CarDark sequence is challenging because the target object undergoes the motion blur in a night environment with low contrast and strong reflection interference. Due to the strong reflection interference, the MIL tracker drifts a little from the target at frame 122 and lost the target at frame 202 and then regards the other car as the tracked target. The IVT, L1, and VTD trackers drift away from the target at frame 277. The OAB tracker performs well on this sequence and yields the second best results. The proposed method can accurately track the target object in the whole sequence with small center position error and high overlap rate.Experiment 2 (occlusion, scale change, and rotation of nonrigid object).
The target object in FaceOcc2 sequence undergoes the drastic occlusion and in-plane rotation. As shown in Figure3(a), when there is a small occlusion with a book at frames 128~185 and frames 245~279, all methods perform well. But when partial occlusion and in-plane rotation occur together at frames 392~510, most of trackers have poor performances. When the target almost fully occludes by a book and a hat at frames 688~740, all methods except this paper method drift away from the target at different degrees. Since the proposed method uses multiple complementary features to build transient and stable appearance models and update the template library online, it can effectively handle the occlusion and give satisfactory tracking results.Figure 3
The tracking results of object undergoing occlusion, scale change, and rotation.
(a)
FaceOcc2
(b)
Freeman1
(c)
GirlThe Freeman1 sequence is challenging because the interested man’s face undergoes large scale changes and view variations. Due to the large scale changes, MIL drifts away from the target at frame 32. As the view changes from the left to right, the L1 and OAB trackers totally lost the target at frames 131 and 176, respectively. The tracking methods like IVT, VTD, and the proposed method perform well on this challenging sequence and can track the target accurately.The Girl sequence has drastic appearance changes because of the out of plane rotation and similar target occlusion. When out of plane rotation occurs at frames 90~122 and 169~260, IVT tracker totally fails to track the Girl’s face; other trackers drift at different degrees. The OAB and VTD trackers fail to track the target object and track the similar target when the Girl’s face is occluded by a man’s face at frames 420~470. MIL tracker can successfully track the target except some errors like frames 303 and 433. L1 tracker and the proposed method perform well on this sequence.Experiment 3 (illumination, scale change, and occlusion of nonrigid object).
The track target in shaking sequence undergoes drastic illumination and poses changes on the whole video sequence. It brings more challenges to accurately track the target because the color of object appearance is similar to the stage lighting. IVT and OAB almost fail to track the target at frame 23 and cannot recover at the rest frames. The MIL tracker drifts a little at frame 61 due to the drastic illumination changes. Although the stage lights change drastically accompanied with the serious head shaking, the L1 and VTD trackers perform well except some errors. The proposed method can effectively adapt to the severe object appearance changes when those variations occur together and achieve satisfactory results.The Woman sequence is very challenging because the target object undergoes large scale changes, view variations, and occlusions simultaneously. As shown in Figure4(b), all trackers merely perform well except the proposed method. The results of the proposed method show slight drift; other methods like IVT, L1, MIL, OAB, and VTD trackers totally lose the target when heavy occlusion occurs at frame 130 and never recover to track the target in the subsequent video sequence except the OAB tracker. The OAB tracker recaptures the target at frame 337 and keeps to track the target until the end of sequence with a little drift.Figure 4
The tracking results of object undergoing severe occlusion, illumination, and pose changes.
(a)
Shaking
(b)
Woman
(c)
JoggingThe Jogging sequence is more challenging to track because the tracked target is fully occluded by a stem and undergoes large scale change and fast motion simultaneously. The IVT, L1, MIL, OAB, and VTD trackers completely fail to track the target when the target is fully occluded by a stem at frames 50~62, and the OAB tracker recaptures the target at frame 106. The proposed method can accurately track the target on the entire sequence.From some sampled tracking results of the proposed method and other five methods on 9 image sequences, we can conclude that the algorithm in this paper can accurately and robustly track the target under the environment with illumination variation, scale change, and motion blur.


### 4.3. Quantitative Comparison

Two metrics are used to evaluate the proposed tracker with reference trackers in gray-scale videos. The first is the center position error, which is applied to evaluate the distance between the ground-truth centerR
g and tracked object center R
t (in pixels) at each frame by the Euclidean distance. The other metric is the overlap rate [28], which is defined as score
=
area
(
R
t
∩
R
g
)
/
a
r
e
a
(
R
t
∪
R
g
), where R
t denotes the bounding box generated by a tracking method and R
g is the ground-truth bounding box.Table2 and Table 3 show the average center position errors and the average overlap rate for all trackers. Figures 5 and 6 show the center position error curve and overlap rate evaluation curve of different trackers on 9 video sequences at each frame. It can be seen that the proposed algorithm has the optimal or suboptimal performance in terms of average center position errors and average overlap rate in most test video sequences compared with other methods. Most competing tracking methods do not give a satisfactory result; the center position error is larger and the overlap rate is lower. The average position error of this paper at 9 videos is only 5.53 pixels, which is far less than other trackers; the average central position errors of other trackers like IVT, L1, MIL, OAB, and VTD trackers are 59.58 pixels, 76.14 pixels, 51.91 pixels, 46.01 pixels, and 41.54 pixels, respectively. The average overlap rate of the proposed method reaches 73.6%, which is higher than other trackers; the average overlap rates of IVT, L1, MIL, OAB, and VTD trackers are 39.9%, 34%, 40%, 33.7%, and 32.7%, respectively. They both highlight the advantages of the algorithm in this paper. Overall, the effectiveness of the proposed tracker method is verified.Table 2
The average position errors (pixels).
Video
IVT
L1
MIL
OAB
VTD
Our method
Car4
16.31
209.22
55.10
103.49
34.46
2.97
CarScale
42.64
53.52
43.63
15.66
29.35
8.44
CarDark
22.03
19.85
48.41
3.90
23.29
1.92
FaceOcc2
16.21
15.54
18.17
24.32
10.69
8.42
Freeman1
6.85
61.74
18.98
25.44
10.71
7.44
Girl
29.37
3.76
14.11
8.55
11.70
3.63
Shaking
87.90
21.05
9.59
144.66
13.42
7.11
Woman
186.70
151.79
128.08
63.76
119.01
3.27
Jogging
128.29
148.81
131.14
24.29
121.31
6.54
Average
59.58
76.14
51.91
46.01
41.54
5.53
Note: the optimal result is shown as bold and the suboptimal one as italic.Table 3
The average overlap rate for all of trackers.
Video
IVT
L1
MIL
OAB
VTD
Our method
Car4
0.570
0.196
0.257
0.151
0.357
0.873
CarScale
0.469
0.515
0.421
0.451
0.408
0.768
CarDark
0.480
0.513
0.153
0.741
0.461
0.868
FaceOcc2
0.572
0.641
0.611
0.549
0.70
0.723
Freeman1
0.493
0.234
0.241
0.288
0.343
0.587
Girl
0.147
0.676
0.404
0.577
0.521
0.710
Shaking
0.032
0.50
0.663
0.017
0.613
0.735
Woman
0.129
0.145
0.147
0.217
0.144
0.654
Jogging
0.139
0.142
0.127
0.60
0.134
0.711
Average
0.337
0.396
0.336
0.399
0.409
0.736
Note: the optimal result is shown as bold and the suboptimal one as italic.Figure 5
The center position errors of different trackers on 9 video sequences.Figure 6
The overlap rate evaluation of different trackers on 9 video sequences in each frame.


## 4.1. Parameters Setting

For all test video sequences, we manually select the initial target location. Each image patch is normalized to32
∗
32 pixels and sparsity parameters τ
1
=
τ
2
=
0.001 and γ
=
0.1 and the dimensions of intensity features, HOG features, and Haar-like feature is 1024, 1296, and 1760, respectively. The number of particles is N
=
600, and the number of template samples is L
=
16. Table 1 lists the characteristics of the evaluated sequences used in the experiments of this paper.Table 1
Tracking sequence in the experiments.
Sequence
Frame
Main challenge
(1) Car4
659
Occlusion and illumination change
(2) CarScale
252
Occlusion and scale change
(3) CarDark
393
Illumination change, occlusion, and motion blur
(4) FaceOcc2
812
Occlusion and in-plane rotation
(5) Freeman1
326
Scale change and view change
(6) Girl
500
Out of plane rotation, occlusion, and scale and pose change
(7) Shaking
365
Illumination change and pose change
(8) Woman
569
Scale change, view variation, and occlusion
(9) Jogging
307
Occlusion, abrupt motion, and scale change


## 4.2. Qualitative Comparison

Experiment 1 (illumination variation, occlusion, scale change, and fast motion of rigid object).
The sequence of Car4 is to track a car in an open road with illumination variation and partial occlusion as shown in Figure2(a). At frame 86, the OAB tracker appears to slightly drift due to the trees and bridge occlusion and fails to track the car at frame 233. The L1, MIL, and VTD trackers start to drift away from the target when drastic illumination changes occur at frame 195 and fail to track the target at frame 255. The IVT and the proposed method can successfully track the target because they dynamically updated the template, which is beneficial for coping with the occlusion and illumination changes. However, the result of IVT is less satisfied because the tracking box is larger than the target object from frame 195 to the end sequences.Figure 2
The tracking results of the rigid object undergoing severe occlusion, illumination, and scale change.
(a)
Car4
(b)
CarScale
(c)
CarDarkIn the CarScale sequence, the tracking target is a fast motion car in an open road. Compared with the Car4 sequence, this sequence is more challenging because the tracked car undergoes large scale changes and fast motion on the entire sequence. Due to the fast motion accompanied with the tree’s occlusion, IVT, L1, MIL, and VTD trackers drift with different degree at frame 164 and finally lost the target at frame 171. The proposed method gives the best results followed by the OAB tracker.The CarDark sequence is challenging because the target object undergoes the motion blur in a night environment with low contrast and strong reflection interference. Due to the strong reflection interference, the MIL tracker drifts a little from the target at frame 122 and lost the target at frame 202 and then regards the other car as the tracked target. The IVT, L1, and VTD trackers drift away from the target at frame 277. The OAB tracker performs well on this sequence and yields the second best results. The proposed method can accurately track the target object in the whole sequence with small center position error and high overlap rate.Experiment 2 (occlusion, scale change, and rotation of nonrigid object).
The target object in FaceOcc2 sequence undergoes the drastic occlusion and in-plane rotation. As shown in Figure3(a), when there is a small occlusion with a book at frames 128~185 and frames 245~279, all methods perform well. But when partial occlusion and in-plane rotation occur together at frames 392~510, most of trackers have poor performances. When the target almost fully occludes by a book and a hat at frames 688~740, all methods except this paper method drift away from the target at different degrees. Since the proposed method uses multiple complementary features to build transient and stable appearance models and update the template library online, it can effectively handle the occlusion and give satisfactory tracking results.Figure 3
The tracking results of object undergoing occlusion, scale change, and rotation.
(a)
FaceOcc2
(b)
Freeman1
(c)
GirlThe Freeman1 sequence is challenging because the interested man’s face undergoes large scale changes and view variations. Due to the large scale changes, MIL drifts away from the target at frame 32. As the view changes from the left to right, the L1 and OAB trackers totally lost the target at frames 131 and 176, respectively. The tracking methods like IVT, VTD, and the proposed method perform well on this challenging sequence and can track the target accurately.The Girl sequence has drastic appearance changes because of the out of plane rotation and similar target occlusion. When out of plane rotation occurs at frames 90~122 and 169~260, IVT tracker totally fails to track the Girl’s face; other trackers drift at different degrees. The OAB and VTD trackers fail to track the target object and track the similar target when the Girl’s face is occluded by a man’s face at frames 420~470. MIL tracker can successfully track the target except some errors like frames 303 and 433. L1 tracker and the proposed method perform well on this sequence.Experiment 3 (illumination, scale change, and occlusion of nonrigid object).
The track target in shaking sequence undergoes drastic illumination and poses changes on the whole video sequence. It brings more challenges to accurately track the target because the color of object appearance is similar to the stage lighting. IVT and OAB almost fail to track the target at frame 23 and cannot recover at the rest frames. The MIL tracker drifts a little at frame 61 due to the drastic illumination changes. Although the stage lights change drastically accompanied with the serious head shaking, the L1 and VTD trackers perform well except some errors. The proposed method can effectively adapt to the severe object appearance changes when those variations occur together and achieve satisfactory results.The Woman sequence is very challenging because the target object undergoes large scale changes, view variations, and occlusions simultaneously. As shown in Figure4(b), all trackers merely perform well except the proposed method. The results of the proposed method show slight drift; other methods like IVT, L1, MIL, OAB, and VTD trackers totally lose the target when heavy occlusion occurs at frame 130 and never recover to track the target in the subsequent video sequence except the OAB tracker. The OAB tracker recaptures the target at frame 337 and keeps to track the target until the end of sequence with a little drift.Figure 4
The tracking results of object undergoing severe occlusion, illumination, and pose changes.
(a)
Shaking
(b)
Woman
(c)
JoggingThe Jogging sequence is more challenging to track because the tracked target is fully occluded by a stem and undergoes large scale change and fast motion simultaneously. The IVT, L1, MIL, OAB, and VTD trackers completely fail to track the target when the target is fully occluded by a stem at frames 50~62, and the OAB tracker recaptures the target at frame 106. The proposed method can accurately track the target on the entire sequence.From some sampled tracking results of the proposed method and other five methods on 9 image sequences, we can conclude that the algorithm in this paper can accurately and robustly track the target under the environment with illumination variation, scale change, and motion blur.


## 4.3. Quantitative Comparison

Two metrics are used to evaluate the proposed tracker with reference trackers in gray-scale videos. The first is the center position error, which is applied to evaluate the distance between the ground-truth centerR
g and tracked object center R
t (in pixels) at each frame by the Euclidean distance. The other metric is the overlap rate [28], which is defined as score
=
area
(
R
t
∩
R
g
)
/
a
r
e
a
(
R
t
∪
R
g
), where R
t denotes the bounding box generated by a tracking method and R
g is the ground-truth bounding box.Table2 and Table 3 show the average center position errors and the average overlap rate for all trackers. Figures 5 and 6 show the center position error curve and overlap rate evaluation curve of different trackers on 9 video sequences at each frame. It can be seen that the proposed algorithm has the optimal or suboptimal performance in terms of average center position errors and average overlap rate in most test video sequences compared with other methods. Most competing tracking methods do not give a satisfactory result; the center position error is larger and the overlap rate is lower. The average position error of this paper at 9 videos is only 5.53 pixels, which is far less than other trackers; the average central position errors of other trackers like IVT, L1, MIL, OAB, and VTD trackers are 59.58 pixels, 76.14 pixels, 51.91 pixels, 46.01 pixels, and 41.54 pixels, respectively. The average overlap rate of the proposed method reaches 73.6%, which is higher than other trackers; the average overlap rates of IVT, L1, MIL, OAB, and VTD trackers are 39.9%, 34%, 40%, 33.7%, and 32.7%, respectively. They both highlight the advantages of the algorithm in this paper. Overall, the effectiveness of the proposed tracker method is verified.Table 2
The average position errors (pixels).
Video
IVT
L1
MIL
OAB
VTD
Our method
Car4
16.31
209.22
55.10
103.49
34.46
2.97
CarScale
42.64
53.52
43.63
15.66
29.35
8.44
CarDark
22.03
19.85
48.41
3.90
23.29
1.92
FaceOcc2
16.21
15.54
18.17
24.32
10.69
8.42
Freeman1
6.85
61.74
18.98
25.44
10.71
7.44
Girl
29.37
3.76
14.11
8.55
11.70
3.63
Shaking
87.90
21.05
9.59
144.66
13.42
7.11
Woman
186.70
151.79
128.08
63.76
119.01
3.27
Jogging
128.29
148.81
131.14
24.29
121.31
6.54
Average
59.58
76.14
51.91
46.01
41.54
5.53
Note: the optimal result is shown as bold and the suboptimal one as italic.Table 3
The average overlap rate for all of trackers.
Video
IVT
L1
MIL
OAB
VTD
Our method
Car4
0.570
0.196
0.257
0.151
0.357
0.873
CarScale
0.469
0.515
0.421
0.451
0.408
0.768
CarDark
0.480
0.513
0.153
0.741
0.461
0.868
FaceOcc2
0.572
0.641
0.611
0.549
0.70
0.723
Freeman1
0.493
0.234
0.241
0.288
0.343
0.587
Girl
0.147
0.676
0.404
0.577
0.521
0.710
Shaking
0.032
0.50
0.663
0.017
0.613
0.735
Woman
0.129
0.145
0.147
0.217
0.144
0.654
Jogging
0.139
0.142
0.127
0.60
0.134
0.711
Average
0.337
0.396
0.336
0.399
0.409
0.736
Note: the optimal result is shown as bold and the suboptimal one as italic.Figure 5
The center position errors of different trackers on 9 video sequences.Figure 6
The overlap rate evaluation of different trackers on 9 video sequences in each frame.


## 5. Conclusion

In this paper, we propose a robust tracking algorithm that leverages hierarchical sparse coding to optimize the image representation from multifeature. We compare our tracking method with other five state-of-the-art trackers on nine sequences to validate the robustness and accurateness of the proposed method. The experiment results show that our method can effectively and robustly handle the challenging scenes where the target object undergoes drastic variation in pose, scale, rotation, occlusion, and illumination. The success of our method can be attributed to constructing multiple observation models that form the multifeature by hierarchical sparse coding, which takes the spatial relationship of neighborhood features into consideration and solves the sparse approximation problem byl
_
RLS. The appearance model constructed by instantaneous and stable appearance features with two-stage sparse representation coding is more robust to cope with appearance change in complex environment and more effective to select a set of discriminative features to separate the target from its background. In the proposed method, we compute the reliability of each tracker by the tracker likelihood function that accounts for transient and reconstructed appearance model and select the most reliable one among multiple trackers. The training set and the template library are both incrementally online updated. All of these are beneficial to cope with the appearance change and can improve the tracking performance in dynamic environments.However, the limitation of the proposed method is mainly focused on the following.(
1
) The tracking system is not effective enough for real-time tracking because multifeatures are calculated at the same time for test video sequences, which is time-consuming. In addition, it cannot be adapted to extract the feature according to the video attribute. (
2
) The ability of each feature to describe the target cannot be effectively measured. (
3
) It cannot successfully track the target when the object leaves out the scene but reappears in subsequent frames.In the future, we will improve the proposed method in some aspects.(
1
) We will improve the algorithm in real-time by proposing a method to adaptively extract the multifeatures according to the video attribute, which can reduce the time and computation load of the feature extraction. (
2
) We would improve the tracking performance by introducing the occlusion mechanism and drift mechanism, which can alleviate updating the template with wrong samples when the target object is occluded or drifted. Both strategies are useful to deal with appearance changes and beneficial to robustly track the target in complex environments.



---

*Source: 5894639-2016-08-18.xml*