Comparative accuracy of radiostereometric and optical tracking systems

Development and Assessment of a Microcomputed Tomography Compatible Five Degrees-of-Freedom Knee Joint Motion Simulator

Currently available knee joint kinematic tracking systems fail to nondestructively capture the subtle variation in joint and soft tissue kinematics that occur in native, injured, and reconstructed joint states. Microcomputed tomography (CT) imaging has the potential as a noninvasive, high-resolution kinematic tracking system, but no dynamic simulators exist to take advantage of this. The purpose of this work was to develop and assess a novel micro-CT compatible knee joint simulator to quantify the knee joint's kinematic and kinetic response to clinically (e.g., pivot shift test) and functionally (e.g., gait) relevant loading. The simulator applies closed-loop, load control over four degrees-of-freedom (DOF) (internal/external rotation, varus/valgus rotation, anterior/posterior translation, and compression/distraction), and static control over a fifth degree-of-freedom (flexion/extension). Simulator accuracy (e.g., load error) and repeatability (e.g., coefficient of variation) were assessed with a cylindrical rubber tubing structure and a human cadaveric knee joint by applying clinically and functionally relevant loads along all active axes. Micro-CT images acquired of the joint at a loaded state were then used to calculate joint kinematics. The simulator loaded both the rubber tubing and the cadaveric specimen to within 0.1% of the load target, with an intertrial coefficient of variation below 0.1% for all clinically relevant loading protocols. The resultant kinematics calculated from the acquired images agreed with previously published values, and produced errors of 1.66 mm, 0.90 mm, 4.41 deg, and 1.60 deg with respect to anterior translation, compression, internal rotation, and valgus rotation, respectively. All images were free of artifacts and showed knee joint displacements in response to clinically and functionally loading with isotropic CT image voxel spacing of 0.15 mm. The results of this study demonstrate that the joint-motion simulator is capable of applying accurate, clinically and functionally relevant loads to cadaveric knee joints, concurrent with micro-CT imaging. Nondestructive tracking of bony landmarks allows for the precise calculation of joint kinematics with less error than traditional optical tracking systems.

Validation of functional calibration and strap-down joint drift correction for computing 3D joint angles of knee, hip, and trunk in alpine skiing

To obtain valid 3D joint angles with inertial sensors careful sensor-to-segment calibration (i.e. functional or anatomical calibration) is required and measured angular velocity at each sensor needs to be integrated to obtain segment and joint orientation (i.e. joint angles). Existing functional and anatomical calibration procedures were optimized for gait analysis and calibration movements were impractical to perform in outdoor settings. Thus, the aims of this study were 1) to propose and validate a set of calibration movements that were optimized for alpine skiing and could be performed outdoors and 2) to validate the 3D joint angles of the knee, hip, and trunk during alpine skiing. The proposed functional calibration movements consisted of squats, trunk rotations, hip ad/abductions, and upright standing. The joint drift correction previously proposed for alpine ski racing was improved by adding a second step to reduce separately azimuth drift. The system was validated indoors on a skiing carpet at the maximum belt speed of 21 km/h and for measurement durations of 120 seconds. Calibration repeatability was on average <2.7° (i.e. 3D joint angles changed on average <2.7° for two repeated sets of calibration movements) and all movements could be executed wearing ski-boots. Joint angle precision was <4.9° for all angles and accuracy ranged from -10.7° to 4.2° where the presence of an athlete-specific bias was observed especially for the flexion angle. The improved joint drift correction reduced azimuth drift from over 25° to less than 5°. In conclusion, the system was valid for measuring 3D joint angles during alpine skiing and could be used outdoors. Errors were similar to the values reported in other studies for gait. The system may be well suited for within-athlete analysis but care should be taken for between-athlete analysis because of a possible athlete-specific joint angle bias.

Probabilistic sensitivity analysis of in-vehicle reach tasks for digital human models considering anthropometric measurement uncertainty

For design using digital human models, human anthropometry data are required as input and are extracted from measurements. There is inherent error associated with these measurements which impacts the output of the simulation. Current techniques in digital human modeling applications primarily employ deterministic methods which are not well suited for handling variability in anthropometric measurement. An alternative to deterministic methods is probabilistic/sensitivity analysis. This study presents a probabilistic sensitivity approach to gain insights into how uncertainty in anthropometric measurements can affect the results of a digital human model with the specific application of vehicle-related reach tasks. Sensitivity levels are found to determine the importance of variability in each joint angle and link length to the final reach. A55-degree of freedom (DOF) digital human model is introduced to demonstrate the sensitivity approach for reach tasks. Seven right-hand reach target points and two left-hand reach target points (creating a total of 14 reach tasks) within a vehicle are used to compare the sensitivities in the joint angles and link lengths resulting from measurement uncertainty. The results show that the importance of each joint angle or link length is dependent on the characteristics of the reach task and sensitivities for joint angles, and link lengths are different for each reach task.

Cervical spine intervertebral kinematics with respect to the head are different during flexion and extension motions

Previous dynamic imaging studies of the cervical spine have focused entirely on intervertebral kinematics while neglecting to investigate the relationship between head motion and intervertebral motion. Specifically, it is unknown if the relationship between head and intervertebral kinematics is affected by movement direction. We tested the hypothesis that there would be no difference in sagittal plane intervertebral angles at identical head orientations during the flexion and extension movements. Nineteen asymptomatic subjects performed continuous head flexion-extension movements while biplane radiographs were collected at 30 images per second. A previously validated model-based volumetric tracking process determined three-dimensional vertebral position with sub-millimeter accuracy throughout the flexion-extension motion. Head movement was recorded at 60 Hz using conventional motion analysis and reflective markers. Intervertebral angles were determined at identical head orientations during the flexion and extension movements. Cervical motion segments were in a more extended orientation during flexion and in a more flexed orientation during extension for any given head orientation. The results suggest that static radiographs cannot accurately represent vertebral orientation during dynamic motion. Further, data should be collected during both flexion and extension movements when investigating intervertebral kinematics with respect to global head orientation. Also, in vitro protocols that use intervertebral total range of motion as validation criteria may be improved by assessing model fidelity using continuous intervertebral kinematics in flexion and in extension. Finally, musculoskeletal models of the head and cervical spine should account for the direction of head motion when determining muscle moment arms because vertebral orientations (and therefore muscle attachment sites) are dependent on the direction of head motion.

Effect of human link length determination on posture reconstruction

Motion capture experiment results are often used as a means of validation for digital human simulations. Motion capture results are marker positions and joint centers in Cartesian space. However, joint angles are more intuitive and easy to understand compared to marker or joint center positions. Posture reconstruction algorithms are used to map Cartesian space to joint space by re-creating experimental postures with simulation models. This allows for direct comparison between the experimental results and digital human simulations. Besides the inherent experimental errors from motion capture system, one source of simulation error is the determination of the link lengths to be used in the simulation model. The link length errors can propagate through all simulation results. Therefore, it is critical to eliminate the link length errors. The objective of this study is to determine the best method of determining link lengths for the simulation model to best match the model to the experiment results containing errors. Specifically, the way that the link lengths are calculated in the posture reconstruction process from motion capture data has a significant effect on the recreated posture for the simulation model. Three link length calculation methods (experimental-average method, trial-specific method, and T-pose method) are developed and compared to a benchmark method (frame-specific method) for calculating link lengths. The results indicate that using the trial-specific method is the most accurate method when referring to calculating frame-specific link lengths.

The accuracy and repeatability of an automatic 2D-3D fluoroscopic image-model registration technique for determining shoulder joint kinematics

Fluoroscopic imaging, using single plane or dual plane images, has grown in popularity to measure dynamic in vivo human shoulder joint kinematics. However, no study has quantified the difference in spatial positional accuracy between single and dual plane image-model registration applied to the shoulder joint. In this paper, an automatic 2D-3D image-model registration technique was validated for accuracy and repeatability with single and dual plane fluoroscopic images. Accuracy was assessed in a cadaver model, kinematics found using the automatic registration technique were compared to those found using radiostereometric analysis. The in vivo repeatability of the automatic registration technique was assessed during the dynamic abduction motion of four human subjects. The in vitro data indicated that the error in spatial positional accuracy of the humerus and the scapula was less than 0.30mm in translation and less than 0.58° in rotation using dual plane images. Single plane accuracy was satisfactory for in-plane motion variables, but out-of-plane motion variables on average were approximately 8 times less accurate. The in vivo test indicated that the repeatability of the automatic 2D-3D image-model registration was 0.50mm in translation and 1.04° in rotation using dual images. For a single plane technique, the repeatability was 3.31mm in translation and 2.46° in rotation for measuring shoulder joint kinematics. The data demonstrate that accurate and repeatable shoulder joint kinematics can be obtained using dual plane fluoroscopic images with an automatic 2D-3D image-model registration technique; and that out-of-plane motion variables are less accurate than in-plane motion variables using a single plane technique.

Image intensifier distortion correction for fluoroscopic RSA: the need for independent accuracy assessment

Fluoroscopic images suffer from multiple modes of image distortion. Therefore, the purpose of this study was to compare the effects of correction using a range of two-dimensional polynomials and a global approach. The primary measure of interest was the average error in the distances between four beads of an accuracy phantom, as measured using RSA. Secondary measures of interest were the root mean squared errors of the fit of the chosen polynomial to the grid of beads used for correction, and the errors in the corrected distances between the points of the grid in a second position. Based upon the two-dimensional measures, a polynomial of order three in the axis of correction and two in the perpendicular axis was preferred. However, based upon the RSA reconstruction, a polynomial of order three in the axis of correction and one in the perpendicular axis was preferred. The use of a calibration frame for these three-dimensional applications most likely tempers the effects of distortion. This study suggests that distortion correction should be validated for each of its applications with an independent "gold standard" phantom.

The effect of decreasing computed tomography dosage on radiostereometric analysis (RSA) accuracy at the glenohumeral joint

Standard, beaded radiostereometric analysis (RSA) and markerless RSA often use computed tomography (CT) scans to create three-dimensional (3D) bone models. However, ethical concerns exist due to risks associated with CT radiation exposure. Therefore, the aim of this study was to investigate the effect of decreasing CT dosage on RSA accuracy. Four cadaveric shoulder specimens were scanned using a normal-dose CT protocol and two low-dose protocols, where the dosage was decreased by 89% and 98%. 3D computer models of the humerus and scapula were created using each CT protocol. Bi-planar fluoroscopy was used to image five different static glenohumeral positions and two dynamic glenohumeral movements, of which a total of five static and four dynamic poses were selected for analysis. For standard RSA, negligible differences were found in bead (0.21±0.31mm) and bony landmark (2.31±1.90mm) locations when the CT dosage was decreased by 98% (p-values>0.167). For markerless RSA kinematic results, excellent agreement was found between the normal-dose and lowest-dose protocol, with all Spearman rank correlation coefficients greater than 0.95. Average root mean squared errors of 1.04±0.68mm and 2.42±0.81° were also found at this reduced dosage for static positions. In summary, CT dosage can be markedly reduced when performing shoulder RSA to minimize the risks of radiation exposure. Standard RSA accuracy was negligibly affected by the 98% CT dose reduction and for markerless RSA, the benefits of decreasing CT dosage to the subject outweigh the introduced errors.

An alternative definition of the scapular coordinate system for use with RSA

When performing radiostereometric analysis (RSA), computed tomography scans are often taken to obtain the landmarks used to create anatomical coordinate systems (CSs) for quantifying joint kinematics. Different conventions for defining CSs lead to an inability to compare results among studies. The International Society of Biomechanics (ISB) has proposed a set of CSs; however, the landmarks needed to create the recommended scapular CS require the entire scapula to be scanned, thereby also exposing breast and other tissues to radiation. The main purpose of this work was to investigate an alternate definition of the CS that has repeatably attainable landmarks and axes as close as possible to those recommended by the ISB, while limiting the portion of the scapula requiring scanning. Intra- and inter-investigator variabilities of landmark digitization were quantified in one model of a scapula and one cadaveric specimen. Based on the variability of the digitizations, an alternative CS was defined. The differences between the ISB and alternative CSs were evaluated on 11 cadaveric specimens. Beaded biplanar RSA was performed on the glenohumeral joint model in 15 different configurations and the resulting kinematics were calculated for each set of landmark digitizations using both sets of coordinate systems. While the kinematic angles obtained using the alternative CS were statistically different from those obtained using the ISB standard, these differences were small (on the order of 5 degrees) and therefore considered to be of little clinical significance. In all likelihood, the benefits of decreasing radiation exposure outweigh these differences in angles.