A comparison of calibration methods for stereo fluoroscopic imaging systems

Hindfoot joint kinematics analysis after the resection of talocalcaneal coalition

BACKGROUND

The subtalar joint movement between the talus and calcaneus is restricted in patients with talocalcaneal coalition (TCC). When the motion of the subtalar joint is restricted, shock absorption in the foot decreases, leading to pain during walking. Resection methods to maintain subtalar motion by removing abnormal unions have been proposed. The purpose of this study was to analyze the joint kinematics of patients who underwent TCC resection and to quantitatively evaluate the results of the surgery based on the measured kinematics.

METHODS

Joint kinematics of five patients with TCC were obtained using a biplane fluoroscopic imaging system and an intensity-based two-/three-dimensional registration method. The joint kinematics of the tibiotalar and subtalar joints and the tibiocalcaneal motion during the stance phase of walking were obtained. From the kinematics of the hindfoot joints, the inversion/eversion range of motion (ROM) of the patients before and after resection was statistically analyzed using the Wilcoxon signed-rank test to test whether TCC resection improved the ROM.

RESULTS

During the loading response period, the eversion ROM of the subtalar joint and tibiocalcaneal motion significantly increased postoperatively. In addition, a significant postoperative increase was observed in the subtalar and tibiocalcaneal inversion ROM during the pre-swing period.

SIGNIFICANCE

TCC resection surgery increased the ROM of the subtalar joint, which in turn contributed to the increase in tibiocalcaneal ROM. Increased subtalar and tibiocalcaneal ROM could result in increased shock attenuation and may be a contributing factor to pain relief during walking.

Knee kinematics are not different between asymmetrical and symmetrical tibial baseplates in total knee arthroplasty: A fluoroscopic analysis of step-up and lunge motions

PURPOSE

This clinical fluoroscopy study investigated knee kinematics of two different cemented fixed-bearing, posterior-stabilised (PS) total knee arthroplasty (TKA) designs: an asymmetric tibial component including an asymmetric insert designed to optimise personalised balance and fit and its precursor symmetrical design with symmetric insert.

METHODS

A consecutive series of patients (16 TKAs from each treatment group) participating in a randomised controlled trial comparing TKA migration was included. The exclusion criterion was the use of walking aids. Flat-panel fluoroscopic recordings of step-up and lunge motions were acquired 1-year postoperatively. Medial and lateral contact points (CPs) were determined to calculate CP displacement, femoral axial rotation and pivot position. Using linear mixed-effects modelling techniques, kinematics between TKA designs were compared.

RESULTS

During knee extension between 20° flexion and full extension, the CPs moved anteriorly combined with a small internal femoral rotation (a screw-home mechanism). Whereas CP movement was reversed: femoral rollback, external femoral rotation while flexing the knee between full extension and 20° knee flexion, At larger flexion angles, femoral axial rotation (FAR) occurred around a lateral pivot point both during step-up and lunge. The symmetric design had a 2.3° larger range of FAR compared to the asymmetric design during lunge (p = 0.02). All other kinematics were comparable.

CONCLUSION

Despite the differences in design, this study showed that the asymmetric and symmetric PS TKA designs had mostly comparable knee kinematics during step-up and lunge motions. It is therefore expected that the functionality of the successor TKA design is similar to that of its precursor design.

LEVEL OF EVIDENCE

Level II, prospective comparative study.

Development, construction, and validation of a thinner uniplanar calibration cage for radiostereometry

Radiostereometric analysis (RSA) is an accurate and precise radiographic method that can be used to measure micromotion of implants and study joint kinematics in vivo. A calibration cage with radiopaque markers is used to calibrate the RSA images; however, the thickness (250 mm) of the calibration cage restricts the available area for the patient and equipment during RSA recordings. A thinner calibration cage would increase the recording area, facilitate handling of the cage, and ease integration of the cage with the RSA system. We developed a thinner calibration cage without compromise of accuracy and precision. First, we performed numerical simulations of an RSA system, and showed that the calibration cage thickness could be decreased to 140 mm maintaining accuracy and precision using 40 fiducial and 30 control markers. Second, we constructed a new calibration cage (NRT cage) according to the simulation results. Third, we validated the new calibration cage against two state-of-the-art calibration cages (Umeaa cage and Leiden cage) in a phantom study. All cages performed similar for marker-based analysis, except for y-rotation, where the Umeaa cage (SD = 0.064 mm) was less precise compared to the NRT (SD = 0.038 mm) and Leiden cages (0.042 mm) (p = .01). For model-based analysis the NRT cage had superior precision for translations (SD ≤ 0.054 mm) over the Leiden cage (SD ≤ 0.118 mm) and Umeaa cage (SD ≤ 0.093 mm) (p < .01). The combined study confirmed that the new and thinner calibration cage maintained accuracy and precision at the level of existing thicker calibration cages.

Three-Dimensional Subject-Specific Knee Shape Reconstruction with Asynchronous Fluoroscopy Images Using Statistical Shape Modeling

Background and objectives: Statistical shape modeling (SSM) based on computerized tomography (CT) datasets has enabled reasonably accurate reconstructions of subject-specific 3D bone morphology from one or two synchronous radiographs for clinical applications. Increasing the number of radiographic images may increase the reconstruction accuracy, but errors related to the temporal and spatial asynchronization of clinical alternating bi-plane fluoroscopy may also increase. The current study aimed to develop a new approach for subject-specific 3D knee shape reconstruction from multiple asynchronous fluoroscopy images from 2, 4, and 6 X-ray detector views using a CT-based SSM model; and to determine the optimum number of planar images for best accuracy via computer simulations and in vivo experiments.

Methods: A CT-based SSM model of the knee was established from 60 training models in a healthy young Chinese male population. A new two-phase optimization approach for 3D subject-specific model reconstruction from multiple asynchronous clinical fluoroscopy images using the SSM was developed, and its performance was evaluated via computer simulation and in vivo experiments using one, two and three image pairs from an alternating bi-plane fluoroscope.

Results: The computer simulation showed that subject-specific 3D shape reconstruction using three image pairs had the best accuracy with RMSE of 0.52 ± 0.09 and 0.63 ± 0.085 mm for the femur and tibia, respectively. The corresponding values for the in vivo study were 0.64 ± 0.084 and 0.69 ± 0.069 mm, respectively, which was significantly better than those using one image pair (0.81 ± 0.126 and 0.83 ± 0.108 mm). No significant differences existed between using two and three image pairs.

Conclusion: A new two-phase optimization approach was developed for SSM-based 3D subject-specific knee model reconstructions using more than one asynchronous fluoroscopy image pair from widely available alternating bi-plane fluoroscopy systems in clinical settings. A CT-based SSM model of the knee was also developed for a healthy young Chinese male population. The new approach was found to have high mode reconstruction accuracy, and those for both two and three image pairs were much better than for a single image pair. Thus, two image pairs may be used when considering computational costs and radiation dosage. The new approach will be useful for generating patient-specific knee models for clinical applications using multiple asynchronous images from alternating bi-plane fluoroscopy widely available in clinical settings. The current SSM model will serve as a basis for further inclusion of training models with a wider range of sizes and morphological features for broader applications.

Robust Self-Supervised Learning of Deterministic Errors in Single-Plane (Monoplanar) and Dual-Plane (Biplanar) X-Ray Fluoroscopy

Fluoroscopic imaging that captures X-ray images at video framerates is advantageous for guiding catheter insertions by vascular surgeons and interventional radiologists. Visualizing the dynamical movements non-invasively allows complex surgical procedures to be performed with less trauma to the patient. To improve surgical precision, endovascular procedures can benefit from more accurate fluoroscopy data via calibration. This paper presents a robust self-calibration algorithm suitable for single-plane and dual-plane fluoroscopy. A three-dimensional (3D) target field was imaged by the fluoroscope in a strong geometric network configuration. The unknown 3D positions of targets and the fluoroscope pose were estimated simultaneously by maximizing the likelihood of the Student-t probability distribution function. A smoothed k-nearest-neighbour (kNN) regression is then used to model the deterministic component of the image reprojection error of the robust bundle adjustment. The Maximum Likelihood Estimation step and the kNN regression step are then repeated iteratively until convergence. Four different error modeling schemes were compared while varying the quantity of training images. It was found that using a smoothed kNN regression can automatically model the systematic errors in fluoroscopy with similar accuracy as a human expert using a small training dataset. When all training images were used, the 3D mapping error was reduced from 0.61-0.83 mm to 0.04 mm post-calibration (94.2-95.7% improvement), and the 2D reprojection error was reduced from 1.17-1.31 to 0.20-0.21 pixels (83.2-83.8% improvement). When using biplanar fluoroscopy, the 3D measurement accuracy of the system improved from 0.60 mm to 0.32 mm (47.2% improvement).

Toward an End-to-End Calibration for Mobile C-Arm in Combination with a Depth Sensor for Surgical Augmented Reality Applications

C-arm X-ray imaging is commonly applied in operating rooms for guiding orthopedic surgeries. Augmented Reality (AR) with C-arm X-ray images during surgery is an efficient way to facilitate procedures for surgeons. However, the accurate calibration process for surgical AR based on C-arm is essential and still challenging due to the limitations of C-arm imaging systems, such as instability of C-arm calibration parameters and the narrow field of view. We extend existing methods using a depth camera and propose a new calibration procedure consisting of calibration of the C-arm imaging system, and 3D/2D calibration of an RGB-D camera and C-arm system with a new method to achieve reliable data and promising accuracy and, at the same time, consistent with standard surgical protocols. For the calibration procedure, we apply bundle adjustment equations with a 3D designed Lego multi-modal phantom, in contrast to the previous methods in which planar calibration phantoms were applied. By using our method, the visualization of the X-ray image upon the 3D data was done, and the achieved mean overlay error was 1.03 mm. The evaluations showed that the proposed calibration procedure provided promising accuracy for AR surgeries and it improved the flexibility and robustness of existing C-arm calibration methods for surgical augmented reality (using C-arm and RGB-D sensor). Moreover, the results showed the efficiency of our method to compensate for the effects of the C-arm movement on calibration parameters. It was shown that the obtained overlay error was improved for the non-zero rotation movement of C-arm by using a virtual detector.

Mechanical characterization of the ligaments in subject-specific models of the patellofemoral joint using in vivo laxity tests

BACKGROUND

The purpose of this study was to propose a methodology for mechanical characterization of the ligaments in subject-specific models of the patellofemoral joint (PFJ) of living individuals.

METHOD

PFJ laxity tests were performed on a healthy volunteer using a specially designed loading apparatus under biplane fluoroscopy. A three-dimensional (3D) parametric model of the PFJ was developed in the framework of the rigid body spring model using the geometrical data acquired from the subject's computed tomography and magnetic resonance images. The stiffness and pre-strains of the medial and lateral PFJ ligaments were characterized using a two-step optimization procedure which minimized the deviation between the model predictions and the calibration test results. Sensitivity analyses were performed to investigate the effect of mechanical properties of the non-characterized model components on the characterization procedure and its results.

RESULTS

The overall findings indicate that the proposed methodology is applicable and can improve the model predictions effectively. For the subject under study, ligament characterization reduced the root mean square of the deviations between the patellar shift and tilt predicted by the model and obtained experimentally for the validation laxity test (from 6.2 mm to 0.5 mm, and from 8.4° to 1.5°, respectively) and passive knee flexion test (from 1.4 mm to 0.3 mm, and from 2.3° to 1.3°, respectively). The non-characterized mechanical properties were found to have a minimal effect on the characterization procedure and its results.

CONCLUSION

The proposed methodology can help in developing truly patient-specific models of the PFJ, to be used for personalized preplanning of the clinical interventions.

A model-based tracking method for measuring 3D dynamic joint motion using an alternating biplane x-ray imaging system

PURPOSES

To propose a new model-based tracking method for measuring three-dimensional (3D) dynamic joint kinematics using a clinical alternating biplane x-ray imaging system; and to quantify in vitro its errors in measuring ankle and knee motions at different motion speeds.

METHODS

A new model-based tracking method based on motion component partition and interpolation (MCPI) was developed for measuring 3D dynamic joint kinematics based on a clinical alternating biplane x-ray imaging system. Two detectors of the biplane imaging system placed perpendicular to each other were operated to collect alternating fluoroscopic images of the targeted joint during tasks. The CT data of the joint were also acquired for the reconstruction of volumetric and surface models of each of the associated bones. The CT-based models of the bones were first registered to the alternating images using a model-to-single-plane fluoroscopic image registration method, and the resulting bone poses were then refined using a two-level optimization with motion component partition and model vertex trajectory interpolation. The MCPI method was evaluated in vitro for measurement errors for an ankle and a knee specimen moving at different speeds against a standard reference provided by a highly accurate motion capture system. The positional and rotational errors of the measured bone poses were quantified in terms of the bias, precision, and root-mean-squared errors (RMSE), as well as the mean target registration error (mTRE), a final mTRE less than 2.5 mm indicating a successful registration.

RESULTS

The new method was found to have RMSE of bone pose measurements of less than 0.18 mm for translations and 0.72° for rotations for the ankle, and 0.33 mm and 0.74° for the knee with a high successful registration rate (>97%), and did not appear to be affected by joint motion speeds. Given the same alternating fluoroscopic images, the MCPI method outperformed the typical biplane analysis method assuming zero time offset between the two fluoroscopic views. The differences in performance between the methods were increased with increased joint motion speed. With the accurate bone pose data, the new method enabled talocrural, subtalar, and tibiofemoral kinematics measurements with submillimeter and subdegree accuracy, except for an RMSE of 1.04° for the internal/external rotation of the talocrural joint.

CONCLUSIONS

A new model-based tracking method based on MCPI has been developed for measuring dynamic joint motions using an alternating biplane x-ray imaging system widely available in medical centers. The MCPI method has been demonstrated in vitro to be highly accurate in determining the 3D kinematics of the bones of both the ankle joint complex and the knee. The current results suggest that the MCPI method would be an effective approach for measuring in vivo 3D kinematics of dynamic joint motion in a clinical setting equipped with an alternating biplane x-ray imaging system.

Soft tissue artefacts of skin markers on the lower limb during cycling: Effects of joint angles and pedal resistance

Soft tissue artefacts (STA) are a major error source in skin marker-based measurement of human movement, and are difficult to eliminate non-invasively. The current study quantified in vivo the STA of skin markers on the thigh and shank during cycling, and studied the effects of knee angles and pedal resistance by using integrated 3D fluoroscopy and stereophotogrammetry. Fifteen young healthy adults performed stationary cycling with and without pedal resistance, while the marker data were measured using a motion capture system, and the motions of the femur and tibia/fibula were recorded using a bi-plane fluoroscopy-to-CT registration method. The STAs with respect to crank and knee angles over the pedaling cycle, as well as the within-cycle variations, were obtained and compared between resistance conditions. The thigh markers showed greater STA than the shank ones, the latter varying linearly with adjacent joint angles, the former non-linearly with greater within-cycle variability. Both STA magnitudes and within-cycle variability were significantly affected by pedal resistance (p<0.05). The STAs appeared to be composed of one component providing the stable and consistent STA patterns and another causing their variations. Mid-segment markers experienced smaller STA ranges than those closer to a joint, but tended to have greater variations primarily associated with pedal resistance and muscle contractions. The current data will be helpful for a better choice of marker positions for data collection, and for developing methods to compensate for both stable and variation components of the STA.

Mobile Biplane X-Ray Imaging System for Measuring 3D Dynamic Joint Motion During Overground Gait

Most X-ray fluoroscopy systems are stationary and impose restrictions on the measurement of dynamic joint motion; for example, knee-joint kinematics during gait is usually measured with the subject ambulating on a treadmill. We developed a computer-controlled, mobile, biplane, X-ray fluoroscopy system to track human body movement for high-speed imaging of 3D joint motion during overground gait. A robotic gantry mechanism translates the two X-ray units alongside the subject, tracking and imaging the joint of interest as the subject moves. The main aim of the present study was to determine the accuracy with which the mobile imaging system measures 3D knee-joint kinematics during walking. In vitro experiments were performed to measure the relative positions of the tibia and femur in an intact human cadaver knee and of the tibial and femoral components of a total knee arthroplasty (TKA) implant during simulated overground gait. Accuracy was determined by calculating mean, standard deviation and root-mean-squared errors from differences between kinematic measurements obtained using volumetric models of the bones and TKA components and reference measurements obtained from metal beads embedded in the bones. Measurement accuracy was enhanced by the ability to track and image the joint concurrently. Maximum root-mean-squared errors were 0.33 mm and 0.65° for translations and rotations of the TKA knee and 0.78 mm and 0.77° for translations and rotations of the intact knee, which are comparable to results reported for treadmill walking using stationary biplane systems. System capability for in vivo joint motion measurement was also demonstrated for overground gait.

Performance of local optimization in single-plane fluoroscopic analysis for total knee arthroplasty

Fluoroscopy-derived joint kinematics plays an important role in the evaluation of knee prostheses. Fluoroscopic analysis requires estimation of the 3D prosthesis pose from its 2D silhouette in the fluoroscopic image, by optimizing a dissimilarity measure. Currently, extensive user-interaction is needed, which makes analysis labor-intensive and operator-dependent. The aim of this study was to review five optimization methods for 3D pose estimation and to assess their performance in finding the correct solution. Two derivative-free optimizers (DHSAnn and IIPM) and three gradient-based optimizers (LevMar, DoNLP2 and IpOpt) were evaluated. For the latter three optimizers two different implementations were evaluated: one with a numerically approximated gradient and one with an analytically derived gradient for computational efficiency. On phantom data, all methods were able to find the 3D pose within 1mm and 1° in more than 85% of cases. IpOpt had the highest success-rate: 97%. On clinical data, the success rates were higher than 85% for the in-plane positions, but not for the rotations. IpOpt was the most expensive method and the application of an analytically derived gradients accelerated the gradient-based methods by a factor 3-4 without any differences in success rate. In conclusion, 85% of the frames can be analyzed automatically in clinical data and only 15% of the frames require manual supervision. The optimal success-rate on phantom data (97% with IpOpt) on phantom data indicates that even less supervision may become feasible.

Plantar-flexion of the ankle joint complex in terminal stance is initiated by subtalar plantar-flexion: A bi-planar fluoroscopy study

Gross motion of the ankle joint complex (AJC) is a summation of the ankle and subtalar joints. Although AJC kinematics have been widely used to evaluate the function of the AJC, the coordinated movements of the ankle and subtalar joints are not well understood. The purpose of this study was to accurately quantify the individual kinematics of the ankle and subtalar joints in the intact foot during ground walking by using a bi-planar fluoroscopic system. Bi-planar fluoroscopic images of the foot and ankle during walking and standing were acquired from 10 healthy subjects. The three-dimensional movements of the tibia, talus, and calcaneus were calculated with a three-dimensional/two-dimensional registration method. The skeletal kinematics were quantified from 9% to 86% of the full stance phase because of the limited camera speed of the X-ray system. At the beginning of terminal stance, plantar-flexion of the AJC was initiated in the subtalar joint on average at 75% ranging from 62% to 76% of the stance phase, and plantar-flexion of the ankle joint did not start until 86% of the stance phase. The earlier change to plantar-flexion in the AJC than the ankle joint due to the early plantar-flexion in the subtalar joint was observed in 8 of the 10 subjects. This phenomenon could be explained by the absence of direct muscle insertion on the talus. Preceding subtalar plantar-flexion could contribute to efficient and stable ankle plantar-flexion by locking the midtarsal joint, but this explanation needs further investigation.

A Reconfigurable High-Speed Stereo-Radiography System for Sub-Millimeter Measurement of In Vivo Joint Kinematics

Relative motions within normal and pathological joints of the human body can occur on the sub-millimeter and sub-degree scale. Dynamic radiography can be used to create a rapid sequence of images from which measurements of bone motion can be extracted, but available systems have limited speed and accuracy, limit normal subject movement, and do not easily integrate into existing traditional motion capture laboratories. A high-speed stereo radiography (HSSR) system is described that addresses these limitations. The custom radiography system was placed on a standalone reconfigurable gantry structure designed to allow freedom of subject movement while integrating into an existing motion capture laboratory. Validation of the system and measurement of knee kinematics of subjects during gait confirmed the ability to record joint motion with high accuracy and high-speed.

Rigorous geometric self-calibrating bundle adjustment for a dual fluoroscopic imaging system

High-speed dual fluoroscopy is a noninvasive imaging technology for three-dimensional skeletal kinematics analysis that finds numerous biomechanical applications. Accurate reconstruction of bone translations and rotations from dual-fluoroscopic data requires accurate calibration of the imaging geometry and the many imaging distortions that corrupt the data. Direct linear transformation methods are commonly applied for performing calibration using a two-step process that suffers from a number of potential shortcomings including that each X-ray source and corresponding camera must be calibrated separately. Consequently, the true imaging set-up and the constraints it presents are not incorporated during calibration. A method to overcome such drawbacks is the single-step self-calibrating bundle adjustment method. This procedure, based on the collinearity principle augmented with imaging distortion models and geometric constraints, has been developed and is reported herein. Its efficacy is shown with a carefully controlled experiment comprising 300 image pairs with 48 507 image points. Application of all geometric constraints and a 31 parameter distortion model resulted in up to 91% improvement in terms of precision (model fit) and up to 71% improvement in terms of 3-D point reconstruction accuracy (0.3-0.4 mm). The accuracy of distance reconstruction was improved from 0.3±2.0 mm to 0.2 ±1.1 mm and angle reconstruction accuracy was improved from -0.03±0.55(°) to 0.01±0.06(°). Such positioning accuracy will allow for the accurate quantification of in vivo arthrokinematics crucial for skeletal biomechanics investigations.

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.

Effect of calibration error on bone tracking accuracy with fluoroscopy

Model-based 3D-fluoroscopy can quantify joint kinematics with 1 mm and 1 deg accuracy level. A calibration based on the acquisition of devices of known geometry is usually applied to size the system. This study aimed at quantifying the sensitivity of the fluoroscopic pose estimation accuracy specifically to errors in the calibration process, excluding other sources of error. X-ray focus calibration error was quantified for different calibration setups, and its propagation to the pose estimation was characterized in-silico. Focus reference position influenced the calibration error dispersion, while calibration cage pose affected its bias. In the worst-case scenario, the estimation error of the principal point and of the focus distance was lower than 1 mm and 2 mm, respectively. The consequent estimation of joint angles was scarcely influenced by calibration errors. A linear trend was highlighted for joint translations, with a sensitivity proportional to the distance between the model and the image plane, resulting in a submillimeter error for realistic calibration errors. The biased component of the error is compensated when computing relative joint kinematics between two segments.

Intervertebral anticollision constraints improve out-of-plane translation accuracy of a single-plane fluoroscopy-to-CT registration method for measuring spinal motion

PURPOSE

The study aimed to propose a new single-plane fluoroscopy-to-CT registration method integrated with intervertebral anticollision constraints for measuring three-dimensional (3D) intervertebral kinematics of the spine; and to evaluate the performance of the method without anticollision and with three variations of the anticollision constraints via an in vitro experiment.

METHODS

The proposed fluoroscopy-to-CT registration approach, called the weighted edge-matching with anticollision (WEMAC) method, was based on the integration of geometrical anticollision constraints for adjacent vertebrae and the weighted edge-matching score (WEMS) method that matched the digitally reconstructed radiographs of the CT models of the vertebrae and the measured single-plane fluoroscopy images. Three variations of the anticollision constraints, namely, T-DOF, R-DOF, and A-DOF methods, were proposed. An in vitro experiment using four porcine cervical spines in different postures was performed to evaluate the performance of the WEMS and the WEMAC methods.

RESULTS

The WEMS method gave high precision and small bias in all components for both vertebral pose and intervertebral pose measurements, except for relatively large errors for the out-of-plane translation component. The WEMAC method successfully reduced the out-of-plane translation errors for intervertebral kinematic measurements while keeping the measurement accuracies for the other five degrees of freedom (DOF) more or less unaltered. The means (standard deviations) of the out-of-plane translational errors were less than -0.5 (0.6) and -0.3 (0.8) mm for the T-DOF method and the R-DOF method, respectively.

CONCLUSIONS

The proposed single-plane fluoroscopy-to-CT registration method reduced the out-of-plane translation errors for intervertebral kinematic measurements while keeping the measurement accuracies for the other five DOF more or less unaltered. With the submillimeter and subdegree accuracy, the WEMAC method was considered accurate for measuring 3D intervertebral kinematics during various functional activities for research and clinical applications.

Comparisons of surface vs. volumetric model-based registration methods using single-plane vs. bi-plane fluoroscopy in measuring spinal kinematics

Several 2D-to-3D image registration methods are available for measuring 3D vertebral motion but their performance has not been evaluated under the same experimental protocol. In this study, four major types of fluoroscopy-to-CT registration methods, with different use of surface vs. volumetric models, and single-plane vs. bi-plane fluoroscopy, were evaluated: STS (surface, single-plane), VTS (volumetric, single-plane), STB (surface, bi-plane) and VTB (volumetric, bi-plane). Two similarity measures were used: 'Contour Difference' for STS and STB and 'Weighted Edge-Matching Score' for VTS and VTB. Two cadaveric porcine cervical spines positioned in a box filled with paraffin and embedded with four radiopaque markers were CT scanned to obtain vertebral models and marker coordinates, and imaged at ten static positions using bi-plane fluoroscopy for subsequent registrations using different methods. The registered vertebral poses were compared to the gold standard poses defined by the marker positions determined using CT and Roentgen stereophotogrammetry analysis. The VTB was found to have the highest precision (translation: 0.4mm; rotation: 0.3°), comparable with the VTS in rotations (0.3°), and the STB in translations (0.6mm). The STS had the lowest precision (translation: 4.1mm; rotation: 2.1°).

Coracohumeral Distances and Correlation to Arm Rotation: An In Vivo 3-Dimensional Biplane Fluoroscopy Study

Background:

Reduced coracohumeral distances have been reported to be associated with anterior shoulder disorders such as subscapularis tears, biceps tendon injuries, and leading edge supraspinatus tears.

Purpose:

To determine the variability in coracohumeral distance as a function of arm rotation in healthy male subjects. The null hypothesis was that no differences in coracohumeral distance would exist with respect to arm rotation.

Study Design:

Descriptive laboratory study.

Methods:

A total of 9 male participants who had full range of motion, strength, and no prior surgery or symptoms in their tested shoulders were enrolled in this institutional review board–approved study. Computed tomography scans of the shoulder were obtained for each subject. A dynamic biplane fluoroscopy system recorded internal and external shoulder rotation with the arm held in the neutral position. Three-dimensional reconstructions of each motion were generated, and the coracohumeral distance and coracoid index (lateral extension of the coracoid) were measured.

Results:

The mean coracohumeral distance in neutral rotation was 12.7 ± 2.1 mm. A significantly shorter minimum coracohumeral distance of 10.6 ± 1.8 mm was achieved (P = .001) at a mean glenohumeral joint internal rotation angle of 36.6° ± 19.2°. This corresponded to a reduction in coracohumeral distance of 16.4% (range, 6.6%-29.8%). The mean coracoid index was 14.2 ± 6.8 mm. A moderate correlation (R = −0.75) existed between the coracohumeral distance and coracoid index.

Conclusion:

Coracohumeral distance was reduced during internal rotation. Decreased coracohumeral distance was correlated with larger coracoid indices.

Clinical Relevance:

This study provides a reference value for coracohumeral distance in the healthy male population. Knowledge of how coracohumeral distance varies over the range of arm internal-external rotation may improve the clinical diagnosis and treatment plan for patients with anterior shoulder pathology, specifically subcoracoid impingement. Imaging of the coracohumeral distance during internal rotation with the hand at approximately midline should be considered to assess patients with anterior shoulder pain.

Effect of plane of arm elevation on glenohumeral kinematics: a normative biplane fluoroscopy study

BACKGROUND

Understanding glenohumeral motion in normal and pathologic states requires the precise measurement of shoulder kinematics. The effect of the plane of arm elevation on glenohumeral translations and rotations remains largely unknown. The purpose of this study was to measure the three-dimensional glenohumeral translations and rotations during arm elevation in healthy subjects.

METHODS

Eight male subjects performed scaption and forward flexion, and five subjects (three men and two women) performed abduction, inside a dynamic biplane fluoroscopy system. Bone geometries were extracted from computed tomography images and used to determine the three-dimensional position and orientation of the humerus and scapula in individual frames. Descriptive statistics were determined for glenohumeral joint rotations and translations, and linear regressions were performed to calculate the scapulohumeral rhythm ratio.

RESULTS

The scapulohumeral rhythm ratio was 2.0 ± 0.4:1 for abduction, 1.6 ± 0.5:1 for scaption, and 1.1 ± 0.3:1 for forward flexion, with the ratio for forward flexion being significantly lower than that for abduction (p = 0.002). Humeral head excursion was largest in abduction (5.1 ± 1.1 mm) and smallest in scaption (2.4 ± 0.6 mm) (p < 0.001). The direction of translation, as determined by the linear regression slope, was more inferior during abduction (-2.1 ± 1.8 mm/90°) compared with forward flexion (0.1 ± 10.9 mm/90°) (p = 0.024).

CONCLUSIONS

Scapulohumeral rhythm significantly decreased as the plane of arm elevation moved in an anterior arc from abduction to forward flexion. The amount of physiologic glenohumeral excursion varied significantly with the plane of elevation, was smallest for scaption, and showed inconsistent patterns across subjects with the exception of consistent inferior translation during abduction.

High knee valgus in female subjects does not yield higher knee translations during drop landings: a biplane fluoroscopic study

The goal of this study was to determine the effects of peak knee valgus angle and peak knee abductor moment on the anterior, medial, and lateral tibial translations (ATT, MTT, LTT) in the "at risk" female knee during drop landing. Fifteen female subjects performed drop landings from 40 cm. Three-dimension knee motion was simultaneously recorded using a high speed, biplane fluoroscopy system, and a video-based motion analysis system. Valgus knee angles and knee abduction moments were stratified into low, intermediate, and high groups and peak ATT, MTT, and LTT were compared between these groups with ANOVA (α = 0.05). Significant differences were observed between stratified groups in peak knee valgus angle (p < 0.0001) and peak knee abduction moment (p < 0.0001). However, no corresponding differences in peak ATT, LTT, and MTT between groups exhibiting low to high-peak knee valgus angles (ATT: p = 0.80; LTT: p = 0.25; MTT: p = 0.72); or, in peak ATT (p = 0.61), LTT (p = 0.26) and MTT (p = 0.96) translations when stratified according to low to high knee abduction moments, were found. We conclude that the healthy female knee is tightly regulated with regard to translations even when motion analysis derived knee valgus angles and abduction moments are high.

The effects of arm elevation on the 3-dimensional acromiohumeral distance: a biplane fluoroscopy study with normative data

HYPOTHESIS AND BACKGROUND

Narrowing of the subacromial space has been implicated in several shoulder pathologies. However, the location of the minimum distance points during clinical testing has not been defined. We sought to measure the in vivo minimum distance and location of the minimum distance points on the acromion and proximal humerus during arm elevation.

METHODS

Eight healthy male subjects (mean age, 30 years) underwent a dynamic in vivo biplane fluoroscopy assessment of scaption and forward elevation. For each frame, the 3-dimensional position and orientation of the humerus and scapula were determined, and the acromiohumeral distance (AHD) was measured as the shortest distance between the acromion and proximal humerus.

RESULTS

The minimum AHD was 2.6 ± 0.8 mm during scaption and 1.8 ± 1.2 mm during forward flexion at elevation angles of 83° ± 13° and 97° ± 23°, respectively. The minimum distance point was located on the articular surface of the humeral head from the neutral arm position until 34° ± 8° for scaption and 36° ± 6° for forward flexion. Upon further elevation, the minimum distance point was located within the footprint of the supraspinatus muscle until 72° ± 12° for scaption and 65° ± 8° for forward flexion. At greater elevation angles, the minimum distance points were between the acromion and the proximal humeral shaft, distal from the greater tuberosity.

CONCLUSIONS

The shortest AHD was at approximately 90° of arm elevation. The AHD was no longer measured intra-articularly or within the supraspinatus footprint above approximately 70° of arm elevation.

Accuracy of a contour-based biplane fluoroscopy technique for tracking knee joint kinematics of different speeds

While measuring knee motion in all six degrees of freedom is important for understanding and treating orthopaedic knee pathologies, traditional motion capture techniques lack the required accuracy. A variety of model-based biplane fluoroscopy techniques have been developed with sub-millimeter accuracy. However, no studies have statistically evaluated the consistency of the accuracy across motions of varying intensity or between degrees of freedom. Therefore, this study evaluated the bias and precision of a contour-based tracking technique by comparing it to a marker-based method (gold standard) during three movements with increasing intensity. Six cadaveric knees with implanted tantalum markers were used to simulate knee extension, walking and drop landings, while motion was recorded by a custom biplane fluoroscopy system. The 3D geometries of the bones were reconstructed from CT scans and anatomical coordinate systems were assigned. The position and orientation of the bone and marker models were determined for an average of 27 frames for each trial and knee joint kinematics were compared. The average bias and precision was 0.01 ± 0.65° for rotations and 0.01 ± 0.59 mm for joint translations. Rotational precision was affected by motion (p=0.04) and depended on the axis of rotation (p=0.02). However, the difference in average precision among motions or axes was small (≤ 0.13°) and not likely of consequence for kinematic measurements. No other differences were found. The contour-based technique demonstrated sub-millimeter and sub-degree accuracy, indicating it is a highly accurate tool for measuring complex three dimensional knee movements of any intensity.

Statistical shape model-based femur kinematics from biplane fluoroscopy

Studying joint kinematics is of interest to improve prosthesis design and to characterize postoperative motion. State of the art techniques register bones segmented from prior computed tomography or magnetic resonance scans with X-ray fluoroscopic sequences. Elimination of the prior 3D acquisition could potentially lower costs and radiation dose. Therefore, we propose to substitute the segmented bone surface with a statistical shape model based estimate. A dedicated dynamic reconstruction and tracking algorithm was developed estimating the shape based on all frames, and pose per frame. The algorithm minimizes the difference between the projected bone contour and image edges. To increase robustness, we employ a dynamic prior, image features, and prior knowledge about bone edge appearances. This enables tracking and reconstruction from a single initial pose per sequence. We evaluated our method on the distal femur using eight biplane fluoroscopic drop-landing sequences. The proposed dynamic prior and features increased the convergence rate of the reconstruction from 71% to 91%, using a convergence limit of 3 mm. The achieved root mean square point-to-surface accuracy at the converged frames was 1.48 ± 0.41 mm. The resulting tracking precision was 1-1.5 mm, with the largest errors occurring in the rotation around the femoral shaft (about 2.5° precision).

The long head of the biceps tendon has minimal effect on in vivo glenohumeral kinematics: a biplane fluoroscopy study

BACKGROUND

The in vivo stabilizing role of the long head of the biceps tendon (LHB) is poorly understood. While cadaveric studies report that the loaded LHB constrains translations in all directions, clinical data suggest that there is no clinically demonstrable alteration in glenohumeral position after LHB tenodesis or tenotomy. The purpose of this study was to investigate potential alterations in glenohumeral kinematics after LHB tenodesis during 3 dynamic in vivo motions using a biplane fluoroscopy system.

HYPOTHESIS

Our hypothesis was that there would be no difference in glenohumeral translations greater than 1.0 mm between shoulders after biceps tenodesis and healthy contralateral shoulders.

STUDY DESIGN

Controlled laboratory study.

METHODS

Five patients who underwent unilateral, open subpectoral tenodesis performed abduction, a simulated late cocking phase of a throw, and simulated lifting with both their tenodesed shoulder and their contralateral healthy shoulder inside a biplane fluoroscopy system. Dynamic 3-dimensional glenohumeral positions and electromyography activity of the biceps brachii muscle were determined and compared.

RESULTS

Significant glenohumeral translations occurred in both shoulders for abduction (3.4 mm inferiorly; P < .01) and simulated late cocking (2.6 mm anteriorly; P < .01). The mean difference for each motion in glenohumeral position between the tenodesed and the contralateral healthy shoulders was always less than 1.0 mm. The tenodesed shoulders were more anterior (centered) during abduction (0.7 mm; P < .01) and for the eccentric phase of the simulated late cocking motion (0.9 mm; P < .02). No significant differences were found during the simulated lifting motion and in the superior-inferior direction.

CONCLUSION

The effect of biceps tenodesis on glenohumeral position during the motions studied in vivo was minimal compared with physiological translations and interpatient variability.

CLINICAL RELEVANCE

Our findings demonstrated that LHB tenodesis does not dramatically alter glenohumeral position during dynamic motions, suggesting the risk for clinically significant alterations in glenohumeral kinematics after tenodesis is low in otherwise intact shoulders.