MRI vs CT-based 2D-3D auto-registration accuracy for quantifying shoulder motion using biplane video-radiography

Biplane 2D-3D registration approaches have been used for measuring 3D, in vivo glenohumeral (GH) joint kinematics. Computed tomography (CT) has become the gold standard for reconstructing 3D bone models, as it provides high geometric accuracy and similar tissue contrast to video-radiography. Alternatively, magnetic resonance imaging (MRI) would not expose subjects to radiation and provides the ability to add cartilage and other soft tissues to the models. However, the accuracy of MRI-based 2D-3D registration for quantifying glenohumeral kinematics is unknown. We developed an automatic 2D-3D registration program that works with both CT- and MRI-based image volumes for quantifying joint motions. The purpose of this study was to use the proposed 2D-3D auto-registration algorithm to describe the humerus and scapula tracking accuracy of CT- and MRI-based registration relative to radiostereometric analysis (RSA) during dynamic biplanar video-radiography. The GH kinematic accuracy (RMS error) was 0.6–1.0 mm and 0.6–2.2° for the CT-based registration and 1.4–2.2 mm and 1.2–2.6° for MRI-based registration. Higher kinematic accuracy of CT-based registration was expected as MRI provides lower spatial resolution and bone contrast as compared to CT and suffers from spatial distortions. However, the MRI-based registration is within an acceptable accuracy for many clinical research questions.     

The precision of radiostereometric measurements. Manual vs. digital measurements.

The precision of digital vs. manual radiostereometric measurements in total hip arthroplasty was evaluated using repeated stereoradiographic exposures with an interval of 10-15 min. Ten Lubinus SP2 stems cemented into bone specimens and 12 patients with the same stem design were used to evaluate the precision of stem translations and rotations. The precision of translations and rotations of the cup and femoral head penetration was studied in 12 patients with whole polyethylene cups.The use of a measurement method based on digitised radiographs improved the precision for some of the motion parameters, whereas many of them did not change. A corresponding pattern was observed for both the intra- and interobserver error. Of the wear parameters, the most pronounced improvements were the 3D wear and in the proximal-distal direction, although the anterior-posterior precision was also improved. The mean errors of rigid body and elliptic fitting decreased in all evaluations but one, consistent with a more reproducible identification of the markers centres and the edge of the femoral head.Increased precision of radiostereometric measurements may be used to increase the statistical power of future randomised studies and to study new fields in orthopaedics requiring higher precision than has been available with RSA based on manual measurements.     

A new model-based RSA method validated using CAD models and models from reversed engineering

Roentgen stereophotogrammetric analysis (RSA) was developed to measure micromotion of an orthopaedic implant with respect to its surrounding bone. A disadvantage of conventional RSA is that it requires the implant to be marked with tantalum beads. This disadvantage can potentially be resolved with model-based RSA, whereby a 3D model of the implant is used for matching with the actual images and the assessment of position and rotation of the implant. In this study, a model-based RSA algorithm is presented and validated in phantom experiments. To investigate the influence of the accuracy of the implant models that were used for model-based RSA, we studied both computer aided design (CAD) models as well as models obtained by means of reversed engineering (RE) of the actual implant. The results demonstrate that the RE models provide more accurate results than the CAD models. If these RE models are derived from the very same implant, it is possible to achieve a maximum standard deviation of the error in the migration calculation of 0.06 mm for translations in x- and y-direction and 0.14 mm for the out of plane z-direction, respectively. For rotations about the y-axis, the standard deviation was about 0.1 degrees and for rotations about the x- and z-axis 0.05 degrees. Studies with clinical RSA-radiographs must prove that these results can also be reached in a clinical setting, making model-based RSA a possible alternative for marker-based RSA.     

In-vivo measurement of dynamic joint motion using high speed biplane radiography and CT: application to canine ACL deficiency

Dynamic assessment of three-dimensional (3D) skeletal kinematics is essential for understanding normal joint function as well as the effects of injury or disease. This paper presents a novel technique for measuring in-vivo skeletal kinematics that combines data collected from high-speed biplane radiography and static computed tomography (CT). The goals of the present study were to demonstrate that highly precise measurements can be obtained during dynamic movement studies employing high frame-rate biplane video-radiography, to develop a method for expressing joint kinematics in an anatomically relevant coordinate system and to demonstrate the application of this technique by calculating canine tibio-femoral kinematics during dynamic motion. The method consists of four components: the generation and acquisition of high frame rate biplane radiographs, identification and 3D tracking of implanted bone markers, CT-based coordinate system determination, and kinematic analysis routines for determining joint motion in anatomically based coordinates. Results from dynamic tracking of markers inserted in a phantom object showed the system bias was insignificant (-0.02 mm). The average precision in tracking implanted markers in-vivo was 0.064 mm for the distance between markers and 0.31 degree for the angles between markers. Across-trial standard deviations for tibio-femoral translations were similar for all three motion directions, averaging 0.14 mm (range 0.08 to 0.20 mm). Variability in tibio-femoral rotations was more dependent on rotation axis, with across-trial standard deviations averaging 1.71 degrees for flexion/extension, 0.90 degree for internal/external rotation, and 0.40 degree for varus/valgus rotation. Advantages of this technique over traditional motion analysis methods include the elimination of skin motion artifacts, improved tracking precision and the ability to present results in a consistent anatomical reference frame.     

Comparative accuracy of radiostereometric and optical tracking systems

This study aims to quantify and compare the accuracy of traditional radiostereometric analysis (RSA), fluoroscopic RSA (fRSA), and optical tracking systems. Three phantoms were constructed, each having three stainless steel spheres and three reflective markers. One phantom was mounted to the base of a precision cross-slide table, one to the base of a precision rotation table, and the third was mounted to each moveable tabletop. Two dial-gauges, rigidly mounted to the cross-slide table and rotation table, quantified translations and rotations. Two fluoroscopy units placed orthogonally tracked the steel spheres while a four-camera optical motion capture system tracked the reflective markers in three-dimensional space. RSA was performed with both digital radiography and fluoroscopy. Three axes of translation were tested: parallel to one fluoroscopy image, parallel to the other fluoroscopy image, and at approximately 45° to each image. One axis of rotation was tested. Intraclass correlation coefficients indicated excellent agreement between the actual (dial-gauge) and measured translations for all modalities (ICCs>0.99) and excellent agreement between actual and measured rotations for RSA and fRSA (ICCs>0.99). Standard errors of measurement ranged from 0.032 mm and 0.121° for RSA, to 0.040 mm and 0.229° for fRSA, and to 0.109 mm and 0.613° for optical tracking. Differences between actual and measured translations along the 45° axis were significantly smaller than the two parallel axes. These findings suggest that under ideal conditions, accuracy of fRSA is comparable to traditional RSA, and superior to optical tracking. Accuracy is highest when measured at 45° to the fluoroscopy units.     

Digital automated RSA compared to manually operated RSA

The accuracy of digital Roentgen stereophotogrammetric analysis (RSA) was compared to the accuracy of a manually operated RSA system. For this purpose, we used radiographs of a phantom and radiographs of patients. The radiographs of the patients consisted of double examinations of 12 patients that had a tibial osteotomy and of double examinations of 12 patients that received a total hip prosthesis. First, the radiographs were measured manually with an accurate measurement table. Subsequently, the images were digitized by a film scanner at 150 DPI and 300 DPI resolutions and analyzed with the RSA-CMS software. In the phantom experiment, the manually operated system produced significantly better results than the digital system, although the maximum difference between the median values of the manually operated system and the digital system was as low as 0.013 mm for translations and 0.033° for rotations. In the radiographs of the patients, the manually operated system and the digital system produced equally accurate results: no significant differences in translations and rotations were found. We conclude that digital RSA is an accurate, fast, and user friendly alternative for manually operated RSA. Currently, digital RSA systems are being used in a growing number of clinical RSA-studies.     

Recalage 2D-3D dans le domaine fréquentiel pour l’étude du mouvement en condition quasi-statique : étude pilote sur l’articulation du genou

ObjectivesIn order to make kinematics studies, an estimation of the bones motions is necessary. In this paper, we develop a new method to make 2D-3D registration in order to estimate the bones motion of the knee (femur and tibia). Our registration is made in the frequency domain using different pairs of radiographs and 3D reconstruction of the bone in the initial position.Material and methodWe use the knee bones radiographs of two healthy persons. These acquisitions are made by a new low dose radiographic system called EOS. Acquisitions are made in full extension, in 30° and 60° of knee flexion. A 3D reconstruction of the bones is made in the extension position and a new 2D-3D registration method based on the frequency domain is applied to make the motion estimation for the flexion positions.ResultsOur registration method is made in the frequency domain and does not need to make radiograph simulations as it is commonly made in 2D-3D registration methods. We compare our results to a manual registration of the data. The precision we get is similar to the previous works precision.DiscussionFor acquisitions, we use a low dose radiographic system which limits the X-ray irradiation for the patients.     

Implementation and validation of an implant-based coordinate system for RSA migration calculation

An in vitro radiostereometric analysis (RSA) phantom study of a total knee replacement was carried out to evaluate the effect of implementing two new modifications to the conventional RSA procedure: (i) adding a landmark of the tibial component as an implant marker and (ii) defining an implant-based coordinate system constructed from implant landmarks for the calculation of migration results. The motivation for these two modifications were (i) to improve the representation of the implant by the markers by including the stem tip marker which increases the marker distribution (ii) to recover clinical RSA study cases with insufficient numbers of markers visible in the implant polyethylene and (iii) to eliminate errors in migration calculations due to misalignment of the anatomical axes with the RSA global coordinate system. The translational and rotational phantom studies showed no loss of accuracy with the two new measurement methods. The RSA system employing these methods has a precision of better than 0.05 mm for translations and 0.03° for rotations, and an accuracy of 0.05 mm for translations and 0.15° for rotations. These results indicate that the new methods to improve the interpretability, relevance, and standardization of the results do not compromise precision and accuracy, and are suitable for application to clinical data.     

Accuracy of a CT-based bone contour registration method to measure relative bone motions in the hindfoot

For measuring the in-vivo range of motion of the hindfoot, a CT-based bone contour registration method (CT-BCM) was developed to determine the three-dimensional position and orientation of bones. To validate this technique, we hypothesized that the range of motion in the hindfoot is equally, accurately measured by roentgen stereophotogrammetric analysis (RSA) as by the CT-BCM technique.Tantalum bone markers were placed in the distal tibia, talus and calcaneus of one cadaver specimen. With a fixed lower leg, the cadaveric foot was held in neutral and subsequently loaded in eight extreme positions. Immediately after acquiring a CT-scan with the foot in a position, RSA radiographs were made. Bone contour registration and RSA was performed. Helical axis parameters were calculated for talocrural and subtalar joint motion from neutral to extreme positions and between opposite extreme positions. Differences between CT-BCM and RSA were calculated.Compared with RSA, the CT-BCM data registered an overall root mean square difference (RMSd) of 0.21° for rotation about the helical axis, and 0.20 mm translation along the helical axis for the talocrural and subtalar joint and for all motions combined. The RMSd of the position and direction of the helical axes was 3.3 mm and 2.4°, respectively. The latter errors were larger with smaller helical rotations.The differences are similar to those reported for validated RSA and thus are not clinically relevant. Concluding, CT-BCM is an accurate and accessible alternative for studying joint motion, as it does not have the risk of infection and overlapping bone markers.     

使用Radon变换进行二维MRI图像配准

使用了一种基于Radon变换的技术来进行二维的MRI图像配准。MRI的图像配准一般使用灰度配准,而Radon变换一般用于CT图像的重建,虽然现已经存在使用Radon变换进行图像配准,但是比较繁琐,我们对这一配准算法进行了简化。     

Imaging knee position using MRI, RSA/CT and 3D digitisation

The purpose of this study was to compare 3 methods of imaging knee position. Three fresh cadaver knees were imaged at 6 flexion angles between 0 degrees and 120 degrees by MRI, a combination of RSA and CT and 3D digitisation (in two knees). Virtual models of all 42 positions were created using suitable computer software. Each virtual model was aligned to a newly defined anatomically based Cartesian coordinate system. The angular rotations around the 3 coordinate system axes were calculated directly from the aligned virtual models using rigid body kinematics and found to be equally accurate for the 3 methods. The 3 rotations in each knee could be depicted using anatomy-based diagrams for all 3 methods. We conclude that the 3 methods of data acquisition are equally and adequately accurate in vitro. MRI may be the most useful in vivo.     

Estimation of 3D rotation of femur in 2D hip radiographs

Femoral radiographs are affected by the degree of rotation of the femur with respect to the plane of projection. We aimed to determine the 3D rotation of the proximal femur in 2D radiographs. A 3D Statistical Appearance Model (SAM), which was built from CT images of cadaver proximal femurs (n=33) was randomly sampled to form a training set of 500 bones. Nineteen clinical CT images were collected for testing. All CT images were rotated to ±20° in 2° division around the shaft axis, ±10° around medial-lateral axis, and by simultaneous rotation of both axes (±16° and ±8° around shaft and medial-lateral axes). In each orientation, a 2D projection was recorded for generating a 2D SAM. The outcome parameters of the 2D SAM were used as input for a linear regression model and an artificial neural network to predict the rotation. The artificial neural network estimated the rotation more accurately than the linear regression. For artificial neural networks the mean errors were 4.0° and 2.0° around the shaft and medial-lateral axes, respectively. For an individual radiograph, the confidence interval of estimation was still relatively large. However, this method has high potential to differentiate the amount of rotations in two image sets.     

Results of automatic image registration are dependent on initial manual registration

Measurement of static alignment of articulating joints is of clinical benefit and can be determined using image-based registration. We propose a method that could potentially improve the outcome of image-based registration by using initial manual registration. Magnetic resonance images of two wrist specimens were acquired in the relaxed position and during simulated grasp. Transformations were determined from voxel-based image registration between the two volumes. The volumes were manually aligned to match as closely as possible before auto-registration, from which standard transformations were obtained. Then, translation/rotation perturbations were applied to the manual registration to obtain altered initial positions, from which altered auto-registration transformations were obtained. Models of the radiolunate joint were also constructed from the images to simulate joint contact mechanics. We compared the sensitivity of transformations (translations and rotations) and contact mechanics to altering the initial registration condition from the defined standard. We observed that with increasing perturbation, transformation errors appeared to increase and values for contact force and contact area appeared to decrease. Based on these preliminary findings, it appears that the final registration outcome is sensitive to the initial registration.     

Variable angle stereo imaging for rapid patient position correction in an in-house real-time position monitoring system

PurposeTo develop and validate a variable angle stereo image based position correction methodology in an X-ray based in-house online position monitoring system.Materials and methodsA stereo imaging module that enables 3D position determination and couch correction of the patient based on images acquired at any arbitrary angle and arbitrary angular separation was developed and incorporated to the in-house SeedTracker real-time position monitoring system. The accuracy of the developed system was studied by imaging an anthropomorphic phantom implanted with radiopaque markers set to known offset positions from its reference position in an Elekta linear accelerator (LA) and associated XVI imaging system. The accuracy of the system was further validated using CBCT data set from 10 prostate SBRT patients. The time gains achieved with the stereo image based position correction was compared with the manual matching of seed positions in Digitally Reconstructed Radiographs (DRRs) and kV images in the Mosaiq record and verify system.ResultsBased on phantom and patient CBCT dataset study stereo imaging module implemented in the SeedTracker shown to have an accuracy of 0.1(σ = 0.5) mm in detecting the 3D position offset. The time comparison study showed that stereo image based methodology implemented in SeedTracker was a minimum of 80(4) s faster than the manual method implemented in Mosaiq R&V system with a maximum time saving of 146(6) s.ConclusionThe variable angle stereo image based position correction method was shown to be accurate and faster than the standard manual DRR–kV image based correction approach, leading to more efficient treatment.     

Theoretical accuracy of model-based shape matching for measuring natural knee kinematics with single-plane fluoroscopy

Quantification of knee motion under dynamic, in vivo loaded conditions is necessary to understand how knee kinematics influence joint injury, disease, and rehabilitation. Though recent studies have measured three-dimensional knee kinematics by matching geometric bone models to single-plane fluoroscopic images, factors limiting the accuracy of this approach have not been thoroughly investigated. This study used a three-step computational approach to evaluate theoretical accuracy limitations due to the shape matching process alone. First, cortical bone models of the femur tibia/fibula, and patella were created from CT data. Next, synthetic (i.e., computer generated) fluoroscopic images were created by ray tracing the bone models in known poses. Finally, an automated matching algorithm utilizing edge detection methods was developed to align flat-shaded bone models to the synthetic images. Accuracy of the recovered pose parameters was assessed in terms of measurement bias and precision. Under these ideal conditions where other sources of error were eliminated, tibiofemoral poses were within 2 mm for sagittal plane translations and 1.5 deg for all rotations while patellofemoral poses were within 2 mm and 3 deg. However, statistically significant bias was found in most relative pose parameters. Bias disappeared and precision improved by a factor of two when the synthetic images were regenerated using flat shading (i.e., sharp bone edges) instead of ray tracing (i.e., attenuated bone edges). Analysis of absolute pose parameter errors revealed that the automated matching algorithm systematically pushed the flat-shaded bone models too far into the image plane to match the attenuated edges of the synthetic ray-traced images. These results suggest that biased edge detection is the primary factor limiting the theoretical accuracy of this single-plane shape matching procedure.     

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.     

A simplified radiographic method for measuring bone-component motion in total knees.

Measurement of the relative motion between an implant and the surrounding bone over different time periods is valuable for assessing and comparing the component stability and predicting the potential future outcome. The RSA method, where small beads are implanted in the bone adjacent to the component, can measure implant-bone position to an accuracy of about 0.1 mm. However, the method involves special radiographic views and analytical software, not readily available. For purposes of component assessment on a more routine basis and for multi-centre trials, a method was developed where standard A-P and M-L radiographs were used. Computer software was written which estimated the out-of-plane rotations of the component relative to the plane of the film, and then carried out corrections to enhance the accuracy of calculation of the bead heights relative to the component. The theoretical errors were shown to be less than 0.04 mm for the expected range of out-of-plane rotations. When radiographs of components in simulated bones were taken at a range of rotations, the 95% confidence limits for axial displacement were found to be less than +0.3 mm, and for rotation in the plane to be 0.6 degrees. This indicated that the method was useful for studying knee components where the sinkage could reach a range of 0.5-2 mm in a 2 year period.     

Applications of surface-surface matching algorithms for determination of orthodontic tooth movements

Orthodontic tooth movements are described as the differences between initial and final tooth positions. A computer based method for determination of tooth movements for different treatment methods was developed. A total of 20 casts of the upper jaw of patients treated with tooth positioners or fixed appliances were used as a basis for this study. Tooth movement was analysed on casts before (Ci) and after treatment (Cf). The casts were digitized either with a COMT or 3D laser scanning systems. After digitization, the models were superimposed in the palate by using a surface-surface matching algorithm. Tooth surfaces of the orthodontically moved teeth were segmented and determination of tooth movement was accomplished by matching the moved teeth from Ci to Cf. The resulting transformations delivered three dimensional information on translations and rotations. An accuracy of 0.2 mm in translations and 1 degree in rotations could be demonstrated, showing the different efficiency of treatment schemes.     

Fusion of radiostereometric analysis data into computed tomography space: application to the elbow joint

Improvement of joint prostheses is dependent upon information concerning the biomechanical properties of the joint. Radiostereometric analysis (RSA) and electromagnetic techniques have been applied in previous cadaver and in vivo studies on the elbow joint to provide valuable information concerning joint motion axes. However, such information is limited to mathematically calculated positions of the axes according to an orthogonal coordinate system and is difficult to relate to individual skeletal anatomy. The aim of this study was to evaluate the in vivo application of a new fusion method to provide three-dimensional (3D) visualization of flexion axes according to bony landmarks. In vivo RSA data of the elbow joint's flexion axes was combined with data obtained by 3D computed tomography (CT). Results were obtained from five healthy subjects after one was excluded due to an instable RSA marker. The median error between imported and transformed RSA marker coordinates and those obtained in the CT volume was 0.22 mm. Median maximal rotation error after transformation of the rigid RSA body to the CT volume was 0.003 degrees . Points of interception with a plane calculated in the RSA orthogonal coordinate system were imported into the CT volume, facilitating the 3D visualization of the flexion axes. This study demonstrates a successful fusion of RSA and CT data, without significant loss of RSA accuracy. The method could be used for relating individual motion axes to a 3D representation of relevant joint anatomy, thus providing important information for clinical applications such as the development of joint prostheses.     

Regulation of protein secretion by crinophagy in perfused rat liver.

We examined the secretion of three serum proteins, albumin (RSA), alpha 2 mu-globulin (alpha 2 mu G), and transferrin (Trf), in the isolated perfused liver. Within 4 h of perfusion, only 20 to 35% of previously synthesized proteins were secreted by the liver into the recirculating medium. Low temperature inhibited the secretion of alpha 2 mu G and Trf, but not RSA. The amount of RSA secreted by the liver increased twofold in the presence of leupeptin, a proteinase inhibitor, or primaquine, a weak base capable of neutralizing acidic compartments. Neither drug affected Trf secretion, while the release of alpha 2 mu G was enhanced threefold by primaquine treatment. Only 55 to 70% of the total amount of these serum proteins present in the liver at the onset of perfusion could be accounted for after 4 h of perfusion. Our evidence suggests that these losses are due to protein degradation. The degradation of RSA and alpha 2 mu G was inhibited at 15 degrees C and by both leupeptin and primaquine. Contrary, RSA degradation was not altered when livers were perfused at 20 degrees C. Morphological techniques combined with immunological probes were utilized to identify possible intracellular sites of RSA degradation. RSA and cathepsin L were colocalized to large vacuoles found near the cell periphery. Entry of RSA into these vacuoles occurred at 20 degrees C but not at 15 degrees C. Our results using perfused rat livers suggest that as much as 40% of hepatic serum proteins are degraded via fusion of secretory vesicles with lysosomes (e.g., crinophagy).