object of interest from any direction, it is obvious
to expect its movement to follow the shortest path.
This typical scenario is illustrated in Fig. 1(a). The
intermediate space is unconstrained and the target can
move freely while heading towards the static object
of interest. A block based geometry, as depicted in
Fig. 1(b), is used to illustrate the different motion
configurations possible within such an environment.
That is, let us consider that the environment is
partitioned into a rectangular grid of a predefined
size. Therefore, as a target moves towards the static
object of interest, the path is through the intermediate
blocks (C
1
,C
2
,C
3
) considering its initial position
is in one of the outer layer blocks, e.g. B
1
, ..., B
5
.
There can be different ways a target can reach to a
particular inner block from a set of outer blocks. For
example, if the initial location of the target is in (B
1
),
and if the shortest path assumption holds, then the
target must go through (C
1
) as depicted in Fig. 1(c).
Other possibilities are described in the successive
images. Under these conditions, a theoretical model
can be used to describe the movement. Suppose, the
probability of a target being in one of the outermost
blocks is represented as P. It can easily be verified
from the images shown in Figs. 1(c-f) that, a target
can reach to one inner block from an outer block
in three possible ways. Therefore, probability of
reaching any of the inner block becomes three times
the probability of an outer block. Since the static
object of interest is at the center and surrounded by
eight equip-probable blocks, the possible ways of
reaching the target is eight times higher than of its
immediate neighboring blocks. Figs. 1(g-h) illustrate
the method for a three layer geometry. However, this
can easily be extended for any desired number of
layers.
Figure 2: (a) Division of a space assuming the object of in-
terest is located at the center and corresponding frequency
of a block being used while a target approaches the object
of interest. (b) pdf of the importance of a block in terms
of number of times a block is accessed while a target ap-
proaches the object of interest located at center.
The Hypothesis. Suppose, x is a random variable that
denotes the probability of a target visiting a particular
block while approaching an object of interest. If we
plot the different number of ways an inner block
can be reached from an outer block, the normalized
probability computed shall represent a distribution of
the importance of the blocks. Now, we hypothesize
that, given a scenario where sufficient number of
targets approach towards an object of interest, if the
surveillance space is divided into rectangular blocks
as shown in Fig. 2(a), then the target motion model
will usually follow the theoretical distribution shown
in Fig. 2(b).
Proposed Trajectory Analysis. The proposed
trajectory analysis method has been designed with
the knowledge of the theoretical model described
earlier. We have used the target detection and
tracking algorithm proposed in (Dinh and Medioni,
2011) to extract the trajectories of moving targets.
A spatial domain heuristic, where-in a point on
the trajectory is removed if it deviates abruptly
from its usual path, has been applied to remove
noise/outliers from the trajectories. Further, the
importance of a block is estimated from the cleaned
trajectories using the steps detailed below. During
pre-processing, the surveillance scene is divided
into a rectangular grid of uniform dimension as
shown in Fig. 1(b) where the total number of blocks
is denoted by M. Assuming that N trajectories of
multiple targets are available for analysis; it can be
represented as a set, say T = {t
1
,t
2
...,t
N
} such that
t
i
= (< x
1
, y
1
>, < x
2
, y
2
>, ......, < x
m
i
, y
m
i
>
) is a
trajectory of length m
i
. The importance of a block is
estimated as:
• Step I. The average velocity (V
O
i
avg
) of a target
(O
i
) is estimated from the uniformly sampled seg-
ments of its trajectory. First, the minimum and
maximum values of the velocity of the targets are
calculated using (1) and (2) such that p
j
and p
j+1
denote successive points on the trajectory t
i
that is
bounded by 0 < j <
|
t
i
|
.
V
O
i
min
= min |p
j
− p
j+1
| (1)
V
O
i
max
= max |p
j
− p
j+1
| (2)
Next, range of the velocity [V
O
i
max
− V
O
i
min
] is di-
vided into R uniform segments and a histogram of
the instantaneous velocity is generated. Finally,
the mean of all instantaneous velocities under the
largest bin is taken as the average velocity of the
target. This will remove any bias that may occur
due to fast moving segments within a given trajec-
tory.
• Step II. Next, the total number of times a block is
visited by various targets is computed. We call it
global visit (G
M
k
) where M
k
is the block. Initially,
global visit parameter for all the blocks are set to
zero. However, when a target enters into a new
block, we increment its global visit value by one.
VISAPP2015-InternationalConferenceonComputerVisionTheoryandApplications
480