In vivo measurement of the dynamic 3-D kinematics of the ovine stifle joint

The ovine stifle joint is a promising animal model for investigation of joint mechanobiology. A method for in vivo measurement of dynamic 3-D kinematics of the ovine stifle joint is described (accuracy: 0.36 +/- 0.39 mm). Inter-subject variability in kinematics is greater than both intra-subject and inter-session variability. For future studies in which joint kinematics are measured prior to and following controlled orthopaedic interventions, pooling of data should be avoided and each subject should act as its own control.     

Technical issues in using robots to reproduce joint specific gait

Reproduction of the in vivo motions of joints has become possible with improvements in robot technology and in vivo measuring techniques. A motion analysis system has been used to measure the motions of the tibia and femur of the ovine stifle joint during normal gait. These in vivo motions are then reproduced with a parallel robot. To ensure that the motion of the joint is accurately reproduced and that the resulting data are reliable, the testing frame, the data acquisition system, and the effects of limitations of the testing platform need to be considered. Of the latter, the stiffness of the robot and the ability of the control system to process sequential points on the path of motion in a timely fashion for repeatable path accuracy are of particular importance. Use of the system developed will lead to a better understanding of the mechanical environment of joints and ligaments in vivo.     

Design and validation of an unconstrained loading system to measure the envelope of motion in the rabbit knee joint

An unconstrained loading system was developed to measure the passive envelope of joint motion in an animal model commonly used to study ligament healing and joint arthritis. The design of the five-degree-of-freedom system allowed for unconstrained knee joint loading throughout flexion with repeated removal and reapplication of the device to a specimen. Seven New Zealand White rabbit knees were subjected to varus, valgus, internal and external loads, and the resulting envelopes of motion were recorded using an electromagnetic tracking device. Intra-specimen reproducibility was excellent when measured in one specimen, with maximal rotational differences of 0.6 and 0.3 deg between the fourth and fifth testing cycles for the varus (VR) and valgus (VL) envelopes, respectively. Similarly, the maximal internal (INT) and external (EXT) envelope differences were 0.5 and 0.4 deg, respectively, between the fourth and fifth cycles. Good inter-animal envelope reproducibility was also observed with consistent motion pathways for each loading condition. A maximal VR-VL laxity of 17.9 +/- 2.3 deg was recorded at 95 deg flexion for the seven knees tested. The maximal INT-EXT laxity of 75.2 +/- 4.8 deg occurred at 50 deg flexion. Studies on measurement reproducibility of re-applying individual testing components demonstrated a maximal error of 1.2 +/- 0.7 deg. Serial removal and re-application (test-retest) of the complete measuring system to one cadaveric knee demonstrated maximal envelope differences of less than 0.7 deg for VR-VL rotation and 2.1 deg for INT-EXT rotation. Our results demonstrate that the measuring system is reproducible and capable of accurate evaluation of knee joint motion. Baseline in vitro data were generated on normal joint kinematics for future in-vivo studies with this system, evaluating ligament healing and disease progression in arthritis models.     

Validation of a new model-based tracking technique for measuring three-dimensional, in vivo glenohumeral joint kinematics

Shoulder motion is complex and significant research efforts have focused on measuring glenohumeral joint motion. Unfortunately, conventional motion measurement techniques are unable to measure glenohumeral joint kinematics during dynamic shoulder motion to clinically significant levels of accuracy. The purpose of this study was to validate the accuracy of a new model-based tracking technique for measuring three-dimensional, in vivo glenohumeral joint kinematics. We have developed a model-based tracking technique for accurately measuring in vivo joint motion from biplane radiographic images that tracks the position of bones based on their three-dimensional shape and texture. To validate this technique, we implanted tantalum beads into the humerus and scapula of both shoulders from three cadaver specimens and then recorded biplane radiographic images of the shoulder while manually moving each specimen's arm. The position of the humerus and scapula were measured using the model-based tracking system and with a previously validated dynamic radiostereometric analysis (RSA) technique. Accuracy was reported in terms of measurement bias, measurement precision, and overall dynamic accuracy by comparing the model-based tracking results to the dynamic RSA results. The model-based tracking technique produced results that were in excellent agreement with the RSA technique. Measurement bias ranged from -0.126 to 0.199 mm for the scapula and ranged from -0.022 to 0.079 mm for the humerus. Dynamic measurement precision was better than 0.130 mm for the scapula and 0.095 mm for the humerus. Overall dynamic accuracy indicated that rms errors in any one direction were less than 0.385 mm for the scapula and less than 0.374 mm for the humerus. These errors correspond to rotational inaccuracies of approximately 0.25 deg for the scapula and 0.47 deg for the humerus. This new model-based tracking approach represents a non-invasive technique for accurately measuring dynamic glenohumeral joint motion under in vivo conditions. The model-based technique achieves accuracy levels that far surpass all previously reported non-invasive techniques for measuring in vivo glenohumeral joint motion. This technique is supported by a rigorous validation study that provides a realistic simulation of in vivo conditions and we fully expect to achieve these levels of accuracy with in vivo human testing. Future research will use this technique to analyze shoulder motion under a variety of testing conditions and to investigate the effects of conservative and surgical treatment of rotator cuff tears on dynamic joint stability.     

A comparison of passive flexion-extension to normal gait in the ovine stifle joint

Obtaining accurate values of joint tissue loads in human subjects and animals in vivo requires exact 3D-reproduction of joint kinematics and comparisons of in vivo motions between subjects and animals, and also necessitates an accurate reference position. For the knee, passive flexion-extension of isolated joints by hand has been assumed to produce bony motions similar to those of normal gait. We hypothesized that passive flexion-extension kinematics would not accurately reproduce in vivo gait, and, further, that such kinematics would vary significantly between testers. In vivo gait motions of four ovine stifle joints were measured in six degrees of freedom, as were passive flexion-extension motions after sacrifice. Passive flexion-extension motions were performed by three testers on the same stifle joints used in vitro. Results showed statistically significant differences in all degrees of freedom, with the largest differences in the proximal-distal and internal-external directions. Differences induced by muscle loads and kinetic factors in vivo were most evident during stance and hoof-off phases of gait. The in vitro passive paths generated by hand created motions with large variability both between and within individual testers. The user dependence and "area" of motion of passive flexion-extension indicates that passive flexion-extension is contained in a volume of motion, rather than constrained to a unique path. The assumption that the passive path has relevance to precise bone positions during normal in vivo gait is not supported by these results. Thus, using passive flexion-extension as a reference between joints may introduce large motion variability in the observed outcome, and large potential errors in determining joint tissue loads.     

Development and validation of a multi-body model of the canine stifle joint

Multi-body musculoskeletal models that can be used concurrently to predict joint contact pressures and muscle forces would be extremely valuable in studying the mechanics of joint injury. The purpose of this study was to develop an anatomically correct canine stifle joint model and validate it against experimental data. A cadaver pelvic limb from one adult dog was used in this study. The femoral head was subjected to axial motion in a mechanical tester. Kinematic and force data were used to validate the computational model. The maximum RMS error between the predicted and measured kinematics during the complete testing cycle was 11.9 mm translational motion between the tibia and the femur and 4.3° rotation between patella and femur. This model is the first step in the development of a musculoskeletal model of the hind limb with anatomically correct joints to study cartilage loading under dynamic conditions.     

An adaptive system identification model of the biomechanical response of the human trunk during sudden loading

Sudden loading injuries to the low back are a concern. Current models are limited in their ability to quantify the time-varying nature of the sudden loading event. The method of approach used six males who were subjected to sudden loads. Response data (EMG and kinematics) were input into a system identification model to yield time-varying torso stiffness estimates. The results show estimates of system stiffness in good agreement with values in the literature. The average root mean square error of the model's predictions of sagittal motion was equal to 0.1 deg. In conclusion, system identification can be implemented with minimal error and used to gain more insight into the time-dependent trunk response to sudden loads.     

The effects of orientation, temperature, and displacement magnitude changes on the sonometrics system accuracy

The accuracy and reliability of a sonomicrometry system (Sonometrics Corporation, Ontario, Canada) was evaluated for its potential use in measuring 3-D in vivo joint kinematics. Distances between different sets of piezoelectric crystals were measured through a salt solution using ultrasound technology. We evaluated crystal-to-crystal distance under simulated in vivo conditions of changing crystal orientation and displacement magnitude. Crystal-to-crystal distance was also evaluated under changing solution temperature, since the crystals may be used at different temperatures. The 2 mm round and peg crystals were accurate to within 0.5mm for 0 through 180 degrees rotations, but the 2mm round suture loop crystals were only reliable at 0 degrees rotation. The speed of sound through a salt solution (and hence the distance between crystals) versus temperature was fit using a second order polynomial, C=1421.1+3.9808T-3.09x10(-2)T2, with an R2 value of 0.9998. The translational error was less than 0.072 mm for crystal displacements of 0.012, 0.2, 1.0, and 5.0 mm. The system was also accurate under dynamic conditions with translational errors that were less than 0.045 mm under 0.65 Hz motion. These results suggest that the Sonometrics crystals possess attributes (translational accuracy and rotational independence) that could provide the basis for a system capable of measuring joint kinematics.     

Simultaneous measurement of three-dimensional joint kinematics and ligament strains with optical methods

The objective of this study was to assess the precision and accuracy of a nonproprietary, optical three-dimensional (3D) motion analysis system for the simultaneous measurement of soft tissue strains and joint kinematics. The system consisted of two high-resolution digital cameras and software for calculating the 3D coordinates of contrast markers. System precision was assessed by examining the variation in the coordinates of static markers over time. Three-dimensional strain measurement accuracy was assessed by moving contrast markers fixed distances in the field of view and calculating the error in predicted strain. Three-dimensional accuracy for kinematic measurements was assessed by simulating the measurements that are required for recording knee kinematics. The field of view (190 mm) was chosen to allow simultaneous recording of markers for soft tissue strain measurement and knee joint kinematics. Average system precision was between +/-0.004 mm and +/-0.035 mm, depending on marker size and camera angle. Absolute error in strain measurement varied from a minimum of +/-0.025% to a maximum of +/-0.142%, depending on the angle between cameras and the direction of strain with respect to the camera axes. Kinematic accuracy for translations was between +/-0.008 mm and +/-0.034 mm, while rotational accuracy was +/-0.082 deg to +/-0.160 deg. These results demonstrate that simultaneous optical measurement of 3D soft tissue strain and 3D joint kinematics can be performed while achieving excellent accuracy for both sets of measurements.     

Design and demonstration of a new instrumented spatial linkage for use in a dynamic environment: application to measurement of ankle rotations during snowboarding

Joint injuries during sporting activities might be reduced by understanding the extent of the dynamic motion of joints prone to injury during maneuvers performed in the field. Because instrumented spatial linkages (ISLs) have been widely used to measure joint motion, it would be useful to extend the functionality of an ISL to measure joint motion in a dynamic environment. The objectives of the work reported by this paper were to (i) design and construct an ISL that will measure dynamic joint motion in a field environment, (ii) calibrate the ISL and quantify its static measurement error, (iii) quantify dynamic measurement error due to external acceleration, and (iv) measure ankle joint complex rotation during snowboarding maneuvers performed on a snow slope. An "elbow-type" ISL was designed to measure ankle joint complex rotation throughout its range (+/-30 deg for flexion/extension, +/-15 deg for internal/external rotation, and +/-15 deg for inversion/eversion). The ISL was calibrated with a custom six degree-of-freedom calibration device generally useful for calibrating ISLs, and static measurement errors of the ISL also were evaluated. Root-mean-squared errors (RMSEs) were 0.59 deg for orientation (1.7% full scale) and 1.00 mm for position (1.7% full scale). A custom dynamic fixture allowed external accelerations (5 g, 0-50 Hz) to be applied to the ISL in each of three linear directions. Maximum measurement deviations due to external acceleration were 0.05 deg in orientation and 0.10 mm in position, which were negligible in comparison to the static errors. The full functionality of the ISL for measuring joint motion in a field environment was demonstrated by measuring rotations of the ankle joint complex during snowboarding maneuvers performed on a snow slope.     

In situ comparison of A-mode ultrasound tracking system and skin-mounted markers for measuring kinematics of the lower extremity

Skin-mounted marker based motion capture systems are widely used in measuring the movement of human joints. Kinematic measurements associated with skin-mounted markers are subject to soft tissue artifacts (STA), since the markers follow skin movement, thus generating errors when used to represent motions of underlying bone segments. We present a novel ultrasound tracking system that is capable of directly measuring tibial and femoral bone surfaces during dynamic motions, and subsequently measuring six-degree-of-freedom (6-DOF) tibiofemoral kinematics. The aim of this study is to quantitatively compare the accuracy of tibiofemoral kinematics estimated by the ultrasound tracking system and by a conventional skin-mounted marker based motion capture system in a cadaveric experimental scenario. Two typical tibiofemoral joint models (spherical and hinge models) were used to derive relevant kinematic outcomes. Intra-cortical bone pins equipped with optical markers were inserted in the tibial and femoral bones to serve as a reference to provide ground truth kinematics. The ultrasound tracking system resulted in lower kinematic errors than the skin-mounted markers (the ultrasound tracking system: maximum root-mean-square (RMS) error 3.44° for rotations and 4.88 mm for translations, skin-mounted markers with the spherical joint model: 6.32° and 6.26 mm, the hinge model: 6.38° and 6.52 mm). Our proposed ultrasound tracking system has the potential of measuring direct bone kinematics, thereby mitigating the influence and propagation of STA. Consequently, this technique could be considered as an alternative method for measuring 6-DOF tibiofemoral kinematics, which may be adopted in gait analysis and clinical practice.     

An MRI-based method to align the compressive loading axis for human cadaveric knees

There is a need to align the mechanical axis of the tibia with the axis of loading for studies involving tibiofemoral compression to interpret results and to ensure repeatability of loading within and among specimens. Therefore, the objectives of this study were (1) to develop a magnetic resonance imaging (MRI)-based alignment method for use with apparatuses applying tibiofemoral joint compression, (2) to demonstrate the usefulness of the method by aligning cadaveric knees in an apparatus that could apply tibiofemoral joint compression, and (3) to quantify the error associated with the alignment method. A four degree-of-freedom adjustable device was constructed to allow determination and alignment of the mechanical axis of the tibia of cadaveric knee joints with the axis of loading of an apparatus applying tibiofemoral joint compression. MRI was used to determine the locations of bony landmarks in three dimensions defining the mechanical axis of the tibia relative to an initial orientation of the four degree-of-freedom device. Adjustment values of the device were then computed and applied to the device to align the mechanical axis of the tibia with the axis of a compressive loading apparatus. To demonstrate the usefulness of the method, four cadaveric knees were aligned in the compressive loading apparatus. The vectors describing the mechanical axis of the tibia and the loading axis of the apparatus before and after adjustment of the four degree-of-freedom device were computed for each cadaveric knee. After adjustment of the four degree-of-freedom device, the mechanical axis of the tibia was collinear with the loading axis of the apparatus for each cadaveric knee. The errors in the adjustment values introduced by inaccuracies in the MR images were quantified using the Monte Carlo technique. The precisions in the translational and rotational adjustments were 1.20 mm and 0.90 deg respectively. The MR-based alignment method will allow consistent interpretation of results obtained during tibiofemoral compressive studies conducted using the apparatus described in this paper by providing a well-defined loading axis. The alignment method can also be adapted for use with other apparatuses applying tibiofemoral compression.     

Development and validation of a computed tomography-based methodology to measure carpal kinematics

Motion of the wrist bones is complicated and difficult to measure. Noninvasive measurement of carpal kinematics using medical images has become popular This technique is difficult and most investigators employ custom software. The objective of this paper is to describe a validated methodology for measuring carpal kinematics from computed tomography (CT) scans using commercial software. Four cadaveric wrists were CT imaged in neutral, full flexion, and full extension. A registration block was attached to the distal radius and used to align the data sets from each position. From the CT data, triangulated surface models of the radius, lunate, and capitate bones were generated using commercial software. The surface models from each wrist position were read into engineering design software that was used to calculate the centroid (position) and principal mass moments of inertia (orientation) of (1) the capitate and lunate relative to the fixed radius and (2) the capitate relative to the lunate. These data were used to calculate the helical axis kinematics for the motions from neutral to extension and neutral to flexion. The kinematics were plotted in three dimensions using a data visualization software package. The accuracy of the method was quantified in a separate set of experiments in which an isolated capitate bone was subjected to two different known rotation/translation motions for ten trials each. For comparison to in vivo techniques, the error in distal radius surface matching was determined using the block technique as a gold standard. The motion that the lunate and capitate underwent was half that of the overall wrist flexion-extension range of motion. Individually, the capitate relative to the lunate and the lunate relative to the radius generally flexed or extended about 30 deg, while the entire wrist (capitate relative to radius) typically flexed or extended about 60 deg. Helical axis translations were small, ranging from 0.6 mm to 1.8 mm across all motions. The accuracy of the method was found to be within 1.4 mm and 0.5 deg (95% confidence intervals). The mean error in distal radius surface matching was 2.4 mm and 1.2 deg compared to the use of a registration block. Carpal kinematics measured using the described methodology were accurate, reproducible, and similar to findings of previous investigators. The use of commercially available software should broaden the access of researchers interested in measuring carpal kinematics using medical imaging.     

Static and dynamic error of a biplanar videoradiography system using marker-based and markerless tracking techniques

The use of biplanar videoradiography technology has become increasingly popular for evaluating joint function in vivo. Two fundamentally different methods are currently employed to reconstruct 3D bone motions captured using this technology. Marker-based tracking requires at least three radio-opaque markers to be implanted in the bone of interest. Markerless tracking makes use of algorithms designed to match 3D bone shapes to biplanar videoradiography data. In order to reliably quantify in vivo bone motion, the systematic error of these tracking techniques should be evaluated. Herein, we present new markerless tracking software that makes use of modern GPU technology, describe a versatile method for quantifying the systematic error of a biplanar videoradiography motion capture system using independent gold standard instrumentation, and evaluate the systematic error of the W.M. Keck XROMM Facility's biplanar videoradiography system using both marker-based and markerless tracking algorithms under static and dynamic motion conditions. A polycarbonate flag embedded with 12 radio-opaque markers was used to evaluate the systematic error of the marker-based tracking algorithm. Three human cadaveric bones (distal femur, distal radius, and distal ulna) were used to evaluate the systematic error of the markerless tracking algorithm. The systematic error was evaluated by comparing motions to independent gold standard instrumentation. Static motions were compared to high accuracy linear and rotary stages while dynamic motions were compared to a high accuracy angular displacement transducer. Marker-based tracking was shown to effectively track motion to within 0.1?mm and 0.1 deg under static and dynamic conditions. Furthermore, the presented results indicate that markerless tracking can be used to effectively track rapid bone motions to within 0.15 deg for the distal aspects of the femur, radius, and ulna. Both marker-based and markerless tracking techniques were in excellent agreement with the gold standard instrumentation for both static and dynamic testing protocols. Future research will employ these techniques to quantify in vivo joint motion for high-speed upper and lower extremity impacts such as jumping, landing, and hammering.     

ACL/MCL transection affects knee ligament insertion distance of healing and intact ligaments during gait in the Ovine model

The objective of this study was to assess the impact of combined transection of the anterior cruciate and medial collateral ligaments on the intact and healing ligaments in the ovine stifle joint. In vivo 3D stifle joint kinematics were measured in eight sheep during treadmill walking (accuracy: 0.4±0.4 mm, 0.4±0.4°). Kinematics were measured with the joint intact and at 2, 4, 8, 12, 16 and 20 weeks after either surgical ligament transection (n=5) or sham surgery without transection (n=3). After sacrifice at 20 weeks, the 3D subject-specific bone and ligament geometry were digitized, and the 3D distances between insertions (DBI) of ligaments during the dynamic in vivo motion were calculated. Anterior cruciate ligament/medial collateral ligament (ACL/MCL) transection resulted in changes in the DBI of not only the transected ACL, but also the intact lateral collateral ligament (LCL) and posterior cruciate ligament (PCL), while the DBI of the transected MCL was not significantly changed. Increases in the maximal ACL DBI (2 week: +4.2 mm, 20 week: +5.7 mm) caused increases in the range of ACL DBI (2 week: 3.6 mm, 20 week: +3.8 mm) and the ACL apparent strain (2 week: +18.9%, 20 week: +24.0%). Decreases in the minimal PCL DBI (2 week: −3.2 mm, 20 week: −4.3 mm) resulted in increases in the range of PCL DBI (2 week: +2.7 mm, 20 week: +3.2 mm). Decreases in the maximal LCL DBI (2 week: −1.0 mm, 20 week: −2.0 mm) caused decreased LCL apparent strain (2 week: −3.4%, 20 week: −6.9%). Changes in the mechanical environment of these ligaments may play a significant role in the biological changes observed in these ligaments.     

Design of an MRI-compatible robotic stereotactic device for minimally invasive interventions in the breast

The objective of this work was to develop a robotic device to perform biopsy and therapeutic interventions in the breast with real-time magnetic resonance imaging (MRI) guidance. The device was designed to allow for (i) stabilization of the breast by compression, (ii) definition of the interventional probe trajectory by setting the height and pitch of a probe insertion apparatus, and (iii) positioning of an interventional probe by setting the depth of insertion. The apparatus is fitted with five computer-controlled degrees of freedom for delivering an interventional procedure. The entire device is constructed of MR compatible materials, i.e. nonmagnetic and non-conductive, to eliminate artifacts and distortion of the MR images. The apparatus is remotely controlled by means of ultrasonic motors and a graphical user interface, providing real-time MR-guided planning and monitoring of the operation. Joint motion measurements found probe placement in less than 50 s and sub-millimeter repeatability of the probe tip for same-direction point-to-point movements. However, backlash in the rotation joint may incur probe tip positional errors of up to 5 mm at a distance of 40 mm from the rotation axis, which may occur for women with large breasts. The imprecision caused by this backlash becomes negligible as the probe tip nears the rotation axis. Real-time MR-guidance will allow the physician to correct this error Compatibility of the device within the MR environment was successfully tested on a 4 Tesla MR human scanner     

Suture augmentation following ACL injury to restore the function of the ACL, MCL, and medial meniscus in the goat stifle joint

Functional tissue engineering (FTE) approaches have shown promise in healing an injured anterior cruciate ligament (ACL) of the knee. Nevertheless, additional mechanical augmentation is needed to maintain joint stability and appropriate loading of the joint while the ACL heals. The objective of this study was to quantitatively evaluate how mechanical augmentation using sutures restores the joint kinematics as well as the distribution of loading among the ACL, medial collateral ligament, and medial meniscus (MM) in response to externally applied loads. Eight goat stifle joints were tested on a robotic/universal force-moment sensor testing system under two loading conditions: (1) a 67N anterior tibial load (ATL) and (2) a 67N ATL with 100N axial compression. For each joint, four experimental conditions were tested at 30°, 60°, and 90° of flexion: the (1) intact and (2) ACL-deficient joint, as well as following (3) suture repair of the transected ACL, and (4) augmentation using sutures passed from the femur to the tibia. Under the 67N ATL, suture augmentation could restore the anterior tibial translation (ATT) to within 3mm of the intact joint (p>0.05), representing a 54-76% improvement over suture repair (p<0.05). With the additional axial compression, the ATT and in-situ forces of the sutures following suture augmentation remained 2-3 times closer to normal (p<0.05). Also, the in-situ forces in the MM were 58-73% lower (p<0.05). Thus, suture augmentation may be helpful in combination with FTE approaches for ACL healing by providing the needed initial joint stability while lowering the loads on the MM.     

Increased musculoskeletal stiffness during load carriage at increasing walking speeds maintains constant vertical excursion of the body center of mass

The primary objective of this research was to determine changes in body and joint stiffness parameters and kinematics of the knee and body center of mass (COM), that result from wearing a backpack (BP) with a 40% body weight load at increasing speeds of walking. It was hypothesized that there would be speed and load-related increases in stiffness that would prevent significant deviations in the COM trajectory and in lower-extremity joint angles. Three independent biomechanical models employing kinematic data were used to estimate global lower-extremity stiffness, vertical stiffness and knee joint rotational stiffness in the sagittal plane during walking on a treadmill at speeds of 0.6-1.6 ms(-1) in 0.2 ms(-1) increments in BP and no backpack conditions. Kinematic data were collected using an Optotrak, three-dimensional motion analysis system. Knee angles and vertical excursion of the COM during the compression (loading phase) increased as a function of speed but not load. All three estimates of stiffness showed significant increases as a function of both speed and load. Significant interaction effects indicated a convergence of load-related stiffness values at lower speeds. Results suggested that increases in muscle-mediated stiffness are used to maintain a constant vertical excursion of the COM under load across the speeds tested, and thereby limit increases in metabolic cost that would occur if the COM would travel through greater vertical range of motion.     

A method for measurement of joint kinematics in vivo by registration of 3-D geometric models with cine phase contrast magnetic resonance imaging data

A new method is presented for measuring joint kinematics by optimally matching modeled trajectories of geometric surface models of bones with cine phase contrast (cine-PC) magnetic resonance imaging data. The incorporation of the geometric bone models (GBMs) allows computation of kinematics based on coordinate systems placed relative to full 3-D anatomy, as well as quantification of changes in articular contact locations and relative velocities during dynamic motion. These capabilities are additional to those of cine-PC based techniques that have been used previously to measure joint kinematics during activity. Cine-PC magnitude and velocity data are collected on a fixed image plane prescribed through a repetitively moved skeletal joint. The intersection of each GBM with a simulated image plane is calculated as the model moves along a computed trajectory, and cine-PC velocity data are sampled from the regions of the velocity images within the area of this intersection. From the sampled velocity data, the instantaneous linear and angular velocities of a coordinate system fixed to the GBM are estimated, and integration of the linear and angular velocities is used to predict updated trajectories. A moving validation phantom that produces motions and velocity data similar to those observed in an experiment on human knee kinematics was designed. This phantom was used to assess cine-PC rigid body tracking performance by comparing the kinematics of the phantom measured by this method to similar measurements made using a magnetic tracking system. Average differences between the two methods were measured as 2.82 mm rms for anterior/posterior tibial position, and 2.63 deg rms for axial rotation. An intertrial repeatability study of human knee kinematics using the new method produced rms differences in anterior/posterior tibial position and axial rotation of 1.44 mm and 2.35 deg. The performance of the method is concluded to be sufficient for the effective study of kinematic changes caused to knees by soft tissue injuries.