Bone orientation and position estimation errors using Cosserat point elements and least squares methods: Application to gait

Conclusion or Illusion: Quantifying Uncertainty in Inverse Analyses From Marker-Based Motion Capture due to Errors in Marker Registration and Model Scaling

Estimating kinematics from optical motion capture with skin-mounted markers, referred to as an inverse kinematic (IK) calculation, is the most common experimental technique in human motion analysis. Kinematics are often used to diagnose movement disorders and plan treatment strategies. In many such applications, small differences in joint angles can be clinically significant. Kinematics are also used to estimate joint powers, muscle forces, and other quantities of interest that cannot typically be measured directly. Thus, the accuracy and reproducibility of IK calculations are critical. In this work, we isolate and quantify the uncertainty in joint angles, moments, and powers due to two sources of error during IK analyses: errors in the placement of markers on the model (marker registration) and errors in the dimensions of the model’s body segments (model scaling). We demonstrate that IK solutions are best presented as a distribution of equally probable trajectories when these sources of modeling uncertainty are considered. Notably, a substantial amount of uncertainty exists in the computed kinematics and kinetics even if low marker tracking errors are achieved. For example, considering only 2 cm of marker registration uncertainty, peak ankle plantarflexion angle varied by 15.9°, peak ankle plantarflexion moment varied by 26.6 N⋅m, and peak ankle power at push off varied by 75.9 W during healthy gait. This uncertainty can directly impact the classification of patient movements and the evaluation of training or device effectiveness, such as calculations of push-off power. We provide scripts in OpenSim so that others can reproduce our results and quantify the effect of modeling uncertainty in their own studies.

Estimating soft tissue artifact of the thigh in high knee flexion tasks using optical motion Capture: Implications for marker cluster placement

Soft tissue artifact in motion capture is widely accepted as a significant source of error in kinematic and kinetic measurements. Non-invasive methods of estimating soft tissue artifact, those requiring only motion capture, provide a feasible method to evaluate marker placement on a segment and enable recommendations for marker configurations which can minimize soft tissue artifact. The purpose of this study was to investigate the effect of thigh marker cluster location on soft tissue artifact during high knee flexion (>120 deg) as unique deformation of soft tissue occurs in this range (e.g. thigh-calf contact). Motion of the pelvis and lower limbs were recorded during squatting and kneeling in fifty participants. Six rigid marker clusters were affixed to the skin on the anterior, lateral, and anterolateral aspect, at the distal and middle third of the thighs. To estimate soft tissue artifact, the functional hip joint center was reconstructed relative to the pelvis cluster and each of the six thigh clusters throughout motion. The difference in the position of these two points was input into Bland-Altman analyses and compared between the thigh clusters. Across the tasks, the total mean difference ranged from 2.81 to 8.95 cm while the lower and upper limits of agreement ranged from -0.79 to 2.54 cm and 5.04 to 17.65 cm, respectively. Using this non-invasive method, the mid-anterolateral cluster was least susceptible to soft tissue artifact and thus would be recommended, while the lateral clusters were most susceptible and should avoided in high knee flexion and similar tasks.

A joint kinematics driven model of the pelvic soft tissue artefact

When skin-markers trajectories are used in human movement analysis, compensating for their relative movement with respect to the underlying bone (soft tissue artefact, STA) is essential for accurate bone-pose estimation; information about the artefact is required in the form of a mathematical model. Such model, not available for pelvic artefacts, could allow pelvic STA compensation in routine gait analysis by embedding it in skeletal kinematics estimators and developing ad-hoc optimization problems for the estimate of subject-specific model parameters. It was developed as driven by adjacent body segment kinematics. Model architecture feasibility was tested; its compensation effectiveness was assessed evaluating the error in pelvic orientation after removing the modelled artefact from the measured one. Five volunteers with a wide body mass range (BMI: 22-37) underwent MRI scans to reconstruct subject-specific pelvic digital bone models. Multiple anatomical calibrations performed in different static postures, as occurring during walking and star-arc movements, registering the bone-models with points digitized through stereophotogrammetry over pelvic bony prominences, allowed to define the relevant poses of a pelvis-embedded anatomical coordinate system. Such approach allowed to measure STAs over several pelvic anatomical landmarks, for each posture and subject. Model parameters were estimated by minimizing the least squares difference between measured and modelled STAs. The measured STAs were appropriately modelled with subject-specific calibrations, both in terms of shape (correlation coefficient: median [inter-quartile-range]: 0.72 [0.36]) and amplitude (root mean square residual: 3.0 [3.2] mm). Consequently, the overall error in pelvic orientation vector (5.1 [4.4] deg) was reduced after removing the modelled artefacts (2.5 [1.9] deg).

A Framework for Measuring the Time-Varying Shape and Full-Field Deformation of Residual Limbs Using 3-D Digital Image Correlation

Effective prosthetic socket design following lower limb amputation depends upon the accurate characterization of the shape of the residual limb as well as its volume and shape fluctuations.

OBJECTIVE

This study proposes a novel framework for the measurement and analysis of residual limb shape and deformation, using a high-resolution and low-cost system.

METHODS

A multi-camera system was designed to capture sets of simultaneous images of the entire residuum surface. The images were analyzed using a specially developed open-source three-dimensional digital image correlation (3D-DIC) toolbox, to obtain the accurate time-varying shapes as well as the full-field deformation and strain maps on the residuum skin surface. Measurements on a transtibial amputee residuum were obtained during knee flexions, muscle contractions, and swelling upon socket removal.

RESULTS

It was demonstrated that 3D-DIC can be employed to quantify with high resolution time-varying residuum shapes, deformations, and strains. Additionally, the enclosed volumes and cross-sectional areas were computed and analyzed.

CONCLUSION

This novel low-cost framework provides a promising solution for the in vivo evaluation of residuum shapes and strains, as well as has the potential for characterizing the mechanical properties of the underlying soft tissues.

SIGNIFICANCE

These data may be used to inform data-driven computational algorithms for the design of prosthetic sockets, as well as of other wearable technologies mechanically interfacing with the skin.

Chest Wall Kinematics Using Triangular Cosserat Point Elements in Healthy and Neuromuscular Subjects

Optoelectronic plethysmography (OEP) is a noninvasive method for assessing lung volume variations and the contributions of different anatomical compartments of the chest wall (CW) through measurements of the motion of markers attached to the CW surface. The present study proposes a new method for analyzing the local CW kinematics from OEP measurements based on the kinematics of triangular Cosserat point elements (TCPEs). 52 reflective markers were placed on the anterior CW to create a mesh of 78 triangles according to an anatomical model. Each triangle was characterized by a TCPE and its kinematics was described using four time-variant scalar TCPE parameters. The total CW volume ([Formula: see text]) and the contributions of its six compartments were also estimated, using the same markers. The method was evaluated using measurements of ten healthy subjects, nine patients with Pompe disease, and ten patients with Duchenne muscular dystrophy (DMD), during spontaneous breathing (SB) and vital capacity maneuvers (VC) in the supine position. TCPE parameters and compartmental volumes were compared with [Formula: see text] by computing the phase angles [Formula: see text] (for SB) and the correlation r (for VC) between them. Analysis of [Formula: see text] and r of the outward translation parameter [Formula: see text] of each TCPE revealed that for healthy subjects it provided similar results to those obtained by compartmental volumes, whereas for the neuromuscular patients the TCPE method was capable of detecting local asynchronous and paradoxical movements also in cases where they were undistinguished by volumes. Therefore, the TCPE approach provides additional information to OEP that may enhance its clinical evaluation capabilities.