Can magnetic resonance imaging-derived bone models be used for accurate motion measurement with single-plane three-dimensional shape registration?

The influence of handheld weight on the scapulohumeral rhythm

Scapulohumeral rhythm (SHR) provides insight to neuromuscular control and fundamental biomechanics of the shoulder. This rhythm often is disrupted in pathologic shoulders. As the first step, we sought to quantify SHR in healthy subjects for diagnostic assessment of shoulder function. Ten healthy shoulders were studied. Three-dimensional models of the humerus and scapula were created from computed tomography scans. Dynamic shoulder motion was recorded by use of single-plane fluoroscopy during arm abduction with 0-kg and 3-kg handheld loads. Shoulder kinematics were quantified by use of model-based 3-dimensional-to-2-dimensional registration techniques. SHR decreased (more scapular motion) with increasing abduction. With a 3-kg load, scapulothoracic motion was significantly reduced through the range of 35 degrees to 45 degrees of glenohumeral motion. Muscular stabilization of the scapula increased with external loading, as shown by decreased SHR during early lifting. Dynamic scapular stabilization provides a critical platform for upper extremity activity.

Sensitivity of knee replacement contact calculations to kinematic measurement errors

The ability to measure in vivo knee kinematics accurately makes it tempting to calculate in vivo contact forces, pressures, and areas directly from kinematic data. However, the sensitivity of contact calculations to kinematic measurement errors has not been adequately investigated. To address this issue, we developed a series of sensitivity analyses derived from a validated in vivo computational simulation of gait. The simulation used an elastic foundation contact model to reproduce in vivo contact force, center of pressure, and fluoroscopic motion data collected from an instrumented knee replacement. Treating each degree of freedom (DOF) in the simulation as motion controlled, we first quantified how errors in measured relative pose of the implant components affected contact calculations. Pose variations of +/-0.1 mm or degree over the entire gait cycle changed maximum contact force, pressure, and area by 204, 100, and 117%, respectively. Larger variations of +/-0.5 mm or degree changed these same quantities by 1157, 108, and 578%, respectively. In both cases, the largest sensitivities were to errors in superior-inferior translation and varus-valgus rotation, with loss of contact occurring on one or both sides. We then quantified how switching the sensitive DOFs from motion to load control affected the sensitivity results. Pose variations of +/-0.5 mm or degree in the remaining DOFs changed maximum contact quantities by at most 3%. These results suggest that accuracy on the order of microns and milliradians is needed to estimate contact forces, pressures, and areas directly from in vivo kinematic measurements, and that use of load rather than motion control for the sensitive DOFs may improve the accuracy of in vivo contact calculations.

Dynamic activity dependence of in vivo normal knee kinematics

Dynamic knee kinematics were analyzed for normal knees in three activities, including two different types of maximum knee flexion. Continuous X-ray images of kneel, squat, and stair climb motions were taken using a large flat panel detector. CT-derived bone models were used for model registration-based 3D kinematic measurement. Three-dimensional joint kinematics and contact locations were determined using three methods: bone-fixed coordinate systems, interrogation of CT-based bone model surfaces, and interrogation of MR-based articular cartilage model surfaces. The femur exhibited gradual external rotation throughout the flexion range. Tibiofemoral contact exhibited external rotation, with contact locations translating posterior while maintaining 15 degrees to 20 degrees external rotation from 20 degrees to 80 degrees of flexion. From 80 degrees to maximum flexion, contact locations showed a medial pivot pattern. Kinematics based on bone-fixed coordinate systems differed from kinematics based on interrogation of CT and MR surfaces. Knee kinematics varied significantly by activity, especially in deep flexion. No posterior subluxation occurred for either femoral condyle in maximum knee flexion. Normal knees accommodate a range of motions during various activities while maintaining geometric joint congruency.

Image-based RSA: Roentgen stereophotogrammetric analysis based on 2D–3D image registration

Image-based Roentgen stereophotogrammetric analysis (IBRSA) integrates 2D-3D image registration and conventional RSA. Instead of radiopaque RSA bone markers, IBRSA uses 3D CT data, from which digitally reconstructed radiographs (DRRs) are generated. Using 2D-3D image registration, the 3D pose of the CT is iteratively adjusted such that the generated DRRs resemble the 2D RSA images as closely as possible, according to an image matching metric. Effectively, by registering all 2D follow-up moments to the same 3D CT, the CT volume functions as common ground. In two experiments, using RSA and using a micromanipulator as gold standard, IBRSA has been validated on cadaveric and sawbone scapula radiographs, and good matching results have been achieved. The accuracy was: |mu |< 0.083 mm for translations and |mu| < 0.023 degrees for rotations. The precision sigma in x-, y-, and z-direction was 0.090, 0.077, and 0.220 mm for translations and 0.155 degrees , 0.243 degrees , and 0.074 degrees for rotations. Our results show that the accuracy and precision of in vitro IBRSA, performed under ideal laboratory conditions, are lower than in vitro standard RSA but higher than in vivo standard RSA. Because IBRSA does not require radiopaque markers, it adds functionality to the RSA method by opening new directions and possibilities for research, such as dynamic analyses using fluoroscopy on subjects without markers and computer navigation applications.

Determination of in vivo glenohumeral translation using fluoroscopy and shape-matching techniques

The purpose of this study was to investigate glenohumeral translation in-vivo during active shoulder abduction in the scapular plane. Three-dimensional (3D) models of 9 shoulders were created from CT scans. Fluoroscopic views aligned to the plane of the scapula were recorded during active arm abduction with neutral rotation. 3D motions were determined using model-based 3D-to-two-dimensional (2D) registration. Humeral translation was referenced to the glenoid center in the superior/inferior direction. The humerus moved an average of 1.7 mm superior with arm abduction, from an inferior location to the glenoid center. The humeral head was centered within 1 mm from the glenoid center above 80 degrees abduction. Variability in glenohumeral translation between shoulders decreased significantly from initial to final arm abduction. Our findings agree with some authors' observations of inferior-to-central translation of the humerus and behavior as a congruent ball and socket. We believe this information will help improve the understanding of shoulder function.