Next Article in Journal
Fusion of Smartphone Motion Sensors for Physical Activity Recognition
Next Article in Special Issue
Robust Pedestrian Tracking and Recognition from FLIR Video: A Unified Approach via Sparse Coding
Previous Article in Journal
Affine Transform to Reform Pixel Coordinates of EOG Signals for Controlling Robot Manipulators Using Gaze Motions
Previous Article in Special Issue
Sparse Representation for Infrared Dim Target Detection via a Discriminative Over-Complete Dictionary Learned Online
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Joint Target Tracking, Recognition and Segmentation for Infrared Imagery Using a Shape Manifold-Based Level Set

1
School of Mechatronical Engineering, Beijing Institute of Technology, No. 5, Zhongguancun South Street, Haidian District, Beijing 100081, China
2
School of Electrical and Computer Engineering, Oklahoma State University, 202 Engineering South, Stillwater, OK 74078, USA
3
School of Electrical and Computer Engineering, University of Oklahoma, 110 West Boyd, DEH 150Norman, OK 73019, USA
*
Author to whom correspondence should be addressed.
Sensors 2014, 14(6), 10124-10145; https://0-doi-org.brum.beds.ac.uk/10.3390/s140610124
Submission received: 24 March 2014 / Revised: 22 May 2014 / Accepted: 23 May 2014 / Published: 10 June 2014

Abstract

: We propose a new integrated target tracking, recognition and segmentation algorithm, called ATR-Seg, for infrared imagery. ATR-Seg is formulated in a probabilistic shape-aware level set framework that incorporates a joint view-identity manifold (JVIM) for target shape modeling. As a shape generative model, JVIM features a unified manifold structure in the latent space that is embedded with one view-independent identity manifold and infinite identity-dependent view manifolds. In the ATR-Seg algorithm, the ATR problem formulated as a sequential level-set optimization process over the latent space of JVIM, so that tracking and recognition can be jointly optimized via implicit shape matching where target segmentation is achieved as a by-product without any pre-processing or feature extraction. Experimental results on the recently released SENSIAC ATR database demonstrate the advantages and effectiveness of ATR-Seg over two recent ATR algorithms that involve explicit shape matching.

1. Introduction

As a challenging problem in pattern recognition and machine learning for decades, automatic target tracking and recognition (ATR) has been an important topic for many military and civilian applications. Infrared (IR) ATR is a more challenging problem due to two main reasons. First, an IR target's appearance may change dramatically under different working conditions and ambient environment. Second, the IR imagery usually has poor quality compared with the visible one. There are two important and related research issues in ATR research, appearance representation and motion modeling [1]. The former one focuses on capturing distinct and salient features (e.g., edge, shape, texture) of a target, and the latter one tries to predict the target's state (e.g., position, pose, velocity) during sequential estimation. They could play a complementary role in an ATR process [2].

Shape is a simple yet robust, feature for target representation in many ATR applications. There are three commonly used ways of shape representation: a 3D mesh model [3], 2D shape templates [4,5] and a manifold-based shape generative model learned from 2D snapshots [68]. When a 3D model was used, a 3D-to-2D projection is needed to get the 2D shapes according to the camera model and the target's position. Using a 3D model for shape modeling usually needs more memory and expensive computational resources. In [5], a 2D shape template was used to represent the target's appearance, and an online learning was used to update this shape model under different views. Manifold learning methods have proven to be powerful for shape modeling by providing a variety of meaningful shape prior to assist or constrain the shape matching process. In [8], a couplet of view and identity manifolds (CVIM) was proposed for multi-view and multi-target shape modeling, where target pre-segmentation was implemented via background subtraction and the ATR inference involves explicit shape matching between segmented targets and shapes hypothesis generated by CVIM.

In this work, we propose a new particle filter-based ATR-Seg (segmentation) algorithm that integrates JVIM (joint view-identity manifold) with a shape-aware level set energy function which leads to a joint tracking, recognition and segmentation framework. JVIM encapsulates two shape variables, identity and view, in a unified latent space, which is embedded with one view-independent identity manifold and infinite identity-dependent view manifolds. Unlike CVIM obtained via nonlinear tensor decomposition, JVIM is learned via a modified Gaussian process latent variable model [9] which leads to a probabilistic shape model. Also, a stochastic gradient descent method [10] is developed to speed up JVIM learning, and a local approximate method is used for fast shape interpolation and efficient shape inference. Furthermore, we integrate JVIM with a level set energy function that is able to evaluate how likely a shape synthesized by JVIM can segment out a valid target from an image. This energy function is adopted as the likelihood function in the particle filter where a general motion model is used for handling highly maneuverable targets. The performance of ATR-Seg was evaluated using the SENSIAC (Military Sensing Information Analysis Center) IR dataset [11], which demonstrated the advantage of the proposed method over several methods that involve target pre-segmentation and explicit shape matching.

The remainder of this paper is organized as follow. In Section 2, we review some related works on shape manifold learning and shape matching. In Section 3, we use a graphical model to develop a probabilistic framework of our ATR-Seg algorithm. In Section 4, we introduce JVIM for general shape modeling. In Section 5, we present a shape-aware level set energy function for implicit shape matching.

In Section 6, we present a particle filter-based sequential inference method for ATR-Seg. In Section 7, we evaluate the proposed ATR-Seg algorithm in two aspects, i.e., JVIM-based shape modeling and implicit shape matching which are involved in the likelihood function of the particle filter. We conclude our paper in Section 8.

2. Related Works

ATR itself is a broad field involving diverse topics. Due to the fact that shape modeling is the key issue in our ATR research, our review below will be focused on two shape-related topics, manifold-based shape modeling and shape matching.

2.1. Manifold-Based Shape Modeling

A manifold-based shape model can be learned from a set of exemplar shapes and is able to interpolate new shapes from the low-dimensional latent space. Roughly speaking, there are three manifold learning approaches for shape modeling, geometrically-inspired methods, latent variable models, and hybrid models. The first approach seeks to preserve the geometric relationships among the high-dimensional data in the low-dimensional space, e.g., IsoMap [12], Local Linear Embedding (LLE) [13], Diffusion Maps [14] and Laplacian Eigenmaps [15]. These methods focus on how to explore the geometric structure among the high-dimensional data and how to maintain this structure in the low dimensional embedding space. However, the mapping relationship from the latent space and the data space is not available and has to be learned separately. The second approach represents the shape data by a few latent variables along with a mapping from the latent space to the data space, such as PCA [16], PPCA [17], KPCA [18], Gaussian Process Latent Variable Models (GPLVM) [19] and tensor decomposition [20], et al. GPLVM is a probabilistic manifold learning method which employs the Gaussian process as the nonlinear mapping function. Above approaches are data driven shape modeling methods without involving prior knowledge in the latent space, and as a result, the shape-based inference process may be less intuitive due to the lack of a physically meaningful manifold structure.

To support a more meaningful and manageable manifold structure while preserving the mapping function, there is a trend to combine the first two approaches along with some topology prior for manifold learning [21]. In [9], the local linear GPLVM (LL-GPLVM) was proposed for complex motion modeling, which incorporates a LLE-based topology prior in the latent space. Specifically, a circular-shaped manifold prior is used to jointly model both “walking” and “running” motion data in a unified cylinder-shaped manifold. In [8], CVIM was proposed for shape modeling via nonlinear tensor decomposition where two independent manifolds, an identity manifold and a view manifold, were involved. Specifically, the view manifold was assumed to be a hemisphere that represents all possible viewing angles for a ground target, and the identity manifold was learned from the tensor coefficient space that was used to interpolate “intermediate” or “unknown” target types from known ones. A key issue about the identity manifold is the determination of manifold topology, i.e., the ordering relationship across all different target types. Sharing a similar spirit of IsoMap, the shortest-closed-path is used to find the optimal manifold topology that allows targets with similar shapes to stay closer and those with dissimilar shapes far away. This arrangement ensures the best local smoothness and global continuity that are important for valid shape interpolation along the identity manifold.

2.2. Shape Matching

In shape-based tracking algorithms, there are two ways to measure shape similarity: explicit shape matching and implicit shape matching. The former one involves a direct spatial comparison between two shapes, an observed one and a hypothesized one, by using a certain distance metric. In such a case, pre-processing or feature extraction, e.g., background subtraction in [8], is needed prior to tracking and recognition, which is relatively manageable for a stationary sensor platform and may need additional computational load in a general case. Moreover, the overall ATR performance could be sensitive to the pre-processing results. The latter one represents a shape implicitly by a level set embedding function which can be used to evaluate the segmentation quality of a given shape in an image. For example, a shape-constrained energy function was used in [6,7] to evaluate how likely the given shape can segment out a valid object, where a gradient descent method was used to optimize this energy function to achieve tracking and segmentation jointly. Therefore, implicit shape matching does not involve any pre-processing or feature extraction beforehand, however, due to the lack of dynamic modeling in level set optimization, it is still hard to track highly maneuverable targets by the traditional data-driven gradient descent optimization method. As pointed in [22], motion/dynamic modeling is an important step for most ATR applications. This motivates our research to augment a motion model in implicit shape matching for maneuverable target tracking.

3. ATR-Seg Problem Formulation

We list all symbols used in this paper in Table 1. Given the observed video sequence It, with t = 1,…, T , where T is the total number of image frames, the objective of ATR-Seg is to (1) find the 3D position of a target in the camera coordinate p (tracking) or 2D image coordinate, (2) to identify the target type α (recognition), along with the view angle φ (pose estimation), and (3) to segment the target-of-interest that best explains the observation data Φ (segmentation). The 2D shape of a target can be determined by the target type α, and view angle φ, so we define Θ = [α, φ] to represent two shape related variables. The conditional dependency among all variables is shown in Figure 1.

According to Figure 1, we define the objective function of ATR-Seg from the joint distribution p(pt, Θt, Φ, It) which can be written (t is omitted for simplicity) as:

p ( p , Θ , Φ , I ) = p ( I | p , Φ ) p ( Φ | Θ ) p ( Θ ) p ( p )

By using the Bayesian theorem, we can get the posterior as:

p ( p , Θ , Φ , I ) = p ( I | p , Φ ) p ( Φ | Θ ) p ( Θ ) p ( p )
which encapsulates three major components in the proposed ATR-Seg algorithm, as shown below:
  • Shape manifold learning provides a mapping from Θ to Φ, i.e., p(Φ|Θ). In Section 4, JVIM is proposed for multi-view and multi-target shape modeling, which features a novel manifold structure with one view-independent identity manifold and infinite identity-dependent view manifolds to impose a conditional dependency between the two shape-related factors, view and identity, in a unified latent space.

  • Shape-aware level set p(I|p, Φ) measures how likely Φ can segment a valid target at position p in image I. In Section 5, a shape-aware level set energy function is proposed for implicit shape matching, which evaluates the segmentation quality.

  • Pose/position priors Θ and p, i.e., p(Θ) and p(p), i.e., are important to track highly manoeuverable targets in a sequential manner. In Section 6, sequential shape inference method is presented that involve dynamic priors for Θ and p using a 3D motion model.

The flowchart for ATR-Seg is shown in Figure 2 where four steps are involved sequentially and recursively. First, state prediction will draw a set of samples to predict all state variables (position/angle/identity). Second, a series of shape hypotheses are created via JVIM in some hypothesized locations according to predicted state information. Third, a level-set energy function is used as the likelihood function to weight each hypothesized shape/location that quantifies how well that shape can segment a valid target in that particular location. Fourth, state estimation at the current frame is obtained by the conditional mean of all weighted samples and will be used for state prediction in the next frame.

4. Joint View-Identity Manifold (JVIM)

JVIM is learned from a set of 2D shape exemplars Y generated from a set of 3D CAD models. The latent space X can be represented by two variables, identity α and view φ (including aspect angle θ and elevation angle ϕ), which are defined along their respective manifolds. Considering the fact that all targets have different 3D structures, leading to different view manifolds, and they keep the same identity under different views, we impose a conditional dependency between α and φ in JVIM that encapsulates one view-independent identity manifold and infinite identity-dependent view manifolds. Specifically, the identity manifold represents the view-independent shape variability across different target types, and an identity-specific view manifold captures the shape variability of a target under different views. Motivated by [8,23], the identity manifold is simplified to have a circular-shaped topology prior, which facilitates manifold learning and shape inference. Intuitively, a hemispherical-shaped topology prior is assumed for identity-specific view manifold, which represents all possible aspect and elevation angles for ground vehicle. All topology priors are encoded by LLE and incorporated into the GPLVM-based learning framework, as shown in Figure 3.

The objective of JVIM learning is to find X and β by maximizing p(Y|X, β, w), where β is the mapping parameter and w represents the LLE-based topology prior in the latent space. The Gaussian process (GP) is used as the nonlinear mapping function from the latent space to the shape space (XY), and the objective function of JVIM learning is written as:

p ( Y | X , β , w ) = p ( Y | X , β ) p ( X | w )
where:
p ( Y | X , β ) 1 ( 2 π ) N d | K Y | d exp ( 1 2 t r ( K Y 1 YY T ) )
where d is the dimension of the shape space and β denotes the kernel hyper-parameters used in the covariance matrix, KY. It is worth noting that Equation (4) is similar to the objective function of GPLVM [19], and:
p ( X | w ) = 1 Z exp ( 1 σ 2 i = 1 N X i j = 1 N w i j X j 2 )
where w is the set of LLE weights to reconstruct each latent point from its local neighboring points by minimizing f ( w ) = i = 1 N X i j = 1 N w i j X j 2, Z is a normalization constant, σ2 represents a global scaling of the prior and N the number of training samples. Furthermore, the negative log operation is used to simplify the objective function as:
L JVIM = log p ( Y | X , β ) p ( X | w ) = L D + L T + C
where C is a constant.

JVIM learning involves a gradient descent method to minimize the objective function defined in Equation (3) with respect to X and β. With an O(N3) operation required at each iteration, it is computationally prohibitive for a large training data set. The stochastic gradient descent proposed in [10] is adapted to be a local updating according to the unique structure of JVIM to approximate the gradients locally. At each iteration, the reference point, xr, is chosen randomly, and the derivatives w.r.t XR and β are calculated as:

L D X R ( K R 1 Y R Y R T K R 1 d K R 1 ) K R X R
L JVIM β ( K R 1 Y R Y R T K R 1 d K R 1 ) K R β
where XR is the neighborhood for a reference point, xr, of size M1, YR is the corresponding shape data and KR (M1 × M1) is the kernel matrix of XR. The neighborhood for each training data can be pre-assigned according to the topology structure, and the gradients are estimated stochastically, locally and efficiently.

As a generative model, given an arbitrary latent point in X, JVIM can generate the corresponding shape via Gaussian Process (GP) mapping. For real-time applications, shape interpolation must be carried out efficiently, which is difficult for a large training data set with high dimensionality. Inspired by [25], a GP can be approximated by a set of local GPs, in JVIM-based shape interpolation, the kernel matrix is computed locally from a set of training data that are close to the given point. Given x′, we first find its closest training point, which has a pre-assigned neighborhood, X′, of size M2; then, X′ and the corresponding shape data Y′ are used to approximate the mean and variance of GP mapping as:

μ ^ x = k x X T K Y 1 Y
σ ^ x 2 = k ( x , x ) k x X T K Y 1 k x X
where kx′X′ is a vector made of k(x′, xi) (xiX′) and KY′ (M2 × M2) is the local covariance matrix computed from X′. More detail about JVIM learning and inference can be found in our previous work [26], where explicit shape matching is involved. In the following section, we will introduce implicit shape matching by incorporating a shape-aware level set for target tracking and recognition, where target segmentation becomes a by-product.

5. Shape-Aware Level Set

JVIM is used to provide a useful shape prior that can be further combined with the level set to define an energy function for implicit shape matching. This is called the shape-aware level set, which does not involve feature extraction or target pre-segmentation. The shape-aware level set in this work is distinct from that in [6,7,27] primarily in two aspects. Firstly, the shape generated model in [6,7], which was less structured with little semantic meaning and, was limited to object recognition/segmentation under the same view or human pose estimation for the same person along the same walking path. JVIM is a multi-view and multi-target shape model that has a well-defined semantic structure, which supports robust ATR for different targets under arbitrary view angles. Secondly, a gradient decent method was used for level set optimization in [6,7,27], which does not involve a motion model and makes it hard to track highly maneuverable targets. In this work, a 3D motion model is used to combine the position/pose priors into a sequential optimization model to improve the robustness and accuracy of ATR-Seg.

As shown in Figure 4a, we represent an image by I = {xi, yi}, where 1 ≤ in, n is the number of pixels in I and x and y are the pixel 2D location and pixel intensity value, respectively. We introduce a parameter, M, to represent the foreground/background models M = {Mf, Mb}; then, the original graphical model of ATR-Seg in Figure 1 will become the one in Figure 4c. which defines a joint distribution of all parameters for each pixel (xi, yi) as

p ( x i , y i , p , Θ , Φ , M ) = p ( x i | p , Φ , M ) p ( y i | M ) p ( Φ | Θ ) p ( M ) p ( Θ ) p ( p )
where Φ is a shape represented by the level set embedding function shown in Figure 4b and p(Φ|Θ) corresponds to JVIM-based shape interpolation via GP mapping. A histogram is used for foreground/background appearance model p(yi|M), and the number of bins is dependent on the size of the target and gray scale. In order to get the posterior, p(p, Θ, Φ, M|xi, yi), which will be used to develop the objective function for ATR-Seg, we take the same strategy as in [27]. First, divide Equation (11) by p ( y i ) = j { f , b } p ( y i | M j ) p ( M j ):
p ( x i , p , Θ , Φ , M | y i ) = p ( x i | p , Φ , M ) p ( M | y i ) p ( Φ | Θ ) p ( Θ ) p ( p )
where p(M|yi) is given by:
p ( M j | y i ) = p ( y i | M j ) p ( M j ) k { f , b } p ( y i | M k ) p ( M k ) , j { f , b }

Upon dividing Equation (12) by p(xi) = 1/n and marginalizing over the models, M, we obtain:

p ( p , Θ , Φ | x i , y i ) = n j { f , b } p ( x i | p , Φ , M j ) p ( M j | y i ) p ( Φ | Θ ) p ( Θ ) p ( p )

Assuming all pixels are independent, the posterior for all pixels in a frame is then given by:

p ( x i | p , Φ , M f ) = H ϵ [ Φ ( x i ) ] / η f p ( x i | p , Φ , M b ) = { 1 H ϵ [ Φ ( x i ) ] } / η b
where Hϵ [·] is the smoothed Heaviside step function, η f = i = 1 n H ϵ [ Φ ( x i ) ], η b = i = 1 n { 1 H ϵ [ Φ ( x i ) ] } and p(Mj) = ηj/n for j ∈ {f, b}.

Then, from Equations (2), (14) and (15), we have:

p ( I | p , Φ ) i = 1 n j { f , b } p ( x i | p , Φ , M j ) p ( M j | y i )
which evaluates how likely shape Φ can segment a valid target from I at position p. The objective function in Equation (2) can be optimized through a gradient descent method similar to the one in [7], which is illustrated in Figure 5. As shown in Figure 5, JVIM is firstly used to generate a shape hypothesis, Φ0, given initial identity and view angle Θ0; then, Φ0 is used to initialize the objective function, p(p, Θ, I), for initial position, p0. We take the derivative of p(p, Θ, Φ|I) with respect to Θ and p to get p ( p , Θ , Φ | I ) Θ and p ( p , Θ , Φ | I ) p, which will be used to update Θ and p until the objective function converges. When p(p, Θ, Φ|I) is maximized, we output the updated target's 2D position, p*, target identity and view angle Θ*, as well as the updated shape Φ* that can best segment the image.

This method works well on a single image when a good initialization is given in the latent space of JVIM. However, it may fail quickly when dealing with an image sequence with a highly maneuverable target, due to four possible cases of shape ambiguity, as shown in Figure 6, which makes data-driven optimization not reliable in practice. (1) The first is due to the possible shape mismatch between the CAD models and real targets, even for the same target type (Figure 6a). (2) The second is due to the symmetry property of a target's shape (Figure 6b), which means a target may present a similar shape at different (e.g., supplement) aspect angles, especially when the elevation angle is zero (Figure 6b). (3) The third is due to the ambiguity of the front/rear views when a target looks very similar (Figure 6c). (4) The fourth is similar to the previous one in which many targets look alike at the front/rear views (Figure 6d). These factors make the gradient-based approach not effective at dealing with a maneuvering target. A possible remedy is to introduce a dynamic motion model to support robust sequential shape inference based on JVIM, as to be discussed below.

6. Sequential Shape Inference

Essentially, the objective of ATR-Seg is to perform sequential shape inference from an image sequence by maximizing the posterior of p(pt, Θt, Φt|It). According to Figure 1 in Section 3, Φ is only dependent on Θ, so the objective function can be rewritten as:

p ( p t , Θ t , Φ t | I t ) = p ( p t , Θ t | I t ) p ( Φ t | Θ t )
where p(Φt|Θt) is JVIM-based shape interpolation via GP mapping. Since p(Φt|Θt) is not related to the observation, so the main computational load is the maximization of p(pt, Θt|It). For sequential ATR-Seg, the optimization of p(pt, Θt|It) has two stages: prediction and update. In the first stage (prediction), we use a motion model to predict p(pt, Θt|It−1) from the previous result p(pt−1, Θt−1|It−1) as:
p ( p t , Θ t | I t 1 ) = p ( p t 1 , Θ t 1 | I t 1 ) p ( p t | p t 1 ) p ( Θ t | Θ t 1 ) d Θ t 1 d p t 1
where p(pt|pt−1) and p(Θt|Θt−1) are used to predict the position and identity/view of a moving target. They are related a motion model that characterizes the target's dynamics and kinematics. In the second stage (update stage), we use the Bayes' rule to compute the posterior as:
p ( p t , Θ t , Φ t | I t ) = p ( p t , Θ t , Φ t | I t 1 ) p ( I t | p t , Φ t ) p ( I t | I t 1 )
where p(pt, Θt, Φt|It−1) = p(pt,Θt|It−1)p(Φt|Θt) and we have p(It|pt, Φt, Θt) = p(It|pt, Φt). Hence, the objective function of the sequential ATR-Seg algorithm can be further rewritten as:
p ( p t , Θ t , Φ t | I t ) p ( I t | p t , Φ t ) p ( Φ t | Θ t ) p ( Θ t | Θ t 1 ) p ( p t 1 , Θ t 1 , Φ t 1 | I t 1 ) p ( p t | p t 1 ) d Θ t 1 d p t 1

Due to the nonlinear nature of Equation (20), we resort to a particle filter-based inference framework [29] for sequential optimization, as represented by the graphic model in Figure 7 (left). Thanks to the compact and continuous nature of JVIM, we can draw samples from its latent space for efficient shape interpolation. In the inference process, the state vector is defined as Z t = [ p t T , υ t , ψ t , α t ] T where p t = [ p x t , p y t , p z t ] T represents the target's 3D position, with the xyz axes denoting the horizon (x), elevation (y) and range (z) directions, respectively (as shown in Figure 7 (right)); υt is the velocity along the heading direction, ψt. A 3D-2D camera projection, W (p), is needed to project a 3D position to a 2D position in an image that is assumed to be unchanging for a stationary sensor platform. It is worth noting that we can compute θt (the aspect angle) from ψt (the heading direction) or vice versa. As a matter of fact, the two angles are similar for distant targets when the angle between the lineof sight and the optical axis along the range direction (z) is very small. Because the target is a ground vehicle and to keep it general, a white noise acceleration model is used to represent the dynamics of Zt, where a simple random walk is applied on the heading direction, ψt, to represent arbitrary maneuvering. Moreover, we define the dynamics of αt (target identity) to be a simple random walk along the identity manifold by which the estimated identity value normally quickly converges to the correct one.

{ ψ t = ψ t 1 + ζ t ψ , υ t = υ t 1 + ζ t υ , p x t = p x t 1 + υ t 1 sin ( ψ t 1 ) Δ t + ζ t x , p y t = p y t 1 + ζ t y , p z t = p z t 1 + υ t 1 cos ( ψ t 1 ) Δ t + ζ t z , α t = α t 1 + ζ t α ,
where Δt is the time interval between two adjacent frames. The process noises associated with the target kinematics, ζ t ψ, ζ t υ, ζ t x, ζ t y, ζ t z, and ζ t α, are usually assumed to be a zero-mean Gaussian.

In a particle filter-based inference algorithm, samples were first drawn according to the dynamics of the state vector and the previous state value, and then, the implicit shape matching defined in Equation (16) was performed to assign a weight for each particle. The mean estimation of weighted samples produces the solution in the present frame. The pseudo-code for the ATR-Seg algorithm is given in Table 2. Thanks to the unique structure of JVIM, we can capture the continuous and smooth shape evolution during target tracking and recognition, where segmentation Φt is also archived as a by-product via the shape-aware level set. We expect that the proposed ATR-Seg algorithm has some advantages over other methods that require pre-processing or feature extraction prior to ATR inference [8,26].

7. Experimental Results

This experimental section provides a detailed evaluation of the ATR-Seg algorithm in six parts. First, we briefly talk about training data collection for JVIM learning with some qualitative results of shape interpolation. Second, we introduce the infrared ATR database used in this work and how different shape models are to be evaluated collectively and fairly. Third, we present the results of the particle filter-based infrared ATR algorithm, where four shape models (JVIM, CVIM, LL-GPLVM, nearest neighbor (NN)) are compared in the case of explicit shape matching. Fourth, we discuss the results of the proposed ATR-Seg algorithm, which involves JVIM-based implicit shape matching and is compared with the algorithms using explicit shape matching (with JVIM and CVIM). Fifth, we discuss the target segmentation results, which are the by-product of the ATR-Seg algorithm. We will also discuss some limitation of ATR-Seg along with some failed cases.

7.1. Training Data Collection

In our work, we considered six target classes as [8], i.e., SUVs, mini-vans, cars, pick-ups, tanks and armored personnel carriers (APCs), each of which has six sub-classes, resulting in a total of 36 targets, as shown in Figure 8. These 36 targets were ordered along the view-independent identity manifold according to a unique topology optimized by the class-constrained shortest-closed-path method proposed in [8] (before training). We considered aspect and elevation angles in the ranges 0 ≤ θ < 2π and 0 ≤ ϕ < π/4, which are digitized in the interval of π/15 and π/18 rad, respectively. A total of 150 training viewpoints were used for each target; all training data are generated by their 3D CAD models. In order to reduce the data dimension, the DCT-based shape descriptor proposed in [7] was used here to represent all training shapes for manifold learning. We first detect the contour of a 2D shape (120 × 80) and then apply the signed distance transform to the contour image, followed by the 2D DCT. Only about 10% DCT coefficients are used to represent a shape, which are sufficient for nearly lossless shape reconstruction. Another advantage of this shape descriptor is that we can do zero-padding prior to inverse DCT to accommodate an arbitrary scaling factor without additional zooming or shrinking operations.

JVIM-based shape interpolation is demonstrated in Figure 9, which manifests its capability of handling a variety of target shapes with respect to viewpoint changes for a known target, as well as the generalization to previously unseen target types. In Figure 9a, we pick one target type from each of the six classes. For each target type, we can obtain an identity-specific view manifold from JVIM along which we can interpolate new target shapes of intermediate views (in black) between two training view-points. A smooth shape transition is observed across all interpolated shapes, despite the strong nonlinearity of training shapes. Figure 9b shows the shape interpolation results (in black) along the view-independent identity manifold for the same side view. Although the interpolated shapes are not as smooth as previous ones, most of them are still meaningful, with a mixed nature of two adjacent training target types along the identity manifold. Compared to CVIM in [8], which assumes that the identity and view manifolds are independent, JVIM shows better shape interpolation results by imposing a conditional dependency between the two manifolds and is also more computationally efficient due to local inference. A detailed comparison can be found in [26], where JVIM is found to be advantageous over CVIM and several GPLVM-based shape models, both qualitatively and quantitatively.

7.2. Infrared ATR Database and Shape Models

We have obtained a set of mid-wave IR sequences from the SENSIAC ATR database [11], which includes IR imagery of civilian and military ground vehicles maneuvering around a closed-circular path at ranges from 1–3 km. Forty sequences from eight different target types at ranges of 1.0 km, 1.5 km, 2.0 km, 2.5 km and 3 km were selected for this work. For each sequence, tracking was performed on 500 frames. Background subtraction [30] was applied to each frame for clutter rejection, which is needed for two competing algorithms involving explicit shape matching. For each tracking method, the particle filter was initialized with the ground-truth in the first frame. Similar to [8], the process noise of the heading direction ζ t ψ is assumed to be a non-zero mean Gaussian to accommodate the circular moving trajectory which is necessary due to the ill-posed nature of image-based 3D tracking. This assumption can be relaxed if 3D pose estimation is not needed. Using the metadata provided with the database and a calibrated camera model, we computed the 3D ground-truth of position and aspect angle (in the sensor-centered coordinate system) for each frame. We refer the readers to [26] for more details about the ATR database.

In the following infrared ATR evaluation, we compare JVIM with LL-GPLVM [9] and CVIM [8], as well as the traditional nearest neighbor shape interpolation (NN). Both JVIM and CVIM treat shape factors (view and identity) continuously. To make a fair comparison, we learned a set of target-specific view manifolds by using LL-GPLVM, which involves a hemisphere as the topology constraint for manifold-based shape modeling. Then, we augment a “virtual” circular-shaped identity manifold (similar to that in JVIM and CVIM) for LL-GPLVM, where a NN method is used to “interpolate” arbitrary target types via training ones. Likewise, two “virtual manifolds” are introduced for the NN-based shape model, where we use the nearest neighbor to find the best matched training shapes. Thus, the two shape variables for four shape models can be inferred in a similar continuous way during ATR inference.

7.3. ATR with Explicit Shape Matching

We adopted the particle filter-based ATR algorithm used in [8], where JVIM, CVIM, LL-GPLVM and NN are evaluated in the case of explicit shape matching. In the CVIM-based ATR algorithm, two independent dynamical models are used. In JVIM-based tracking, the dynamic model is a two-stage one, where the first stage is along the view-independent identity manifold, while the second stage along the identity-dependent view manifold. For the LL-GPLVM-based ATR algorithm, one dynamic model is defined on each target-specific view manifold and one on the virtual identity manifold, where NN is used for identity interpolation. For the NN-based ATR algorithm, we employ two dynamic models on two virtual manifolds, like those in CVIM, where shape interpolation is done via NN (i.e., just using the training shapes).

The ATR performance of four shape models was evaluated with respect to three figures of merit: (1) px (horizontal) position error (in meters); (2) pz (slant range) position error (in meters); and (3) heading direction error ψ (in rads). Quantitative tracking performance results are reported in Figure 10, which give the horizontal, slant range and heading direction tracking errors, respectively, averaged over the eight target types for each range. It is shown that JVIM gains 9%, 10% and 35% improvements over CVIM, LL-GPLVM and NN along the horizontal direction, respectively, 35%, 31% and 72% along the slant range, respectively, and 5%, 13% and 62% along the heading direction, respectively. The results demonstrate that JVIM delivers better tracking performance with respect to all three figures of merit, with the advantage over CVIM, LL-GPLVM and NN being particularly significant for the range estimation.

7.4. ATR-Seg with Implicit Shape Matching

The proposed ATR-Seg algorithm (noted as Method I in the following) was tested against 15 SENSIAC sequences of five targets (SUV, BMP2, BTR70, T72 and ZSU23) under three ranges (1 km, 1.5 km and 2 km). Two more methods, Method II (JVIM with explicit shape matching, [23,26]) and Method III (CVIM [8]), were considered for comparison. All methods share a similar inference algorithm shown in Figure 7. Both Methods II and III involve explicit shape matching, and JVIM was used for both Methods I and II, while CVIM was used for method III. The tracking results are shown in Table 3. Results for tanks were averaged over T72 and ZSU23, and those for APCs averaged over BTR70 and BMP2. It is shown that Method I outperformed Methods II and III by providing lower tracking errors. More importantly, unlike Methods II and III, which require target pre-segmentation, Method I accomplishes target segmentation along with tracking and recognition as a by-product.

During tracking, the target identity is also estimated frame-by-frame by three methods, and the recognition accuracy is calculated as the percentage of frames where the target types were correctly classified in terms of the six target classes. The overall recognition results of three methods are shown in Table 4, where all methods perform well, and Method I (ATR-Seg) still slightly and moderately outperforms Methods II and III, respectively. Especially, when the range is large, e.g., 2 km, the advantage of Method I over Method III is more significant. This is mainly due to the fact that target segmentation is less reliable when the target is small.

The tracking, recognition and segmentation results of Method I (ATR-Seg) against five 1.5-km sequences were shown in Figure 11, where the two best matched target types are presented to show sub-class target recognition. As shown in Figure 11 (the forth tracking result of ZSU23), part of ZSU23 is missing during tracking; the proposed method still can give an accurate segmentation and tracking result. ATR-Seg uses the intensity information from the present frame to build the energy term in Equation (20) that reduces the error accumulation over time and then evaluates how likely a hypothesized shape created by JVIM can segment a valid target at the predicted position. On the other hand, Method III in [8] uses the background subtraction results and involves an explicit shape comparison for evaluation, so the tracking and recognition results highly depend on the pre-segmentation results.

7.5. ATR-Seg Segmentation Results

We evaluated the segmentation performance of ATR-Seg using the metric of the overlap ratio. The ground-truth segmentation results were generated manually for five randomly selected frames in each of 15 sequences. For a comparison, we also computed the overlap ratios for background subtraction results, which are averaged around 81%. While those of ATR-Seg are averaged around 85%. It is worth noting that the segmentation results of ATR-Seg are essentially constrained by the training shapes created from the CAD models, and the training models may have some shape discrepancy with the observed targets in the SENSIAC data. Another source of segmentation errors is due to tracking errors. Some segmentation results of five targets at 1.5 km were shown in Figure 12, where we overlaid the ATR-Seg results (contours) over the ground-truth ones. Background subtraction is not easy for a moving platform and is susceptive to the occlusion problem, while ATR-Seg is more flexible and robust to the sensor ego-motion and has great potential for occlusion handling, due to the shape prior involved [6,7].

7.6. Limitation and Discussion

There are two limitations of ATR-Seg due to the unsupervised nature of the level set, where no prior is used for foreground/background pixel intensities, and the mismatching between the training targets and the test ones. Thus, when a major part of a target is occluded or part of a target is similar to the background, the shape-aware set will lose the sensitivity for segmentation evaluation, leading to tracking failure, as shown in Figure 13 (first row), which shows the failed results for the pick-up sequence at 1.5 km. The mismatching and the low-quality data are the main reasons for the tracking failure of 2S3 at a range of 1.5 km (second row in Figure 13). One possible remedy is to incorporate some pixel priors of background and foreground into the level set energy function. However, an online learning scheme may be needed to update the pixel priors that are usually necessary for a long infrared sequence [31]. It is worth emphasizing that the goal of this work is to test the potential of a “model-based” approach that only uses CAD models for training. It is a natural extension to incorporate real infrared data for training that is likely to improve the algorithm robustness and applicability significantly.

8. Conclusion

A new algorithm, called ATR-Seg, is proposed for joint target tracking, recognition and segmentation in infrared imagery, which has three major technical components. First is a novel GPLVM-based shape generative model, the joint view-identity manifold (JVIM), which unifies one view-independent identity manifold and infinite identity-dependent view manifolds jointly in a semantically meaningful latent space. Second is the incorporation of a shape-aware level set energy function that evaluates how likely a valid target can be segmented by a shape synthesized by JVIM. Third, a particle filter-based sequential inference algorithm is developed to jointly accomplish target tracking, recognition and segmentation. Specifically, the level set energy function is used as the likelihood function in the particle filter that performs implicit shape matching, and a general motion model is involved to accommodate a highly maneuvering target. Experimental results on the recent SENSIAC ATR database manifest the advantage of ATR-Seg over two existing methods using explicit shape matching. This work is mainly focused on a shape-based approach. One possible future research issue is to involve other visual cues, such as pixel intensities or textures, to enhance the sensitivity and discriminability of the shape-aware level set energy function, which could mitigate the limitations of the ATR-Seg algorithm.

Author Contributions

This work was supported in part by the U.S. Army Research Laboratory (ARL) and U.S. Army Research Office (ARO) under grant W911NF-04-1-0221 and and the Oklahoma Center for the Advancement of Science and Technology (OCAST) under grant HR12-30.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Yilmaz, A.; Javed, O.; Shah, M. Object tracking: A survey. ACM Comput. Surv. 2006, 38, 1–45. [Google Scholar]
  2. Fan, X.; Fan, G.; Havilcek, J. Generative Models for Maneuvering Target Tracking. IEEE Trans. Aerospace Electron. Syst. 2010, 46, 635–655. [Google Scholar]
  3. Srivastava, A. Bayesian filtering for tracking pose and location of rigid targets. Proc. SPIE 2000, 4052, 160–171. [Google Scholar]
  4. Shaik, J.; Iftekharuddin, K. Automated tracking and classification of infrared images. Proceedings of International Joint Conference on Neural Networks, Portland, OR, USA, 20–24 July 2003; Volume 2, pp. 1201–1206.
  5. Khan, Z.; Gu, I.H. Tracking visual and infrared objects using joint Riemannian manifold appearance and affine shape modeling. Proceedings of IEEE International Conference on Computer Vision Workshops (ICCVW), Barcelona, Spain, 6–13 November 2011.
  6. Prisacariu, V.; Reid, I. Nonlinear shape manifolds as shape priors in level set segmentation and tracking. Proceedings of IEEE International Conference on Computer Vision and Pattern Recognition (CVPR), Colorado Springs, CO, USA, 20–25 June 2011; pp. 2185–2192.
  7. Prisacariu, V.; Reid, I. Shared shape spaces. Proceedings of IEEE International Conference on Computer Vision (ICCV), Barcelona, Spain, 6–13 November 2011; pp. 2587–2594.
  8. Venkataraman, V.; Fan, G.; Yu, L.; Zhang, X.; Liu, W.; Havlicek, J.P. Automated Target Tracking and Recognition using Coupled View and Identity Manifolds for Shape Representation. EURASIP J. Adv. Signal Process. 2011, 124, 1–17. [Google Scholar]
  9. Urtasun, R.; Fleet, D.J.; Geiger, A.; Popović, J.; Darrell, T.J.; Lawrence, N.D. Topologically-constrained latent variable models. Proceedings of International Conference on Machine learning (ICML), Helsinki, Finland, 5–9 July 2008; pp. 1080–1087.
  10. Yao, A.; Gall, J.; Gool, L.; Urtasun, R. Learning Probabilistic Non-Linear Latent Variable Models for Tracking Complex Activities. Proceedings of Annual Conference on Neural Information Processing Systems (NIPS), Granada, Spain, 12–14 December 2011; pp. 1–9.
  11. Military Sensing Information Analysis Center (SENSIAC). Available online: https://www.sensiac.org/ (accessed on 12 December 2012).
  12. Yankov, D.; Keogh, E. Manifold Clustering of Shapes. Proceedings of International Conference on Data Mining, Hong Kong, China, 18–22 December 2006.
  13. Roweis, S.T.; Saul, L.K. Nonlinear dimensionality Reduction by Local Linear Embedding. Science 2000, 290, 2323–2326. [Google Scholar]
  14. Etyngier, P.; Segonne, F.; Keriven, R. Shape Priors using Manifold Learning Techniques. Proceedings of IEEE International Conference on Computer Vision (ICCV), Rio de Janeiro, Brazil, 14–20 October 2007.
  15. Etyngier, P.; Keriven, R.; Segonne, F. Projection onto a Shape Manifold for Image Segmentation with Prior. Proceedings of IEEE International Conference on Image Processing (ICIP), San Antonio, TX, USA, 16–19 October 2007; Volume 4, pp. 361–364.
  16. He, R.; Lei, Z.; Yuan, X.; Li, S. Regularized active shape model for shape alignment. Proceedings of IEEE International Conference on Automatic Face Gesture Recognition (FG), Amsterdam, the Netherlands, 17–19 September 2008.
  17. Lüthi, M.; Albrecht, T.; Vetter, T. Probabilistic Modeling and Visualization of the Flexibility in Morphable Models. Proceedings of IMA International Conference on Mathematics of Surfaces XIII, York, UK, 7–9 September 2009; Volume 5654, pp. 251–264.
  18. Dambreville, S.; Rathi, Y.; Tannenbaum, A. A Framework for Image Segmentation Using Shape Models and Kernel Space Shape Priors. IEEE Trans. Pattern Anal. Mach. Intell. 2008, 30, 1385–1399. [Google Scholar]
  19. Lawrence, N. Probabilistic Non-linear Principal Component Analysis with Gaussian Process Latent Variable Models. J. Mach. Learn. Res. 2005, 6, 1783–1816. [Google Scholar]
  20. Elgammal, A.; Lee, C.S. Separating style and content on a nonlinear manifold. Proceedings of IEEE International Conference on Computer Vision and Pattern Recognition (CVPR), Washington, DC, USA, 27 June– 2 July 2004; Volume 1, pp. 478–485.
  21. Lee, C.; Elgammal, A. Modeling View and Posture Manifolds for Tracking. Proceedings of IEEE International Conference on Computer Vision, Rio de Janeiro, Brazil, 14–21 October 2007.
  22. Li, X.R.; Jilkov, V. Survey of maneuvering target tracking. Part I. Dynamic models. IEEE Trans. Aerospace Electron. Syst. 2003, 39, 1333–1364. [Google Scholar]
  23. Gong, J.; Fan, G.; Yu, L.; Havlicek, J.; Chen, D. Joint view-identity manifold for target tracking and recognition. Proceedings of 2012 19th IEEE International Conference on Image Processing (ICIP), Orlando, FL, USA, 30 September– 3 October 2012; pp. 1357–1360.
  24. Gong, J.; Fan, G.; Yu, L.; Havlicek, J.P.; Chen, D.; Fan, N. Joint view-identity manifold for infrared target tracking and recognition. Comput. Vis. Image Underst. 2014, 118, 211–224. [Google Scholar]
  25. Urtasun, R.; Darrell, T. Sparse probabilistic regression for activity-independent human pose inference. Proceedings of IEEE International Conference on Computer Vision and Pattern Recognition (CVPR), Anchorage, AK, USA, 23–28 June 2008.
  26. Gong, J.; Fan, G.; Yu, L.; Havlicek, J.P.; Chen, D.; Fan, N. Joint View-Identity Manifold for Infrared Target Tracking and Recognition. Comput. Vis. Image Underst. 2014, 118, 211–224. [Google Scholar]
  27. Bibby, C.; Reid, I. Robust Real-Time Visual Tracking Using Pixel-Wise Posteriors. Proceedings of the 10th European Conference on Computer Vision: Part II, Marseille, France, 12–18 October 2008; pp. 831–844.
  28. Gong, J.; Fan, G.; Havlicek, J.P.; Fan, N.; Chen, D. Infrared Target Tracking Recognition Segmentation using Shape-Aware Level Set. Proceedings of IEEE International Conference on Image Processing (ICIP), Melbourne, Australia, 15–18 September 2013.
  29. Arulampalam, S.; Maskell, S.; Gordon, N.; Clapp, T. A Tutorial on Particle Filters for Online Non-linear/Non-Gaussian Bayesian Tracking. IEEE Trans. Signal Process. 2002, 50, 174–188. [Google Scholar]
  30. Zivkovic, Z.; van der Heijden, F. Efficient adaptive density estimation per image pixel for the task of background subtraction. Pattern Recognit. Lett. 2006, 27, 773–780. [Google Scholar]
  31. Venkataraman, V.; Fan, G.; Havlicek, J.; Fan, X.; Zhai, Y.; Yeary, M. Adaptive Kalman Filtering for Histogram-based Appearance Learning in Infrared Imagery. IEEE Trans. Image Process. 2012, 21, 4622–4635. [Google Scholar]
Figure 1. Graphical modeling for the proposed ATR-Seg algorithm, where It represents an image frame, p 3D target position, Φ target segmentation, and Θ the set of shape variables.
Figure 1. Graphical modeling for the proposed ATR-Seg algorithm, where It represents an image frame, p 3D target position, Φ target segmentation, and Θ the set of shape variables.
Sensors 14 10124f1 1024
Figure 2. Flowchart for ATR-Seg.
Figure 2. Flowchart for ATR-Seg.
Sensors 14 10124f2 1024
Figure 3. JVIM learning and shape interpolation, where one view-independent identity manifold and six identity-dependent view manifolds are color-coded according to the uncertainty of GP mapping. (Adapted from [24], with permission from Elsevier.)
Figure 3. JVIM learning and shape interpolation, where one view-independent identity manifold and six identity-dependent view manifolds are color-coded according to the uncertainty of GP mapping. (Adapted from [24], with permission from Elsevier.)
Sensors 14 10124f3 1024
Figure 4. Shape-aware level set model for implicit shape matching. (a) Illustration of a target in an infrared image: foreground Ωf and background Ωb, foreground/background intensity models M, and the 3D-2D camera projection W(x,p). (b) The shape embedding function Φ. (c) The graphical model for shape-aware level set, where p is the target 3D location of a ground-vehicle, and Θ is the shape parameter in JVIM.
Figure 4. Shape-aware level set model for implicit shape matching. (a) Illustration of a target in an infrared image: foreground Ωf and background Ωb, foreground/background intensity models M, and the 3D-2D camera projection W(x,p). (b) The shape embedding function Φ. (c) The graphical model for shape-aware level set, where p is the target 3D location of a ground-vehicle, and Θ is the shape parameter in JVIM.
Sensors 14 10124f4 1024
Figure 5. Optimization of ATR-Seg by a gradient descent method.
Figure 5. Optimization of ATR-Seg by a gradient descent method.
Sensors 14 10124f5 1024
Figure 6. Possible reasons for the failure of the gradient descent method.
Figure 6. Possible reasons for the failure of the gradient descent method.
Sensors 14 10124f6 1024
Figure 7. The graphical model representation of ATR-Seg and the 3D camera coordinate. (Reprint from [28] with permission from IEEE).
Figure 7. The graphical model representation of ATR-Seg and the 3D camera coordinate. (Reprint from [28] with permission from IEEE).
Sensors 14 10124f7 1024
Figure 8. All 36 CAD models used in this work, which are ordered according to the class-constrained shortest-closed-path [8]. (Reprint from [24], with permission from Elsevier.)
Figure 8. All 36 CAD models used in this work, which are ordered according to the class-constrained shortest-closed-path [8]. (Reprint from [24], with permission from Elsevier.)
Sensors 14 10124f8 1024
Figure 9. Qualitative analysis of JVIM shape interpolation: (a) along six identity-specific view manifolds. (b) along the view-independent identity manifold between two training target types. (Reprint from [24], with permission from Elsevier.)
Figure 9. Qualitative analysis of JVIM shape interpolation: (a) along six identity-specific view manifolds. (b) along the view-independent identity manifold between two training target types. (Reprint from [24], with permission from Elsevier.)
Sensors 14 10124f9 1024
Figure 10. Comparison of the tracking errors of the horizontal position, slant range and heading direction. In each plot, the results for each method averaged over eight target types for each range. From left to right, the plot gives the results for JVIM (first, red), couplet of view and identity manifolds (CVIM) (second, blue), local linear (LL)-Gaussian process latent variable model (GPLVM) (third, cyan) and nearest neighbor (NN) (forth, green). (Reprint from [24], with permission from Elsevier).
Figure 10. Comparison of the tracking errors of the horizontal position, slant range and heading direction. In each plot, the results for each method averaged over eight target types for each range. From left to right, the plot gives the results for JVIM (first, red), couplet of view and identity manifolds (CVIM) (second, blue), local linear (LL)-Gaussian process latent variable model (GPLVM) (third, cyan) and nearest neighbor (NN) (forth, green). (Reprint from [24], with permission from Elsevier).
Sensors 14 10124f10 1024
Figure 11. ATR-Seg results for five IR sequences. Column 1: truth target types. Columns 2–5: selected IR frames overlayed with the segmentation results. Columns 6–7: the two best matched training targets along the identity manifold. (Reprint from [28], with permission from IEEE.)
Figure 11. ATR-Seg results for five IR sequences. Column 1: truth target types. Columns 2–5: selected IR frames overlayed with the segmentation results. Columns 6–7: the two best matched training targets along the identity manifold. (Reprint from [28], with permission from IEEE.)
Sensors 14 10124f11 1024
Figure 12. Segmentation results of five targets at the range of 1.5 km.
Figure 12. Segmentation results of five targets at the range of 1.5 km.
Sensors 14 10124f12 1024
Figure 13. Tracking failure for the pick-up and 2S3 sequences at 1.5 km.
Figure 13. Tracking failure for the pick-up and 2S3 sequences at 1.5 km.
Sensors 14 10124f13 1024
Table 1. Descriptions of all mathematical symbols.
Table 1. Descriptions of all mathematical symbols.
Symbols used in Problem Formulation (Section 3)

pthe target's 3D position p = [px, py, pz]T
αthe target type
φthe view angle
Θshape related variables ([α, φ])
Φa target shape segmentation
Itan observed image frame at time t

Symbols used in JVIM-based shape modeling (Section 4)

YJVIM training data
XJVIM latent space
θthe aspect angle of a target
ϕthe elevation angle of a target
βthe kernel hyper-parameters of JVIM
dthe dimension of the shape space
wthe LLE coefficients for local topology encoding
LJVIMthe JVIM objective function
LDthe data term in LJVIM
LTthe topology term in LJVIM
KYthe covariance matrix of JVIM learning
xra reference latent point in JVIM learning
XRthe neighborhood of xr for local learning
M1the size of XR (the range of local learning)
YRthe corresponding shape for XR
Nthe size of training data
x′a new latent point for JVIM-based shape interpolation
X′the neighborhood of x′ for local inferencing
M2the size of X′ (the range of local inferencing)
Y′the corresponding shape data for X′
k(x1, x2)a RBF kernel function in JVIM
μ̂x′an interpolated shape at x′ via JVIM
σ ^ x 2uncertainty of shape interpolation at x′

Symbols used in shape-aware level set (Section 5)

xa 2D pixel location in an image frame
ya pixel intensity value
Mforeground/background models M = {Mf, Mb}
Hϵ [·]the smoothed Heaviside step function

Symbols used in sequential inference (Section 6)

ψtthe heading direction of a ground vehicle in frame t
υtthe target velocity along ψt in frame t
Δtthe time interval of two adjacent frames
Ztthe state vector in frame t ( Z t = [ p t T , υ t , ψ t , α t ] T)
Table 2. Pseudo-code for ATR-Seg algorithm.
Table 2. Pseudo-code for ATR-Seg algorithm.
• Initialization: Initialize the target position, p0, type α0, heading direction ψ0 and speed υ0 according to the ground-truth and get the initial state, Z0. Draw Z 0 j N ( Z 0 , 1 ) , j { 1 , , N p }, Np is the number of particles.
• For t = 1,…, T (number of frames)
 1. For j = 1, …, Np
   1.1 Draw samples Z t j p ( Z t j | Z t 1 j ) as in Equation (21).
   1.2 Generate the target shape according to the target state using Equations (9) and (10).
   1.3 Compute weights w t j = p ( z t | α t j , Z t j ) using Equation (16).
  End
 2. Normalize the weights, such that j = 1 N p w t j = 1.
 3. Compute the mean estimates of the target state, Z ^ t = j = 1 N p w t j Z t j
 4. Set Z t j = resample ( Z k j , w k j ) to increase the effective number of particles [29].
• End
Table 3. Tracking errors for three ATR methods (Method I/Method II/Method III). (Reprint from [28], with permission from IEEE).
Table 3. Tracking errors for three ATR methods (Method I/Method II/Method III). (Reprint from [28], with permission from IEEE).
RangeError inTankAPCSUVTotal
1 kmpx (m)0.22/0.25/0.220.22/0.25/0.180.16/0.17/0.190.21/0.23/0.20
pz (m)5.06/8.67/7.534.19/4.48/5.149.03/8.36/10.955.51/6.93/7.26
ψ (rad)0.13/0.17/0.180.15/0.32/0.150.11/0.24/0.220.13/0.24/0.18
1.5 kmpx (m)0.24/0.19/0.180.15/0.21/0.200.16/0.56/0.600.27/0.27/0.27
pz (m)4.40/7.20/7.284.70/5.88/5.96– – NA – –4.55/6.54/6.26
ψ (rad)0.16/0.22/0.240.18/0.53/0.510.11/0.32/0.350.20/0.36/0.37
2 kmpx (m)0.31/0.27/0.280.23/0.19/0.360.13/0.17/0.350.24/0.22/0.32
pz (m)8.68/10.6/8.588.95/9.28/7.955.35/8.09/14.258.19/9.55/9.46
ψ (rad)0.19/0.38/0.260.41/0.18/0.380.08/0.41/0.310.26/0.31/0.32
Table 4. Recognition accuracy (%) for Methods I, II and III. (Reprint from [28], with permission from IEEE).
Table 4. Recognition accuracy (%) for Methods I, II and III. (Reprint from [28], with permission from IEEE).
TargetsTanksAPCsSUVTotal
1 km100/96/96100/100/94100/100/100100/98/96
1.5 km98/96/9499/100/89100/100/10099/98/93
2 km98/92/86100/100/85100/100/9899/98/88

Share and Cite

MDPI and ACS Style

Gong, J.; Fan, G.; Yu, L.; Havlicek, J.P.; Chen, D.; Fan, N. Joint Target Tracking, Recognition and Segmentation for Infrared Imagery Using a Shape Manifold-Based Level Set. Sensors 2014, 14, 10124-10145. https://0-doi-org.brum.beds.ac.uk/10.3390/s140610124

AMA Style

Gong J, Fan G, Yu L, Havlicek JP, Chen D, Fan N. Joint Target Tracking, Recognition and Segmentation for Infrared Imagery Using a Shape Manifold-Based Level Set. Sensors. 2014; 14(6):10124-10145. https://0-doi-org.brum.beds.ac.uk/10.3390/s140610124

Chicago/Turabian Style

Gong, Jiulu, Guoliang Fan, Liangjiang Yu, Joseph P. Havlicek, Derong Chen, and Ningjun Fan. 2014. "Joint Target Tracking, Recognition and Segmentation for Infrared Imagery Using a Shape Manifold-Based Level Set" Sensors 14, no. 6: 10124-10145. https://0-doi-org.brum.beds.ac.uk/10.3390/s140610124

Article Metrics

Back to TopTop