Accuracy and precision of gait events derived from motion capture in horses during walk and trot

This study aimed to create an evidence base for detection of stance-phase timings from motion capture in horses. The objective was to compare the accuracy (bias) and precision (SD) for five published algorithms for the detection of hoof-on and hoof-off using force plates as the reference standard.     

Predicting muscle forces in gait from EMG signals and musculotendon kinematics

An EMG-driven muscle model for determining muscle force-time histories during gait is presented. The model, based on Hill's equation (1938), incorporates morphological data and accounts for changes in musculotendon length, velocity, and the level of muscle excitation for both concentric and eccentric contractions. Musculotendon kinematics were calculated using three-dimensional cinematography with a model of the musculoskeletal system. Muscle force-length-EMG relations were established from slow isokinetic calibrations. Walking muscle force-time histories were determined for two subjects. Joint moments calculated from the predicted muscle forces were compared with moments calculated using a linked segment, inverse dynamics approach. Moment curve correlations ranged from r = 0.72 to R = 0.97 and the root mean square (RMS) differences were from 10 to 20 Nm. Expressed as a relative RMS, the moment differences ranged from a low of 23% at the ankle to a high of 72% at the hip. No single reason for the differences between the two moment curves could be identified. Possible explanations discussed include the linear EMG-to-force assumption and how well the EMG-to-force calibration represented excitation for the whole muscle during gait, assumptions incorporated in the muscle modeling procedure, and errors inherent in validating joint moments predicted from the model to moments calculated using linked segment, inverse dynamics. The closeness with which the joint moment curves matched in the present study supports using the modeling approach proposed to determine muscle forces in gait.     

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.     

Correlates between kinematics and baropodometric measurements for an integrated in-vivo assessment of the segmental foot function in gait

Baropodometry and multi-segmental foot kinematics are frequently employed to obtain insight into the mechanics of the foot-ground interaction in both basic research and clinical settings. However, nothing hitherto has been reported on the full integration of kinematics with baropodometric parameters, and only a few studies have addressed the association between intersegmental kinematics and plantar loading within specific foot regions. The aim of this study was to understanding the relationships between foot joint mobility and plantar loading by focusing on the correlation between these two measures.     

In vivo tibiofemoral cartilage deformation during the stance phase of gait

The knowledge of articular cartilage contact biomechanics in the knee joint is important for understanding the joint function and cartilage pathology. However, the in vivo tibiofemoral articular cartilage contact biomechanics during gait remains unknown. The objective of this study was to determine the in vivo tibiofemoral cartilage contact biomechanics during the stance phase of treadmill gait. Eight healthy knees were magnetic resonance (MR) scanned and imaged with a dual fluoroscopic system during gait on a treadmill. The tibia, femur and associated cartilage were constructed from the MR images and combined with the dual fluoroscopic images to determine in vivo cartilage contact deformation during the stance phase of gait. Throughout the stance phase of gait, the magnitude of peak compartmental contact deformation ranged between 7% and 23% of the resting cartilage thickness and occurred at regions with thicker cartilage. Its excursions in the anteroposterior direction were greater in the medial tibiofemoral compartment as compared to those in the lateral compartment. The contact areas throughout the stance phase were greater in the medial compartment than in the lateral compartment. The information on in vivo tibiofemoral cartilage contact biomechanics during gait could be used to provide physiological boundaries for in vitro testing of cartilage. Also, the data on location and magnitude of deformation among non-diseased knees during gait could identify where loading and later injury might occur in diseased knees.     

Talocrural and talocalcaneal joint kinematics and kinetics during the stance phase of walking

The purpose of this investigation was to study the kinematics and kinetics of the joints between the leg and calcaneus during the stance phase of walking. The talocrural and talocalcaneal joints were each assumed to act as monocentric single degree of freedom hinge joints. Motion at one joint was defined by the relative rotation of a point on the opposing joint. The results, based upon the gait of three subjects, showed that the hinge joint assumption may be reasonable. A discrepancy in the kinematics was shown between the talocrural joint rotation and its commonly assumed sagittal plane representation, especially during initial flatfoot. This discrepancy is due to the fact that the sagittal plane rotation is created by the combined rotations of the talocrural and talocalcaneal joints. The talocalcaneal joint showed a peak 25-30 Nm supinatory moment at 80% of stance. The talocrural joint moment was qualitatively similar to the commonly measured sagittal plane moment, but the present results show that the sagittal plane moment overpredicted the true moment by 6-22% due to the two-dimensional assumption.     

Effects of bilateral medial knee osteoarthritis on intra- and inter-limb contributions to body support during gait

Patients with knee OA show altered gait patterns, affecting their quality of living. The current study aimed to quantify the effects of bilateral knee OA on the intra-limb and inter-limb sharing of the support of the body during gait. Fifteen patients with mild, 15 with severe bilateral knee OA, and 15 healthy controls walked along a walkway while the kinematic and kinetic data were measured. Compared with the controls, the patients significantly reduced their knee extensor moments and the corresponding contributions to the total support moment in the sagittal plane (p<0.05). For compensation, the mild OA group significantly increased the hip extensor moments (p<0.05) to maintain close-to-normal support and a more symmetrical inter-limb load-sharing during double-limb support. The severe OA group involved compensatory actions of both the ankle and hip, but did not succeed in maintaining a normal sagittal total support moment during late stance, nor a symmetrical inter-limb load-sharing during double-limb support. In the frontal plane, the knee abductor moments and the corresponding contributions to the total support moment were not affected by the changes in the other joints, regardless of the severity of the disease. The observed compensatory changes suggest that strengthening of weak hip muscles is essential for body support during gait in patients with knee OA, but that training of weak ankle muscles may also be needed for patients with severe knee OA.     

Kinematic evaluation of cruciate-retaining total knee replacement patients during level walking: A comparison with the displacement-controlled ISO standard

Differences between wear-scar features of simulator-tested and retrieved tibial total knee replacement (TKR) liners have been reported. This disagreement may result from differences between in vivo kinematic profiles and those defined by the International Organization for Standardization (ISO). The purpose of this study was to determine the knee kinematics of a TKR subject group during level walking and compare them with the motion profiles defined by the ISO standard for a displacement-controlled knee wear testing simulator. Twenty-nine patients with a posterior cruciate ligament-retaining TKR design were gait tested using the point cluster technique to obtain flexion–extension (FE) rotation, anterior–posterior (AP) translation and internal–external (IE) rotation knee motions during a complete cycle of level walking. Relative ranges of motion and timing of key points within the in vivo motion data were compared against the same ranges and same key points from the input profiles of the displacement-controlled wear testing standard ISO14243-3. The subjects exhibited a FE pattern similar to ISO, with an insignificant difference in range of FE rotation from midstance to terminal stance. However, the subjects had a significantly higher range of knee flexion from terminal stance into swing. The subjects also exhibited a phase delay for the entire gait cycle. For AP translation, the standard profile had statistically significant lower magnitudes than seen in vivo. Opposite pattern of AP motion was also apparent from midstance and swing. Similarly, ISO specified a smaller IE total range of rotation with a motion pattern in complete opposition to that seen in vivo. In conclusion, significant differences were found in both the magnitudes and pattern of in vivo motion compared with ISO.     

Bilateral knee osteoarthritis does not affect inter-joint coordination in older adults with gait deviations during obstacle-crossing

Fifteen elderly subjects with bilateral medial knee osteoarthritis (OA) and 15 healthy elderly subjects walked and crossed obstacles with heights of 10%, 20%, and 30% of their leg lengths while sagittal angles and angular velocities of each joint were measured and their phase angles () calculated. Continuous relative phase (CRP) were also obtained, i.e., _hip−knee and _knee−ankle. The standard deviations of the CRP curve points were averaged to obtain deviation phase (DP) values for the stance and swing phases. Significant differences between the OA and control groups were found in several of the peak and crossing angles, and angular velocities at the knee and ankle. Both groups had similar CRP patterns, and the DP values of the hip–knee and knee–ankle CRP curves were not significantly different between the two groups. Despite significant changes in the joint kinematics, knee OA did not significantly change the way the motions of the lower limb joints are coordinated during obstacle-crossing. It appears that the OA groups adopted a particular biomechanical strategy among all possible strategies that can accommodate the OA-induced changes of the knee mechanics using unaltered inter-joint coordination control. This enabled the OA subjects to accommodate reliably the mechanical demands related to bilateral knee OA in the sagittal plane during obstacle-crossing. Maintaining normal and reliable inter-joint coordination may be considered a goal of therapeutic intervention, and the patterns and variability of inter-joint coordination can be used for the evaluation of treatment effects.     

Verification of accuracy and validity of gait phase detection system using motion sensors for applying walking assistive FES

In this study, we have analysed heel strike (HS) and toe off (TO) of normal individuals and hemiplegic patients, taking advantage of output curves acquired from various sensors, and verified the validity of sensor detection methods and their effectiveness when they were used for hemiplegic gaits. Gait phase detections using three different motion sensors were valid, since they all had reliabilities more than 95%, when compared with foot velocity algorithm. Results showed that the tilt sensor and the gyrosensor could detect gait phase more accurately in normal individuals. Vertical acceleration could detect HS most accurately in hemiplegic patient group A. The gyrosensor could detect HS and TO most accurately in hemiplegic patient groups A and B. The detection of TO from all sensor signals was valid in both the patient groups A and B. However, the vertical acceleration detected HS validly in patient group A and the gyrosensor detected HS validly in patient group B.     

In-vivo 6 degrees-of-freedom kinematics of metal-on-polyethylene total hip arthroplasty during gait

Knowledge of accurate in-vivo 6 degree-of-freedom (6-DOF) kinematics of total hip arthroplasty (THA) during daily activities is critical for improvement of longevity of the components. Previous studies assessed in-vivo THA kinematics using skin marker-based motion analysis. However, skin markers are prone to move with respect to the underlying bones. A non-invasive dual fluoroscopic imaging system (DFIS) based tracking technique has been used to avoid skin artifacts and provide accurate 6-DOF kinematic measurement. This study aimed to quantify in-vivo 6-DOF THA kinematics during gait using DFIS. Twenty eight well-functioning THAs were evaluated during treadmill gait under DFIS surveillance. The maximum translations of the femoral head were 0.46±0.10 mm and 0.45±0.10 mm during the stance and swing phases (p=0.57), respectively. The range of hip flexion was from 8.7° to 47.6°, adduction from 3.0° to 12.5° and external rotation from 19.2° to 29.7°. The THA was flexed, externally rotated and adducted throughout the gait. The magnitudes of the femoral head translations were found to be within the manufacture tolerance of the components, suggesting that in-vivo hip “pistoning” during gait cycle may be minimal in well-functioning THAs. The 6-DOF kinematics could be used as the baseline knowledge for further improvement of wear-testing of hip implant, implants manufacturing and implant positioning during surgery.     

In vivo biomechanical behavior of the human heel pad during the stance phase of gait

A technique is introduced for simultaneous measurements of the heel pad tissue deformation and the heel–ground contact stresses developing during the stance phase of gait. Subjects walked upon a gait platform integrating the contact pressure display optical method for plantar pressure measurements and a digital radiographic fluoroscopy system for skeletal and soft tissue motion recording. Clear images of the posterior-plantar aspect of the calcaneus and enveloping soft tissues were obtained simultaneously with the pressure distribution under the heel region throughout the stance phase of gait. The heel pad was shown to undergo a rapid compression during initial contact and heel strike, reaching a strain of 0.39±0.05 in about 150 ms. The stress–strain relation of the heel pad was shown to be highly non-linear, with a compression modulus of 105±11 kPa initially and 306±16 kPa at 30% strain. The energy dissipation during heel strike was evaluated to be 17.8±0.8%. The present technique is useful for biomechanical as well as clinical evaluation of the stress–strain and energy absorption characteristics of the heel pad in vivo, during natural gait.     

Increased hip internal abduction moment and reduced speed are the gait strategies used by women with knee osteoarthritis

The purpose of this study was to identify the gait strategies in women with mild and moderate knee osteoarthritis (OA). Forty women diagnosed with OA of the knee and 40 healthy women participated in the study. Toe-out progression angle, trunk lateral lean, hip internal abduction moment and gait speed were measured using Qualisys ProReflex System and two force plates. Principal component analysis was applied to extract features from the gait waveforms data that characterized the waveforms main modes of temporal variation. Discriminant analysis with a stepwise model was conducted to determine which strategies could best discriminate groups. According to the discriminant model, the PC2 of the internal abduction moment of the hip and the gait speed were the most discriminatory variables between the groups. The OA group showed decreased gait speed, decreased hip internal abduction moment during the loading response phase, and increased hip internal abduction moment during the mid and terminal stance phases. Interventions that may increase hip internal abduction moment, such as the strengthening of the hip abductors muscles, may benefit women with knee OA. Training slower than normal gait speeds must be considered in light of potential adverse implications on overall physical function, daily tasks, and safety.     

Influence of gait speed on the control of mediolateral dynamic stability during gait initiation

This study investigated the influence of gait speed on the control of mediolateral dynamic stability during gait initiation. Thirteen healthy young adults initiated gait at three self-selected speeds: Slow, Normal and Fast. The results indicated that the duration of anticipatory postural adjustments (APA) decreased from Slow to Fast, i.e. the time allocated to propel the centre of mass (COM) towards the stance-leg side was shortened. Likely as an attempt at compensation, the peak of the anticipatory centre of pressure (COP) shift increased. However, COP compensation was not fully efficient since the results indicated that the mediolateral COM shift towards the stance-leg side at swing foot-off decreased with gait speed. Consequently, the COM shift towards the swing-leg side at swing heel-contact increased from Slow to Fast, indicating that the mediolateral COM fall during step execution increased as gait speed rose. However, this increased COM fall was compensated by greater step width so that the margin of stability (the distance between the base-of-support boundary and the mediolateral component of the “extrapolated centre of mass”) at heel-contact remained unchanged across the speed conditions. Furthermore, a positive correlation between the mediolateral extrapolated COM position at heel-contact and step width was found, indicating that the greater the mediolateral COM fall, the greater the step width. Globally, these results suggest that mediolateral APA and step width are modulated with gait speed so as to maintain equivalent mediolateral dynamical stability at the time of swing heel-contact.     

Prediction of ground reaction forces during gait based on kinematics and a neural network model

Kinetic information during human gait can be estimated with inverse dynamics, which is based on anthropometric, kinematic, and ground reaction data. While collecting ground reaction data with a force plate is useful, it is costly and requires regulated space. The goal of this study was to propose a new, accurate methodology for predicting ground reaction forces (GRFs) during level walking without the help of a force plate. To predict GRFs without a force plate, the traditional method of Newtonian mechanics was used for the single support phase. In addition, an artificial neural network (ANN) model was applied for the double support phase to solve statically indeterminate structure problems. The input variables of the ANN model, which were selected to have both dependency and independency, were limited to the trajectory, velocity, and acceleration of the whole segment's mass centre to minimise errors. The predicted GRFs were validated with actual GRFs through a ten-fold cross-validation method, and the correlation coefficients (R) for the ground forces were 0.918 in the medial–lateral axis, 0.985 in the anterior–posterior axis, and 0.991 in the vertical axis during gait. The ground moments were 0.987 in the sagittal plane, 0.841 in the frontal plane, and 0.868 in the transverse plane during gait. The high correlation coefficients(R) are due to the improvement of the prediction rate in the double support phase. This study also proved the possibility of calculating joint forces and moments based on the GRFs predicted with the proposed new hybrid method. Data generated with the proposed method may thus be used instead of raw GRF data in gait analysis and in calculating joint dynamic data using inverse dynamics.     

Muscular coordination of knee motion during the terminal-swing phase of normal gait

Children with cerebral palsy often walk with diminished knee extension during the terminal-swing phase, resulting in a troublesome "crouched" posture at initial contact and a shortened stride. Treatment of this gait abnormality is challenging because the factors that extend the knee during normal walking are not well understood, and because the potential of individual muscles to limit terminal-swing knee extension is unknown. This study analyzed a series of three-dimensional, muscle-driven dynamic simulations to quantify the angular accelerations of the knee induced by muscles and other factors during swing. Simulations were generated that reproduced the measured gait dynamics and muscle excitation patterns of six typically developing children walking at self-selected speeds. The knee was accelerated toward extension in the simulations by velocity-related forces (i.e., Coriolis and centrifugal forces) and by a number of muscles, notably the vasti in mid-swing (passive), the hip extensors in terminal swing, and the stance-limb hip abductors, which accelerated the pelvis upward. Knee extension was slowed in terminal swing by the stance-limb hip flexors, which accelerated the pelvis backward. The hamstrings decelerated the forward motion of the swing-limb shank, but did not contribute substantially to angular motions of the knee. Based on these data, we hypothesize that the diminished knee extension in terminal swing exhibited by children with cerebral palsy may, in part, be caused by weak hip extensors or by impaired hip muscles on the stance limb that result in abnormal accelerations of the pelvis.     

In vivo three-dimensional intervertebral kinematics of the subaxial cervical spine during seated axial rotation and lateral bending via a fluoroscopy-to-CT registration approach

Accurate measurement of the coupled intervertebral motions is helpful for understanding the etiology and diagnosis of relevant diseases, and for assessing the subsequent treatment. No study has reported the in vivo, dynamic and three-dimensional (3D) intervertebral motion of the cervical spine during active axial rotation (AR) and lateral bending (LB) in the sitting position. The current study fills the gap by measuring the coupled intervertebral motions of the subaxial cervical spine in ten asymptomatic young adults in an upright sitting position during active head LB and AR using a volumetric model-based 2D-to-3D registration method via biplane fluoroscopy. Subject-specific models of the individual vertebrae were derived from each subject’s CT data and were registered to the fluoroscopic images for determining the 3D poses of the subaxial vertebrae that were used to obtain the intervertebral kinematics. The averaged ranges of motion to one side (ROM) during AR at C3/C4, C4/C5, C5/C6, and C6/C7 were 4.2°, 4.6°, 3.0° and 1.3°, respectively. The corresponding values were 6.4°, 5.2°, 6.1° and 6.1° during LB. Intervertebral LB (ILB) played an important role in both AR and LB tasks of the cervical spine, experiencing greater ROM than intervertebral AR (IAR) (ratio of coupled motion (IAR/ILB): 0.23–0.75 in LB, 0.34–0.95 in AR). Compared to the AR task, the ranges of ILB during the LB task were significantly greater at C5/6 (p=0.008) and C6/7 (p=0.001) but the range of IAR was significantly smaller at C4/5 (p=0.02), leading to significantly smaller ratios of coupled motions at C4/5 (p=0.0013), C5/6 (p<0.001) and C6/7 (p=0.0037). The observed coupling characteristics of the intervertebral kinematics were different from those in previous studies under discrete static conditions in a supine position without weight-bearing, suggesting that the testing conditions likely affect the kinematics of the subaxial cervical spine. While C1 and C2 were not included owing to technical limitations, the current results nonetheless provide baseline data of the intervertebral motion of the subaxial cervical spine in asymptomatic young subjects under physiological conditions, which may be helpful for further investigations into spine biomechanics.     

Kinematic motion of the anterior cruciate ligament deficient knee during functionally high and low demanding tasks

The purpose of this study was to determine whether mechanical adaptations were present in patients with anterior cruciate ligament (ACL)-deficient knees during high-demand activities. Twenty-two subjects with unilateral ACL deficiency (11 males and 11 females, 19.6 months after injury) performed five different activities at a comfortable speed (level walking, ascending and descending steps, jogging, jogging to a 90-degree side cutting toward the opposite direction of the tested side). Three-dimensional knee kinematics for the ACL-deficient knees and uninjured contralateral knees were evaluated using the Point Cluster Technique. There was no significant difference in knee flexion angle, but an offset toward the knee in less valgus and more external tibial rotation was observed in the ACL-deficient knee. The tendency was more obvious in high demand motions, and a significant difference was clearly observed in the side cutting motions. These motion patterns, with the knee in less valgus and more external tibial rotation, are proposed to be an adaptive movement to avoid pivot shift dynamically, and reveal evidence in support of a dynamic adaptive motion occurring in ACL-deficient knees.     

Comparison of hierarchical and six degrees-of-freedom marker sets in analyzing gait kinematics

The objective of this study was to determine how marker spacing, noise, and joint translations affect joint angle calculations using both a hierarchical and a six degrees-of-freedom (6DoF) marker set. A simple two-segment model demonstrates that a hierarchical marker set produces biased joint rotation estimates when sagittal joint translations occur whereas a 6DoF marker set mitigates these bias errors with precision improving with increased marker spacing. These effects were evident in gait simulations where the 6DoF marker set was shown to be more accurate at tracking axial rotation angles at the hip, knee, and ankle.