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 kinematic model of the shoulder complex to evaluate the arm-reachable workspace

Upper-arm evaluation including shoulder motion in physiotherapy has no three-dimensional tool for an arm-functioning evaluation, which hampers an uniform, objective comparison. Human shoulder complex models suffer from lack of shoulder girdle kinematic data. A kinematic shoulder-complex model with six degrees of freedom is proposed as the composition of the inner joint representing the shoulder-girdle joints and outer joint representing the glenohumeral joint. The outer shoulder joint has three perpendicular rotations: adduction/abduction, retroflexion/flexion and internal/external rotation of the humerus. The inner shoulder joint has two rotations, depression/elevation and retraction/protraction, and one translation, which are all dependent on the elevation angle of the humerus. The human arm-reachable workspace that represents the area within reach of the wrist is calculated on the basis of the shoulder-complex model and the additional elbow-joint direct kinematics. It was demonstrated that cross-sections of the calculated workspace are in agreement with the measured arm-reachable workspace in all three anatomical planes. The arm-reachable workspace volume and graphics were calculated and a comparison of the arm's workspaces during a patient's shoulder treatment was made. The obtained numerical and graphical arm-reachable workspaces can be used for arm-functioning evaluations in rehabilitation and ergonomics.     

Non-invasive determination of coupled motion of the scapula and humerus--an in-vitro validation

Measuring the motion of the scapula and humerus with sub-millimeter levels of accuracy in six-degrees-of-freedom (6-DOF) is a challenging problem. The current methods to measure shoulder joint motion via the skin do not produce clinically significant levels of accuracy. Thus, the purpose of this study was to validate a non-invasive markerless dual fluoroscopic imaging system (DFIS) model-based tracking technique for measuring dynamic in-vivo shoulder kinematics. Our DFIS tracks the positions of bones based on their projected silhouettes to contours on recorded pairs of fluoroscopic images. For this study, we compared markerlessly tracking the bones of the scapula and humerus to track them with implanted titanium spheres using a radiostereometric analysis (RSA) while manually manipulating a cadaver specimen's arms. Additionally, we report the repeatability of the DFIS to track the scapula and humerus during dynamic shoulder motion. The difference between the markerless model-based tracking technique and the RSA was ±0.3 mm in translation and ±0.5° in rotation. Furthermore, the repeatability of the markerless DFIS model-based tracking technique for the scapula and humerus was ±0.2 mm and ±0.4°, respectively. The model-based tracking technique achieves an accuracy that is similar to an invasive RSA tracking technique and is highly suited for non-invasively studying the in-vivo motion of the shoulder. This technique could be used to investigate the scapular and humeral biomechanics in both healthy individuals and in patients with various pathologies under a variety of dynamic shoulder motions encountered during the activities of daily living.     

Effect of shoulder model complexity in upper-body kinematics analysis of the golf swing

The golf swing is a complex full body movement during which the spine and shoulders are highly involved. In order to determine shoulder kinematics during this movement, multibody kinematics optimization (MKO) can be recommended to limit the effect of the soft tissue artifact and to avoid joint dislocations or bone penetration in reconstructed kinematics. Classically, in golf biomechanics research, the shoulder is represented by a 3 degrees-of-freedom model representing the glenohumeral joint. More complex and physiological models are already provided in the scientific literature. Particularly, the model used in this study was a full body model and also described motions of clavicles and scapulae. This study aimed at quantifying the effect of utilizing a more complex and physiological shoulder model when studying the golf swing. Results obtained on 20 golfers showed that a more complex and physiologically-accurate model can more efficiently track experimental markers, which resulted in differences in joint kinematics. Hence, the model with 3 degrees-of-freedom between the humerus and the thorax may be inadequate when combined with MKO and a more physiological model would be beneficial. Finally, results would also be improved through a subject-specific approach for the determination of the segment lengths.     

Shoulder joint kinematics during elevation measured by ultrasound-based measuring system.

In order to analyze shoulder joint movements, the authors use a ZEBRIS CMS-HS ultrasound-based movement analysis system. In essence, the measurement involves the determination of the spatial position of the 16 anatomical points, which are specified on the basis of the coordinates of ultrasound-based triplets positioned on the upper limb, the scapula, and the thorax; their spatial position is measured in the course of motion. Kinematic characteristics of 74 shoulder joints of 50 healthy persons were identified during elevation in the plane of the scapula. Kinematic characteristics of motion were identified by scapulothoracic, glenohumeral, and humeral elevation angles; range of angles; scapulothoracis and glenohumeral rhythm; scapulothoracic, glenohumeral, and scapuloglenoid ratios; and the relative displacement between the rotation centers of the humerus and the scapula. Motion of the humerus and the scapula relative to each other was characterized by their rotation as well as the relative displacement between the rotation centers of scapula and humerus. The biomechanical model of the shoulder joint during elevation can be described by analyzing the results of the measurements performed.     

On the kinematic modelling and the parameter estimation of the human shoulder.

This paper presents some results on the modelling and the corresponding parameter estimation of the human shoulder. This system consists of the clavicle, the scapula, the humerus and the various joints between these bodies and the trunk through the sternum; it will be represented as a succession of a rotational joint between the sternum and the clavicle and a constant distance joint, representing the scapula between, the clavicle and the humerus head. The parameters of this system are the components of the position vectors of the joint characteristic points (the corresponding centres of the rotations). Experimental results are presented as well as a validation of the proposed model.     

Towards a model for force predictions in the human shoulder.

In this paper the concept of a three-dimensional biomechanical model of the human shoulder is introduced. This model is used to analyze static load sharing between the muscles, the bones and the ligaments. The model consists of all shoulder structures, which means that different positions and different load situations may be analyzed using the same model. Solutions can be found for the complete range of shoulder motion. However, this article focuses only on elevation in the scapular plane and on forces in structures attached to the humerus. The intention is to expand the model in future studies to also involve the forces acting on the other shoulder bones: the scapula and the clavicle. The musculoskeletal forces in the shoulder complex are predicted utilizing the optimization technique with the sum of squared muscle stresses as an objective function. Numerical results predict that among the muscles crossing the glenohumeral joint parts of the deltoideus, the infraspinatus, the supraspinatus, the subscapularis, the pectoralis major, the coracobrachialis and the biceps are the muscles most activated during this sort of abduction. Muscle-force levels reached values of 150 N when the hand load was 1 kg. The results from the model seem to be qualitatively accurate, but it is concluded that in the future development of the model the direction of the contact force in the glenohumeral joint must be constrained.     

Effect of axis alignment on in vivo shoulder kinematics

BACKGROUND. To describe 3D shoulder joint movements, the International Society of Biomechanics (ISB) recommends using segment coordinate systems (SCSs) on the humerus, scapula and thorax, and joint coordinate systems (JCSs) on the shoulder. However, one of the remaining problems is how to define the zero angles when the arm is in an initial reference position. The aim of this paper is to compare various methods of determining the JCSs of the shoulder that make it possible to define the zero angles of the arm in the resting position. METHODS. Able-bodied subjects performed elevation movements in the scapular plane, specifically neutral, internal and external rotations of the humerus. The initial humerus position (at the beginning of the arm movement) and range of motion were analysed for the purpose of clinical interpretation of arm attitude and movement. The following four different JCSs were explored: (1) the standard JCS, defined as recommended by the ISB, (2) a first aligned JCS, where the humerus SCS is initially aligned with the scapula SCS, (3) a second aligned JCS, where the opposite operation is performed and 4) a third aligned JCS, where both the humerus and the scapular SCS are initially aligned with the thorax SCS. FINDINGS. The second aligned JCS was the only method that did not produce any exaggerated range of movement in either anatomical plane. INTERPRETATION. Mathematical JCS alignment allows clearer clinical interpretation of arm attitude and movement.     

Measuring dynamic in-vivo glenohumeral joint kinematics: technique and preliminary results

Rotator cuff tears are a common injury that affect a significant percentage of the population over age 60. Although it is widely believed that the rotator cuff's primary function is to stabilize the humerus against the glenoid during shoulder motion, accurately measuring the three-dimensional (3D) motion of the shoulder's glenohumeral joint under in-vivo conditions has been a challenging endeavor. In particular, conventional motion measurement techniques have frequently been limited to static or two-dimensional (2D) analyses, and have suffered from limited or unknown in-vivo accuracy. We have recently developed and validated a new model-based tracking technique that is capable of accurately measuring the 3D position and orientation of the scapula and humerus from biplane X-ray images. Herein we demonstrate the in-vivo application of this technique for accurately measuring glenohumeral joint translations during shoulder motion in the repaired and contralateral shoulders of patients following rotator cuff repair. Five male subjects were tested at 3-4 months following arthroscopic rotator cuff repair. Superior-inferior humeral translation was measured during elevation, and anterior-posterior humeral translation was measured during external rotation in both the repaired and contralateral shoulders. The data failed to detect statistically significant differences between the repaired and contralateral shoulders in superior-inferior translation (p=0.74) or anterior-posterior translation (p=0.77). The measurement technique overcomes the limitations of conventional motion measurement techniques by providing accurate, 3D, in-vivo measures of glenohumeral joint motion during dynamic activities. On-going research is using this technique to assess the effects of conservative and surgical treatment of rotator cuff tears.     

Comparison of an EMG-based and a stress-based method to predict shoulder muscle forces

The estimation of muscle forces in musculoskeletal shoulder models is still controversial. Two different methods are widely used to solve the indeterminacy of the system: electromyography (EMG)-based methods and stress-based methods. The goal of this work was to evaluate the influence of these two methods on the prediction of muscle forces, glenohumeral load and joint stability after total shoulder arthroplasty. An EMG-based and a stress-based method were implemented into the same musculoskeletal shoulder model. The model replicated the glenohumeral joint after total shoulder arthroplasty. It contained the scapula, the humerus, the joint prosthesis, the rotator cuff muscles supraspinatus, subscapularis and infraspinatus and the middle, anterior and posterior deltoid muscles. A movement of abduction was simulated in the plane of the scapula. The EMG-based method replicated muscular activity of experimentally measured EMG. The stress-based method minimised a cost function based on muscle stresses. We compared muscle forces, joint reaction force, articular contact pressure and translation of the humeral head. The stress-based method predicted a lower force of the rotator cuff muscles. This was partly counter-balanced by a higher force of the middle part of the deltoid muscle. As a consequence, the stress-based method predicted a lower joint load (16% reduced) and a higher superior–inferior translation of the humeral head (increased by 1.2 mm). The EMG-based method has the advantage of replicating the observed cocontraction of stabilising muscles of the rotator cuff. This method is, however, limited to available EMG measurements. The stress-based method has thus an advantage of flexibility, but may overestimate glenohumeral subluxation.     

Evaluation of regression-based 3-D shoulder rhythms

The movements of the humerus, the clavicle, and the scapula are not completely independent. The coupled pattern of movement of these bones is called the shoulder rhythm. To date, multiple studies have focused on providing regression-based 3-D shoulder rhythms, in which the orientations of the clavicle and the scapula are estimated by the orientation of the humerus. In this study, six existing regression-based shoulder rhythms were evaluated by an independent dataset in terms of their predictability. The datasets include the measured orientations of the humerus, the clavicle, and the scapula of 14 participants over 118 different upper arm postures. The predicted orientations of the clavicle and the scapula were derived from applying those regression-based shoulder rhythms to the humerus orientation. The results indicated that none of those regression-based shoulder rhythms provides consistently more accurate results than the others. For all the joint angles and all the shoulder rhythms, the RMSE are all greater than 5°. Among those shoulder rhythms, the scapula lateral/medial rotation has the strongest correlation between the predicted and the measured angles, while the other thoracoclavicular and thoracoscapular bone orientation angles only showed a weak to moderate correlation. Since the regression-based shoulder rhythm has been adopted for shoulder biomechanical models to estimate shoulder muscle activities and structure loads, there needs to be further investigation on how the predicted error from the shoulder rhythm affects the output of the biomechanical model.     

Morphometric affinities of the human shoulder

To analyze differences between apes and monkeys and the affinities of man, we have studied the shoulder girdle of 327 specimens of anthropoid primates. The scapula, clavicle and humerus are viewed as an integrated functional complex on the basis of 18 measurements. Several varieties of multivariate analysis show that man is clearly closer to other hominoids than to the included monkey taxa (whether terrestrial or arboreal, Old World or New World). The marked shoulder differences between apes and monkeys and similarities between apes and man correlate with the muscular anatomy, which in hominoids allows the motions involved in their locomotion and feeding behavior. As the hominid-pongid correspondence in shoulder morphology is especially detailed regarding the functionally important joint surfaces, it is consistent with a fairly recent period of common ancestry and behavior. No hypothetical evolutionary pathway or ancestral form of the human shoulder need look far beyond the model afforded by extant pongids. In contrast with previous studies on the primate shoulder, these results agree with information accumulating from other systems—comparative anatomy, primate behavior, and molecular biology — in suggesting very close relationship between man and extant African pongids.     

Effect of axis alignment on in vivo shoulder kinematics

Background. To describe 3D shoulder joint movements, the International Society of Biomechanics (ISB) recommends using segment coordinate systems (SCSs) on the humerus, scapula and thorax, and joint coordinate systems (JCSs) on the shoulder. However, one of the remaining problems is how to define the zero angles when the arm is in an initial reference position. The aim of this paper is to compare various methods of determining the JCSs of the shoulder that make it possible to define the zero angles of the arm in the resting position.

Methods. Able-bodied subjects performed elevation movements in the scapular plane, specifically neutral, internal and external rotations of the humerus. The initial humerus position (at the beginning of the arm movement) and range of motion were analysed for the purpose of clinical interpretation of arm attitude and movement. The following four different JCSs were explored: (1) the standard JCS, defined as recommended by the ISB, (2) a first aligned JCS, where the humerus SCS is initially aligned with the scapula SCS, (3) a second aligned JCS, where the opposite operation is performed and 4) a third aligned JCS, where both the humerus and the scapular SCS are initially aligned with the thorax SCS.

Findings. The second aligned JCS was the only method that did not produce any exaggerated range of movement in either anatomical plane.

Interpretation. Mathematical JCS alignment allows clearer clinical interpretation of arm attitude and movement.     

Methodology and sensitivity studies for finite element modeling of the inferior glenohumeral ligament complex

The objectives of this research were to develop a methodology for three-dimensional finite element (FE) modeling of the inferior glenohumeral ligament complex (IGHL complex) as a continuous structure, to determine optimal mesh density for FE simulations, to examine strains and forces in the IGHL complex in clinically relevant joint positions, and to perform sensitivity studies to assess the effects of assumed material properties. A simple translation test in the anterior direction was performed on a cadaveric shoulder, with the humerus oriented at 60 degrees of glenohumeral abduction and 0 degrees of flexion/extension, at 0 degrees , 30 degrees and 60 degrees of humeral external rotation. The geometries of the relevant structures were extracted from volumetric CT data to create a FE model. Experimentally measured kinematics were applied to the FE model to simulate the simple translation test. First principal strains, insertion site forces and contact forces were analyzed. At maximum anterior humeral translation, strains in the IGHL complex were highly inhomogeneous for all external rotation angles. The motion of the humerus with respect to the glenoid during the simple translation test produced a tangential load at the proximal and distal edges of the IGHL complex. This loading was primarily in the plane of the inferior glenohumeral ligament complex, producing an in-plane shear-loading pattern. There was a significant increase in strain with increasing angle of external rotation. The largest insertion site forces occurred at the axillary pouch insertion to the humerus (36.7N at 60 degrees of external rotation) and the highest contact forces were between the anterior band of the IGHL complex and the humeral cartilage (7.3N at 60 degrees of external rotation). Strain predictions were highly sensitive to changes in the ratio of bulk to shear modulus of the IGHL complex, while predictions were moderately sensitive to changes in elastic modulus of the IGHL complex. Changes to the material properties of the humeral cartilage had little effect on predicted strains. The methodologies developed in this research and the results of the mesh convergence and sensitivity studies provide a basis for the subject-specific modeling of the mechanics of the IGHL complex.     

Predicting the shapes of bones at a joint: application to the shoulder

This paper presents a novel method to explore the intrinsic morphological correlation between the bones of a shoulder joint (humerus and scapula). To model this correlation, canonical correlation analysis (CCA) is used. We also propose a technique to predict a three-dimensional (3D) bone shape from its adjoining segment at a joint based on partial least squares regression (PLS). The high dimensional 3D surface information of a bone is represented by a few variables using principal component analysis, which also captures the pattern of variability of the shapes in our datasets. Our results show that the humerus set and scapula set have highly linear morphological relationship and that the correlation information can be used as a classifier. In this study, primate shoulder bone datasets were categorised into two clusters: great apes (including humans) and monkeys. A leave one out experiment was performed to test the robustness of this prediction method. The prediction behaviour using this method shows statistically significantly better results than using the mean shape from the training set.     

A Biomechanical Model of the Scapulothoracic Joint to Accurately Capture Scapular Kinematics during Shoulder Movements

The complexity of shoulder mechanics combined with the movement of skin relative to the scapula makes it difficult to measure shoulder kinematics with sufficient accuracy to distinguish between symptomatic and asymptomatic individuals. Multibody skeletal models can improve motion capture accuracy by reducing the space of possible joint movements, and models are used widely to improve measurement of lower limb kinematics. In this study, we developed a rigid-body model of a scapulothoracic joint to describe the kinematics of the scapula relative to the thorax. This model describes scapular kinematics with four degrees of freedom: 1) elevation and 2) abduction of the scapula on an ellipsoidal thoracic surface, 3) upward rotation of the scapula normal to the thoracic surface, and 4) internal rotation of the scapula to lift the medial border of the scapula off the surface of the thorax. The surface dimensions and joint axes can be customized to match an individual’s anthropometry. We compared the model to “gold standard” bone-pin kinematics collected during three shoulder tasks and found modeled scapular kinematics to be accurate to within 2mm root-mean-squared error for individual bone-pin markers across all markers and movement tasks. As an additional test, we added random and systematic noise to the bone-pin marker data and found that the model reduced kinematic variability due to noise by 65% compared to Euler angles computed without the model. Our scapulothoracic joint model can be used for inverse and forward dynamics analyses and to compute joint reaction loads. The computational performance of the scapulothoracic joint model is well suited for real-time applications; it is freely available for use with OpenSim 3.2, and is customizable and usable with other OpenSim models.     

Main component of soft tissue artifact of the upper-limbs with respect to different functional,daily life and sports movements

Soft tissue artifact (STA) is the main source of error in kinematic estimation of human movements based on skin markers. Our objective was to determine the components of marker displacements that best describe STA of the shoulder and arm (i.e. clavicle, scapula and humerus). Four participants performed arm flexion and rotation, a daily-life and a sports movement. Three pins with reflective markers were inserted into the clavicle, scapula and humerus. In addition, up to seven skin markers were stuck on each segment. STA was described with a modal approach: individual marker displacements or marker-cluster (i.e. translations, rotations, homotheties and stretches) relative to the local segment coordinate system defined by markers secured to the pins. The modes were then ranked according to the percentage of total STA energy that they explained. Both individual skin marker displacements and marker-cluster geometrical transformations were task-, location-, segment- and subject-specific. However, 85% of the total STA energy was systematically explained by the rigid transformations (i.e. translations and rotations of the marker-cluster). In conclusion, large joint dislocations and limited efficiency of least squares bone pose estimators are expected for the computation of upper limb joint kinematics from skin markers. Future developments shall consider the rigid transformations of marker-clusters in the implementation of an STA model to reduce its effects on kinematics estimation.     

Comparative evaluation of scapular and humeral coordinate systems based on biomedical images of the glenohumeral joint

Bio-imaging techniques represent a powerful tool for shoulder joint biomechanical analysis. However, the restricted field of view may prevent the acquisition of complete scapula and humerus bone models and hence limiting the applicability of standardized anatomical coordinate system (ACS) definitions. The aim of this study was to propose ACS definitions for both scapula and humerus which can be implemented when limited portions of the relevant bones are available. Magnetic resonance (MR) images of twenty right humeri and scapulae were acquired. The proposed ACSs were assessed in terms of (1) sensitivity to bone morphological variation, (2) intra – and inter – operator repeatability and (3) consistency with the anatomical cardinal directions. A comparison with alternative ACS definitions was also performed. Overall, our ACS scapular proposal and that presented in Kedgley and Dunning (2010) were found to be the least sensitive to the morphometric variability (mean angular absolute deviation lower than 8.3 deg) and they were characterized by a high intra – and inter – operator repeatability (mean angular absolute deviation lower than 1.5 deg). The humeral ACS proposal showed a morphometric variability similar to Amadi et al. (2009b) (mean angular absolute deviation lower than 8.3 deg) but a higher reproducibility. The scapular and humeral ACS mean angular deviation from the reference anatomical cardinal directions were smaller than 15 deg and 8.6 deg, respectively. The proposed scapular and humeral ACS definitions are therefore suitable to be applied when a limited portion of the glenohumeral joint is available as it may occur in standard shoulder clinical exams.