Research Article  Open Access
Articulated Human Motion Tracking Using Sequential Immune Genetic Algorithm
Abstract
We formulate human motion tracking as a highdimensional constrained optimization problem. A novel generative method is proposed for human motion tracking in the framework of evolutionary computation. The main contribution is that we introduce immune genetic algorithm (IGA) for pose optimization in latent space of human motion. Firstly, we perform human motion analysis in the learnt latent space of human motion. As the latent space is low dimensional and contents the prior knowledge of human motion, it makes pose analysis more efficient and accurate. Then, in the search strategy, we apply IGA for pose optimization. Compared with genetic algorithm and other evolutionary methods, its main advantage is the ability to use the prior knowledge of human motion. We design an IGAbased method to estimate human pose from static images for initialization of motion tracking. And we propose a sequential IGA (SIGA) algorithm for motion tracking by incorporating the temporal continuity information into the traditional IGA. Experimental results on different videos of different motion types show that our IGAbased pose estimation method can be used for initialization of motion tracking. The SIGAbased motion tracking method can achieve accurate and stable tracking of 3D human motion.
1. Introduction
Tracking articulated 3D human motion from video is an important problem in computer vision which has many potential applications, such as virtual character animation, human computer interface, intelligent visual surveillance, and biometrics. Despite having been attacked by many researchers, this challenging problem is still long standing because of the difficulties conduced mainly by the complicated nature of 3D human motion, selfocclusions, and highdimensional search space.
In the previous work, two main classes of motion tracking approaches can be identified: discriminative approaches and generative approaches [1]. Discriminative methods attempt to learn a direct mapping from image features to 3D pose using training data. The mapping is often approximated using nearest neighbor [2], regression models [3] or mixture of regressors [4]. Discriminative approaches are effective and fast. However, they need a large training database and are limited to fixed classes of motion. Moreover, the inherent onetomany mapping from 2D images to 3D poses is difficult to learn accurately. In contrast, generative methods exploit the fact that although the mapping from visual features to poses is complex and multimodal, the reverse mapping is often well posed. Therefore, pose recovery is tackled by optimizing an object function that encodes the posefeature correspondence [5] or by sampling posterior pose probabilities [6]. Compared with discriminative methods, generative methods are usually more accurate. However, generative methods are generally computationally expensive because one has to perform complex search over the highdimensional pose state space in order to locate the peaks of the observation likelihood. Moreover, prediction model and initialization are also the bottlenecks of the approach in the tracking scenario. In this work, we focus on recovering 3D human pose within the generative framework.
In general, highdimensional state space and search strategy are two main problems in generative approaches. Highdimensional pose state space makes pose analysis computationally expensive or even infeasible. Despite the high dimensionality of the configuration space, many human motion activities lie intrinsically on lowdimensional latent space [7, 8]. Motivated by this observation, we use ISOMAP, a nonlinear dimensionality reduction method, to learn the lowdimensional latent space of pose state, by which the aim of both reducing dimensionality and extracting the prior knowledge of human motion are achieved simultaneously. On the other hand, search strategy, in general how to track in the lowdimensional latent space, is another important problem. The search strategy should suit for the characteristics of the subspace and be global, optimal, and convergent. Although considerable work has already been done, a more effective search strategy is still intensively needed for robust visual tracking. In our opinion, motion prior knowledge has great influence on the search strategy, which can aid in performing more stable tracking. Compared with the previous methods, extracting the prior knowledge and introducing it in the designing of search strategy are of particular interests to us.
In this paper, we propose a novel generative approach in the framework of evolutionary computation, by which we try to widen the bottlenecks mentioned above with effective search strategy embedded in the extracted state subspace. The framework of our approach is illustrated in Figure 1. Firstly, we use ISOMAP to learn this latent space. Then we propose a manifold reconstruction method to establish the inverse mapping, which enables pose analysis in this latent space. As the latent space is low dimensional and contents the prior knowledge of human motion, it makes pose analysis more efficient and accurate. In the search strategy we introduce immune genetic algorithm (IGA) for pose optimization. Details of the implementations, such as encoding and initialization, computation of affinity, and genetic and immunity operators, are designed. We propose an IGAbased method for pose estimation, which can be used for initialization of motion tracking. In order to make IGA suitable for human motion tracking, a sequential IGA (SIGA) framework is proposed by incorporating the temporal continuity information into the traditional IGA. Experimental results on different motion types and different image sequences demonstrate our methods.
The rest of the paper is organized as follows. Section 2 gives an introduction to the related works. Section 3 gives a description of how the latent space is learnt. In Section 4, we give a detailed description of how we apply IGA for pose optimization in the latent space. We then show how to apply IGAbased pose optimization algorithm for pose estimation and tracking in Section 5. Section 6 contains experimental results and comparison with other tracking algorithms. The conclusions and possible extension for future work are given in Section 7.
2. Related Works
There has been a great deal of prior works on human motion analysis from video [1, 8, 9]. Here we focus our survey on the most related research on generative methods. In generative human motion tracking methods, the highdimensional pose state space is the most significant problem. There are several possible strategies for reducing the dimensionality of the configuration space, including using motion models [10], hierarchical search [11], and dimensionality reduction [5, 12]. Motion models are often derived from training data of a single class of movement. Although they can aid in performing more stable tracking, this comes at the cost of putting a strong restriction on the poses that can be recovered. Another way to constrain the configuration space is to perform a hierarchical search. For example, John et al. [11] proposed a hierarchical particle swarm optimization method to search the best pose hierarchically. An inherent problem with this approach is the need to estimate accurately the position and orientation of the initial body segment (typically the torso), as a wrong pose estimate for the initial segment can distort the pose estimates for subsequent limbs. Nowadays, dimensionality reduction has become the most widely used methods. For example, Urtasun et al. [12] construct a differentiable objective function based on the Principle Component Analysis (PCA) of motion capture data and then find the poses of all frames simultaneous by optimizing a function in lowdimensional space. However, this method needs many example sequences of data to perform PCA, and all of these sequences must keep the same length and same phase by interpolating and aligning. Zhao and Liu [5] use PCA to learn the lowdimensional state space of human pose and perform pose analysis in the latent space. However since the mapping between the original pose space and the latent space is in general nonlinear, linear PCA is inadequate. Nonlinear dimensionality reduction methods have also been used. For example, Sminchisescu and Jepson [13] use spectral embedding to learn the embedding which is modeled as a Gaussian mixture model. Radial Basis Functions (RBFs) are learned for inverse mapping. A linear dynamical model is used for tracking. Elgammal and Lee [14] learn viewbased representations of activity manifolds using nonlinear dimensionality reduction method (LLE). Then, the nonlinear mappings from the embedding space into both visual input space and 3D pose space are learnt using the generalized radial basis function. Although nonlinear dimensionality reduction methods can learn this nonlinear mapping, they are not invertible. The smooth inverse mapping is still a not wellsolved problem. In this paper, we use ISOMAP [15], a nonlinear dimensionality reduction method, to learn the low dimensionality subspace of a specific activity. And then, based on the intrinsic executive mechanism of ISOMAP, a manifold reconstruction method is proposed to generate smooth mappings between the subspace and the original space. This enables us to perform human motion tracking in the learned subspace.
Search strategy is another key research problem of pose tracking in the generative framework. They are typically tackled using either deterministic methods or stochastic methods. Deterministic methods usually involve a gradient descent search to minimize a cost function [16]. Although these methods are usually computationally efficient, they easily become trapped in local minima. In contrast, stochastic methods introduce some stochastic factors into the searching process in order to have a higher probability of reaching the global optimum of the cost function. Particle filter [6] is the most wildly studied stochastic method which is based on Monte Carlo sampling. Although, in theory, particle filter is very suitable for tracking, it needs a large number of particles to approximate the posterior distributions, and it tends to suffer from sample impoverishment, so that the final particle sets cannot represent the true distributions. Therefore, many improvements on the traditional particle filter have been proposed. For example, Deutscher et al. [17] introduced the annealed particle filter which combines a deterministic annealing approach with stochastic sampling to reduce the number of samples required. At each time step the particle set is refined through a series of annealing cycles with decreasing temperature to approximate the local maxima in the fitness function. In Krzeszowski et al. [18], a particle swarm optimization algorithm is utilized in the particle filter to shift the particles toward more promising configurations of the human model. Compared with the deterministic counterparts, stochastic methods are usually more robust, but they suffer a large computational load, especially in highdimensional state space. In recent years, evolutional computing methods, such as genetic algorithm [5] and particle swarm optimization [11, 18, 19], have received increasing attention. For example, Zhao and Liu [5] proposed an annealed genetic algorithm to track human motion in compact base space, where the base space is learned using PCA. John et al. [11] proposed a hierarchical particle swarm optimization (HPSO) algorithm for articulated human tracking. Their comparative experimental results show that HPSO is more accurate than particle filter. However, based on our experimental results, due to the highdimensional pose state space and imperfect image observations, HPSO may deviate from the pose state space and result in inaccurate tracking. Evolutionary algorithms are all good searching algorithms with an iterative process of generation and test. Two operators, crossover and mutation, give each individual the chance of optimization and ensure the evolutionary tendency with the select mechanism of survival of the fittest. However, the two operators change individuals randomly and indirectly under some conditions. Therefore, they not only give individuals the evolutionary chance but also cause certain degeneracy. Recently, immune algorithms have been another hotspot succeeding genetic algorithm and particle swarm optimization for its success in solving pattern recognition and optimization problems. Its main advantage, compared with GA and PSO, is it has the ability to use the prior knowledge of problem by vaccination and immune selection [20]. In this paper, we apply immune genetic algorithm (IGA) [20], a novel immune method, for pose optimization. We propose an IGAbased method for pose estimation from monocular images. In order to make IGA suitable for pose tracking, we propose a sequential IGA (SIGA) algorithm by incorporating the temporal continuity information into the traditional IGA. To the best of our knowledge, the proposed algorithm is new in the human motion tracking literature.
3. Learning the Latent Space of Human Motion
Tracking in a lowdimensional latent space requires three components [8]. First, a mapping between original pose space and lowdimensional subspace must be learned. Second, an inverse mapping must be defined. Third, how tracking within the lowdimensional space occurs must be defined. In this section, we first learn the lowdimensional subspace using ISOMAP [15]. Then, we propose a manifold reconstruction method to establish the mappings between high and lowdimensional states.
3.1. ISOMAPBased Latent Space Learning
We describe the human body as a kinematic tree consisting of rigid limbs that are linked by joints. Every joint contains a number of degrees of freedom (DOF), indicating in how many directions the joint can move. All DOF in the body model together form the pose representation. In this paper, the pose is described by a 66D vector , where 3D vector represents the root joint rotations and 63D vector represents the body joints rotations. Apart from the kinematic structure, the human shape is also modeled. Each rigid limb of the body is fleshed out using conic sections with elliptical crosssections (see Figure 2). Human shape will be used to compute the likelihood function (see Section 4.2).
(a)
(b)
Since the mapping between the original pose space and latent space is in general nonlinear, linear PCA is inadequate. So we use ISOMAP to learn the nonlinear mapping. We extract the subspace using motion capture data obtained from the CMU database [21].
As for a special activity, such as walking, running, jumping, and so forth, the original pose state space has no relation with the global motion. Different from the previous methods of learning different manifolds for the same activity (such as walking) of different views, we filter out the rotations of root joint () and represent the pose using the rotations of body joints () only. Assuming is a given sequence of motion capture data corresponding to one motion type, where , is the frame index, is the number of total frames, and is the original pose state space, the subspace is extracted by ISOMAP as follows.(1) Construct Neighborhood Graph. Define the graph over all data points (in our method the data point is one frame in motion sequence) by connecting point and if . Set edge length to be . Moreover, āwhere is the dimensionality of pose state space; here.(2) Compute Shortest Paths. Initialize if and are linked by an edge; otherwise. Then for each value of in turn, replace all entries by . The matrix of final values will contain the shortest path distances between all pairs of points in .(3) ConstructāāDimensional Embedding. Let be the th eigenvalue (in decreasing order) of the matrix and the th component of the th eigenvector. Then set the th component of the dimensional coordinate vector to be equal to .
The subspace learned by ISOMAP is shown in Figure 3. Actually, similar lowdimensional subspace can be extracted from the training sequences that belong to the same type of motions but performed by different subjects. And the training sequences corresponding to different types of motions produce different subspaces. For example, experiments demonstrate that different walking sequences generate similar manifolds in the 3D subspace, which is different from that of running motion.
(a)
(b)
ISOMAP cannot only reduce the dimensionality of highdimensional input space, but also find meaningful lowdim structures hidden behind their highdim observations. In doing so, infeasible solutions, namely, the absurd poses, can be avoided naturally during optimization, which will make pose tracking in this subspace more efficient and accurate.
3.2. Mapping between High and LowDimensional States
Traditional ISOMAP can only learn the mapping from the original pose space to the latent space but not the inverse mapping. However, in order to track human motion in the lowdimensional manifold, the inverse mapping is required. Based on the intrinsic executive mechanism of ISOMAP, we proposed an ISOMAPbased manifold reconstruction method to establish the mapping between high and lowdimensional states.
Suppose the pose state space to be and the lowdim state space to be . Denote the mapping as ,āā and ,āā, where , are the high and lowdimensional vectors, respectively. The set of input instances is , and their corresponding points in the embedding space learned by ISOMAP are . Assume are the neighbors of point , where is the number of neighbors. And their corresponding points in the embedding space are . Then our ISOMAPbased manifold reconstruction method can be described as in Algorithmāā1.

Using the ISOMAPbased manifold reconstruction method, we can generate smooth mapping between the original pose space and the latent space, which enables us to track human pose in the latent space. In the following section, we will show how tracking within the latent space occurs.
4. Immune Genetic Algorithm for Pose Optimization
We formulate pose estimation as a constrained optimization problem and solve it using immune genetic algorithm. In this section, we first give a brief introduction to IGA. Then, we design an IGAbased method for pose optimization.
4.1. Immune Genetic Algorithm
In IGA, the idea of immunity is mainly realized through two steps based on reasonably selecting vaccines, that is, a vaccination and an immune selection, of which the former is used for raising fitness and the latter is for preventing the deterioration. A very clear overview of IGA, from immunology and engineering points of view, is presented in [20].
4.1.1. Immunological Terms
In order to describe the IGA clearly, some immunological terms will be given first [22].
Antigen. In immunology, an antigen is any substance that causes immune system to produce antibodies against it. In this paper, IGA is used for optimization: where , is the feasible region, is the number of problematic parameter, and the antigen is defined as the objective function .
Antibody and Antibody Population. In this paper, an antibody is a representation of a candidate solution of an antigen. The antibody is the coding of variable , denoted by , and is called the decoding of antibody , expressed as . The representation of antibody varies with antigen and can be binary string, real number sequence, symbolic sequence, and characteristic sequence. In this study, we adopt realcoded representation, that is, .
An antibody population, is an dimensional group of antibody , where the positive integer is the size of antibody population .
Affinity. In immunology, affinity is the fitness measurement for an antibody. For the optimization problem, the affinity, , is a mapping of the objective function for a given antibody .
4.1.2. Description of IGA
In this paper, we use the IGA for optimization task. The flow chart of IGA is shown in Figure 4. The main steps of our modified IGA can be summarized as follows.
Step 1. Initialization: randomly generate the initial antibody population ; set .
Step 2. Vaccine construction: abstract vaccines according to the prior knowledge.
Step 3. Evaluation: calculate the affinities of all antibodies in .
Step 4. Termination test: if termination criteriion is satisfied, export the antibody having the highest affinity in as the output of the algorithm and stop the algorithm; otherwise, continues.
Step 5. Genetic operators: perform genetic operators on the th parent āāand obtain the results .
Step 6. Vaccination: perform vaccination on and obtain the results .
Step 7. Immune selection: perform immune selection on and obtain the next parent , then go to Step 3.
In general, the IGA algorithm is to be implemented as the following evolvement process: where , , and are the antibody populations during different periods in a single evolution generation, is the iterative step. , , and are the genetic, vaccination, and immune selection operators, respectively.
4.2. Apply IGA for Pose Optimization
In this section, we apply IGA for human pose optimization. Some details of our implementations are discussed below.
4.2.1. Encoding and Initialization
In IGA, each antibody represents a potential solution in the search space. For our problem, we perform human motion analysis in the latent space. So an antibody is corresponding to a pose vector in the latent space. In this paper, we represent the full 3D pose vector as , where 3D vector represents the root joint rotations, corresponds to the pose vector in latent space; we set here. So is a 9dimensional vector. We use real encodings. We represent the antibody population as , whereāā is the size of population.
4.2.2. Computation of Affinity
For each antibody, an affinity measure needs to be computed to estimate how well a given antibody (pose) matches the observed images. Here we use the bidirectional likelihood proposed by [23]. Let represent the binary silhouette map for the body model and the image foreground. We seek to minimize the nonoverlapping regions, red and blue, therefore maximizing the Yellow region (see Figure 5). The size of each region can be computed by summing over all image pixels using
(a)
(b)
(c)
(d)
Then the objective function of candidate pose with regard to image can be calculated as where is the weight. We set in this paper.
Affinity is the fitness measurement for an antibody. As defined above, the affinity, , for a given antibody is a linear mapping of the objective function . Therefore, we define the affinity of antibody as where is a positive constant; we set in this paper.
4.2.3. Genetic Operators
We design five genetic operators, which are executed orderly in IGA. We introduce the operators by evolving an example antibody . The new antibody generated by the operators is denoted as . Assuming the positions generated randomly are numbers 2 and 6 or 3 (for point mutation operator) of , for example, the five operators are illustrated in Table 1. The application order of the genetic operators in the algorithm is just as that listed in Table 1.

The genetic operators were represented as . We perform genetic operator on the th parent and obtain the results .
4.2.4. Vaccine Construction and Immune Operators
Genetic operators give each antibody the chance of optimization and ensure the evolutionary tendency with the select mechanism of survival of the fittest. However, it changes individuals randomly and indirectly under some conditions. Therefore, they not only give individuals the evolutionary chance but also cause certain degeneracy. In IGA, the idea of immunity is mainly realized through two steps based on reasonably selecting vaccines, that is, a vaccination and an immune selection, of which the former is used for raising fitness and the latter is for preventing the deterioration.
In this section, we extract the prior knowledge of human motion and construct two vaccines. Then we design the vaccination and immune selection operators.(1)Vaccine Construction. A vaccine is abstracted from the prior knowledge of the pending problem. Human pose in subspace is located on a manifold structure but not the whole subspace. Actually, pose subspace is a compact space. We constrained the subspace of human motion and construct two vaccines for our human pose estimation problem. (i) Vaccine 1. Every dimensionality of subspace pose should be distributed in a scope as , where the bound and are learned from the motion training data.(ii) Vaccine 2. Vaccine 2 is motivated by the fact that every generated pose should locate on the manifold. Based on the consistency of human motion, we partition the manifolds into different subparts with means clustering, where the number of class is 5 in this paper (see Figure 6). For each class , , we assume the poses in it is of Gaussian distribution, described as follows: āwhere, is the mean vector, is the covariance matrix, is the dimensionality of the pose subspace. Then the vaccine 2 can be described as for all , such that , where .
(a)
(b)
4.3. IGABased Human Pose Optimization
Based on the description above, the IGAbased pose optimization algorithm can be described as in Algorithmāā2.

We will show how to apply IGAbased pose optimization algorithm for pose estimation and tracking in the next section.
5. Sequential Immune Genetic Algorithm for Pose Tracking
In tracking applications, the data is typically a time sequence, and hence the task is essentially a dynamic optimization problem which distinguishes it from traditional optimization problems. In tracking situation, the previous estimation results can be used to cut the current search space. From the Bayesā view, we can formulate the pose tracking problem as where and represent temporal states and observations, respectively. How to determine the conditional distribution effectively is the core problem for 3D human pose tracking.
In this paper, we proposed a sequential IGA (SIGA)based framework for human motion tracking. The flowchart of the SIGA framework is shown in Figure 7. First, we perform human pose estimation on the first frame of the video as initialization for tracking. Then, the previous converged antibodies at time are randomly propagated as initial antibodies for the next time (frame) . Finally, we perform IGAbased pose optimization on current antibodies. The individual with best affinity is used to approximate the tracking result of time , and the converged antibodies are used to initial the next frame. There are three major stages in the SIGA framework: automatically initialization, next frame propagation, and IGAbased optimization.
5.1. Pose Estimation for Initialization of Motion Tracking
Initialization is an important problem of human motion tracking. How to begin the tracking process from a good starting point sometimes is an intractable problem. We achieve the automatic initialization by determining the pose of the first frame using the IGAbased human pose estimation algorithm, which can be described as follows.
Pose estimation is the process to estimate articulated human pose from a single image which can be formulated as an optimization process. We apply IGA for pose estimation. For clarity, we redefine the full 3D pose vector as , where is the global motion of human body with respect to the camera and is the pose vector in state subspace. We perform the state posterior inference by optimizing the affinity function. The optimal pose can be represented as
We maximize the search efficiency by embedding the global search capability of IGA into the local conditions of state subspace.
The global motion of human body is very important for its visual appearance in an image and is also critical in disambiguating the leftright confusion. Determining this motion accurately makes our method viewpoint invariant. With the aim of both cutting the search space and determining the motion direction roughly, we incorporated the global motion process step [5] into the framework of IGA. The global motion process can be summarized as follows. (1) In state vector , the global motion include the rotation of the full body about the coordinate axes , , and , respectively. In the first round of state evolution (), we only actually search the optimal solutions of global motion. Other state components () are taken as one of the clustering centers ,āā, randomly. The variance domain and of is computed by storing the best antibodies. is determined empirically according to the threshold value of affinity. (2) In the rest rounds of state evolution, the antibody is evolved normally as described in Section 4. In doing so, we can get the coarse scopes of global motion in the first round of state evolution, and the fine tuning of these parameters can be achieved in the followed evolution rounds.
Based on the proposed IGA pose optimization algorithm, the antibody with the highest affinity in population will be selected to be the optimal pose. Figure 8 is the process of pose estimation, where (a) is one frame of input video, (b) is the initialized poses, (c), and (d) are results with 10, 40 times of iteration, respectively. We can see that the poses generated by our initialization method can cover the whole walking pose state space, and the poses become convergent with times of iteration increase.
5.2. NextFrame Propagation
Nextframe propagation is the key stage in the SIGA framework which aims to find the dynamic model . In this paper, we design a randomly propagation method. The randomly propagation method is actually a firstorder GaussMarkov dynamical model. Given the converged antibodies at frame , the antibodies in the next frame are initialized by sampling a Gaussian distribution centered in the current best antibodies.
Consider, where are the initial antibodies at time , are the converged antibodies at time , , and is the covariance matrix of Gaussian distribution. Low value will promote temporal consistency but is likely to lose the diversity. We set it empirically according to the motion type and speed. SIGA propagates only a minimal amount of information between frames and does not incorporate any motion model. Although randomly propagation is simple, it is sufficient because it is only used to produce an initial value for a subsequent search for the optimal state.
We do not incorporate any learnt constant motion model here, which is motivated by two considerations.(1)Generality: many prior motion models are derived from training data. A possible weakness of these motion models is that the ability to accurately represent the space of realizable human movements generally depends significantly on the amount of available training data. This comes as a cost of putting a strong restriction on the poses that can be recovered. Therefore, we do not use any constant learnt motion models here. (2)The effectiveness of our IGA pose optimization algorithm, which can explore efficiently large portions of the search space starting from the initial distribution of antibodies.
Actually, the SIGA framework is a āsampleandrefineā search strategy. Firstly, the initial antibodies are sampled for the transition distribution as . Then the antibodies are updated according to the newest observations in each IGA iteration. Through the IGA iteration, the antibodies are moved towards the region where the likelihood of observation has larger values and are finally relocated to the dominant modes of the likelihood. And in a Bayesian inference view, the IGA iterations are essentially a multilayer importance sampling strategy which incorporates the new observations into a sampling stage and thus avoids the sample impoverishment problem suffered by the particle filter [6].
5.3. Sequential Immune Genetic AlgorithmBased Pose Tracking
Based on the designing above, we can formulate our sequential IGA for pose tracking as in Algorithmāā3.

6. Experimental Results
6.1. Experimental Data and Evaluation Measures
Experimental Data
The data for latent space training is from CMU Database [21]. We quantitatively evaluate our method on synthesized image sequences as in [3]. We also give experimental results on real image sequences from [24], CMU Database [21], and HumanEva [23].
Evaluation Measures
In this paper, we use the evaluation measures proposed in [23]. The average error over all joint angles (in degrees) is defined as where and are the ground truth pose and the estimated pose date, respectively. For the sequence of frames, the average performance and the standard deviation of the performance are computed using the following:
6.2. The Convergence of IGA
It is understood that the number of antibodies and iteration times will affect the convergence. We take pose estimation experiment on a single image and report the affinities of the best antibody during iteration. Figure 9 demonstrates the convergence process. Different lines represent different numbers of antibodies used. The axis is the times of iteration while the axis is the affinity value. As shown in Figure 9, the affinities will converge as the times of iteration increase. The experimental results demonstrate that our IGAbased pose optimization is convergent.
We have ascertained experimentally that higher numbers of and will achieve better results. However, in order to deal with the tradeoff of computational time and accuracy, we set ,āā.
6.3. IGABased Pose Estimation Results
We test our IGAbased pose estimation method on three image sequences, including one straight walk sequence [24], one turning walk sequence [23], and one run sequence [21]. The purpose is to test the capability of the method to cope with limb occlusion, leftright ambiguity, and viewpoint problems, which are the main challenges that a pose estimation method has to deal with. As mentioned in Section 3, we first learn the subspace of walking and running. To extract the motion subspace of walking, a data set consisting of motion capture data of a single subject was used. The total number of frame is 316. For running subspace learning, a data set with 186 frames was used. It was found that the different subjects and different frame numbers can produce generally identical subspace. So the learned subspaces are also used in the tracking experiments.
For pose estimation on a single image, the parameters of IGA are set as and to deal with the tradeoff of computational time and accuracy. We test our IGAbased pose estimation method on 100 frames of images for all three types of motions, and the mean errors of joint angle are reported, which are shown in Figure 10. From Figure 10 we can see that, except for some particular joints, the mean errors of most joints for three sequences are less than 5 degrees. The mean errors of some joint angles are larger than others because they have wider range of variation or less observability related to 2D image features. Our results are competitive with others reported in the related literatures.
(a)
(b)
(c)
Table 2 shows the ground truth and estimated values of some joint angles in an example frame. Three values in each cell are the rotation angles of the joints around ,āā,āāandāā axes, respectively. The values come from a frame on the level of average error. Actually, other frames show generally the similar comparison results. From Table 2 we can see that estimated joint angles are close to the ground truth data. The experiment results demonstrate that our IGAbased pose estimation method is effective to analyze articulated human pose from a single image.

The results on real images are shown in Figure 11. From the above experiment results, we can see that, on most of the frames, the occlusion and leftright confusion problems are tackled by searching the optimal pose in the extracted subspace because the prior knowledge about motions is contained in this subspace. And the pose estimator is view invariant, mainly because of the viewpointindependent manifold learning and special step for searching the global motion. In addition, the experiment results on walking and running sequence demonstrate that our algorithm is efficient for different types of motions. Actually, our method can be generalized to any other types of motions as long as the corresponding subspace can be properly extracted from training data.
6.4. SIGABased Pose Tracking Results
We demonstrate our tracking algorithm on walking and running image sequences. And then we compare SIGA quantitatively with other tracking methods and include particle filter (PF) method [6], particle swarm optimization (PSO) [11], and pose tracking in linear subspace using annealing genetic algorithm (PCA + GA) [5].
As suggested in [6], for a human model with DOF between 6 and 12, PF needs about 1000 particle to run. And in [17], PF used 4000 particles for a 29 DOF human model. While in [11], 7200 particles are used for a 31 DOF human model. In this paper, the human model in the original space is with 66 DOF; we set the particles size to be 12000 for PF. While in IGA, the quantitative results of experiments show that IGA with 40 antibodies yields results, under similar testing conditions, more accurate than PF available to us. For motion tracking, the iteration time is set to be . Thus, the number of likelihood evaluations for a single image would be 800 at most, which is much less than 4000 for GA [5] (size of population is 100, iteration time is 40), 7200 for PSO [13], and 12000 for PF.
We first use IGAbased pose estimation method to analyze human pose on the first image of the video for initialization, where the parameters are set as , for careful search of the state space in initialization. While on the following frames, we set the iteration times to be . It is mainly because our nextframe propagation strategy can produce a compact antibodies population for optimization. And in our experiment, we set for straight walking sequences and for running sequences.
The mean errors of different methods over all joint angles of the test sequences are shown in Figure 12. And Table 3 is the statistics of the average errors and the standard deviations. From Figure 12 and Table 3, we can see that our method achieve better results. The average errors and the standard deviations over all frames are near 3Ā° and 1Ā°, respectively, in general. It also can be found that the change of mean error of our method in whole sequence is small, which indicates that our method can achieve stable tracking of 3D human pose.

(a)
(b)
Figure 13 is the tracking results on walking and running image sequences, respectively. From the above experimental results we can see that our IGAbased pose estimation method can successfully be used for initialization of tracking. Acutely, our IGAbased pose estimation method is also used for initialization of PF in our experiments. Experimental results on different types of motion sequence show that SIGA has good performance even without any learnt constant motion models, which demonstrate our nextframe propagation strategy is effective to generate initial distribution of antibodies for the next frame.
(a)
(b)
(c)
Experimental results demonstrate that our SIGAbased tracking method can achieve accurate and stable tracking of 3D human motion. However, our method has some drawbacks as discussed below. Firstly, though pose optimization in the latent space makes our method more effective and accurate, it makes our method not suitable for more complicated motion analysis. So in our future work, we will extend our algorithm to cover a wider class of human motions and explore switch mechanism between different subspaces. Secondly, in generative tracking approaches, the time taken by an algorithm depends mostly on the number of likelihood evaluations. In our IGA pose optimization method, the time complexity would be , which makes it cannot work for real time applications. In addition, our method is dependent on the silhouette detection from video. But human silhouette detection from video is difficult especially in uncontrolled environment. More robust human silhouette detection method and more sophisticated image likelihood function will be considered in our future work.
Recently, Gaussian Process Latent Variable Models (GPLVM [25]) has been another widely studied latent space learning method for human motion tracking. Compared with manifold learning method (ISOMAP), GPLVM could building the inverse mapping easily. However, GPLVM cannot work well on small training dataset and highdimensional data. So in our future work, we will study how to apply GPLVM for motion tracking effectively. And more, studies on motion tracking using evolutional computing methods are still limited. In our future work, we will consider to apply other evolutional computing methods for motion tracking.
7. Conclusions
We presented a novel generative approach to reconstruct 3D human pose from a single monocular image as well as from monocular image sequences. The main contribution is to optimize human pose in learnt latent human motion space. Pose analysis in the latent space learnt using ISOMAP happens to be more efficient and accurate. In the search strategy, we apply the immune genetic algorithm for pose estimation. A sequential IGA framework is proposed for pose tracking by incorporating the temporal continuity information into the traditional IGA. Compared with GA and PSO, IGA has the ability to use the prior knowledge of human motion. Experiment results on different motion types and image sequences demonstrated that our IGAbased method for pose estimation is effective to deal with occlusion, leftright ambiguity, and the viewpoint problem. The sequential IGA method can achieve stable and accurate pose tracking. Quantitative experiments compared with other stateofart methods show that our methods achieve better results.
In the future work, we will extend our algorithm to cover a wider class of human motions and explore switch mechanism between different subspaces. In addition, we will also consider more sophisticated image likelihood and how to reduce the computation time. How to apply other evolutionary computation methods for human motion tracking will also be considered in our future work.
Acknowledgments
This work is supported by the National High Technology Research and Development Program of China (2007AA01Z334), National Natural Science Foundation of China (61272219, 61021062, and 61100110), Program for New Century Excellent Talents in University of China (NCET0404605), Natural Science Foundation of Jiangsu Province (BK2010375), and Key Technology R&D Program of Jiangsu Province (BY2012190, BE2010072, and BE2011058).
References
 C. Sminchisescu, ā3D human motion analysis in monocular video, techniques and challenges,ā in Human Motion Understanding, Modeling, Capture and Animation, R. Kleete, D. Metaxas, and B. Rosenhahn, Eds., Springer, New York, NY, USA, 2007. View at: Google Scholar
 G. Mori and J. Malik, āRecovering 3D human body configurations using shape contexts,ā IEEE Transactions on Pattern Analysis and Machine Intelligence, vol. 28, no. 7, pp. 1052ā1062, 2006. View at: Publisher Site  Google Scholar
 A. Agarwal and B. Triggs, āRecovering 3D human pose from monocular images,ā IEEE Transactions on Pattern Analysis and Machine Intelligence, vol. 28, no. 1, pp. 44ā58, 2006. View at: Publisher Site  Google Scholar
 A. Agarwal and B. Triggs, āMonocular human motion capture with a mixture of regressors,ā in Proceedings of the IEEE Computer Society Conference on Computer Vision and Pattern Recognition (CVPRW '05), p. 72, San Diego, Calif, USA, June 2005. View at: Google Scholar
 X. Zhao and Y. Liu, āGenerative tracking of 3D human motion by hierarchical annealed genetic algorithm,ā Pattern Recognition, vol. 41, no. 8, pp. 2470ā2483, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 M. Isard and A. Blake, āCONDENSATION—conditional density propagation for visual tracking,ā International Journal of Computer Vision, vol. 29, no. 1, pp. 5ā28, 1998. View at: Publisher Site  Google Scholar
 L. Raskin, M. Rudzsky, and E. Rivlin, āDimensionality reduction using a Gaussian process annealed particle filter for tracking and classification of articulated body motions,ā Computer Vision and Image Understanding, vol. 115, no. 4, pp. 503ā519, 2011. View at: Publisher Site  Google Scholar
 R. Poppe, āVisionbased human motion analysis: an overview,ā Computer Vision and Image Understanding, vol. 108, no. 12, pp. 4ā18, 2007. View at: Publisher Site  Google Scholar
 T. B. Moeslund, A. Hilton, and V. Krüger, āA survey of advances in visionbased human motion capture and analysis,ā Computer Vision and Image Understanding, vol. 104, no. 23, pp. 90ā126, 2006. View at: Publisher Site  Google Scholar
 C. S. Lee and A. Elgammal, āCoupled visual and kinematic manifold models for tracking,ā International Journal of Computer Vision, vol. 87, no. 12, pp. 118ā139, 2010. View at: Publisher Site  Google Scholar
 V. John, E. Trucco, and S. Ivekovic, āMarkerless human articulated tracking using hierarchical particle swarm optimisation,ā Image and Vision Computing, vol. 28, no. 11, pp. 1530ā1547, 2010. View at: Publisher Site  Google Scholar
 R. Urtasun, D. J. Fleet, and P. Fua, āMonocular 3D tracking of the golf swing,ā in Proceedings of the IEEE Computer Society Conference on Computer Vision and Pattern Recognition (CVPR '05), pp. 932ā938, San Diego, Calif, USA, June 2005. View at: Google Scholar
 C. Sminchisescu and A. Jepson, āGenerative modeling for continuous nonlinearly embedded visual inference,ā in Proceedings of the 21st International Conference on Machine Learning (ICML '04), pp. 759ā766, Alberta, Canada, July 2004. View at: Google Scholar
 A. Elgammal and C. S. Lee, āInferring 3D body pose from silhouettes using activity manifold learning,ā in Proceedings of the IEEE Computer Society Conference on Computer Vision and Pattern Recognition (CVPR '04), pp. 681ā688, Washington, DC, USA, July 2004. View at: Google Scholar
 J. B. Tenenbaum, V. de Silva, and J. C. Langford, āA global geometric framework for nonlinear dimensionality reduction,ā Science, vol. 290, no. 5500, pp. 2319ā2323, 2000. View at: Publisher Site  Google Scholar
 N. R. Howe, M. E. Leventon, and W. T. Freeman, āBayesian Reconstruction of 3D human motion from singlecamera video,ā in Proceedings of Neural Information Processing Systems (NIPS '00), pp. 820ā826, Denver, Colo, USA, 2000. View at: Google Scholar
 J. Deutscher, A. Blake, and I. Reid, āArticulated body motion capture by annealed particle filtering,ā in Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition (CVPR '00), pp. 126ā133, Hilton Head, SC, USA, June 2000. View at: Google Scholar
 T. Krzeszowski, B. Kwolek, and K. Wojciechowski, āArticulated body motion tracking by combined particles swarm optimization and particle filtering,ā in Proceedings of the International Conference on Computer Vision and Graphics (ICCVG '10), Part I, vol. 6374 of Lecture Notes in Computer Science, pp. 147ā154, 2010. View at: Google Scholar
 X. Zhang, W. Hu, S. Maybank, and L. Xi, āSequential particle swarm optimization for visual tracking,ā in Proceedings of the IEEE Computer Society Conference on Computer Vision and Pattern Recognition (CVPR '08), pp. 23ā28, Anchorage, Alaska, USA, 2008. View at: Google Scholar
 W. Lei and J. Licheng, āImmune evolutionary algorithms,ā in Proceedings of the International Conference of Signal Process (ICSP '00), pp. 1655ā1662, Beijing, China, August 2000. View at: Google Scholar
 CMU Database, http://mocap.cs.cmu.edu/.
 M. Gong, L. Jiao, and L. Zhang, āBaldwinian learning in clonal selection algorithm for optimization,ā Information Sciences, vol. 180, no. 8, pp. 1218ā1236, 2010. View at: Publisher Site  Google Scholar
 L. Sigal, A. O. Balan, and M. J. Black, āHumanEva: synchronized video and motion capture dataset and baseline algorithm for evaluation of articulated human motion,ā International Journal of Computer Vision, vol. 87, no. 12, pp. 4ā27, 2010. View at: Publisher Site  Google Scholar
 D. Ormoneit, H. Sidenbladh, M. J. Black, and T. Hastie, āLearning and tracking cyclic human motion,ā in Proceedings of Neural Information Processing Systems (NIPS '01), pp. 894ā900, Vancouver, Canada, December 2001. View at: Google Scholar
 R. Urtasun, D. J. Fleet, A. Hertzmann, and P. Fua, āPriors for people tracking from small training sets,ā in Proceedings of the 10th IEEE International Conference on Computer Vision (ICCV '05), pp. 403ā410, Beijing, China, October 2005. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2013 Yi Li and Zhengxing Sun. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.