Bilateral step length estimation using a single inertial measurement unit attached to the pelvis
 Alper Köse^{1},
 Andrea Cereatti^{1} and
 Ugo Della Croce^{1}Email author
DOI: 10.1186/1743000399
© Köse et al; licensee BioMed Central Ltd. 2012
Received: 16 August 2011
Accepted: 8 February 2012
Published: 8 February 2012
Abstract
Background
The estimation of the spatiotemporal gait parameters is of primary importance in both physical activity monitoring and clinical contexts. A method for estimating step length bilaterally, during level walking, using a single inertial measurement unit (IMU) attached to the pelvis is proposed. In contrast to previous studies, based either on a simplified representation of the human gait mechanics or on a general linear regressive model, the proposed method estimates the step length directly from the integration of the acceleration along the direction of progression.
Methods
The IMU was placed at pelvis level fixed to the subject's belt on the right side. The method was validated using measurements from a stereophotogrammetric system as a gold standard on nine subjects walking ten laps along a closed loop track of about 25 m, varying their speed. For each loop, only the IMU data recorded in a 4 m long portion of the track included in the calibrated volume of the SP system, were used for the analysis. The method takes advantage of the cyclic nature of gait and it requires an accurate determination of the foot contact instances. A combination of a Kalman filter and of an optimally filtered direct and reverse integration applied to the IMU signals formed a single novel method (Kalman and Optimally filtered Step length Estimation  KOSE method). A correction of the IMU displacement due to the pelvic rotation occurring in gait was implemented to estimate the step length and the traversed distance.
Results
The step length was estimated for all subjects with less than 3% error. Traversed distance was assessed with less than 2% error.
Conclusions
The proposed method provided estimates of step length and traversed distance more accurate than any other method applied to measurements obtained from a single IMU that can be found in the literature. In healthy subjects, it is reasonable to expect that, errors in traversed distance estimation during daily monitoring activity would be of the same order of magnitude of those presented.
Keywords
Inertial measurement Gait analysis Gait parameters Accelerometer Step length Gait monitoring Stride length Inertial sensor Wearable.Background
The measurement of temporal and spatial features of gait is essential for the assessment of gait abnormalities and the quantitative evaluation of treatment outcomes [1]. In particular, amplitude, variability and asymmetry of step length (SL) have been shown to be effective outcomes of walking ability. In fact, they are strongly related to the propulsion generation and can be representative of the compensatory mechanisms adopted in pathological walking [2, 3]. Having access to instruments capable of gathering information about the patient walking ability outside the laboratory [4], during daily life with no space limitations and for prolonged periods of time is of paramount importance in numerous clinical applications [5].
Inertial measurement units (IMU) are strong candidates for these applications. Those IMUs, featuring 3axis accelerometers and gyroscopes [6], can be employed to estimate the SL during walking. An estimate of the SL can be obtained by double integrating the IMU acceleration component in the direction of progression (DoP) between the instants of two consecutive heel strikes (HS).
 a)
the identification of the foot contact time instances (gait events);
 b)
 c)
 d)
the estimation of the initial velocity values for the integration of the acceleration along the DoP [12].
The most straightforward solution to determine both right and left SL (rSL and lSL) would be to place an IMU on each foot so that velocity and the orientation of each IMU can be set to zero at the beginning of the integration interval [7, 13, 14].
However, when the focus is on the description of the individual motor capacity ("can do in daily environment") and performance ("does do in daily environment") [15], the requirement of a light, small and minimally obstructive setup is of primary importance. Therefore, it would be desirable to obtain the same information using a single, discomfortfree IMU.
To the authors' knowledge, all the studies in the literature, based upon the use of a single IMU, determined the SL through indirect estimation methodologies [16–18].
Zijlstra and Hof (2003) [16] proposed a method for estimating spatialtemporal parameters of gait from trunk accelerations. The IMU was placed on the back at the level of the S2 vertebrae. The method used zero crossing of the forward accelerations for detecting foot contact instances and an inverted pendulum model to estimate the SL. However, in their method, left and right foot contacts identification failed for 6 out of 15 subjects, 12% of the times and the SLs were underestimated in all subjects and at all speeds.
Gonzales and colleagues (2007) [17] proposed a modified version of the pendulum model for taking into account the single stance and doublestance, separately. The IMU was placed on the back at the L3 vertebral level. Traversed distance estimation ranged between 94.5% to 106.7% of the actual traversed distance. No information was provided about the errors associated to the SL estimates.
Shin and Park [18] determined the SL using a linear combination of the walking frequency and the variance of the accelerometric signals recorded by an IMU attached to the subject's waist belt. Experiments were carried out on a single subject. The accuracy of the SL estimation is not reported, but the accuracy of traversed distance resulted to be of about 96%.
The abovementioned methods for the SL estimate are based either on a simplified representation of the human gait mechanics (inverted pendulum model) or on a general linear regressive model. Hence, the SL estimates provided are expected to be affected by the errors intrinsic to the specific model formulation.
Conversely, this study presents a method for the estimation of rSL and lSL during level walking using a single IMU attached to the pelvis without using gait models. In line with previous approaches, the proposed method requires the identification of the gait events as a preliminary step. To minimize the detrimental drift effects, the IMU acceleration along the DoP is double integrated by means of a combination of direct and reverse integrations [19] of an optimally filtered acceleration.
We hypothesized that the proposed method for the estimation of the SL, would be more accurate than methods found in the literature based on the use of inverted pendulum or regressive models.
The performance of the proposed approach was evaluated on data collected on healthy subjects walking at various speeds, using stereophotogrammetric (SP) measurements as a gold standard.
Methods
Instrumentation
An IMU (FreeSense, Sensorize^{®}) featuring a triaxial accelerometer and two biaxial gyroscopes (acceleration resolution 0.0096 m/s^{2}, angular rate resolution 0.2441 deg/s, unit weight 93 g, unit size 85 mm × 49 mm × 21 mm) sampling data at 100 frames/s was used. A 10camera BTS^{®} SMARTD stereophotogrammetric system acquiring at 100 frames/s, (positional accuracy 0.3 mm) was used for validation purposes.
Step length estimation
The method proposed in this paper is based on the assumption that the pelvic displacement along the DoP, between the instants of contralateral HS and the first following HS, is equal to the ipsilateral step length.
 a)
the identification of gait events;
 b)
the determination of pelvis displacement along the DoP.
Identification of gait events
A gait cycle begins when a foot hits the ground (heel strike  HS) and ends when the next HS of the same foot occurs. When gait cycles can be identified for both sides, then also steps can be identified. A step begins when a contralateral HS occurs and ends at the first following HS.
A waveletbased approach was developed to identify intervals of interest where gait event candidates were searched in the accelerometer signals. The X_{ L } and Y_{ L } accelerometer signals were decomposed using a "Stationary wavelet decomposition" [20]. A Daubechies level 5 ("db5") mother wavelet was chosen given its similarity to the IMU signals in the proximity of HS. The original signals were then decomposed in an approximation curve plus ten levels of detail. Thresholds were applied to the first three detail levels and the other detail levels were discarded. Thresholds for these levels were 1/5, 1/4 and 1/3 of its magnitude for the first, second and third level, respectively. The signals were reconstructed using only the first three levels of detail after thresholds were applied. An interval of interest for the accelerometer signals was defined as the interval of time during which the reconstructed signals differed from zero (Figure 2b and Figure 2c). The intervals of interest allow to identify the regions of the signals with localized high frequency components where rHS and lTO are expected to occur.
Based on the preliminary visual investigation, the right HS (rHS) was detected as the instant of time in the middle between the maximum of the Y_{ L } accelerometer signal and the first minimum of the accelerometer signal along the X_{ L } accelerometer signal in the corresponding intervals of interest. The identification of the left HS (lHS) required the identification of both right and left toe off instances (rTO, lTO). The lTO was detected as the instant of time in the middle between the first maximum of the Y_{ L } accelerometer signal after the rHS and the second minimum of the X_{ L } accelerometer signal in the corresponding intervals of interest. The rTO was found as the time of the minimum negative peak value between two consecutive lTO and rHS. The lHS was found as the first local maximum of the Z_{ L } accelerometer signal before the rTO (Figure 2d) [21].
Estimation of the pelvis displacement along the DoP
 a)
estimation of the IMU acceleration along the DoP;
 b)
integration of the IMU acceleration along the DoP;
 c)
removal of pelvic rotation contribution to IMU displacement.
Estimation of the IMU acceleration along the DoP
To estimate the IMU acceleration along the DoP, the orientation of the LF with respect to the global frame (GF), was estimated using a specifically designed Kalman filter based on both three dimensional acceleration and angular velocity vectors. For a detailed explanation of how the Kalman filter was implemented please refer to Mazzà and colleagues [22].
The GF was defined as follows: the X_{ G } axis coinciding with the direction of gravity, the Y_{ G } axis coinciding with the DoP during level straight walking, and the Z_{ G } axis resulting from the cross product between X_{ G } and Y_{ G } . To define the Y_{ G } axis direction, at the beginning of each acquisition, the IMU Y_{ L } axis was aligned to the DoP of the level straight walking.
The orientation of the LF with respect to the GF at the i^{ th } instant of time was expressed using the orientation quaternion ^{ G }q_{ L }(i).
Integration of the IMU acceleration along the DoP
To obtain the velocity and displacement time series along the DoP, an integration technique, the Optimally Filtered Direct and Reverse Integration (OFDRI), and its simplified version, the Optimally Filtered Integration (OFI) [19], was adapted to manage acceleration signals during gait. Both the OFI and the OFDRI were originally designed for step negotiation motor tasks [19] and require the knowledge of the final value of the integral to set a cut off frequency for the high pass filter employed to reduce the effect of the drift in the accelerometer signals. In this application, the cut off frequency was determined from a gait cycle for which the initial and final velocity of the IMU was assumed to be equal. The resulting cutoff frequency was then applied for filtering the acceleration signal along the DoP (a_{ x } (i)), one gait cycle at a time. Since the trials started with the subject standing, the initial velocity along the DoP was set to zero. The velocity values found for the final HS were used as initial velocity values for the integration of the following gait cycle.
For gait cycles in which the velocity of the final HS was within a tolerance (± ε, ε = 0.3 m/s) of the initial value then it was forced to be exactly equal to it and the OFDRI was applied. The same updated velocity value was used also as initial velocity of the following gait cycle. The value ε = 0.3 m/s was chosen based on a trial and error approach.
The integration of velocity along the DoP (v_{ x } (i)) to obtain displacement, was performed using the OFI. Since, for each cycle, the final value of the integral (the displacement at the final HS) is the unknown quantity, the OFDRI cannot be applied.
where s_{ x } is the displacement along the DoP.
The SL estimation method including the estimation and integration of the IMU acceleration along the DoP was named KOSE (Kalman and Optimal based filtering Step length Estimation).
Removal of pelvic rotation contribution to IMU displacement
When the IMU is located on the right side of the pelvis, the IMU displacement along the DoP, at the end of the right step duration (rT, the time interval between a lHS and the following rHS), is larger than the actual rSL due to the pelvic angular displacement θ, about the vertical axis of the GF, occurring during the same interval of time. Conversely, the IMU displacement at the end of the left step duration (lT, the time interval between the rHS and the lHS) is smaller than the lSL. Therefore, applying the equations (2) without taking into account the pelvic rotation, would result in an overestimate of the rSL and an underestimate of the lSL (vice versa if the IMU is attached to the left side of the pelvis). The amount of the over [under] estimate would be equal to d sin(rθ_{ j } ) [d sin(lθ_{ j } )], where d is the inter ASIS (anterior superior iliac spine) distance, and rθ_{ j } [lθ_{ j } ] represents the pelvic angular displacement about the vertical axis between the beginning and the end of the j^{ th } right [left step]. Both rθ_{ j } and lθ_{ j } are obtained from the IMU orientation provided by the Kalman filter.
when the IMU is on the left side of the pelvis.
Experimental session
A 4 m long portion of the walking track was included in the calibrated volume of the SP system, (Figure 3a) with its reference frame made to coincide with the GF. For each lap, SP data of three consecutive gait cycles were recorded.
Only data recorded by the IMU and SP system within the calibrated volume were used for the method validation. Data recorded during each lap were arranged together to obtain a continuous data set formed by 3 × 10 gait cycles (Figure 3c).
Data analysis
The estimates of the timing of the gait events (rHS, lHS), of the right and left step duration (rT, lT) and of rSL and lSL, obtained from the heel marker coordinates reconstructed with the SP system were used as gold standard measurements when evaluating the accuracy of the estimates obtained using the IMU [23].
The same processing was applied to the left side.
Descriptive statistics of the above errors were computed.
In addition, to compare our results with previous works, the difference between the IMU estimated distance and the actual traversed distance was computed for each subject.
Results
Heel strike detection
Heel strike detection error
[s]  e _{ rHSj }  e _{ lHSj }  

subject  mean  s.d.  mean  s.d 
1  0.019  0.027  0.026  0.019 
2  0.008  0.021  0.017  0.016 
3  0.026  0.017  0.034  0.021 
4  0.012  0.007  0.027  0.025 
5  0.021  0.031  0.029  0.019 
6  0.016  0.013  0.028  0.024 
7  0.024  0.022  0.031  0.017 
8  0.018  0.02  0.028  0.016 
9  0.009  0.019  0.019  0.025 
average  0.017  0.020  0.027  0.020 
Right and left step duration
Step duration error
[s]  e _{ rTj }  e _{ lTj }  

subject  mean  s.d.  %  mean  s.d.  % 
1  0.021  0.031  3.10%  0.024  0.025  2.60% 
2  0.015  0.026  2.80%  0.011  0.017  2.60% 
3  0.014  0.017  2.10%  0.023  0.019  2.70% 
4  0.009  0.013  1.70%  0.017  0.017  2.90% 
5  0.014  0.022  1.90%  0.011  0.023  1.60% 
6  0.011  0.021  1.50%  0.016  0.019  2.00% 
7  0.022  0.027  3.20%  0.019  0.025  2.70% 
8  0.017  0.021  2.50%  0.024  0.025  3.20% 
9  0.01  0.018  1.50%  0.009  0.018  1.40% 
average  0.015  0.022  2.26%  0.017  0.021  2.41% 
Right and left SL estimation
Cutoff frequencies employed in both OFI and OFDRI to highpass filter the original data varied, across the subjects, between the 0.067 Hz and 0.096 Hz. The pelvic rotation about the GF vertical axis was on average across subjects about 5 degrees.
Step length error
[m]  e _{ rSLj }  e _{ lSLj }  

subject  mean  s.d.  %  mean  s.d.  % 
1  0.014  0.016  2.0%  0.010  0.014  1.9% 
2  0.020  0.015  2.9%  0.017  0.013  2.6% 
3  0.006  0.015  0.8%  0.002  0.016  0.3% 
4  0.001  0.021  0.2%  0.002  0.019  0.4% 
5  0.017  0.018  2.1%  0.012  0.018  1.5% 
6  0.019  0.019  2.1%  0.019  0.021  2.2% 
7  0.013  0.015  1.7%  0.008  0.014  1.1% 
8  0.008  0.017  1.2%  0.001  0.019  0.2% 
9  0.003  0.013  0.6%  0.012  0.013  1.4% 
average  0.009  0.017  1.2%  0.008  0.016  1.1% 
Total traversed distance
Total traversed distance
subject [no.]  traversed distance [m]  estimated distance [m]  difference [m]  difference [%] 

1  32.65  33.02  0.37  1.1 
2  30.99  31.46  0.46  1.5 
3  31.17  30.80  0.37  1.2 
4  31.71  32.25  0.54  1.7 
5  30.37  29.91  0.46  1.5 
6  33.44  32.89  0.54  1.6 
7  31.40  31.74  0.34  1.1 
8  32.11  31.78  0.33  1.0 
9  30.63  31.00  0.37  1.2 
average  31.61  31.65  0.04  0.1 
Discussion and Conclusion
A method for determining both right and left SL during level walking using a single IMU to be used either indoor or outdoor, without space limitations and for prolonged periods of time, was presented and validated.
SL estimates were obtained directly from the IMU displacement between two consecutive HS, by employing an original method (the KOSE method) which double integrates the IMU acceleration along the DoP obtained after applying a Kalman filter. The KOSE method includes an adaptation to gait of an optimal integration technique originally proposed by Zok and colleagues [19] to reduce the errors caused by the drift affecting the dynamometric signals.
The method was tested on healthy subjects while they were walking increasing and decreasing their speed, five times along 30 m. The algorithm for the identification of left and right HS never failed in all the subjects analyzed. Both right and left HS were detected with an average time delay corresponding to 1.7 samples and 2.7 samples, respectively (sample frequency 100 frames/s). These errors implied that the right step duration was on average two samples shorter than the left one. Despite the shorter estimated step duration, the rSL resulted, over the different subjects, overestimated by 1.2%. The maximum mean error in estimating the SL was overestimated by 2.9% for the rSL and underestimated by 2.6% for the lSL. The main explanation for the side to side differences, is probably that the effect of the pelvic rotation, associated to the asymmetrical IMU positioning, was not completely compensated for by the correction term d sin(θ_{ j } ). In fact, in this regard, several assumptions were not completely satisfied: being attached to the belt of the subject, the IMU was not rigidly moving with the pelvis and d, which is taken as the interASIS distance may differ from the distance between the LF origin to the pelvis vertical axis of rotation. As expected, the maximum percent error on the total traversed distance was to 1.7% and therefore smaller than the maximum SL percent error. It is worth noticing that the estimate of the traversed distance was never consistently either an overestimate or an underestimate (average over subjects equal to 0.1%).
Several methods to estimate rSL and lSL, using a single IMU, have been proposed in the literature [16–18, 24]. A common feature of these methods is that the SL estimates were derived using indirect methodologies, either based on inverted pendulum models for level human walking representation [25, 26] or based on general regressive equation in which SL is expressed as function of the step frequency and accelerometric signal variations. Unfortunately, no data on the accuracy of the rSL and lSL estimates were provided in abovementioned studies. Therefore, a straightforward comparison among methods is not possible. However, Zijlstra and Hof (2003) [16] reported a consistent SL underestimation. As a consequence, also the mean speed (mean SL divided mean step time) was underestimated and to correct it, they introduced a coefficient of 1.25 heuristically determined. It could then be inferred that the SL estimations were affected by errors of about 2030%. In Gonzales and colleagues (2007) [17], the performance of their method was evaluated by assessing the errors of the distance traversed in each acquisition. Across subjects, errors ranged from 6.5% to +6.0%. Shin and Park [18] tested their method on a single subject walking at a constant speed for a 70 m straight trajectory. Three data set were recorded (slow, normal and fast). Across trials, errors for the traversed distance varied between 1% to 3% with an accuracy error of 3.7% in the worst case.
In our study the errors in estimating the traversed distance ranged from 1.6% to +1.7% across subjects.
The smaller estimation errors compared to previously publish methods, confirm our initial hypothesis that the KOSE method improves the accuracy with respect to indirect methods based on the use of general models, which may not always take into account subject specificity in gait.
In this study, the IMU was attached to the subject's belt on the right side. This location was chosen to reduce subject discomfort and because it is a practical solution for monitoring daily life motor activities. However, it can be hypothesized that attaching the IMU centrally in the back (S2 or L3 vertebral level) similarly to Zijlstra and Hof (2003) [16], Gonzales and colleagues (2007) [17] and Shin and Park [18], the residual errors related to the pelvic rotation correction would be reduced and it would not be necessary to correct for the pelvic rotation effects.
In this study we presented results relative to gait event timings and SL estimations. Additional gait parameters can be easily calculated from the estimated parameters. Therefore, using a single IMU, we were able to provide excellent estimates of both temporal and spatial gait parameters bilaterally.
It is important to stress that the KOSE method was tested on experimental conditions fairly similar to those which can be encountered in the real life: different gait speeds (three gait speed levels), numerous velocity transients (five velocity transients on 30 m). It is reasonable to expect that, in healthy subjects, errors in traversed distance estimation during daily monitoring activity would be of the same order of magnitude of those presented.
The data set, derived from the experiments carried out in the present study, referred to a straight level walking. For this specific method validation, the GF was therefore defined, at the beginning of each acquisition with the Y_{ G } axis coinciding with the DoP. The IMU acceleration component along the DoP, was then calculated by projecting the accelerometric signal, after having removed the gravitational contribution, on the GF. In general, if during the acquisition, the subject varies the DoP more than once, it would be necessary to define for every different DoP, an auxiliary GF with the Y_{ G } axis aligned to the relevant DoP. However, the solution of the latter issue is out of the scope of the present study.
There are some limitations to the current proof of concept which need to be addressed before using the presented methodology to estimate gait spatial and temporal parameters in clinical contexts. One issue is related to the algorithm used to identify the gait events, which is normally based on IMU signal features. Such features may change or even disappear in pathologic subjects. Another issue is related to the assumption that the pelvis displacement along the DoP between HS equals the SL. It is certainly true in a symmetric gait, but in some gait disorders the pelvis displacement may not be in phase with the feet displacement as much as in healthy subjects. Future studies need to be performed to test the method performance on pathological gait [27].
In conclusion, the method proposed allows for an accurate estimate of right and left step length using a minimal invasive experimental set up through an optimal integration procedure of the accelerometric signals.
List of abbreviations
 DoP:

direction of progression
 GF:

global reference frame
 HS:

heel strike
 IMU:

inertial measurement unit
 LF:

local reference frame
 lHS:

left heel strike instant
 lSL:

right step length
 lTO:

left toe off instant
 OFDRI:

optimally filtered direct and reverse integration
 OFI:

optimally filtered integration
 rHS:

right heel strike instant
 rSL:

right step length
 rTO:

right toe off instant
 SL:

step length
 SP:

stereophotogrammetric.
Declarations
Acknowledgements
The authors thank Dr. Marco Donati from the Laboratory of Locomotor Apparatus Bioengineering of the University of Rome "Foro Italico" for providing the code of the Kalman filter.
Authors’ Affiliations
References
 Perry J: Gait Analysis: Normal and pathological function. 1st edition. Thorofare, NJ SLACK Inc.; 1992.
 NanhoeMahabier W, Snijders AH, Delval A, Weerdesteyn V, Duysens J, Overeem S, Bloem BR: Walking patterns in Parkinson's disease with and without freezing of gait. Neuroscience 2011, 182: 217224.View ArticlePubMed
 Allen JL, Kautz SA, Neptune RR: Step length asymmetry is representative of compensatory mechanisms used in poststroke hemiparetic walking. Gait Posture 2011, 33: 538543. 10.1016/j.gaitpost.2011.01.004PubMed CentralView ArticlePubMed
 MoeNilssen , Helbostad : Estimation of gait cycle characteristics by trunk. Journal of Biomechanics 2004,37(1):1216. 10.1016/S00219290(03)002331View ArticlePubMed
 van Dam MS, Kok GJ, Munneke M, Vogelaar FJ, Vliet Vlieland TP, Taminiau AHJ: Measuring physical activity in patients after surgery for a malignant tumour in the leg. The reliability and validity of a continuous ambulatory activity monitor. Bone & Joint Surgery [Br] 2001, 83: 10151019. 10.1302/0301620X.83B7.12195View Article
 Yun X, Bachmann ER, Moore H, Calusdian J: Selfcontained position tracking of human movement using small inertial/magnetic sensor modules. IEEE International Conference on Robotics and Automation 2007, 25262533.
 Sabatini AM, Martelloni C, Scapellato S, Cavallo F: Assessment of walking features from foot inertial sensing. IEEE Transactions on Biomedical Engineering 2005, 52: 486494. 10.1109/TBME.2004.840727View ArticlePubMed
 Schepers HM, Koopman HF, Veltink PH: Ambulatory assessment of ankle and foot dynamics. Transactions on Biomedical Engineering 2007, 54: 895902.View Article
 Woodman OJ: An Introduction to Inertial Navigation. University of Cambridge, Cambridge, UK, InHouse Technical Report 696 2007.
 Djuri Z: Mechanisms of noise sources in microelectromechanical systems. Microelectronic Reliability 2000, 40: 919932. 10.1016/S00262714(00)000044View Article
 Thong YK, Woolfson MS, Crowe JA, HayesGill BR, Jones DA: Numerical double integration of acceleration measurements in noise. Measurement 2004, 36: 7392. 10.1016/j.measurement.2004.04.005View Article
 Peruzzi A, Della Croce U, Cereatti A: Validity of zero velocity update technique for stride length estimation. Journal of Biomechanics 2011, 44: 19914. 10.1016/j.jbiomech.2011.04.035View ArticlePubMed
 Veltink PH, Slycke P, Hemssems J, Buschman R, Bultstra G, Hermens H: Three dimensional inertial sensing of foot movements for automatic tuning of a twochannel implantable dropfoot stimulator. Medical Engineering & Physics 2003, 25: 2128. 10.1016/S13504533(02)000413View Article
 Foxlin E: Pedestrian tracking with SHOEmounted inertial sensors. IEEE Computer Graphics and Applications 2005, 25: 3846. 10.1109/MCG.2005.140View ArticlePubMed
 Holsbeeke L, Ketelaar M, Schoemaker MM, Gorter JW: Capacity, capability, and performance: different constructs or three of a kind? Phys Med Rehabil 2009, 90: 849855. 10.1016/j.apmr.2008.11.015View Article
 Zijlstra W, Hof AL: Assessment of spatiotemporal gait parameters from trunk accelerations during human walking. Gait and Posture 2003, 18: 110.View ArticlePubMed
 González RC, Alvarez D, López AM, Alvarez JC: Modified pendulum model for mean step length estimation. Proc IEEE Eng Med Biol Soc 2007, 13714.
 Shin SH, Park CG: Adaptive step length estimation algorithm using optimal parameters and movement status awareness. Med Eng Phys 2011, 33: 106410671. 10.1016/j.medengphy.2011.04.009View ArticlePubMed
 Zok M, Mazzà C, Della Croce U: Total body centre of mass displacement estimated using ground reactions during transitory motor tasks: application to step ascent. Medical Engineering & Physics 2004, 26: 791798. 10.1016/j.medengphy.2004.07.005View Article
 Pesquet JC, Krim H, Carfatan H: Timeinvariant orthonormal wavelet representations. IEEE Trans. Sign. Proc 1996, 44: 19641970. 10.1109/78.533717View Article
 Kose A, Peruzzi A, Cereatti A, Laudani L, Della Croce U: Detection of heel strikes and toeoffs during gait using a single inertial measurement unit attached to the waist. Proc Second National Congress of Biongineering, Turin, Italy 2010, 233.
 Mazzà C, Donati M, McCamley J, Cappozzo A: Optimization of the Kalman filter parameters for the estimate of trunk angles from inertial sensors data during walking at different speeds. Gait and Posture 2012, 35: 138142. 10.1016/j.gaitpost.2011.08.024View ArticlePubMed
 Hreljac Alan, Marshall RobertN: Algorithms to determine event timing during normal walking using kinematic data. Journal of Biomechanics 2000, 33: 7836. 10.1016/S00219290(00)000142View ArticlePubMed
 Käppi J, Saarinen J, Syrjärinne J: MEMSIMU based pedestrian navigator forhandheld devices. Proc ION GPS, Salt Lake City, UT 2001, 13691373.
 Cavagna GA, Thys H, Zamboni A: The sources of external work in level walking and running. Physiol 1976, 262: 639657.View Article
 Kuo AD, Donelan JM, Ruina A: Energetic consequences of walking like an inverted pendulum: steptostep transitions. Sport Sci Rev 2005, 33: 8897. 10.1097/0000367720050400000006View Article
 Esser P, Dawes H, Collett J, Feltham MG, Howells K: Assessment of spatiotemporal gait parameters using inertial measurement units in neurological populations. Posture 2011, 34: 55860.View Article
Copyright
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.