Gait Correlation Analysis Based Human Identification

Human gait identification aims to identify people by a sequence of walking images. Comparing with fingerprint or iris based identification, the most important advantage of gait identification is that it can be done at a distance. In this paper, silhouette correlation analysis based human identification approach is proposed. By background subtracting algorithm, the moving silhouette figure can be extracted from the walking images sequence. Every pixel in the silhouette has three dimensions: horizontal axis (x), vertical axis (y), and temporal axis (t). By moving every pixel in the silhouette image along these three dimensions, we can get a new silhouette. The correlation result between the original silhouette and the new one can be used as the raw feature of human gait. Discrete Fourier transform is used to extract features from this correlation result. Then, these features are normalized to minimize the affection of noise. Primary component analysis method is used to reduce the features' dimensions. Experiment based on CASIA database shows that this method has an encouraging recognition performance.


Introduction
Biometrics is a technology that makes use of the physiological or behavioral characteristics to authenticate or identify people [1]. The most commonly used biometrics applications are fingerprint and iris based identification.
Human gait was firstly studied in the medical field [2][3][4][5]. Doctors analyzed human gait to find out whether patients had health problem. Later, researchers [5] found that just like fingerprint and iris, almost everyone had his distinctive walking style. So someone believed that gait could also be used as a biological feature to identify a person. Comparing to fingerprint and iris based identification, human gait identification has the following advantages [6]: (1) it does not require the user's interaction and (2) it can be done at a distance, as long as the gait is visible.
Current human gait identification can be divided into two categories: model-based methods and motion-based methods.
Model-based approaches aim to describe human movement using a mathematical model, for example, Cunado et al. [7] used Hough translation to extract arms, legs, and torso and use articulated pendulum to match the moving body parts. Yoo et al. [8] divide the body into head, neck, waist, leg, and arm by image segmentation and then obtained the moving curve of these body parts. Lee and Grimson [9] applied 7 ellipses to model the human body and applied the ellipses' movement features to identify human. Yam et al. [10,11] used dynamically coupled oscillator to describe and analyze the walking and running style of a person. Tafazzoli and Safabakhsh [12] constructed movements model based on anatomical proportions then Fourier transform was used to analysis human walking style. Dupuis et al. [13] created probabilistic based gait modeling to describe the human walking.
Motion-based methods consider the human gait as a sequence of image and extract features from these images. Cheng et al. [14] employed Hidden Markov Models to analyze the relationship among these images. Chen et al. [15] used parallel HMM to describe the features of human gait. Yu et al. [16] overlapped all the images to get Gait Energy Image (GEI) and GEI was used as the features to identify humans. Similar to GEI, Fan et al. [17] took Chrono-Gait Image as the gait features. Kale et al. [18] used "frieze" patterns to  extract features from image sequence and used these features to identify humans. Liu et al. [19] represent one's walking style by outermost contour. The raw features extracted by motion-based methods usually contained several cures or images, so Gabor transform [20], Random Forest algorithm [12], and Fourier transform [21,22] were used to reduce the dimensions of these features.

Overview of Our
Approach. Unlike traditional image correlation based identification [16,23,24], this paper aims to build a human gait identification method not only in spatial domain but also in temporal domain. A median model is used to estimate the background. By background subtraction, we get the binary human silhouette. Every silhouette is correlated with its neighbor at horizontal, vertical, and temporal axis. The size of the correlation result is used to construct the walking feature curves and every feature curve consists of several walking cycles. The feature curve is normalized to build a standardized curve. For every image sequence, we can get totally 13 standardized curves. Primary component analysis is used to extract features from these above curves. The overview of the proposed algorithm is shown in Figure 1.

Our Contribution.
The main purpose and contributions of this paper can be summarized as follows.
(1) We extract 13 feature curves from the human gait images sequence. These 13 curves described the feature of human walking. (2) We proposed a method to construct normalized human gait feature vectors. PCA is used to reduce the dimensions of gait feature vectors.
(3) We discuss the relationship between accuracy of identification and the shift length along horizontal axis and vertical axis and try to find out the best shift length.

Silhouettes Sequence Extraction.
In this paper, the human gait is considered as a binary image. Gauss background model is used to estimate the background. We use the method proposed by Wang et al. [25] to subtract the background from the foreground of the image as follows: ( , ) is the difference between point (pixel in foreground image) and point (pixel in background image).

About the Image Sequence Correlation.
Image correlation is a shift-invariant feature method and was first proposed by Otsu and Kurita [26]. Let ( ) be -dimensions data with = ( 1 , 2 , . . . , ); then the correlation function is defined as where ( ,1 , . . . , , ) is the shift vector. Given a human gait image sequence, ( ) means an image in the sequence whose frame index is . If = ( , , , , , ), ( + ) means we get the frame whose frame index is + , , then shift its , pixels along -axis and , pixels along -axis. The three dimensions of the image are shown in Figure   Given a human walking cycle, images sequence consists of frames, the features of this human gait images sequence can be described as . . .
(3) Figure 4 shows a gait cycle. In the first stance, the person is at rest and the silhouette size is minimum, this corresponds to the valley at Figure 5; then he steps out his leg, the silhouette size will increase with his stepping out, in the sixth stand his legs and arms are completely separated and this corresponds to the maximum at Figure 5. Then, the legs and arms are closed, at stance 11 the size reaches minimum again. In fact, this is not a complete gait cycle but a "half " cycle, but the previous "half " cycle and the remaining "half " cycle are symmetrical, so in this paper we call this "half " gait cycle as one gait cycle. Formula (3) is the matrix expression of the curves shown in Figure 3. Ω can be used as features to identify humans. But these features are not robust. The human silhouettes are prone to be affected by the noise and deformation. We will transform the curve into standard curve to make features more robust by the following steps.

The Normalized Human Gait Image Correlation.
(1) Discrete Fourier transforms: by dividing the walking images into walking cycles as shown in Figure 3, we can find the start (end) of a walking cycle by finding the minimum point in Figure 5.  [ ] is the discrete points and̂[ ] is the Fourier coefficients.
(2) Normalization: Considering that different walking cycles have different frame count ( ), we need to "align" to a fixed value. In this paper, we divide every walking cycle into parts, where is a fixed value (in this paper = 12). That is Given person , the normalized feature matrix for walking cycle can be expressed as 4 The Scientific World Journal Given walking cycles, we use the normalized Ω( , ) as the final features. The final feature matrix for person can be expressed asΩ By inverse Fourier transforms (Formula (5)), we can also get the correlation image size fromΩ( ) inversely; the correlation image size calculated by Fourier transforms is shown in Figure 5.

Reduce Dimensions of Features.
We define Ω as follows: is the human count in the test dataset. That is to say, Ω is the average value of all human gait features. Then, we calculate the covariance matrix COV as follows: Principal component analysis is used to reduce the dimensions ofΩ( ). The reducedΩ( ) can be expressed as where ( 1 , 2 , . . . , ) are the eigenvectors selected by PCA from covariance matrix COV.

Human Gait Identification.
In the scenario of identification, the gait feature distance between human and can be defined as follows: where ‖ ‖ means the Euclidean distance. This probe is assigned to person by using nearest neighbor method. That is In the scenario of verification, the similarity between two feature matrixes is defined as the negative of distance, that is In this paper, the similarity between a probe and in the gallery is defined as -normed similarity [27] as follows:

Experiment and Discussion
We use the CASIA database (dataset B) [28] to evaluate our method. There are total 124 people in dataset B. For every person, there are at least 3 walking sequences captured from 12 view angles. Every walking sequence consisted of 3-4 walking cycles. In this paper, most of the experiments are based on 90 ∘ view angle except with declaration. We randomly select 2 walking sequences as the training data to get the normalized human gait vector. Then we use the remained walking sequence to evaluate performance. In this paper, for the consideration of simplicity, , has only three values: {−Δ , 0, −Δ }, , has three values: {−Δ , 0, −Δ }, and , has two values: {0, Δ }. Then, we can totally get 3 × 3 × 2 = 18 , = 1, . . . , 18. Using these shift vectors , we totally get 18 feature vectors for one human gait images sequence. By using PCA, we reduce these 18 vectors to 3 vectors.

Features Extracted by Image Correlation.
For a human gait image sequence (consisting of about four waking cycles), the 13 features curves are shown in Figure 6 (Δ = 5, Δ = 5, Δ = 1). Figure 6, we can see that the 13 image correlation curves have some similarity. That is to say that there is some redundant information. We use formula (11) to reduce the dimension. We set threshold = 90% and compress the 13 human gait feature vectors into 3 vectors. Figure 7 shows these three vectors for one human walking cycle.

The Result of Dimension Reduction. From
We also compared the recognition accuracy result using different dimensions (that is the in Formula (11)). The experiment result is showed in Figure 8.
From Figure 8, we can see that the recognition result will reach its maximum at = 4.

The Affection of Shift Length.
From the definition of image correlation, we can see that the shift length (Δ , Δ , Δ ) along these three axes affects the correlation result. The shift length along the horizontal and vertical axes should follow below principles: the features distance between two gaits images sequences calculated by these shifts should have the maximum value. That is to say, all samples should have the maximum standard deviation as follows: To find the optimized shift length, we vary shift length along -and -axis from 1 to 10 pixels and calculate the standard deviation. To eliminate the affection of silhouette size, we divide the standard deviation by silhouette size. The calculation result is shown in Figure 9. The accuracy of reorganization is showed in Figure 10. For completeness, we also estimate FAR (False Acceptance Rate) and FRR (False Rejection Rate) in verification mode. The ROC (Receiver Operating Characteristic) curves are shown in Figure 11.
Comparing Figures 9, 10, and 11, we can see that both the accuracy of reorganization and ERR (Equal Error Rate) will reach optimized value at approximately the maximum point of standard derivation in Figure 9. In fact, the 5 pixels length is just the moving length along -axis between two frames and the 2 pixels length is just the average moving pixels along -axis and -axis between two frames.

The Affection of View Angle.
As we know, the view angle may affect the recognition accuracy, so we compare the performance under differential view angle. Table 1 shows the experiment result. From Table 1, we can see that our method has better performance under view angles 0, 90, and 180. That may be because the image token from above the angle is not prone to be affected by deformation.

Discussions and Future
Works. The correlation of two images in temporal domain and spatial domain produces a new image. This image contains much information about the moving object, especially the dynamic features of the human gait. In this paper, only the areas of the correlated images are is used as the features for identification. In the future, we should try to use other features about the correlation images,   for example, the texture features about the correlation images, or features from projecting the correlation images to -oraxis. In this paper, we only talk about two images correlation. Maybe the correlation among three or more images can also be used in the human gait identification. We will also try this in the future.

Conclusions
With the increasing demand of security control, human gait based identification will attract more interest. In the future, gait based human identification might finally become applicable just as fingerprint and iris.
The kernel idea of this paper is based on the analysis of the image correlation. Fourier transform is used to get The Scientific World Journal 7 the features from the image correlation. We discuss the relationship between shift length and the identification accuracy. We also use CASIA gait database to validate our method. The maximum classification accuracy is nearly 90% and EER is approximately 10%.

Conflict of Interests
The author declares that he has no financial and personal relationships with other people or organizations that can inappropriately influence his work. There is no professional or other personal interests of any nature or kind in any product, service, and/or company that could be construed as influencing the position presented in, or the review of, the paper entitled.