Scholarly article on topic 'Gait Characteristic Analysis and Identification Based on theiPhone’s Accelerometer and Gyrometer'

Gait Characteristic Analysis and Identification Based on theiPhone’s Accelerometer and Gyrometer Academic research paper on "Medical engineering"

Share paper
Academic journal
OECD Field of science

Academic research paper on topic "Gait Characteristic Analysis and Identification Based on theiPhone’s Accelerometer and Gyrometer"

Sensors 2014,14, 17037-17054; doi:10.3390/s140917037



ISSN 1424-8220


Gait Characteristic Analysis and Identification Based on the iPhone's Accelerometer and Gyrometer

Bing Sun 1;*, Yang Wang2 and Jacob Banda3

1 School of Electronics and Information Engineering, Beihang University, 37 Xueyuan Road Haidian District, Beijing 100191, China

2 School of Electrical and Electronic Engineering, Nanyang Technological University, 50 Nanyang Avenue, 639798, Singapore; E-Mail:

3 Department of Mathematics, the University of Texas-Pan American, 1201 West University Drive, Edinburg, TX 78539, USA; E-Mail:

* Author to whom correspondence should be addressed; E-Mail:; Tel.: +86-10-8233-8670.

Received: 7 August 2014; in revised form: 29 August 2014 / Accepted: 9 September 2014 / Published: 12 September 2014

Abstract: Gait identification is a valuable approach to identify humans at a distance. In this paper, gait characteristics are analyzed based on an iPhone's accelerometer and gyrometer, and a new approach is proposed for gait identification. Specifically, gait datasets are collected by the triaxial accelerometer and gyrometer embedded in an iPhone. Then, the datasets are processed to extract gait characteristic parameters which include gait frequency, symmetry coefficient, dynamic range and similarity coefficient of characteristic curves. Finally, a weighted voting scheme dependent upon the gait characteristic parameters is proposed for gait identification. Four experiments are implemented to validate the proposed scheme. The attitude and acceleration solutions are verified by simulation. Then the gait characteristics are analyzed by comparing two sets of actual data, and the performance of the weighted voting identification scheme is verified by 40 datasets of 10 subjects.

Keywords: gait; iPhone; accelerometer; identification; weighted voting

1. Introduction

Gait refers to an idiosyncratic feature of a person, that is determined by an individual's weight, limb length, posture combined with characteristic motion [1] and so on. Also gait may reflect the health condition of humans [2,3], such as symptoms of Cerebral Palsy and Parkinson's disease. Hence, gait can be used as a biometric measure to recognize known persons and classify unknown subjects [1]. Gait identification allows a nonintrusive way of data collection from a distance. This procedure can provide supplemental information for the traditional approaches (e.g., fingerprint, iris, face, voice) in security and usability applications, thereby enhancing or complimenting them.

At present, gait analysis can be classified according to two different approaches: non-wearable sensors (NWS) based and wearable sensors (WS) based [4]. NWS systems usually require optical sensors [5-8] or pressure sensors on the floor [9]. Gait features are extracted by processing image or video data. In the security surveillance field, this approach can effectively narrow the range of subjects in gait database and plays an important role. WS systems are based on motion-recording sensors which are attached to moving subjects [10-16]. Motion-record sensors (e.g., accelerometers, gyrometers, force sensors, bend sensors, and so on) are usually attached to various parts of the body such as the ankle, hip and waist [12,16]. Then, gait signatures are extracted by analyzing the recorded signals in both time and frequency domains. With the application of a Microelectromechanical System (MEMS), a device such as an Inertial Measurement Unit (IMU) is easily inserted to shoes, gloves, watches and other wearable attire. The prospect of applying wearable sensors to gait analysis is discussed in [12-15].

Smart mobile phones, nowadays widely equipped with MEMS accelerometers and gyrometers to determine device orientation, provide developers and users with a rich and easily accessible platform for exploration. Plenty of application programs are developed based on smartphone sensors. An iPhone is used as a wireless accelerometer sensor for Parkinson's disease tremor analysis in [17]. The feasibility of gait identification using an iPhone's accelerometer is discussed in [18]. However, current study of gait analysis based on fixed smartphone sensors mostly use data including gravitational acceleration without considering the effects of attitude change. Our goal is mainly to analyze gait characteristics and identification with an unfixed iPhone in pant pockets instead of professional sensors and laboratory conditions. Under this case, attitude of the iPhone should be considered to minimize the effects of attitude change and gravitational acceleration. After extracting some gait information, a voting scheme is investigated for gait identification.

This paper is organized as follows. In Section 2, pre-processing methods for iPhone inertial sensor data are introduced. In Section 3, gait characteristics are analyzed and methods are studied to extract gait parameters from accelerometer and gyrometer data. In Section 4, a weighted voting scheme based on gait signatures is proposed for gait identification. In Section 5, validating experiments are conducted for gait analysis and identification using both simulated data and actually collected data.

2. Pre-Processing of iPhone Measurement Data

Several kinds of sensors are integrated in an iPhone like a GPS sensor, accelerometer, gyrometer, magnetometer, and so on. Theoretically, the motion of an iPhone (including its acceleration, velocity, and position) can be solved based on the sensor data. However, according to experiments, the precision

of iPhone sensors [19] can hardly support accurate solutions of device velocity and position, especially for short-distance measurements. According to the operating principles of iPhone accelerometers, the collected data is in its own inertial coordinate system, and has gravitational acceleration included. During the collection of an individual's gait data, the attitude of an iPhone is changing so it is not an accurate method to directly process the measured acceleration data. In this section, a pre-processing approach is studied to get the linear acceleration, in which the gravitational acceleration is not included. This is achieved by defining a reference coordinate system with its Z-axis corresponding to the opposite direction of gravity.

2.1. Inertial Data Calibration

Although the inertial sensors are calibrated by the manufacturers, there are still some systematic errors for applications. Figure 1 illustrates the direct outputs of the triaxial accelerometer and gyrometer, which are collected when the iPhone is placed on a horizontal surface for 120 s. It is obvious that the values of z-axis acceleration deviate from the theoretical value -1.0 * g0, where g0 = 9.81 m/s2 is the gravitational acceleration constant. A similar situation occurs for the gyrometer. So the calibration of inertial data is needed.

Figure 1. The direct outputs of triaxial (a) accelerometer and (b) gyrometer.

(a) (b)

2.1.1. Static State Judgement

To determine the initial attitude of an iPhone, an approximate static state is required. If the variance of the complex acceleration amplitude is below a threshold in a short period (e.g., 1 s), the iPhone can be regarded as static or at uniform motion in a straight line. Due to ups and downs during walking, uniform motion is almost impossible. In a static state, the measured acceleration data corresponds to the gravity only.

2.1.2. Acceleration Bias

There may be systematic errors in accelerometers because of device differences, device aging, or even occasional damage to the device. The systematic error is coupled with motion error and gravitational

acceleration, and leads to the acceleration bias in measurements. The acceleration bias can be estimated by recording triaxial accelerations under static state and calculating their mean values. Thus in the processing of gait data, acceleration bias can be eliminated by directly subtracting the estimated bias from triaxial acceleration data. In Figure 1a, the main bias comes from z-axis, which is about -0.05 * g0.

2.1.3. Gyrometer Bias and Drift

The theoretical outputs of triaxial gyrometer should be zeros in a static state. We also calculate their mean values to be the systematic biases and subtract these biases from triaxial gyrometer data. Usually the biggest problem for a gyrometer is its drift during long time measurement. However, according to Figure 1b we can see the drift in 120 s is almost zero. Therefore in a short time period the drift can be neglected in gait measurement, which is usually sampled in less than 20 s intervals.

2.2. Attitude during Motion

To study the attitude of the iPhone sensor during motion, a reference coordinate system is defined based on gravity and the initial attitude of the iPhone sensor. In this paper, attitude is described by quaternion, which can be calculated via the combination of acceleration and gyrometer data. Then the coordinate transformation matrix is solved. The specific process is as follows:

(1) To define the reference coordinate system:

The iPhone sensor has its own inertial coordinate system [O,Xi,Yi, Zj], as shown in Figure 2. Define Z-axis of the reference coordinate system opposite to the gravity vector, u as the angle between the Z-axis and Z-axis, and a unit vector C = (c1,c2, c3) which is vertical to both Zj- and Z-axis as the rotation axis. For the initial attitude of iPhone sensor, the triaxial accelerations are A0 = (ax0, ay0, az0)T in iPhone inertial coordinate system, and the gravitational acceleration is B0 = (0,0, — g0)T in reference coordinate system, then u and C can be obtained as u = acos^ \a0\!\bo|), C = . After a rotation

around C by an angle of u, the Zj-axis is rotated to Z-axis, and the Xj-axis and Yj-axis are rotated to the horizontal plane, in which they are now defined as X-axis and Y-axis. [O, X, Y, Z] is the reference coordinate system for further analysis.

According to the quaternion as rotations [20,21], the initial quaternion Q0 = [w0, x0, y0, z0]T is solved as:

( arccosu\

wo = cos ^-2-) (1)

x0 = ci sin u (2)

yo = c2 sin u (3)

z0 = c3 sin u (4)

It means one time rotating around unit vector C=(c1,c2,c3) by u transforms the iPhone inertial coordinate to the defined reference coordinate.

Figure 2. iPhone inertial coordinate system and reference coordinate system.

In a further step, the three Euler angles of the iPhone can be solved based on quaternion. However, the values of Euler angles depend on their definition and rotations [22]. In this paper, quaternion is adopted for iPhone solution of attitude.

(2) Update quaternion:

T / ' \ T

To update quaternion Qi = [wi,Xi,yi, Zi] with combined acceleration data Ai = ialx,ai,alz) and

gyrometer data gi = (дгх,дгу,gl)T using the Fourth-order Runge-Kutta Algorithm [23], the specific process is as follows:

Qi+i = Qi + T (Ki + 2K2 + 2K3 + K4) 6

K = 2 Qi * gi

Qi + § t

Qi + K2 t

K4 = 2 [Qi + K3T ] *gi+i

(8) (9)

where Qi * gi refers to quaternion product,

Qi *gi

—xi -yi -zi

Z' Ш'

у' x' ш'

and 9i+i=2 (gi + gi+i).

(3) Calculate coordinate transformation matrix:

To find the coordinate transformation matrix Tim for t = i, according to real-time quaternion Qi = [wi, xi, yi, zi]T. The coordinate transformation matrix [22] is

1 - 2 (y2 + z2) 2 (xiyi - UiZi) 2 (xiZi + Uiyi)

2 (xiyi + ^iZi) 1 - 2 (x2 + Z2) 2 (yiZi - UiXi) 2 (xiZi - Uiyi) 2 (yiZi + Wixi) 1 - 2 (x2 + y2)

2.3. Linear Acceleration Solution

Gravitational acceleration can be eliminated to get the linear acceleration (llx, ly, lZ)Tusing the coordinate transformation matrix, described as follows

lx 0

ly i T om ay - 0

r _ aZ. _ -g0 _

The first term of the right side of Equation (12) is the total acceleration in reference coordinate system. After subtracting the gravitational acceleration in the same reference coordinate system, we can get the linear acceleration. Since the two components IX and ly correspond to the horizon direction of motion, the following gait characteristic analysis is based mainly on the component of linear acceleration l\.

3. Gait Characteristic Analysis

In this section, we discuss four gait characteristic parameters: gait frequency, symmetry coefficient, dynamic range and similarity degree between characteristic curves. The location-related parameters are not considered because the integration of acceleration will generate accumulation error.

3.1. Gait Frequency

Gait frequency is the most basic and important parameter in gait characteristics analysis. Due to the quasi-periodic characteristic of paces during walking, gait frequency can be extracted by the characteristic of the acceleration signal in the frequency domain. After FFT of the acceleration signal, the maximum frequency, excluding DC component, is not always the gait frequency. An improved method is presented in this paper. The autocorrelation of the acceleration signal is calculated and converted to frequency domain. Then the maximum frequency, excluding zero frequency component, equals to gait frequency. Figure 3 shows the spectrum of acceleration signals and their autocorrelation. Interpolation processing can greatly improve the accuracy of gait frequency while calculating correlation coefficient.

According to the Figure 3, the spectrum lines of the autocorrelation function are more focused than those of the acceleration. Thus the autocorrelation function of acceleration is usually less affected by harmonic components and easier to extract gait frequency.

Figure 3. Acceleration spectrum: (a) FFT of autocorrelation; (b) FFT of acceleration.

3.2. Symmetry Coefficient

Some literature has defined the gait symmetry corresponding their own dataset [24]. In this paper, we define the symmetry coefficient Si of gait according to autocorrelation of acceleration signal, shown in Figure 4.

Si = 1

2 VC C ' °3>

V max m

where Cmax is the maximum of the autocorrelation of acceleration, C and Cr are the maximum values of the points located on the left and right of the peak position half gait period away, respectively.

Figure 4. Schematic definition of gait symmetry. The blue line is the autocorrelation curve of acceleration.

Two possible causes may lead to asymmetry of the acceleration signal: (1) inconsistent gait, such as the gait of patients with unilateral actions and other obstacles; (2) sensors are not strictly centered placed. The first case is not considered in this paper, and all subsequent tests are based on almost symmetrical gait. We will discuss the effect of sensor's positions in Section 5.

3.3. Dynamic Range

The dynamic range of acceleration can reflect a gait characteristic in some aspects. In this paper, the difference between the maximum and minimum values of Z-linear acceleration is defined as dynamic range estimation.

Usually, a few aspects including the speed of walking, pace length, noise and so on will affect the dynamic range. It is reasonable that dynamic range should has less contribution during gait identification.

3.4. Similarity Degree between Characteristic Curves

For typical gait acceleration, the characteristic curve can be extracted according to its fluctuation characteristics, that is the time-amplitude sequence. This sequence can be used as an important basis for the subsequent gait matching and identification. By comparing the characteristic curve with each characteristic curve in database, a high similarity degree indicates the two curves correspond to the same people with a high probability, and it is helpful in gait identification.

To accurately estimate the similarity degree between two characteristic curves, scale transformation and panning of the curves have to be taken into consideration. Even for the same subject, walking speed and magnitude may have subtle variations under different test conditions. The discrete acceleration data cannot be directly used for similarity comparison. Although a one-dimensional invariant matrix method can be used theoretically, the experimental performance is not ideal. In this paper, the similarity degree between characteristic curves is calculated as follows:

1. apply interpolation processing, according to the gait frequencies, to eliminate the effects of scale differences.

2. find the peak positions of two curves as reference, and make a rough alignment of the vectors. The aligned acceleration characteristic curves are signed as li, l2.

3. calculate the similarity coefficient, which is used to describe similarity degree, based on the normalized dot product of two vectors. Moreover, to eliminate noise affect in characteristic curves, we translate one vector within one pace period, which corresponds to N samples, then calculate the normalized dot product of the vectors, and finally select the maximum normal product as the similarity coefficient C.

D = max (/Z) - min (Tz)

E [/i (i + m) • /2 (i)]

C = max

4. Gait Identification

In the previous section, several gait characteristic parameters including gait frequency, symmetry coefficient, dynamic range and similarity coefficient are discussed. Based on these typical parameters, a weighted voting scheme is proposed for gait identification. The specific scheme is as follows:

1. calculate gait frequency Fi, symmetry coefficient Si, dynamic range Di for all the characteristic curves in the database one by one. Suppose the number of the curves in database is M,

i = 1, 2,..., M.

2. calculate gait frequency F, symmetry coefficient S, dynamic range D of the linear acceleration from the measured data.

3. calculate the similarity coefficient Ci between the characteristic curves of linear acceleration and samples in the database.

4. the weighted voting process is as follows:

(a) sort the absolute error between F and Fi in increasing order. We vote the first term 1, the second term 2, and so on. The number of votes of gait frequency is denoted as Vi1.

(b) sort the absolute error between S and Si in increasing order. We vote the first term 1, the second term 2, and so on. The number of votes of symmetry coefficient is denoted as Vi2.

(c) sort the absolute error between D and Di in increasing order. We vote the first term 1, the second term 2, and so on. The number of votes of dynamic range is denoted as Vi3.

(d) sort Ci in decreasing order. We vote the first term 1, the second term 2, and so on. The number of votes of similarity coefficient is denoted as Vi4.

(e) sum the weighted number of votes for each sample in the database.

V = WiV1 + W2V2 + W3V3 + W4V4, i = 1, 2,..., M (16)

where W = [w1, w2, w3, w4] is the weighted coefficient vector.

(f) vote judgment. The sample in the database that has the minimum weighted vote summation matches the current measured data, and we can accomplish gait identification within the given database.

The weighted coefficient can be chosen by further experimental data. The similarity coefficient of characteristic curves is of the most importance for gait identification, followed by gait frequency. The symmetry coefficient has a close relationship with the iPhone's placement. The dynamic range is greatly affected by noise. Therefore, the weighted coefficients are set as W = [2, 2, 1, w4] in this paper, where

' 8 Cmax e [0.9,1]

4 Cmax e O) W4 = < , Cmax = max {Ci} (17)

2 Cmax e [0.7, 0.8) ^^

1 Cmax e [0, 0.7) 5. Experiments

Four experiments are conducted to validate the methods proposed in this paper: (1) simulation experiment to verify the solutions of attitude and linear acceleration; (2) comparison experiment of the same subject walking on cement pavement and on grass; (3) comparison experiment of three different placements of iPhones; (4) gait identification experiment with 10 subjects and 40 sets of data.

In the first experiment, the sinusoidal linear acceleration and attitude data are simulated. The attitude angles (yaw, pitch and roll) are defined by Tait-Bryan angles formalism with rotation sequence z-y-x [25].

Then the triaxial accelerometer and gyrometer data in iPhone's inertial coordinate system is generated as the input of the proposed method in Section 2. After calculation and comparison, the inversion linear acceleration and attitude data, and the errors are shown in Figure 5.

According to simulation results, the inversion errors of attitude are less than 1°, and the inversion errors of acceleration are less than 0.1 m/s2.

In the second experiment, the subject walked on cement pavement and grass, respectively, with the same iPhone placed in the right pant pocket. The characteristic curves of recorded acceleration data are shown in Figure 6.

At least two differences between the two datasets can be seen in Figure 6.

1. The data recorded on cement pavement is nearly 10 cycles, while the data recorded on grass is less than 9 cycles, both recorded in 10 s. They are consistent with the gait periods calculated by gait frequency analysis, which are 1.070 s and 1.174 s, respectively. This is due to the fact that it's more difficult to walk on soft surfaces like grass than on cement pavement.

2. The dynamic range of linear acceleration recorded while walking on cement pavement is a little larger than that on grass, and it can be explained by the buffer caused by grass.

In the third experiment, three iPhones are placed in three different positions: the left pant pocket, the right pant pocket and the center of the back waist area of the subject. The calculated linear accelerations and autocorrelation curves are shown in Figure 7.

According to Figure 7a,c,e, the linear acceleration signal reflects preferable quasi-periodicity. Moreover, the curves in Figure 7a,c are very similar and differ mainly in half gait period in time. The third data shows good symmetry among each single step, and the frequency is twice of the above two curves. According to Figure 7b,d, the accelerations corresponding to the left and right pace are asymmetrical, obviously due to the off-center displacement. Furthermore, the symmetry coefficients are only about 0.4. Figure 7f describes the third case, when the iPhone was placed in the center of the back waist. The left and right paces cause similar linear accelerations, and the symmetry coefficient is more than 0.9. The definition of symmetry coefficient is also validated. Moreover, the dynamic ranges of Figure 7a,b are greater than that in Figure 7c as a result of the up and down leg motion versus the body motion within the same pace period.

In the fourth experiment, a single iPhone is used to record 4 sets of gait data for each of the 10 tested subjects. The phone is placed in their left and right pant pockets twice, respectively. Then one set of data is selected for each subject in order to make the gait database by calculating characteristic parameters and extracting characteristic curves, shown in Table 1 and Figure 8.

Table 1. Indices for the database. D1~D10 represent the 10 samples in the database. The units of gait frequency and dynamic range are Hz and m/s2, respectively.

Indices D1 D2 D3 D4 D5 D6 D7 D8 D9 D10

Gait Frequency Symmetry Coef. Dynamic Range 1.036 0.534 26.480 1.041 0.563 29.301 0.957 0.571 25.300 1.009 0.404 22.379 0.939 0.669 24.891 0.845 0.525 25.618 1.002 0.623 27.051 0.945 0.597 20.674 0.987 0.576 21.567 1.019 0.914 14.736

Figure 5. Inversion results of simulated data. The top, middle and bottom curves in (a) are the input data, inversion data and errors of the roll, pitch and roll attitudes, respectively. The top, middle and bottom curves in (b) are the input data, inversion data and errors of the x-, y- and z- linear accelerations, respectively.

Figure 6. Comparison experiment with different ground conditions. (a) Test conditions and the two test paths. (b) Linear acceleration data of walking on cement pavement along path 1. (c) Linear acceleration data of walking on grass along path 2.


4 5 6 7 8 9 10 11 T£ 13- 14

time(s) (b)

-15 ■

-20 I-'-'-'-'-'-'-'-'-'-

4 5 B 7 E 9 10 11 13 13- 14


Figure 7. Linear acceleration and autocorrelation curves. (a) and (b) correspond to the data collected by the iPhone in left pant pocket. (c) and (d) correspond to the data collected by the iPhone in right pant pocket. (e) and (f) correspond to the data collected by the iPhone in the center of the back waist area.

Figure 8. Characteristics curves of ten tested persons as the database.

time(s) time(s) time(s)

(g) (h) (i)


The remaining 30 sets of recorded data are used as measured data to compare with 10 samples in the database. The voting results of 10 groups of measured data is listed in Table 2. We calculate the weighted votes for measured data M1^M10 according to the method described in Section 4. For example, data M1 yields gait frequency F = 1.039 Hz, symmetry coefficient

Si = 0.554, and dynamic range D = 25.87 m/s2. The similarity coefficients between Ml and D1-D10 are [C1, C2, C3, C4, C5, Cg, C7, C8, C9, Cio]=[0.876, 0.773, 0.667, 0.762, 0.793, 0.498, 0.614, 0.735, 0.794, 0.575]. After comparing the above parameters with Table 1 and sorting the absolute errors in increasing order and similarity coefficients in decreasing order, we can get the votes of gait frequency [V/, V21, V31, V41, V51, Vg1, V1, V81, V91, V110]=[2, 1, 7, 4, 9, 10, 5, 8, 6, 3], the votes of symmetry coefficient [V2, V22, V32, V42, V52, Vg2, V72, V82, V92, V120]=[3, 1, 2, 9, 8, 5, 7, 6, 4, 10], the votes of dynamic range [V3, V23, V33, V43, V53, Vg3, V73, V83, V93, V3,]=[3, 6, 2, 7, 4, 1, 5, 9, 8, 10], and the votes of similarity coefficient [V4, V24, V34, V44, V54, Vg4, V74, V84, V94, V40]=[1, 4, 7, 5, 3, 10, 8, 6, 2, 9]. Following multiplication by W = [2, 2,1,w4], where = 4 is calculated using Equation (17), the summation votes, using Equation (16), equal to [17, 26, 48, 53, 50, 71, 61, 61, 36, 72] and are displayed in the first row in Table 2. The minimum of this row is 17, which corresponds to data D1. Because datasets D1 and M1 are collected from the same person in different time, the identification is true.

Table 2. Voting results for 10 measured data. M1~M10 represent the 10 measured data. D1~D10 represent the 10 samples in the database. The number in the frame is the minimum number in each row. True/False represents the current identification is true or false.

Data Name D1 D2 D3 D4 D5 D6 D7 D8 D9 D10 True/False

Ml M2 M3 M4 M5 M6 M7 M8 M9 M10

46 48 76

66 70 46 58 58 45 61

50 54 29 65 24 53 48

53 47 71 16

52 45 73

65 68 79

71 48 68

72 77 12 87 76 72 90

61 39 55

46 48 22 46 69 66

61 45 31 52 39

43 50 20 48

36 30 45 54 43 41 36 45 22 28

93 82 96 76 72

94 69 22

True True True True True True True True True True

According to Table 2, the weighted voting scheme has a good performance in identification experiments. And all the remaining 20 sets of recorded data are tested and correctly identified.

As a further step, the weighted coefficient vectors are set to W = [0,0,0,1], W = [1,1,1,1], W = [1,1,1, 2], W = [2, 3,1, w4] and W = [2, 2,1, w4], respectively. Then the numbers of successfully recognized datasets are 23, 26, 27, 27 and 30, respectively. It is known that the weighted coefficients have great effects on identification performances. Among the weighted coefficients, the similarity degree of characteristic curves has the highest effect, but the identification cannot depend on this coefficient alone. For experiments conducted in this paper, the weighted coefficient vector W = [2, 2,1, w4] was adopted and resulted in an excellent identification performances, while the error of identification is zero.

6. Conclusions

The accelerometer and gyrometer integrated in smartphones provide users with a convenient platform for gait analysis. In this paper, a new method for gait analysis is proposed based on an iPhone. The

Fourth-order Runge-Kutta algorithm and quaternion are applied to combine inertial data so as to solve linear acceleration and eliminate the errors caused by attitude change and gravitational acceleration. The gait characteristic parameters and characteristic curves are analyzed, then a weighted voting scheme is adopted for gait identification. Simulation and experiment results demonstrate good performance of the proposed scheme. In future study, the gait database will be improved to support gait identification under different scenarios, e.g., different motion status of tested subjects.


The authors would like to thank Ko, Young-woo for providing free iPhone application Sensor Monitor and the volunteers for providing their gait data.

Author Contributions

Bing Sun defined the methodology of the research and the structure of the paper and supervised the entire writing process. Yang Wang performed the main research. Jacob Banda performed the main writing of the paper.

Conflicts of Interest

The authors declare no conflict of interest. References

1. Lee, L.; Grimson, W.E.L. Gait analysis for recognition and classification. In Proceedings of Fifth IEEE International Conference on Automatic Face and Gesture Recognition, Washington, DC, USA, 20-21 May 2002; pp. 148-155.

2. Pogorelc, B.; Bosnic, Z.; Gams, M. Automatic recognition of gait-related health problems in the elderly using machine learning. Multimed. Tools Appl. 2012, 58, 333-354.

3. Azevedo Coste, C.; Sijobert, B.; Pissard-Gibollet, R.; Pasquier, M.; Espiau, B.; Geny, C. Detection of Freezing of Gait in Parkinson Disease: Preliminary Results. Sensors 2014, 14, 6819-6827.

4. Muro-de-la Herran, A.; Garcia-Zapirain, B.; Mendez-Zorrilla, A. Gait Analysis Methods: An Overview of Wearable and Non-Wearable Systems, Highlighting Clinical Applications. Sensors 2014, 14, 3362-3394.

5. Kale, A.; Sundaresan, A.; Rajagopalan, A.; Cuntoor, N.P.; Roy-Chowdhury, A.K.; Kruger, V.; Chellappa, R. Identification of humans using gait. IEEE Trans. Image Process. 2004, 13, 1163-1173.

6. Kale, A.; Cuntoor, N.; Yegnanarayana, B.; Rajagopalan, A.; Chellappa, R. Gait Analysis for Human Identification. In Audio-and Video-Based Biometric Person Authentication; Springer: Berlin Heidelberg, Germany, 2003; pp. 706-714.

7. Tang, J.; Luo, J.; Tjahjadi, T.; Gao, Y. 2.5 D Multi-View Gait Recognition Based on Point Cloud Registration. Sensors 2014,14, 6124-6143.

8. Yoo, J.H.; Hwang, D.; Moon, K.Y.; Nixon, M.S. Automated human recognition by gait using neural network. In Proceedings of First Workshops on Image Processing Theory, Tools and Applications, Sousse, Tunisia, 23-26 November 2008; pp. 1-6.

9. Middleton, L.; Buss, A.A.; Bazin, A.; Nixon, M.S. A floor sensor system for gait recognition. In Proceedings of Fourth IEEE Workshop on Automatic Identification Advanced Technologies, Buffalo, NY, USA, 17-18 October 2005; pp. 171-176.

10. Popovic, M.R.; Keller, T.; Ibrahim, S.; Bueren, G.V.; Morari, M. Gait identification and recognition sensor. In Proceedings of 6th Vienna International Workshop in Functional Electrostimulation, Vienna, Austria, 22-24 September 1998.

11. Sprager, S.; Zazula, D. Gait identification using cumulants of accelerometer data. In Proceedings of the 2nd WSEAS International Conference on Sensors, and Signals and Visualization, Imaging and Simulation and Materials Science, Baltimore, MD, USA, 7-9 November 2009; pp. 94-99.

12. Gafurov, D.; Snekkenes, E. Gait recognition using wearable motion recording sensors. EURASIP J. Adv. Signal Process. 2009, 2009, 1-16.

13. Bamberg, S.J.M.; Benbasat, A.Y.; Scarborough, D.M.; Krebs, D.E.; Paradiso, J.A. Gait analysis using a shoe-integrated wireless sensor system. IEEE Trans. Inf. Technol. Biomed. 2008, 12,413-423.

14. Pan, G.; Zhang, Y.; Wu, Z. Accelerometer-based gait recognition via voting by signature points. Electron. Lett. 2009, 45, 1116-1118.

15. Avvenuti, M.; Casella, A.; Cesarini, D. Using gait symmetry to virtually align a triaxial accelerometer during running and walking. Electron. Lett. 2013,49, 120-121.

16. Gafurov, D.; Snekkkenes, E. Arm swing as a weak biometric for unobtrusive user authentication. In Proceedings of IEEE International Conference on Intelligent Information Hiding and Multimedia Signal Processing, Harbin, China, 15-17 October 2008; pp. 1080-1087.

17. LeMoyne, R.; Mastroianni, T.; Cozza, M.; Coroian, C.; Grundfest, W. Implementation of an iPhone for characterizing Parkinson's disease tremor through a wireless accelerometer application. In Proceedings of 2010 Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC), Buenos Aires, Argentina, 31 August-4 September 2010; pp. 4954-4958.

18. Chan, H.K.; Zheng, H.; Wang, H.; Gawley, R.; Yang, M.; Sterritt, R. Feasibility study on iPhone accelerometer for gait detection. In Proceedings of 2011 5th International Conference on Pervasive Computing Technologies for Healthcare (PervasiveHealth), Dublin, Ireland, 23-26 May 2011; pp. 184-187.

19. Zandbergen, P.A. Accuracy of iPhone locations: A comparison of assisted GPS, WiFi and cellular positioning. Trans. GIS 2009,13,5-25.

20. Shoemake, K. Quaternions, 1994. Available online: quatut.pdf (accessed on 11 September 2014).

21. Haug, E. Computer Aided Analysis and Optimization of Mechanical System Dynamics; NATO ASI Series; Springer-Verlag: Berlin Heidelberg, Germany, 1984.

22. Working, R. Euler Angles, Quaternions, and Transformation Matrices, 1977. Available online: (accessed on 11 September 2014).

23. Zhang, R.; Jia, H.; Chen, T.; Zhang, Y. Attitude solution for strapdown inertial navigation system based on quaternion algorithm. Opt. Precis. Eng. 2008,16, 1963-1970.

24. Sadeghi, H.; Allard, P.; Prince, F.; Labelle, H. Symmetry and limb dominance in able-bodied gait: A review. Gait Posture 2000,12, 34-45.

25. Diebel, J. Representing attitude: Euler angles, unit quaternions, and rotation vectors. Matrix 2006, 58, 15-16.

© 2014 by the authors; licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution license (

Copyright of Sensors (14248220) is the property of MDPI Publishing and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use.