1 Introduction
Spatiotemporal patterns in highdimensional motion data are crucial in many recognition applications. For example, human action is the combination of the body joint movements over a time interval. Facial expression is the result of the facial landmark movements (Figure 1). Understanding the movement trajectories and modeling the underlying spatiotemporal patterns play an important role in recognizing these actions, especially with the recent emergences of reliable algorithms [1, 2] to estimate the positions of body joints and facial landmarks.
In this work, we are interested in developing a probabilistic model to capture the spatiotemporal pattenrs in highdimensional time series for classification purpose. Many recent works develop novel models to capture the spatiotemporal dynamics [3, 4, 5, 6]
. However, most models, such as hidden Markov model (HMM)
[5], dynamic Bayesian network (DBN)
[7] and conditional random field (CRF) [8] are timeslice local models, which assume Markov property and stationary transitions and hence can only capture local dynamics. The local models suffer from two limitations. First, local dynamics may not represent a sequence well because it fails to mpdel the overall dynamics. Second, the stationary transition assumption may not hold for many realworld applications.Compared with timesliced dynamic models, restricted Boltzmann machines (RBMs) has shown strong capability of modeling joint distributions and therefore can capture the global patterns. In literature, RBMs have been successfully applied to separately capture the spatial
[9] or temporal [10] patterns in different types of data. In this work, we propose a variant of RBM that can capture spatial and global temporal patterns simultaneously to comprehensively model the highdimensional motion data.In a typical RBM, since there are no lateral connections among nodes in each layer, input data are independent of each other given the states of hidden layer. This assumption limits RBM’s representation power, since there exist direct dependencies among input data. There are generally two types of data interactions: interactions through latent variable and direct interactions independent of latent variable. For example, as stated in [11], soldiers on a parade follow the commander’s order to march in some direction, but they “form a neat rectangle by interacting with their neighbors”. This example illustrates that soldiers’ behavior are determined by both the commander’s order (latent) and the interactions with their neighbors. In this case, RBM is not effective in modeling the local interaction that is independent of the latent variable. The use of RBM for data representation and classification is further hindered by the difficulty in comparing different RBMs due to the intractable computation of the partition functions.
Allowing interaction among visible units can overcome this shortcoming. We introduce restricted Boltzmann machine with local interactions (LRBM) to capture both the global temporal patterns and local spatial interactions in the input data. Specifically, we add a new pairwise potential term in the learning objective function of RBM to capture the local spatial interactions among components of an input vector. To perform efficient learning for the model, we replace the reconstruction procedure in the conventional Contrastive Divergence
[13] algorithm with a mean field approach.For classification task, RBM is typically used for learning features, as the input to a second stage classifier (e.g., SVM
[14]). Typically one model is trained for all classes. To obtain good features, deep structure is built and back propagation is employed to carefully tune the parameters. In this work, we use RBM as a generative model to capture the spatiotemporal patterns in the data. RBM is used for data representation instead of feature learning. Given an observation, the only output we get from an RBM model is the likelihood of the model. For the recognition purpose, one model is trained for one class of input data. To compare among different models, we employ a method to estimate the relative partition function of a pair of RBMs for binary classification, and a label ranking method is used to extend the binary classification to multiclass classification.To evaluate the performance of LRBM, we apply it to two areas related to complex spatial and temporal patterns: facial expression recognition and human action recognition. Experimental results on benchmark databases demonstrate the effectiveness of the proposed model.
The rest of the paper is structured as follows. Section 2 presents an overview of the related work. Section 3 introduces the LRBM to model the global temporal dynamics and spatial patterns of motion data, as well as the classification method. We will then give the experimental results in 4. The paper is concluded in section 5.
2 Related Work
Capturing and representing spatiotemporal structure in data is important for many recognition and classification tasks. Research for capturing such patterns can be categorized into featurebased and modelbased methods. The most widely used spatiotemporal features include spatiotemporal interest point (STIP) based features [3, 4] and optical flow based features [15]. These features capture local appearance or motion patterns near the interest points or optical flows. Although having been successfully applied to many applications, these features generally focus more on local patterns.
Modelbased methods include probabilistic graphical models such as Hidden Markov Models [5], Dynamic Bayesian Networks [16], Conditional Random Fields [17], and their variants. While capable of simultaneously capturing both spatial and temporal interactions, they can only capture the local spatial and temporal interactions due to the underlying Markov assumption.
Restricted Boltzmann machines (RBMs) have been separately used for modeling spatial correlation or temporal correlation in the data in the last decade. RBM was firstly introduced to learn deep features from handwritings to recognize digits
[18]. In [9], Eslami et al. propose a Deep Belief Network to model the shapes of horses and motorbikes. The samples from the model look realistic and have a good generalization. A more complicated model, proposed by Nair and Hinton
[19], considers the spatial correlation among visible layer using a factored 3way RBM, in which triple potentials are used to model the correlations among pixels in natural images. The intuition is that in natural images, the intensity of each pixel is approximately the average of its neighbors. Wu et al. [20] apply the 3way RBM to facial landmark tracking, and model the relationship between posed faces and frontal faces, under varying facial expressions.For dynamic data modeling, Taylor et al. [21] use a Conditional RBM (CRBM) to model the temporal transitions in human body movements, and reconstruct body movements. Nevertheless, like HMM, CRBM still models local dynamics by assuming ’th order Markov property.
The idea of using RBM to model global pattern is not new. In [22, 23], RBM and CRF are combined for face labeling problem in a simgle image or video sequences. Compared with these works, our work is different for several reasons. First, our goal is to use RBM for data representation for multiclass classification, while the goal of [22, 23] is MAP inference, which is to recover the label for each superpixel. Second, we extend the conventional RBM to capture the local shape and global temporal patterns in a unified model, while in [22, 23] RBM is built on top of the hidden layer of the CRF, as the prior for the labels. Third, our method learns all the parameters simultaneously, while their method performs learning separately.
RBM and its variants have also been used for modeling motion data. For example, Sutskever and Hinton [24] introduce a temporal RBM to model highdimensional sequences. Wang et al. [10] use RBM to capture the global dynamics of finger trace. However, it is limited to model the global dynamics for 1D data only.
To improve the representation power of RBM, semirestricted Boltzmann machine (SRBM) [11] is introduced to model the lateral interactions between visible variables. The main property of SRBM is that given the hidden variables, the visible layer forms a Markov random field. However, for highdimensional motion data, there will be too many parameters if every pair of visible units has an interaction. In this work, we model the dynamic nature of data with fewer parameters than a SRBM, which makes the learning process more efficient.
Besides feature extraction and shape modeling, RBMs have also been used for classification. Larochelle and Bengio
[25] introduce a discriminative RBM as a classifier by including the labels in visible layer, and make predictions by comparing the likelihood of each label vector. In [26] discriminate RBM is introduced to model vector inputs by duplicating the discriminative RBM and adding constraints on the hidden layer. In all RBM related models discussed above, one RBM is trained for all classes. For this work, in contrast, we build one RBM for each class, and perform a multiclass classification task.In this research, we propose to extend the standard RBM to simultaneously capture the spatial and global temporal patterns in highdimensional sequential data and employ such RBM model to different classification tasks in computer vision.
3 Data SpatioTemporal Pattern Modeling
In this section, we will firstly give a brief introduction of restricted Boltzmann machine, and then introduce the proposed RBM with Local Interaction (LRBM) to model multidimensional motion data and its learning method. Finally we present the method to use LRBMs as a set of pairwise classifier to perform multiclass recognition.
3.1 Restricted Boltzmann Machine
A standard RBM is a generative model with two densely connected layers, one visible layer to represent data and one latent layer to extract stochastic binary features from data (Figure 2
). Hidden units are connected to visible nodes using symmetrically weighted connections to model their joint distribution. In our work, we use GaussianBinary RBM, where the hidden units are binary and the visible variables are assumed to follow normal distribution.
The energy function for each pair is parameterized in Equation 1.
(1) 
where is the bias for visible unit ,
is the standard deviation of the Gaussian distribution, which is typically 1 if we normalize the data,
is the bias for the hidden unit , is the weight of the link connecting and .For every possible pair of visible and hidden vector, the network assigns a probability:
(2) 
where is the partition function. With continuous inputs, the is calculated by integrating over all visible nodes and summing over all hidden unit configurations.
3.2 Proposed Model
Using standard RBM to model highdimensional motion data has its limitations. First, the interaction among the data is represented through latent variables, which can be easily represented in direct connections. Second, if a single RBM models the whole sequence, the spatial information in a time slice is treated the same as the temporal information of one dimension. If one RBM models only one time slice (as in [21]), the temporal information remains local.
In this work, we propose the restricted Boltzmann machine with local interaction (LRBM, Figure 3) to overcome such limitations. For sequential data , by allowing local interaction, we have the following energy function,
(3) 
where is a dimension vector representing input vector at time slice , has the dimension of , is the weight vector connecting to a hidden node . and have the same meaning as in Equation 1. is a symmetric matrix with zeros on diagonal, modeling the correlation of each vector .
The proposed energy function models two kinds of data interactions: interaction through latent variables and interaction directly among data. In highdimensional motion data, components in a single time slice (spatial information) is better to be considered differently from components along the timeline (temporal information). Interactions between input data and hidden variables are through weight , which models global pattenrs, because every visible unit is connected to every latent unit. Matrix models direct interactions among input data, as in Figure 3, representing local spatial patterns. This kind of interaction directly affects visible layer without going through latent layer, so it is more effective to model some local spatial patterns. The elements in are pairwise potentials of features in one time slice. Different shapes (spatial pattern) of input will have different contributions to the energy function. Thus, captures spatial patterns among elements of input vector. With temporal information captured via the hidden nodes, our work can capture global spatiotemporal dynamics.
We assume the spatial relationship are constant throughout the whole sequence, so the parameters in are shared in different frames, which means is invariant of . Under this invariance assumption of , the number of parameters is significantly reduced.
From the energy function (Equation 3) and likelihood function (Equation 2), we can derive the probability of an hidden unit to be activated, given an visible layer:
(4) 
As one hidden unit is connected to all visible units, whether it is activated or not depends stochastically on the visible layer. A hidden unit is activated through Equation 4 when it detects some specific pattern in the visible layer. The pattern is captured by the weights connecting each element in to . Therefore the hidden layer is able to capture important global patterns of . Hence, the proposed model can simultaneously capture the global temporal patterns (through ) and local shape patterns (through ).
LRBM is similar to the factored 3way RBM [19] and SRBM [11] with respect to modeling the interactions in the visible layer. However, there exist some significant differences. First, in factored 3way RBM, every pair of visible units has a potential to represent the interaction, while in LRBM, units only interact with neighbors in the time slice. Second, if the number of factors in 3way RBM or the data dimension in SRBM is high, the overall parameters are much more than LRBM. In LRBM, the additional parameters are from matrix , with the parameters, where is the dimension of the vector . This would greatly reduce the computational load. Finally, both 3way RBM and semiRBM are proposed to model the spatial patterns of a single natural image, while LRBM models spatiotemporal patterns of sequential data.
The joint probability of and are the same as in Equation 2. The likelihood of input data is computed by summing over all configurations of hidden units:
(5) 
3.3 Model Learning
In our work, one LRBM is trained for one class. The parameters include bias for visible and hidden units and , weight between two layers , and local interaction matrix . To learn the parameters, we seek to maximize the joint probability of all training data of a class, . Assume the data has been normalized, so can be removed from the energy function. The derivative of the loglikelihood of a training instance with respect to a parameter is given below:
(6)  
(7)  
(8)  
where the angle brackets are used to denote expectations under the distribution specified by the subscript that follows. is the element of at position . are the and components of vector . This leads to a simple learning rule: stochastic steepest ascent in the loglikelihood of training data. Take as an example:
(9) 
where is the learning rate. is easy to get from the data. is much more difficult to compute due to the large dimension of hidden and visible layers. Hinton [27] proposes the CD algorithm to approximate by using a reconstructed sample, which gives the following weight change:
(10) 
Given the visible layer, the hidden units are independent with each other, so sampling the states of hidden units can be performed in parallel using Equation 4.
Given the states of the hidden units, the visible units form a Markov Random Field in which the pairwise interaction weight between and is . They are no longer independent, so sampling cannot be done in parallel. However, we can obtain the conditional probability of each node by fixing its neighbors. The mean field algorithm is used to sample data from the model, and to reconstruct visible layer when learning the parameters of LRBM. Specifically, for each vector , each component is sampled by fixing all its neighbors. Once a component is sampled, the vector is updated with the newly sampled component. This procedure is repeated until all components have been updated, and thus we get a reconstructed vector. The conditional probability of one visible node is given in Equation 11:
(11) 
where are the neighbors of ,
denotes the Gaussian probability density function with mean
and variance
. The superscript means the component of a vector.Compared with distribution of visible units in standard RBM:
(12) 
the mean of the distribution of a visible node is similar in the first two terms, except that in LRBM, it is modified by the pairwise potentials relating one node to all its neighbors.
With the reconstructed data, we are able to calculate the approximate gradient of each parameter, and perform the CD learning procedure. As mentioned in [28]
, the RBMs learn better if more steps of reconstruction are used before collecting the statistics. In practice, we sample 10 times using the mean field method before collecting the reconstructed data in each epoch.
Due to the nonconvex property of the objective function, different parameter initializations in RBM learning can end up with different models. For each class, several candidate models are learned from different initializations, and we select the one that can best discriminate current class from the others.
For example, after learning model for class , given instances from all classes , we expect instance to have greater likelihood on model than all other instances (). In practice, we can sort the likelihoods of instances from all classes, and choose the candidate model that minimizes the rank of the instances from class .
Notice that the likelihood is calculated according to Equation 5, given a single LRBM model, so for the comparison of the likelihoods, we don’t need to estimate the intractable partition function, since it is merely a constant.
3.4 Multiclass Classification
For this work, RBM is used as a generative model for classification. A common strategy for training generative model for classification is to train multiple models for different classes and then evaluate the likelihood of each model during testing. In particular, given a set of properly learned LRBM’s for classes, the basic idea for classification is to find the model that generates the largest likelihood given an instance of data .
(13) 
where is the prediction of our classifier. The likelihood of a data instance is
(14) 
For brevity, we denote as
(15) 
can be computed directly. However the partition functions are intractable for RBMs with large numbers of hidden and visible units. One possible solution is the Annealed Importance Sampling [29] to estimate the partition functions, and directly compare the likelihood. However, such estimation needs too many samples to obtain a good estimation of the partition function. Instead, for binary classification, Schmah et al. [30] propose a method to discriminatively estimate the difference of logpartition functions of two RBMs:
(16) 
We extend this method to multiclass classification. Ideally, . But since the partition function is not directly computed, there is a slight error in the estimation. To address this issue, we perform a label ranking process [31] to make the final decision using a set of pairwise classifiers.
In the simplest case, with all the pairwise classifiers , each prediction is interpreted as a vote for a class, and the class with the highest votes is proposed as a prediction. Alternatively, in confidence estimation, instead of a binary result {0, 1}, a “soft” classifier is employed to map the difference in likelihood into the unit interval [0, 1]:
(17) 
Parameter is searched in a range of to maximize the training accuracy.
The output of such “soft” binary classifier can be interpreted as a confidence value in the classification: the closer the output to 1, the stronger the decision of choosing class is supported. Notice that is not symmetrical. It only gives the preference toward between and . The preference toward is . A valued preference relation matrix is defined for any query instance :
(18) 
In our approach, we evaluate the confidence score by summing up all the confidence value
(19) 
where the index goes through to , meaning the confidence score for one instance on different classes. The label of the model that associates with the highest score is proposed as the final decision.
(20) 
In general, if we have a class classification problem, pairwise classifiers are used to compute the confidence score.
4 Experiments
We evaluate our algorithm on two different areas that involve highdimensional motion data: facial expression recognition and human action recognition. Two benchmark data sets are used in our experiment: the extended CohnKanade data set (CK+) [2] and G3D data set [12]. We will also compare the proposed LRBM model with other related methods.
4.1 Facial Expression Recognition
The Extended CohnKanade data set (CK+) is a complete data set for action units and emotionspecified expressions. In this work, based on a sequence of landmarks on the human face, our goal is to classify it into one of seven expressions: angry, disgust, fear, happy, sadness, surprise and contempt. CK+ data set includes 593 sequences from 123 subjects. From all the sequences, 327 are associated with expression labels. In our experiment, these 327 sequences are used. Landmarks are the positions of 68 facial points from a detection and tracking procedure, which are provided by the database.
As the pose of each subject varies slightly over time, we apply a pose rectification procedure to make each face a frontal one. Then, using the detected eyes and the interocular distance, we perform a geometric normalization by making the size of facial area the same throughout all subjects, and moving the centers of eyes to the same position. A smoothing filter is also used to reduce the noise in the trajectory. To reduce the feature dimension while keeping the useful information, as shown in Figure 4
, we omit the outline landmarks and some other points, and select the more informative points, resulting in a 15point feature for each frame, which is a 30d vector. The features selected are reasonable due to the property of facial symmetry.
An  Di  Fe  Ha  Sa  Su  Co  

An  97.8  0.0  0.0  0.0  0.0  2.2  0.0 
Di  8.5  89.8  0.0  1.7  0.0  0.0  0.0 
Fe  0.0  0.0  84.0  8.0  0.0  8.0  0.0 
Ha  0.0  0.0  0.0  100.0  0.0  0.0  0.0 
Sa  10.7  3.6  3.6  0.0  78.6  3.6  0.0 
Su  0.0  0.0  0.0  0.0  1.2  97.6  1.2 
Co  5.6  0.0  0.0  11.1  11.1  0.0  72.2 
4.1.1 Classification Performance
Intuitively, the spatial patterns are crucial in recognizing an expression. From a single face, human can identify the expression without any difficulty. But in order to explore the underlying facial dynamics as additional features for our system, we choose 10 frames from each sequence, representing the neutral look, intermediate expressions and the peak expression. Since our method captures global temporal patterns, a good alignment of different sequences is important to ensure recognition performance. Linear interpolation is used to get sequences with fixed length. The size of hidden layer is set as 400. To compare with baseline method, we use a leaveonesubjectout crossvalidation configuration. Each time we generate the testing data from sequences of one subject. 25 other subjects form the validation set, and all the other subject form the training set.
Lucey et al. [2]  Wang et al. [32]  LRBM  

An  75.0  91.1  97.8 
Di  94.7  94.0  89.8 
Fe  65.2  83.3  84.0 
Ha  100.0  89.8  100.0 
Sa  68.0  76.0  78.6 
Su  96.0  91.3  97.6 
Co  84.4  78.6  72.2 
Avg.  83.3  86.3  88.6 
The performance of the algorithm is given in Table 1. As we only use the shape features (i.e. facial landmark positions), the comparison is only between methods using the same features. Our method’s hit rates for each expression are: Angry  97.8%, Disgust  89.8%, Fear  84.0%, Happy  100.0%, Sadness 78.6%, Surprise  97.6%, Contempt  72.2%. The average accuracy of all the expressions is 88.6%, which is more than 2% better than stateoftheart method, as reported in [32]. In Table 2, we list the comparison of our model with some recent approaches. Our approach reaches the best accuracy on five expressions (angry, fear, happy, sadness and surprise) among all the methods listed. Although our model does not achieve the best accuracy for every expression, it does not fall behind very much on Disgust, but it fails sometimes on Contempt, because this is a subtle expression that our method needs more features to represent it accurately. In CK+ data set, contempt expression has only 18 samples. However, the average performance of the proposed method has been improved a lot. Another thing to mention is that the method in [2] uses both shape features and appearance features, while we only use the shape features.
If we do not consider the spatial interaction within each frame, the data is simply aligned as a long vector to feed in the conventional RBM. Then the performance is 86.3%, which is less than the LRBM. This also proves the improvement of LRBM over RBM.
We also compare our method with 4 best available results to date, including Timeseries Kernels (Lorincz et al. [33]
), spatiotemporal independent component analysis (Long et al.
[34]), boosted dynamic features (Yang et al. [35]), and nonnegative matrix factorization techniques (Jeni et al. [36]). These methods are different from the baseline method in [2], because contempt emotion is removed from the data set, and binary classification is performed for each expression, so for the classification of one expression, the other five expressions are negative samples. Following the same experiment settings, the result of classification is shown in Figure 5, and detailed comparison is in Table 3. The classification performance of our algorithm is close to, or even better than stateoftheart method, with a better average accuracy. This demonstrate the effectiveness of the proposed LRBM model as both a binary classifier and a multiclass classifier.Method  An  Di  Fe  Ha  Sa  Su  Avg. 

Yang et al. [35]  0.973  0.941  0.916  0.991  0.978  0.998  0.966 
Long et al. [34]  0.933  0.988  0.964  0.993  0.991  0.999  0.978 
Jeni et al. [36]  0.989  0.998  0.977  0.998  0.994  0.994  0.992 
Lorincz et al. [33]  0.991  0.994  0.987  0.999  0.995  0.996  0.994 
LRBM  0.992  0.995  0.995  0.999  0.997  0.995  0.996 
4.1.2 Comparison with Feature Learning Method
To comprehensively evaluate our method, we compare with a feature learning method, since RBMs are typically used for feature learning. A single RBM is trained on all classes of sequences using contrastive divergence. The size of hidden layer is the same as in the generative model. In the test process, given each observation sequence V
, we compute the posterior probability
for each hidden variable according to Equation. 4, as the feature representation. This is reasonable because is the probability of a pattern being activated. A linear multiclass SVM is trained on such features for classification.The confusion matrix of this method is given in Table 4. For most expressions, the featurebased method is not as good as the generative model. The average recognition accuracy is 83.8%, which is approximately 5% below the performance of the generative model. This is reasonable, because the features learned from the model cannot capture the spatial relationship in each frame, which is specifically modeled in our model using the matrix . If we add another latent layer, the classification performance only increases marginally by about 1%.
An  Di  Fe  Ha  Sa  Su  Co  

An  88.9  11.1  0.0  0.0  0.0  0.0  0.0 
Di  11.9  86.4  0.0  1.7  0.0  0.0  0.0 
Fe  0.0  4.0  36.0  12.0  0.0  48.0  0.0 
Ha  0.0  0.0  1.4  97.1  0.0  0.0  1.4 
Sa  17.9  0.0  3.6  0.0  57.1  10.7  10.7 
Su  0.0  2.4  0.0  0.0  0.0  97.6  0.0 
Co  5.6  5.6  5.6  27.8  5.6  5.6  44.4 
4.2 Human Action Recognition
G3D data set is an action data set containing a range of gaming actions captured by Microsoft Kinect. The data set contains 10 subjects performing 20 gaming actions: punch right, punch left, kick right, kick left, defend, golf swing, tennis swing forehand, tennis swing backhand, tennis serve, throw bowling ball, aim and fire gun, walk, run, jump, climb, crouch, steer a car, wave, flap and clap. Synchronized video, depth and skeleton data are available in this data set. We only pick the skeleton data. To reduce the data dimension but keep the useful information, we use the 3D location information of four dominant joints (i.e. two hands and two feet). However the proposed approach can be applied to modeling more joints if we have enough training data.
Before abstracting the features, we perform a normalization procedure to minimize the effect brought by difference in subjects’ body shapes. Specifically, we obtain the average bone lengths from all subjects in the training data, and then change the skeleton tracking results in every frame for each subject with the average shape. Therefore, every subject has the same body shape, and the cross subject distinction is alleviated.
As the size of the visible layer of an RBM is fixed, linear interpolation is performed to convert all sequences into the same length (20 frames for each sequence). The 3D positions of the body joints along all three dimensions (x, y and z) form the 240dimension input for the LRBM model. In the training phase, we set the size of hidden layer as 80.
We use the action segmentation that is provided by the data and the same experiment configuration as in [12]. The data set is split by subjects where the first 4 subjects were used for training, 1 for validation and the remaining 5 subjects for testing.
4.2.1 Classification Performance
The confusion matrix is given in Figure 6. The overall accuracy of LRBM is 90.5%. Our model encounters some failures for the actions of ThrowBowlingBall in Bowling category and Jump. In these actions, occasionally the body parts are occluded from the Kinect sensor, which will lead to estimated positions of joints. Then the positions of the two feet are significantly corrupted. As our method depends fully on the trajectories of joints, corrupted tracking results have a huge influence on the performance of LRBM. If we remove the pairwise potential in the visible layer and make the model a standard RBM, the accuracy drops to 84%, which proves the better performance of LRBM than RBM, as expected, due to modeling the spatial interactions. Comparison with the baseline model [12] is given in Table 5. Other than Bowling, performance on golf action is quite close to the baseline method, and on all the other actions, our approach outperforms the baseline method. The overall accuracy is increased by 15.5% in terms of F1 score.
One reason of the improvement is that the method in [12] is based on 3 frames, hence can only capture local dynamics, while the proposed method is sequence based that can handle global dynamics.
Action  Bloom et al. [12]  LRBM 

Fighting  70.46  97.09 
Golf  83.37  81.82 
Tennis  56.44  84.38 
Bowling  80.78  61.54 
FPS  53.57  95.52 
Driving  84.24  100.00 
Misc  78.21  95.24 
Avg.  72.44  87.94 
To further demonstrate the effectiveness of global dynamic model over local dynamic model, we implement a conditional RBM [21] and a hidden Markov model for action recognition. For the conditional RBM model, each frame depends on 3 previous frames. 20 models are trained for 20 actions. Classification is based on the likelihood of a target sequence on different models. Basically the local dynamic model (CRBM or HMM) models one frame at a time and the global model (RBM or LRBM) models all frames at the same time.
Similar to the expression recognition case, the posterior probability of the hidden variables can be used as the features for classification. We implement the featurebased method using linear SVM as the classifier. Again, it does not perform well compare with other dynamic methods, mainly due to the reason that it cannot capture the local spatial patterns in each frame. Details are given in Table 6.
Method  Accuracy (%) 

HMM  70.3 
1FLSVM  75.9 
2FLSVM  76.7 
CRBM [21]  83.2 
RBM  84.0 
LRBM  90.5 
4.2.2 Handling Noisy and Missing Data
One advantage of generative model is that they can handle noisy or missing data in the input. To demonstrate this point, we design two scenarios with randomly selected noisy or missing data:
(a) Noisy data. For a randomly selected portion of data, we multiply the ground truth data by a Gaussian noise , where is the normal distribution;
(b) Missing data. For a randomly selected portion of data, we assume they are missing, and use the average of their neighbors as the approximation when computing the likelihood.
The data under these scenarios are then fed into our generative model for classification purpose. We vary the portion of noisy or missing data from 0 to 50%, and observe the change of the classification accuracy. The details are given in Figure 7.
In scenario (a), with the portion of noisy data increase, the performance almost decays linearly. In scenario (b), if the percentage of missing data is small (less than 30%), the average of its neighborhood is a decent estimation of the missing value, but with more data missing, the performance decays quite significantly. This is because several consecutive frames are missing at the same time, and the temporal information cannot be recovered effectively. One thing to notice is that with as much as 30% noisy or missing data, our algorithm can achieve 86% accuracy, only 4% below the noiseless situation, which demonstrates the effectiveness of our algorithm to handle noisy or missing data.
4.3 Complexity
In the learning procedure of LRBM, we compute the datadependent expectation and modeldependent expectation using matrix multiplication, so the computational complexity is linear with the size of the weight matrix () and iteration times (). Thus, the computation complexity is . For facial expression recognition task, we set the epoch to 250, and the average running time for training one RBM is 1.77s, compared with 1.30s for standard RBM without local interaction. With 10 candidate models, it takes less than 5 min to train all the models. Experiments are performed on a desktop computer with an Intel i7 3.4GHz CPU and 8GB RAM. For testing, given a target sequence, we need to compute the unnormalized likelihood on each LRBM model using Equation 14, with computational cost . Given the unnormalized likelihood and the relative partition functions, the classification needs comparisons, therefore the overall classification complexity is , which is linear in the dimension of the input sequence. If we linearly increase the length of the sequence by increasing the frame rate, the complexity is also increased linearly.
For multiclass classification with classes, the confidence value matrix is built with complexity . In typical computer vision applications, the number of classes cannot exceed the dimension of the data, which means does not contribute much in the complexity, compared to the first term. To reduce the quadratic term to a linear term , Annealed Importance Sampling [29] can be used to estimate the partition function, thus the likelihood can be directly compared. However, the sampling method requires many samples to get a good estimation of the partition function, and the overall classification complexity is not reduced significantly.
Theoretically, the proposed method can be scaled up to data with thousands of dimensions, but much more training data will be needed since the number of parameters increases linearly with the size of input.
5 Conclusion
In this paper, we study classification problem with multidimensional sequence data. To capture both the global dynamics and the local spatial interactions, we extend the conventional restricted Boltzmann machine by adding pairwise potentials in the energy function. An efficient mean field learning method is proposed to replace the reconstruction procedure in typical Contrastive Divergence learning, to estimate the additional parameters. Also, a novel label ranking approach of estimating the partition function of different LRBMs is presented to perform multiclass classification task. Experiments on two areas that involve multidimensional sequence data are performed to evaluate our approach. Results on two benchmark data sets prove the effectiveness of our algorithm.
Using the generative model for data representation, we do not need to deal with a deep and complicated structure, as in feature learning, and do not need to carefully tune the parameters. This is a major advantage of our method compared with feature learning methods. However, one shortcoming is that is that the data sequences have to be aligned for our global model to function properly, and the classification complexity is quadratic in the number of classes. We will address these issues in future work.
References
References
 [1] J. Shotton, A. Fitzgibbon, M. Cook, T. Sharp, M. Finocchio, R. Moore, A. Kipman, A. Blake, Realtime human pose recognition in parts from single depth images, in: CVPR, IEEE, 2011, pp. 1297–1304.
 [2] P. Lucey, J. F. Cohn, T. Kanade, J. Saragih, Z. Ambadar, I. Matthews, The extended cohnkanade dataset (ck+): A complete dataset for action unit and emotionspecified expression, in: CVPR Workshops, IEEE, 2010, pp. 94–101.
 [3] I. Laptev, On spacetime interest points, International Journal of Computer Vision 64 (23) (2005) 107–123.
 [4] I. Laptev, M. Marszalek, C. Schmid, B. Rozenfeld, Learning realistic human actions from movies, in: CVPR, 2008, pp. 1–8.
 [5] F. Lv, R. Nevatia, Recognition and segmentation of 3d human action using hmm and multiclass adaboost, ECCV (2006) 359–372.
 [6] J. Wang, Z. Liu, Y. Wu, J. Yuan, Mining actionlet ensemble for action recognition with depth cameras, in: CVPR, IEEE, 2012, pp. 1290–1297.
 [7] I. Dagli, M. Brost, G. Breuel, Action recognition and prediction for driver assistance systems using dynamic belief networks, in: Agent Technologies, Infrastructures, Tools, and Applications for EServices, Springer, 2003, pp. 179–194.
 [8] C. Sminchisescu, A. Kanaujia, D. Metaxas, Conditional models for contextual human motion recognition, Computer Vision and Image Understanding 104 (2) (2006) 210–220.
 [9] S. A. Eslami, N. Heess, J. Winn, The shape boltzmann machine: a strong model of object shape, in: CVPR, IEEE, 2012, pp. 406–413.
 [10] Z. Wang, G. Schalk, Q. Ji, Anatomically constrained decoding of finger flexion from electrocorticographic signals, in: NIPS, 2011.
 [11] S. Osindero, G. E. Hinton, Modeling image patches with a directed hierarchy of markov random fields, in: NIPS, 2007, pp. 1121–1128.
 [12] V. Bloom, D. Makris, V. Argyriou, G3d: A gaming action dataset and real time action recognition evaluation framework, in: CVPR Workshops, IEEE, 2012, pp. 7–12.

[13]
M. CarreiraPerpinan, G. Hinton, On contrastive divergence learning, in: Artificial Intelligence and Statistics, Vol. 2005, 2005, p. 17.

[14]
H. Lee, R. Grosse, R. Ranganath, A. Y. Ng, Convolutional deep belief networks for scalable unsupervised learning of hierarchical representations, in: Proceedings of the 26th Annual International Conference on Machine Learning, ACM, 2009, pp. 609–616.
 [15] R. Cutler, M. Turk, Viewbased interpretation of realtime optical flow for gesture recognition, in: Automatic Face and Gesture Recognition, IEEE, 1998, pp. 416–421.
 [16] G. Yang, Y. Lin, P. Bhattacharya, A driver fatigue recognition model based on information fusion and dynamic bayesian network, Information Sciences 180 (10) (2010) 1942–1954.
 [17] B. Packer, K. Saenko, D. Koller, A combined pose, object, and feature model for action understanding, in: CVPR, IEEE, 2012, pp. 1378–1385.

[18]
G. Hinton, R. Salakhutdinov, Reducing the dimensionality of data with neural networks, Science 313 (5786) (2006) 504–507.
 [19] V. Nair, G. Hinton, 3d object recognition with deep belief nets, NIPS 22 (2009) 1339–1347.
 [20] Y. Wu, Z. Wang, Q. Ji, Facial feature tracking under varying facial expressions and face poses based on restricted boltzmann machine, in: CVPR, IEEE, 2013.
 [21] G. Taylor, G. Hinton, S. Roweis, Modeling human motion using binary latent variables, NIPS 19 (2007) 1345.

[22]
A. Kae, K. Sohn, H. Lee, E. LearnedMiller, Augmenting crfs with boltzmann machine shape priors for image labeling, in: Computer Vision and Pattern Recognition (CVPR), 2013 IEEE Conference on, IEEE, 2013, pp. 2019–2026.
 [23] A. Kae, B. Marlin, E. LearnedMiller, The shapetime random field for semantic video labeling, in: Computer Vision and Pattern Recognition (CVPR), 2013 IEEE Conference on, IEEE, 2014.

[24]
I. Sutskever, G. E. Hinton, Learning multilevel distributed representations for highdimensional sequences, in: Artificial Intelligence and Statistics, 2007, pp. 544–551.
 [25] H. Larochelle, Y. Bengio, Classification using discriminative restricted boltzmann machines, in: ICML, ACM, 2008, pp. 536–543.
 [26] J. Louradour, H. Larochelle, Classification of sets using restricted boltzmann machines, in: UAI, IEEE, 2011.
 [27] G. Hinton, Training products of experts by minimizing contrastive divergence, Neural computation 14 (8) (2002) 1771–1800.
 [28] G. Hinton, A practical guide to training restricted boltzmann machines, Momentum 9 (2010) 1.
 [29] R. Salakhutdinov, I. Murray, On the quantitative analysis of deep belief networks, in: Proceedings of the 25th international conference on Machine learning, ACM, 2008, pp. 872–879.
 [30] T. Schmah, G. E. Hinton, S. L. Small, S. Strother, R. S. Zemel, Generative versus discriminative training of rbms for classification of fmri images, in: NIPS, 2008, pp. 1409–1416.
 [31] E. Hüllermeier, J. Fürnkranz, W. Cheng, K. Brinker, Label ranking by learning pairwise preferences, Artificial Intelligence 172 (16) (2008) 1897–1916.
 [32] Z. Wang, S. Wang, Q. Ji, Capturing complex spatiotemporal relations among facial muscles for facial expression recognition, in: CVPR, 2013.
 [33] A. Lorincz, L. A. Jeni, Z. Szabó, J. F. Cohn, T. Kanade, Emotional expression classification using timeseries kernels, in: Computer Vision and Pattern Recognition Workshops (CVPRW), 2013 IEEE Conference on, IEEE, 2013, pp. 889–895.
 [34] F. Long, T. Wu, J. R. Movellan, M. S. Bartlett, G. Littlewort, Learning spatiotemporal features by using independent component analysis with application to facial expression recognition, Neurocomputing 93 (2012) 126–132.
 [35] P. Yang, Q. Liu, D. N. Metaxas, Boosting encoded dynamic features for facial expression recognition, Pattern Recognition Letters 30 (2) (2009) 132–139.
 [36] L. A. Jeni, J. M. Girard, J. F. Cohn, F. De La Torre, Continuous au intensity estimation using localized, sparse facial feature space, in: Automatic Face and Gesture Recognition (FG), 2013 10th IEEE International Conference and Workshops on, IEEE, 2013, pp. 1–7.
Comments
There are no comments yet.