Variability,agreement and reliability of MRI knee landmarks

Surface mesh reconstructions of bones are often required to define landmark-based coordinate systems, regions of interest and morphological features when studying the soft tissues of the knee from MRI scans. This study reports the variability, agreement and reliability of osseous landmarks to better understand their downstream effects. Fifteen landmarks were defined on the distal femur and twelve on the proximal tibia. Surface meshes were created from twenty right knee MRI scans with a mean subject age of 30.9 years. A single observer identified landmarks on all twenty knees, while three observers repeated the observations three times on a subset of eight knees. All observations were aligned to the Procrustes mean shapes. Principal component analysis was used to study inter-subject variability and two-way ANOVA for inter- and intra-observer agreement and reliability. Inter-subject landmark variation ranged from 0.6 to 5.26 mm, while inter- and intra-observer agreement were at most 5.1 and 5.69 mm respectively. Between-observer reliability ranged from 0.07 to 0.98 while within-observer values were between 0.51 and 0.98. Landmarks derived from fitted spheres or circles often performed well, while most others had their poorest agreement or greatest variation limited to only one or two cardinal directions.     

Coupled motions under compressive load in intact and ACL-deficient knees: a cadaveric study

Knowledge of the coupled motions, which develop under compressive loading of the knee, is useful to determine which degrees of freedom should be included in the study of tibiofemoral contact and also to understand the role of the anterior cruciate ligament (ACL) in coupled motions. The objectives of this study were to measure the coupled motions of the intact knee and ACL-deficient knee under compression and to compare the coupled motions of the ACL-deficient knee with those of the intact knee. Ten intact cadaveric knees were tested by applying a 1600 N compressive load and measuring coupled internal-external and varus-valgus rotations and anterior-posterior and medial-lateral translations at 0 deg, 15 deg, and 30 deg of flexion. Compressive loads were applied along the functional axis of axial rotation, which coincides approximately with the mechanical axis of the tibia. The ACL was excised and the knees were tested again. In the intact knee, the peak coupled motions were 3.8 deg internal rotation at 0 deg flexion changing to -4.9 deg external rotation at 30 deg of flexion, 1.4 deg of varus rotation at 0 deg flexion changing to -1.9 deg valgus rotation at 30 deg of flexion, 1.4 mm of medial translation at 0 deg flexion increasing to 2.3 mm at 30 deg of flexion, and 5.3 mm of anterior translation at 0 deg flexion increasing to 10.2 mm at 30 deg of flexion. All changes in the peak coupled motions from 0 deg to 30 deg flexion were statistically significant (p<0.05). In ACL-deficient knees, there was a strong trend (marginally not significant, p=0.07) toward greater anterior translation (12.7 mm) than that in intact knees (8.0 mm), whereas coupled motions in the other degrees of freedom were comparable. Because the coupled motions in all four degrees of freedom in the intact knee and ACL-deficient knee are sufficiently large to substantially affect the tibiofemoral contact area, all degrees of freedom should be included when either developing mathematical models or designing mechanical testing equipment for study of tibiofemoral contact. The increase in coupled anterior translation in ACL-deficient knees indicates the important role played by the ACL in constraining anterior translation during compressive loading.     

Comparison of kinematics in the healthy and ACL injured knee using MRI

Magnetic Resonance Imaging (MRI) was used to examine the characteristics of abnormal motion in the injured knee by mapping tibiofemoral contact. Eleven healthy subjects and 20 subjects with a unilateral ACL injury performed a leg-press against resistance. MRI scans of both knees at 15 degrees intervals from 0 degrees to 90 degrees of flexion were used to record the tibiofemoral contact pattern. The tibiofemoral contact pattern of the injured knees was more posterior on the tibial plateau than the healthy knees, particularly in the lateral compartment. The tibiofemoral contact pattern of the loaded knees did not differ from the unloaded knees. The difference in the tibiofemoral contact pattern in the ACL injured knee was associated with more severe knee symptoms, irrespective of the passive anterior laxity of the knee.     

Rotational flexibility of the human knee due to varus/valgus and axial moments in vivo

Knee ligamentous injuries persist in the sport of Alpine skiing. To better understand the load mechanisms which lead to injury, pure varus/valgus and pure axial moments were applied both singly and in combination to the right knees of six human test subjects. The corresponding relative knee rotations in three degrees of freedom were measured. Knee flexion angles for each test subject were 15 and 60 degrees for the individual moments and 60 degrees for the combination moments. For both knee flexion angles the hip flexion angle was 0 degrees. Leg muscles were quiescent and axial force was minimal during all tests. Tables of data include sample statistics for each of four flexibility parameters in each loading direction. Data were analyzed statistically to test for significant differences in flexibility parameters between the test conditions. In flexing the knee from 15 to 60 degrees, the resulting knee rotations under single moments depended upon flexion angle with varus, valgus, and internal rotations increasing significantly. Also, rotations were different depending on load direction; varus rotation was significantly different and greater than valgus rotation at both flexion angles. Also external rotation was significantly different and greater than internal at 15 degrees flexion, but not at 60 degrees flexion. Coupled rotations under single moments were also observed. Applying pure varus/valgus moments resulted in coupled external/internal rotations which were inconsistent and hence not significant. Applying pure axial moments resulted in consistent and hence significant varus/valgus rotations; an external axial moment induced varus rotation and an internal axial moment induced valgus rotation. For combination moments, varus/valgus rotations decreased significantly from those rotations at similar load levels in the single moment studies. Also, a varus moment significantly increased external rotation and a valgus moment significantly decreased internal rotation. These differences indicate significant interaction between corresponding load combinations. These results suggest that load interaction is a potentially important phenomenon in knee injury mechanics.     

Estimation of ACL forces by reproducing knee kinematics between sets of knees: A novel non-invasive methodology

In situ force in the anterior cruciate ligament (ACL) has been quantified both in vitro in response to relatively simple loads by means of robotic technology, as well as in vivo in response to more complex loads by means of force transducers and computational models. However, a methodology has been suggested to indirectly estimate the in situ forces in the ACL in a non-invasive, non-contact manner by reproducing six-degree of freedom (six-DOF) in vivo kinematics on cadaveric knees using a robotic/UFS testing system. Therefore, the objective of this study was to determine the feasibility of this approach. Kinematics from eight porcine knees (source knees) were collected at 30 degrees , 60 degrees , and 90 degrees of flexion in response to: (1) an anterior load of 100 N and (2) a valgus load of 5 N m. The average of each kinematic data set was reproduced on a separate set of eight knees (target knees). The in situ forces in the ACL were determined for both sets of knees and compared. Significant differences (rho<0.05) were found between the source knees and the target knees for all flexion angles in response to an anterior load. However, in response to valgus loads, there was no significant difference between the source knees and the target knees at 30 degrees and 90 degrees of flexion. It was noted that there was a correlation between anterior knee laxity (the distance along the displacement axis from the origin to the beginning of the linear region of the load-displacement curve) and internal-external rotation. These data suggest that in order to obtain reproducible results one needs to first match knees to knees with comparable anterior knee laxity. Thus, an estimate of the in situ forces in the ACL during in vivo activities might be obtainable using this novel methodology.     

Comparison of 3D dynamic virtual model to link segment model for estimation of net L4/L5 reaction moments during lifting

The purpose of this study was to validate a 3D dynamic virtual model for lifting tasks against a validated link segment model (LSM). A face validation study was conducted by collecting x, y, z coordinate data and using them in both virtual and LSM models. An upper body virtual model was needed to calculate the 3D torques about human joints for use in simulated lifting styles and to estimate the effect of external mechanical devices on human body. Firstly, the model had to be validated to be sure it provided accurate estimates of 3D moments in comparison to a previously validated LSM. Three synchronised Fastrak units with nine sensors were used to record data from one male subject who completed dynamic box lifting under 27 different load conditions (box weights (3), lifting techniques (3) and rotations (3)). The external moments about three axes of L4/L5 were compared for both models. A pressure switch on the box was used to denote the start and end of the lift. An excellent agreement [image omitted] was found between the two models for dynamic lifting tasks, especially for larger moments in flexion and extension. This virtual model was considered valid for use in a complete simulation of the upper body skeletal system. This biomechanical virtual model of the musculoskeletal system can be used by researchers and practitioners to give a better tool to study the causes of LBP and the effect of intervention strategies, by permitting the researcher to see and control a virtual subject's motions.     

A finite element model of the human knee joint for the study of tibio-femoral contact

As a step towards developing a finite element model of the knee that can be used to study how the variables associated with a meniscal replacement affect tibio-femoral contact, the goals of this study were 1) to develop a geometrically accurate three-dimensional solid model of the knee joint with special attention given to the menisci and articular cartilage, 2) to determine to what extent bony deformations affect contact behavior, and 3) to determine whether constraining rotations other than flexion/extension affects the contact behavior of the joint during compressive loading. The model included both the cortical and trabecular bone of the femur and tibia, articular cartilage of the femoral condyles and tibial plateau, both the medial and lateral menisci with their horn attachments, the transverse ligament, the anterior cruciate ligament, and the medial collateral ligament. The solid models for the menisci and articular cartilage were created from surface scans provided by a noncontacting, laser-based, three-dimensional coordinate digitizing system with an root mean squared error (RMSE) of less than 8 microns. Solid models of both the tibia and femur were created from CT images, except for the most proximal surface of the tibia and most distal surface of the femur which were created with the three-dimensional coordinate digitizing system. The constitutive relation of the menisci treated the tissue as transversely isotropic and linearly elastic. Under the application of an 800 N compressive load at 0 degrees of flexion, six contact variables in each compartment (ie., medial and lateral) were computed including maximum pressure, mean pressure, contact area, total contact force, and coordinates of the center of pressure. Convergence of the finite element solution was studied using three mesh sizes ranging from an average element size of 5 mm by 5 mm to 1 mm by 1 mm. The solution was considered converged for an average element size of 2 mm by 2 mm. Using this mesh size, finite element solutions for rigid versus deformable bones indicated that none of the contact variables changed by more than 2% when the femur and tibia were treated as rigid. However, differences in contact variables as large as 19% occurred when rotations other than flexion/extension were constrained. The largest difference was in the maximum pressure. Among the principal conclusions of the study are that accurate finite element solutions of tibio-femoral contact behavior can be obtained by treating the bones as rigid. However, unrealistic constraints on rotations other than flexion/extension can result in relatively large errors in contact variables.     

Motion of the knee after condylar resurfacing--an in vivo study.

The purpose of this study was determine the amounts of rotation and displacement occurring in a relatively unconstrained condylar replacement knee, and to compare the motion with a matched group of normal knees. The motion was measured using an electromagnetic device, the 3-Space Tracker, for various common activities, for 25 normal volunteers and 25 total knee patients. The main variables studied were internal-external rotation of the tibia about its long axis and anterior-posterior displacement of the femoral origin with respect to the tibia, as a function of flexion angle or per cent of gait cycle. The motion of the total knees was very similar to that of the normals. For sitting, standing and free-swing, the knee rotated internally by 5-10 degrees and the femur displaced posteriorly by 9-14 mm, as the knee was flexed from 0 to 90 degrees. For walking and going upstairs and downstairs the absolute values of the rotations and displacements were similar to the above. However, the mean values of rotation were less due to greater variation in the rotation patterns, due probably to the increased laxity of the knee during the swing phases. From these data, taking the mean motions and one standard deviation (S.D.), it is suggested that a knee prosthesis should allow a rotation from minus to plus 12 degrees (a total range of 24 degrees) and an anterior-posterior displacement of 13 mm. It is suggested that these motions be subject to progressive restraint from the neutral position in the manner of the natural knee.     

Computed tomography analysis of knee pose and geometry before and after total knee arthroplasty

Using a three-dimensional (3D) modality to image patients' knees before and after total knee arthroplasty (TKA) allows researchers and clinicians to evaluate causes of pain after TKA, differences in implant design, and changes in the articular geometry as a result of surgery. Computed tomography (CT) has not been fully utilized to date for evaluating the knee after TKA due to metal artifacts obscuring part of the image. We describe an accurate, validated protocol, which has been implemented in vivo, that improves visibility of the patellofemoral joint, matches implant models automatically in 3D, segments preoperative bone semi-automatically, detects and sets coordinate systems automatically, determines the six degrees of freedom of knee pose and geometry, and allows for multiple other measurements that are clinically relevant. Subjects are imaged at 0° and 30° knee flexion, while pushing on a custom-made knee rig to provide partial loadbearing. With some modifications, the protocol can be adopted by any group with access to a CT scanner and image analysis software, allowing for the investigation of numerous clinical and biomechanical questions.     

Statistical shape modelling reveals large and distinct subchondral bony differences in osteoarthritic knees

Knee osteoarthritis (OA) results in changes such as joint-space narrowing and osteophyte formation. Radiographic classification systems group patients by the presence or absence of these gross anatomical features but are poorly correlated to function. Statistical-shape modelling (SSM) can detect subtle differences in 3D-bone geometry, providing an opportunity for accurate predictive models. The aim of this study was to describe and quantify the main modes of shape variation which distinguish end-stage OA from asymptomatic knees. Seventy-six patients with OA and 77 control participants received a CT of their knee. 3D models of the joint were created by manual segmentation. A template mesh was fitted to all meshes and rigidly aligned resulting in a set of correspondent meshes. Principal Component Analysis (PCA) was performed to create the SSM. Logistic regression was performed on the PCA weights to distinguish morphological features of the two groups. The first 7 modes of the SSM captured >90% shape variation with 6 modes best distinguishing between OA and asymptomatic knees. OA knees displayed sub-chondral bone expansion particularly in the condyles and posterior medial tibial plateau of up to 10 mm. The model classified the two groups with 95% accuracy, 96% sensitivity, 94% specificity, and 97% AUC. There were distinct features which differentiated OA from asymptomatic knees. Further research will elucidate how magnitude and location of shape changes in the knee influence clinical and functional outcomes.     

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.     

Non-invasive assessment of soft-tissue artifact and its effect on knee joint kinematics during functional activity

The soft-tissue interface between skin-mounted markers and the underlying bones poses a major limitation to accurate, non-invasive measurement of joint kinematics. The aim of this study was twofold: first, to quantify lower limb soft-tissue artifact in young healthy subjects during functional activity; and second, to determine the effect of soft-tissue artifact on the calculation of knee joint kinematics. Subject-specific bone models generated from magnetic resonance imaging (MRI) were used in conjunction with X-ray images obtained from single-plane fluoroscopy to determine three-dimensional knee joint kinematics for four separate tasks: open-chain knee flexion, hip axial rotation, level walking, and a step-up. Knee joint kinematics was derived using the anatomical frames from the MRI-based, 3D bone models together with the data from video motion capture and X-ray fluoroscopy. Soft-tissue artifact was defined as the degree of movement of each marker in the anteroposterior, proximodistal and mediolateral directions of the corresponding anatomical frame. A number of different skin-marker clusters (total of 180) were used to calculate knee joint rotations, and the results were compared against those obtained from fluoroscopy. Although a consistent pattern of soft-tissue artifact was found for each task across all subjects, the magnitudes of soft-tissue artifact were subject-, task- and location-dependent. Soft-tissue artifact for the thigh markers was substantially greater than that for the shank markers. Markers positioned in the vicinity of the knee joint showed considerable movement, with root mean square errors as high as 29.3 mm. The maximum root mean square errors for calculating knee joint rotations occurred for the open-chain knee flexion task and were 24.3°, 17.8° and 14.5° for flexion, internal–external rotation and abduction–adduction, respectively. The present results on soft-tissue artifact, based on fluoroscopic measurements in healthy adult subjects, may be helpful in developing location- and direction-specific weighting factors for use in global optimization algorithms aimed at minimizing the effects of soft-tissue artifact on calculations of knee joint rotations.     

Correlation of patellar tracking pattern with trochlear and retropatellar surface topographies

The study was aimed to test the hypothesis that in the knee extension range 100 to 30 deg, the patellar "out-of-plane" tracking pattern is controlled by the passive restraint provided by the topographic interaction of the patellofemoral contacting surfaces. The out-of-plane tracking pattern, i.e., the pattern of patellar displacements not in the plane of knee extension/flexion, consists of translation in the medial-lateral direction, and rotations about the anterior-posterior axis (spin) and the proximal-distal axis (tilt). Using 15 fresh-frozen knees subjected to extensor moment magnitudes comparable to those in the "static-lifting" activity (foot-ground reaction = 334 N), the patellar displacements were measured using a calibrated six-degree-of-freedom electromechanical goniometer. The topographies of the trochlear and retropatellar surfaces were then measured using a calibrated traveling dial-gage arrangement and the same coordinate system used for the displacement measurements. Three indices were defined to quantify particular natural features of the three-dimensional topographies that are expected to control the patellar displacements. Correlation of the indices with their corresponding displacements showed that topographic interaction was significant in the control of all three displacements. However, for patellar spin, unlike for the other two displacements, the direction of the active quadriceps tension vector was also a significant controlling factor. Patellar medial-lateral translation was found to be controlled dominantly by the trochlear topography, while retropatellar topography also had a significant role in the control of the other two displacements.     

A reproducible method for studying three-dimensional knee kinematics

The methods used in movement analysis often rely on the definition of joint coordinate systems permitting three-dimensional (3D) kinematics. The first aim of this research project was to present a functional and postural method (FP method) to define a bone-embedded anatomical frame (BAF) on the femur and tibia, and, subsequently, a knee joint coordinate system. The repeatability of the proposed method was also assessed. Using FP method to define the BAFs, 4 kinematic parameters (flexion/extension, abduction/adduction, tibial internal/external rotation, and antero-posterior translation) were computed for 15 subjects walking on a treadmill. The repeatability for all four kinematic parameters was then assessed, using intra- and inter-observer settings. After pooling the results for all observers, the mean repeatability value ranged between 0.4 degrees and 0.8 degrees for rotation angles and between 0.8 and 2.2 mm for translation.     

A multibody knee model with discrete cartilage prediction of tibio-femoral contact mechanics

Combining musculoskeletal simulations with anatomical joint models capable of predicting cartilage contact mechanics would provide a valuable tool for studying the relationships between muscle force and cartilage loading. As a step towards producing multibody musculoskeletal models that include representation of cartilage tissue mechanics, this research developed a subject-specific multibody knee model that represented the tibia plateau cartilage as discrete rigid bodies that interacted with the femur through deformable contacts. Parameters for the compliant contact law were derived using three methods: (1) simplified Hertzian contact theory, (2) simplified elastic foundation contact theory and (3) parameter optimisation from a finite element (FE) solution. The contact parameters and contact friction were evaluated during a simulated walk in a virtual dynamic knee simulator, and the resulting kinematics were compared with measured in vitro kinematics. The effects on predicted contact pressures and cartilage–bone interface shear forces during the simulated walk were also evaluated. The compliant contact stiffness parameters had a statistically significant effect on predicted contact pressures as well as all tibio-femoral motions except flexion–extension. The contact friction was not statistically significant to contact pressures, but was statistically significant to medial–lateral translation and all rotations except flexion–extension. The magnitude of kinematic differences between model formulations was relatively small, but contact pressure predictions were sensitive to model formulation. The developed multibody knee model was computationally efficient and had a computation time 283 times faster than a FE simulation using the same geometries and boundary conditions.     

The effect of coordinate system choice and segment reference on RSA-based knee translation measures

Roentgen stereophotogrammetric analysis (RSA) can be utilized to accurately describe joint kinematics, but even when measuring small displacements within radiographically discernible structures, standardized reference frames are imperative for useful comparison across patients and across studies. In the current paper, accurately controlled laboratory models demonstrated the considerable influence that a mere 1.9-cm offset of the origin of the coordinate system from the rotation axes could exert on translation measures when rotations were occurring. In addition, the use of two different coordinate systems to gauge translation on a radiographic anterior-posterior (A-P) knee laxity exam resulted in a significant correlation (R(2)=0.562) between the two systems; however, differences of up 9.28 mm were found between corresponding measurements. This implies that clinical conclusions can potentially be upheld or refuted, based on the same data set, subject to coordinate system definition. Although the data analyzed presently involved the knee joint, similar issues surround the RSA motion analysis of other joints as well.     

Mandibular kinematics represented by a non-orthogonal floating axis joint coordinate system

There are many methods used to represent joint kinematics (e.g., roll, pitch, and yaw angles; instantaneous center of rotation; kinematic center; helical axis). Often in biomechanics internal landmarks are inferred from external landmarks. This study represents mandibular kinematics using a non-orthogonal floating axis joint coordinate system based on 3-D geometric models with parameters that are "clinician friendly" and mathematically rigorous. Kinematics data for two controls were acquired from passive fiducial markers attached to a custom dental clutch. The geometric models were constructed from MRI data. The superior point along the arc of the long axis of the condyle was used to define the coordinate axes. The kinematic data and geometric models were registered through fiducial markers visible during both protocols. The mean absolute maxima across the subjects for sagittal rotation, coronal rotation, axial rotation, medial-lateral translation, anterior-posterior translation, and inferior-superior translation were 34.10 degrees, 1.82 degrees, 1.14 degrees, 2.31, 21.07, and 6.95 mm, respectively. All the parameters, except for one subject's axial rotation, were reproducible across two motion recording sessions. There was a linear correlation between sagittal rotation and translation, the dominant motion plane, with approximately 1.5 degrees of rotation per millimeter of translation. The novel approach of combining the floating axis system with geometric models succinctly described mandibular kinematics with reproducible and clinician friendly parameters.     

Evaluation of a 3D object registration method for analysis of humeral kinematics

In 3D image-based studies of joint kinematics, 3D registration methods should be automatic, insensitive to segmentation inconsistencies and use coordinate systems that have clinically relevant orientations and locations because this is important for analyzing rotation angles and translation directions. We developed and evaluated a registration method, which is based on the cylindrical geometry of the humerus shaft and an analysis of the inertia moments of the humerus head, in order to consistently and automatically orient the humerus coordinate system according to its anatomy. Registration techniques must be thoroughly evaluated. In this study we used a well-detectable marker as reference, from which coordinate system determination errors of a 3D object could be measured. This allowed us to quantify by means of unique error analysis the translational and rotational errors in terms of how much and about/along which humeral axis errors occurred. The evaluation experiments were performed using virtual rotations of 3D humeral binary image, a humerus model and a 3D image of a volunteer's shoulder. They indicated that the humeral coordinate system determination errors primarily originated from segmentation inconsistencies, which influenced mostly the humeral transverse axes orientation. The error analysis revealed that the developed registration method reduced the effect of manual segmentation inconsistencies on the orientation of the humeral transverse axes up to 37%, in comparison to the commonly used inertia registration.