Validation of a non-invasive fluoroscopic imaging technique for the measurement of dynamic knee joint motion

The accurate measurement of the in vivo knee joint kinematics in six degrees-of-freedom (6DOF) remains a challenge in biomedical engineering. We have adapted a dual fluoroscopic imaging system (DFIS) to investigate the various in vivo dynamic knee joint motions. This paper presents a thorough validation of the accuracy and repeatability of the DFIS system when used to measure 6DOF dynamic knee kinematics. First, the validation utilized standard geometric spheres made from different materials to demonstrate the capability of the DFIS technique to determine the object positions under changing speeds. The translational pose of the spheres could be recreated to less than 0.15±0.09 mm for velocities below 300 mm/s. Next, tantalum beads were inserted into the femur and tibia of two fresh frozen cadaver knees to compare the dynamic kinematics measured by matching knee models to the kinematics from the tantalum bead matching—a technique similar to Roentgen stereophotogrammetric analysis (RSA). Each cadaveric knee was attached to the crosshead of a tensile testing machine and vertically translated at a rate of 16.66 mm/s while images were captured with the DFIS. Subsequently, the tibia was held fixed and the femur manually flexed from full extension to 90° of flexion, as the DFIS acquired images. In vitro translation of the cadaver knee using the tensile testing machine deviated from predicted values by 0.08±0.14 mm for the matched knee models. The difference between matching the knee and tantalum bead models during the dynamic flexion–extension motion of the knee was 0.1±0.65°/s in flexion speed; 0.24±0.16 mm in posterior femoral translation; and 0.16±0.61° in internal–external tibial rotation. Finally, we applied the method to investigate the knee kinematics of a living subject during a step ascent and treadmill gait. High repeatability was demonstrated for the in vivo application. Thus, the DFIS provides an easy and powerful tool for accurately determining 6DOF positions of the knee when performing daily functional activities.     

Bone alignment using the iterative closest point algorithm

Computer assisted surgical interventions and research in joint kinematics rely heavily on the accurate registration of three-dimensional bone surface models reconstructed from various imaging technologies. Anomalous results were seen in a kinematic study of carpal bones using a principal axes alignment approach for the registration. The study was repeated using an iterative closest point algorithm, which is more accurate, but also more demanding to apply. The principal axes method showed errors between 0.35 mm and 0.49 mm for the scaphoid, and between 0.40 mm and 1.22 mm for the pisiform. The iterative closest point method produced errors of less than 0.4 mm. These results show that while the principal axes method approached the accuracy of the iterative closest point algorithm in asymmetrical bones, there were more pronounced errors in bones with some symmetry. Principal axes registration for carpal bones should be avoided.     

In vivo kinematics of two-component total ankle arthroplasty during non-weightbearing and weightbearing dorsiflexion/plantarflexion

Relatively high rates of loosening and implant failure have been reported after total ankle arthroplasty, especially in first and second generation implants. Abnormal kinematics and incongruency of the articular surface may cause increased loads applied to the implant with concomitant polyethylene wear, resulting in loosening and implant failure. The purpose of this study was to measure three-dimensional kinematics of two-component total ankle arthroplasty during non-weightbearing and weightbearing activities, and to investigate incongruency of the articular surfaces during these activities. Forty-seven patients with a mean age of 71 years were enrolled. Radiographs were taken at non-weightbearing maximal dorsiflexion and plantarflexion, and weightbearing maximal dorsiflexion, plantarflexion, and neutral position. 3D-2D model-image registration was performed using the radiographs and the three-dimensional implant models, and three-dimensional joint angles were determined. The implanted ankles showed 18.1±8.6° (mean±standard deviation) of plantarflexion, 0.1±0.7° of inversion, 1.2±2.0° of internal rotation, and 0.8±0.6mm of posterior translation of the talar component in the non-weightbearing activity, and 17.8±7.5° of plantarflexion, 0.4±0.5° of inversion, 1.8±2.0° of internal rotation, and 0.7±0.5mm of posterior translation in the weightbearing activity. There were no significant differences between the non-weightbearing and weightbearing kinematics except for the plantarflexion angle. Incongruency of the articular surface occurred in more than 75% of the ankles. Our observations will provide useful data against which kinematics of other implant designs, such as three-component total ankle arthroplasty, can be compared.     

In vivo determination of the anatomical axes of the ankle joint complex: An optimization approach

This study investigates the feasibility of a subject-specific three-dimensional model of the ankle joint complex for kinematic and dynamic analysis of movement. The ankle joint complex was modelled as a three-segment system, connected by two ideal highe joints: the talocrural and the subtalar joint. A mathematical formulation was developed to express the three-dimensional translation and rotation between the foot and shank segments as a function of the two joint angles, and 12 model parameters describing the locations of the joint axes. An optimization method was used to fit the model parameters to three-dimensional kinematic data of foot and shank markers, obtained during test movements throughout the entire physiological range of motion of the ankle joint. The movement of the talus segment, which cannot be measured non-invasively, is not necessary for the analysis.

This optimization method was used to determine the position and orientation of the joint axes in 14 normal subjects. After optimization, the discrepancy between the best fitting model and actual marker kinematics was between 1 and 3 mm for all subjects. The predicted inclination of the subtalar joint axis from the horizontal plane was 37.4±2.7°, and the medial deviation was 18.0±16.2°. The lateral side of the talucrural axis was directed slightly posteriorly (6.8±8.1°), and inclined downward by 7.0±5.4°. These results are similar to previously reported typical results from anatomical, in vitro, studies. Reproducibility was evaluated by repeated testing of one subject, which resulted in variations of about one-fifth of the standard deviation within the group, the inclination of the subtalar joint axis was significantly correlated to the arch height and a radiographic ‘tarsal index’. It is concluded that this optimization method provides the opportunity to incorporate inter-individual anatomical differences into kinematic and dynamic analysis of the ankle joint complex. This allows a more functional interpretation of kinematic data, and more realistic estimates of internal forces.     

Effect of joint rotation correction when measuring elongation of the gastrocnemius medialis tendon and aponeurosis

It is well known that during maximal plantar flexion contractions the ankle joint rotation overestimates the actual elongation of the tendon and aponeurosis. The aim of this study was to examine the influence of the curve length changes of the Achilles tendon on the joint rotation corrected elongation and strain of the gastrocnemius medialis (GM) tendon and aponeurosis. Nine subjects (age: 29.4 ± 5.7 years, body mass: 78.8 ± 6.8 kg, body height: 178 ± 4 cm) participated in the study. The subjects performed maximal voluntary isometric plantarflexion contractions in the prone position on a Biodex-dynamometer. Ultrasonography (Aloka SSD 4000) was used to visualize the muscle belly of the GM muscle-tendon unit. To calculate the curve length changes of the Achilles tendon its surface contour was reconstructed using a series of small reflective skin markers having a diameter of 2.5 mm. The elongation of the GM tendon and aponeurosis was calculated (a) as the difference of the measured and the passive (due to joint rotation) displacement of the tendon and aponeurosis and (b) as the difference of the measured displacement and the length changes of the reconstructed Achilles tendon surface contour. The absolute difference between the elongation obtained by both methods were 1.2 ± 0.4 mm. These differences were due to the higher changes in length obtained by the reconstruction of the tendon curved surface contour as compared to the changes observed in the passive displacement of the digitised point at the aponeurosis. Without correcting for angle joint rotation, the measured elongation clearly overestimates the actual elongation of the GM tendon and aponeurosis. After the passive displacement correction the calculated elongation still overestimates the actual elongation of the GM tendon and aponeurosis. However, this overestimation has a negligible effect on the examined in vivo strain (0.3%) of the tendon and aponeurosis.     

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.     

In-vivo range of motion of the subtalar joint using computed tomography

Understanding in vivo subtalar joint kinematics is important for evaluation of subtalar joint instability, the design of a subtalar prosthesis and for analysing surgical procedures of the ankle and hindfoot. No accurate data are available on the normal range of subtalar joint motion. The purpose of this study was to introduce a method that enables the quantification of the extremes of the range of motion of the subtalar joint in a loaded state using multidetector computed tomography (CT) imaging. In 20 subjects, an external load was applied to a footplate and forced the otherwise unconstrained foot in eight extreme positions. These extreme positions were foot dorsiflexion, plantarflexion, eversion, inversion and four extreme positions in between the before mentioned positions. CT images were acquired in a neutral foot position and each extreme position separately. After bone segmentation and contour matching of the CT data sets, the helical axes were determined for the motion of the calcaneus relative to the talus between four pairs of opposite extreme foot positions. The helical axis was represented in a coordinate system based on the geometric principal axes of the subjects’ talus. The greatest relative motion between the calcaneus and the talus was calculated for foot motion from extreme eversion to extreme inversion (mean rotation about the helical axis of 37.3±5.9°, mean translation of 2.3±1.1 mm). A consistent pattern of range of subtalar joint motion was found for motion of the foot with a considerable eversion and inversion component.     

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.     

Automatic registration of MRI-based joint models to high-speed biplanar radiographs for precise quantification of in vivo anterior cruciate ligament deformation during gait

Understanding in vivo joint mechanics during dynamic activity is crucial for revealing mechanisms of injury and disease development. To this end, laboratories have utilized computed tomography (CT) to create 3-dimensional (3D) models of bone, which are then registered to high-speed biplanar radiographic data captured during movement in order to measure in vivo joint kinematics. In the present study, we describe a system for measuring dynamic joint mechanics using 3D surface models of the joint created from magnetic resonance imaging (MRI) registered to high-speed biplanar radiographs using a novel automatic registration algorithm. The use of MRI allows for modeling of both bony and soft tissue structures. Specifically, the attachment site footprints of the anterior cruciate ligament (ACL) on the femur and tibia can be modeled, allowing for measurement of dynamic ACL deformation. In the present study, we demonstrate the precision of this system by tracking the motion of a cadaveric porcine knee joint. We then utilize this system to quantify in vivo ACL deformation during gait in four healthy volunteers.     

A methodology to accurately quantify patellofemoral cartilage contact kinematics by combining 3D image shape registration and cine-PC MRI velocity data

Patellofemoral osteoarthritis and its potential precursor patellofemoral pain syndrome (PFPS) are common, costly, and debilitating diseases. PFPS has been shown to be associated with altered patellofemoral joint mechanics; however, an actual variation in joint contact stresses has not been established due to challenges in accurately quantifying in vivo contact kinematics (area and location). This study developed and validated a method for tracking dynamic, in vivo cartilage contact kinematics by combining three magnetic resonance imaging (MRI) techniques, cine-phase contrast (CPC), multi-plane cine (MPC), and 3D high-resolution static imaging. CPC and MPC data were acquired from 12 healthy volunteers while they actively extended/flexed their knee within the MRI scanner. Since no gold standard exists for the quantification of in vivo dynamic cartilage contact kinematics, the accuracy of tracking a single point (patellar origin relative to the femur) represented the accuracy of tracking the kinematics of an entire surface. The accuracy was determined by the average absolute error between the PF kinematics derived through registration of MPC images to a static model and those derived through integration of the CPC velocity data. The accuracy ranged from 0.47 mm to 0.77 mm for the patella and femur and from 0.68 mm to 0.86 mm for the patellofemoral joint. For purely quantifying joint kinematics, CPC remains an analytically simpler and more accurate (accuracy <0.33 mm) technique. However, for application requiring the tracking of an entire surface, such as quantifying cartilage contact kinematics, this combined imaging approach produces accurate results with minimal operator intervention.     

Isolation and characterization of heparin and gelatin binding buffalo seminal plasma proteins and their effect on cauda epididymal spermatozoa

Seventy semen ejaculates were obtained from 14 Murrah buffalo bulls and were subjected to plasma separation immediately after collection by centrifugation at 2000 rpm for 20 min and stored in liquid nitrogen until analysis. In the seminal plasma the total protein concentration were estimated and the heparin and gelatin binding (HB and GB) proteins were isolated using heparin and gelatin affinity column chromatography. The molecular weight of individual isolated HB and GB protein was determined by SDS–PAGE analysis. Buffalo bull spermatozoa was collected from cauda epididymis under aseptic conditions and was used for the in vitro fertility tests (i.e. bovine cervical mucus penetration test (BCMPT) and hypo-osmotic swelling test (HOST)). The heparin and gelatin binding buffalo seminal plasma proteins were used in six concentrations i.e. 10, 20, 30, 40, 50 and 60 μg/ml to test their effect on in vitro fertility assessment of cauda epididymal spermatozoa. The overall mean values of total protein, HB and GB proteins were recorded as 29 ± 2.7, 2.61 and 0.2 mg/ml, respectively. Eighteen total protein bands were observed in the range of 12–127 kDa. Eight major HB proteins were isolated in the range of 13–71 kDa. Seven major GB proteins were isolated in the range of 13–61 kDa in the buffalo seminal plasma. The mean penetration distance (mm) travelled by the buffalo cauda spermatozoa was maximum in HB proteins (26.9 ± 0.6) followed by GB proteins (25.4 ± 0.6) and control (21.2 ± 1.4). The difference in BCMPT values between protein treated and control group was significant (P < 0.05). Almost similar trend in the effect of protein on values of HOST percentage in both HB and GB proteins treated semen samples were recorded (66.4 ± 0.65 and 66.1 ± 0.6, respectively). The difference in HOST values between proteins treated and control group (50.4 ± 2.0) was significant (P < 0.05). The present results indicate that among the isolated proteins, 4 proteins were commonly seen in both the heparin and gelatin–sepharose affinity column chromatography, and the addition of buffalo seminal plasma proteins improved the in vitro sperm functions (40 μg/ml gave best results) of buffalo cauda spermatozoa.     

Aspiration of oocytes from transitional, cycling, and pregnant mares

The aim of this study was to compare the efficacy of three approaches for recovering equine oocytes via transvaginal ultrasound-guided follicular aspiration. Fourteen mares were used as oocyte donors during the spring transition period and physiologic breeding season, and 11 mares were bred for use as oocyte donors during early gestation. In all mares, large (>20 mm) and small (10–20 mm) follicles were aspirated in eight rounds every 10–11 days. In each of the four rounds during the transition period, half the mares received 12.5 mg eFSH once daily for 4 days prior to aspiration. For each of the four rounds during the cycling season, half the mares received 12.5 mg eFSH twice daily for 3 days prior to aspiration. Pregnant mares were aspirated on days 25, 40 and 55 of gestation and received no eFSH. There were more large (>20 mm) follicles in cycling controls (2.25 ± 0.27) and cycling FSH-treated (2.64 ± 0.27) mares than in transitional FSH-treated mares (1.18 ± 0.27). The number of oocytes recovered from small (10–20 mm) follicles varied by mare (P < 0.05) and averaged 1.08 ± 0.22 per aspiration for transitional mares and 1.23 ± 0.22 per aspiration for cycling mares (P > 0.1). The number of oocytes per aspiration from large follicles was greater in cycling FSH-treated mares (0.46 ± 0.09) than in transitional control mares (0.11 ± 0.09). In pregnant mares, more large follicles were present at day 25 than at any other time, and the number of oocytes per aspiration from large follicles was greater at day 25 (0.73 ± 0.16) than at day 55 (0.04 ± 0.18). When compared across all seasons and treatments, the day 25 pregnant mares yielded the greatest number of oocytes per aspiration (2.91 ± 0.66 per mare).     

Influence of glenohumeral mismatch on bone strains and implant displacements in implanted glenoïds. An in vitro experimental study on cadaveric scapulae

In shoulder arthroplasty, there is no consensus about the ideal mismatch between a prosthetic humeral head and a glenoïd component. Thus, investigations into mismatch effects from a biomechanical point of view can be useful. The aim of this in vitro study was to help us understand mismatch influence on bone strains, translational forces in the joint and implant/bone displacements in implanted scapulae.

Five fresh cadaveric scapulae were implanted with a cemented keeled polyethylene implant. The lower part of the scapulae was embedded and the loadings were carried out using five metallic spheres simulating mismatches of 0, 2, 4, 5 and 6 mm. Loadings included a constant compressive preload of 392 N and an anterior, posterior, inferior and superior translation of 2.5 mm. We measured the transversal force necessary to produce the imposed translation, the strains at six locations around the peripheral cortex of the glenoïd using strain gages and the relative implant/bone displacements using CCD cameras.

Generally, the increase of mismatch reduced the translational forces, the strains around the glenoïd and, except for the anterior loading, the relative implant/bone displacements. Few and even no significant differences were observed when the mismatch varied from 0 to 2 mm; the number of significant differences increased when the mismatch varied from 0 to 4 mm and from 0 to 5 mm; the results obtained for a 0–6 mm variation in mismatch were comparable to those obtained for a 0–5 mm variation.

This study underlines that the mismatch has a significant effect on bone strains, relative implant/bone displacements and induced translational forces when a compressive preload and imposed translations were applied on implanted scapulae.     

An instrumented,dynamic test for anterior laxity of the ankle joint complex

Evaluation of anterior laxity of the ankle joint complex is a difficult clinical problem. Currently, the prime determinant for anterolateral ligament function is the subjective manual examination of anterior laxity of the ankle joint complex. An instrumented dynamic test was developed for objective measurement of anterior laxity of the ankle joint complex. The principle of the test was to apply a force-impulse to the calcaneus, within the muscle reflex time, and to measure anterior–posterior and mediolateral rotation. The test was performed on a cadaver specimen and on 15 volunteers of which five subjects suffered from chronic one-sided lateral ankle ligament instability.

In the cadaver test, anterior translation values increased from 5 to 11 mm, after cutting the anterior talofibular ligament and subsequently cutting the calcaneofibular ligament. In the 10 normal subjects, the mean anterior translation value was 6.7 mm (±1.9 mm). The relative variation of the test result within a measurement session was 2.5% (±1.6%). Between the sessions the relative laxity variation was 2.6% (±2.6%). In the ten normal subjects the mean right–left difference was not significantly different from zero. In four out of the five patients it was more than 2 mm. As in the cadaver test in all measurements, the mediolateral rotations were small (<2.5°). The volunteers complained about same pain at the heel after multiple test sessions.

In conclusion the dynamic, functional test appears to be capable of objectively measuring a value for anterior laxity of the ankle joint complex reflecting the functional status of the anterolateral ankle ligaments.     

An in vitro comparison of bone deformation measured with surface and staple mounted strain gauges

Chicken tibiae were chosen as a model for human second metatarsals. Local surface bone deformation in a 4-point bending configuration was measured in vitro by both strain gauge instrumented staples and strain gauges bonded to the bone's cortical surface. A series of staple bridge dimensions (0.5, 0.6, 0.8 and 1.0 mm) was compared to test for staple influence on bone characteristics and greatest measurement validity and reliability. Thicker staple inhibition of bone deformation was the greatest but differences to thinner staples were not statistically significant (p>0.05). All staples except 0.5 mm had maximum deviations from linearity less than 1%. The 1.0 mm staple had an R² value of 0.992±0.006 plotted against the 4-point bending input force and 0.994±0.002 plotted against the surface strain gauge signal. The mean intraclass correlation coefficients (ICC) calculated with four input forces (30, 60, 90 and 120 N) and for loading and unloading conditions for the 0.5, 0.6, 0.8 and 1.0 mm staples were 0.75, 0.83, 0.87 and 0.92, respectively. Finally, the differences in slope of the staple strain gauge signal plotted against surface strain gauge signal between input force loading and unloading conditions (0.32), and between input compression and tension conditions (0.79) was least for the 1.0 mm staple which also resulted in the lowest standard deviations. These results suggested the appropriateness of the 1.0 mm staple for in vivo application.     

Individual variability of the projection line of the common carotid artery

The investigation has been performed on 87 corpses of persons of both sex, that died after 60 years of age from the pathology not connected with any diseases in the neck organs. Individual variability in the common carotid artery projection line has been revealed; it conforms, to a certain extent, with the value of the neck index. When the neck is short and thick, the artery position corresponds to the line that runs across the following points: the superior--0.5 cm forward from the mandibular angle, the inferior--0.5 cm medially from the sternoclavicular joint. When the neck index is within the limits 1.71-1.88, it is expedient to draw the classical projection line. When the neck index is within the limits 1.57-1.69, it is possible to determine the projection zone as an elongated rectangle. At the bottom of every side the border of this zone is a straight line drawn between the sternoclavicular joint and the point situating 0.6 cm laterally from the joint, and at the top--the line connecting the top of the mastoid process with the point 1.0 cm behind the mandibular angle.     

A biplanar reconstruction method based on 2D and 3D contours: application to the distal femur

A three-dimensional (3D) reconstruction algorithm based on contours identification from biplanar radiographs is presented. It requires, as technical prerequisites, a method to calibrate the biplanar radiographic environment and a surface generic object (anatomic atlas model) representing the structure to be reconstructed. The reconstruction steps consist of: the definition of anatomical regions, the identification of 2D contours associated to these regions, the calculation of 3D contours and projection onto the radiographs, the associations between points of the X-rays contours and points of the projected 3D contours, the optimization of the initial solution and the optimized object deformation to minimize the distance between X-rays contours and projected 3D contours. The evaluation was performed on 8 distal femurs comparing the 3D models obtained to CT-scan reconstructions. Mean error for each distal femur was 1 mm.     

Radiographic inspection of porosity in pure titanium dumbbell castings

Radiographic inspection of porosity in pure titanium dumbbell castings Background: Titanium frameworks are frequently indicated for implant supported prostheses; however, voids are usually encountered inside cast titanium. Objective: This study aimed to confirm the efficacy of a radiographic technique for inspection of porosity in commercially pure titanium castings with different diameter. Materials and methods: Sixty dumbbell rods (n = 20) with a central 1.5, 2.0 and 3.5 mm diameter were prepared by lost‐wax casting. Cast specimens were finished and polished and submitted to radiographic examination (90 kV, 15 mA, 0.6 s and 10–13 mm of distance) using periapical film. The radiographs were visually analysed for the presence of porosity in the extension of the dumbbell or in the central portion of the rods. Data were submitted to Pearson Chi‐square test (5%). Results: The tested radiographic method proved to be suitable for the evaluation of cast frameworks. Internal porosities were observed in most of the specimens (91.7%) (p = 0.0005); however, only 20% occurred on the central portion of the rods (p = 0.612). Conclusion: Internal porosities can be visualised through radiographs and occur mostly in small diameter structures. The radiographic evaluation of metal structures can improve the quality of frameworks and thereby potentially increase the longevity of the rehabilitation.     

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.     

Can generic knee joint models improve the measurement of osteoarthritic knee kinematics during squatting activity?

Knee joint kinematics derived from multi-body optimisation (MBO) still requires evaluation. The objective of this study was to corroborate model-derived kinematics of osteoarthritic knees obtained using four generic knee joint models used in musculoskeletal modelling – spherical, hinge, degree-of-freedom coupling curves and parallel mechanism – against reference knee kinematics measured by stereo-radiography. Root mean square errors ranged from 0.7° to 23.4° for knee rotations and from 0.6 to 9.0 mm for knee displacements. Model-derived knee kinematics computed from generic knee joint models was inaccurate. Future developments and experiments should improve the reliability of osteoarthritic knee models in MBO and musculoskeletal modelling.