Kalman smoothing improves the estimation of joint kinematics and kinetics in marker-based human gait analysis.

We developed a Kalman smoothing algorithm to improve estimates of joint kinematics from measured marker trajectories during motion analysis. Kalman smoothing estimates are based on complete marker trajectories. This is an improvement over other techniques, such as the global optimisation method (GOM), Kalman filtering, and local marker estimation (LME), where the estimate at each time instant is only based on part of the marker trajectories. We applied GOM, Kalman filtering, LME, and Kalman smoothing to marker trajectories from both simulated and experimental gait motion, to estimate the joint kinematics of a ten segment biomechanical model, with 21 degrees of freedom. Three simulated marker trajectories were studied: without errors, with instrumental errors, and with soft tissue artefacts (STA). Two modelling errors were studied: increased thigh length and hip centre dislocation. We calculated estimation errors from the known joint kinematics in the simulation study. Compared with other techniques, Kalman smoothing reduced the estimation errors for the joint positions, by more than 50% for the simulated marker trajectories without errors and with instrumental errors. Compared with GOM, Kalman smoothing reduced the estimation errors for the joint moments by more than 35%. Compared with Kalman filtering and LME, Kalman smoothing reduced the estimation errors for the joint accelerations by at least 50%. Our simulation results show that the use of Kalman smoothing substantially improves the estimates of joint kinematics and kinetics compared with previously proposed techniques (GOM, Kalman filtering, and LME) for both simulated, with and without modelling errors, and experimentally measured gait motion.

[1]  G Ferrigno,et al.  Robust recovery of human motion from video using Kalman filters and virtual humans. , 2003, Human movement science.

[2]  J. J. Moré,et al.  Levenberg--Marquardt algorithm: implementation and theory , 1977 .

[3]  D. Sutherland The evolution of clinical gait analysis. Part II kinematics. , 2002, Gait & posture.

[4]  A Leardini,et al.  Position and orientation in space of bones during movement: experimental artefacts. , 1996, Clinical biomechanics.

[5]  T. Başar,et al.  A New Approach to Linear Filtering and Prediction Problems , 2001 .

[6]  J J O'Connor,et al.  Bone position estimation from skin marker co-ordinates using global optimisation with joint constraints. , 1999, Journal of biomechanics.

[7]  F.E. Zajac,et al.  An interactive graphics-based model of the lower extremity to study orthopaedic surgical procedures , 1990, IEEE Transactions on Biomedical Engineering.

[8]  F. Zajac,et al.  A planar model of the knee joint to characterize the knee extensor mechanism. , 1989, Journal of biomechanics.

[9]  J. Miller Numerical Analysis , 1966, Nature.

[10]  Yaakov Bar-Shalom,et al.  Estimation and Tracking: Principles, Techniques, and Software , 1993 .

[11]  A. Cappozzo,et al.  Human movement analysis using stereophotogrammetry. Part 2: instrumental errors. , 2004, Gait & posture.

[12]  Michael O. Kolawole,et al.  Estimation and tracking , 2002 .

[13]  G Ferrigno,et al.  Kinematical models to reduce the effect of skin artifacts on marker-based human motion estimation. , 2005, Journal of biomechanics.

[14]  A. Cappozzo,et al.  Human movement analysis using stereophotogrammetry. Part 3. Soft tissue artifact assessment and compensation. , 2005, Gait & posture.

[15]  C. Spoor,et al.  Rigid body motion calculated from spatial co-ordinates of markers. , 1980, Journal of biomechanics.

[16]  A. Cappozzo,et al.  Human movement analysis using stereophotogrammetry. Part 1: theoretical background. , 2005, Gait & posture.

[17]  G. Ferrigno,et al.  Real-time human motion estimation using biomechanical models and non-linear state-space filters , 2006, Medical and Biological Engineering and Computing.

[18]  C. Striebel,et al.  On the maximum likelihood estimates for linear dynamic systems , 1965 .