Effects of hip joint centre mislocation on gait analysis results

Methods to determine the hip joint centre (HJC) location are necessary in gait analysis. It has been demonstrated that the methods proposed in the literature involve large mislocation errors. The choice should be made according to the extent by which HJC location errors distort the estimates of angles and resultant moments at the hip and knee joints. This study aimed at quantifying how mislocation errors propagate to these gait analysis results. Angles and moments at the hip and knee joint were calculated for five able-bodied subjects during level walking. The nominal position of the HJC was determined as the position of the pivot point of a 3D movement of the thigh relative to the pelvis. Angles and moments were then re-calculated after having added to HJC co-ordinates errors in the range of +/-30 mm. Angles and moments at both hip and knee joints were affected by HJC mislocation. The hip moments showed the largest propagation error: a 30 mm HJC anterior mislocation resulted in a propagated error into flexion/extension component of about -22%. The hip abduction/adduction moment was found the second largest affected quantity: a 30 mm lateral HJC mislocation produced a propagated error of about -15%. Finally, a 30 mm posterior HJC mislocation produced a delay of the flexion-to-extension timing in the order of 25% of the stride duration. HJC estimation methods with minimum antero-posterior error should therefore be preferred.     

In-vitro experimental assessment of a new robust algorithm for hip joint centre estimation

Hip joint centre (HJC) localization is used in several biomedical applications, such as movement analysis and computer-assisted orthopaedic surgery.The purpose of this study was to validate in vitro a new algorithm (MC-pivoting) for HJC computation and to compare its performances with the state-of-the-art (least square approach–LSA). The MC-pivoting algorithm iteratively searches for the 3D coordinates of the point belonging to the femoral bone that, during the circumduction of the femur around the hip joint (pivoting), runs the minimum length trajectory. The algorithm was initialized with a point distribution that can be considered close to a Monte Carlo simulation sampling all around the LSA estimate.The performances of the MC-pivoting algorithm, compared with LSA, were evaluated with tests on cadavers. Dynamic reference frames were applied on both the femur and the pelvis and were tracked by an optical localizer.Results proved the algorithm accuracy (1.7 mm±1.6, 2.3—median value±quartiles), reliability (smaller upper quartiles of the errors distribution with respect to LSA) and robustness (reduction of the errors also in case of large pelvis displacements).     

Subject-specific hip geometry and hip joint centre location affects calculated contact forces at the hip during gait

Hip loading affects the development of hip osteoarthritis, bone remodelling and osseointegration of implants. In this study, we analyzed the effect of subject-specific modelling of hip geometry and hip joint centre (HJC) location on the quantification of hip joint moments, muscle moments and hip contact forces during gait, using musculoskeletal modelling, inverse dynamic analysis and static optimization. For 10 subjects, hip joint moments, muscle moments and hip loading in terms of magnitude and orientation were quantified using three different model types, each including a different amount of subject-specific detail: (1) a generic scaled musculoskeletal model, (2) a generic scaled musculoskeletal model with subject-specific hip geometry (femoral anteversion, neck-length and neck-shaft angle) and (3) a generic scaled musculoskeletal model with subject-specific hip geometry including HJC location. Subject-specific geometry and HJC location were derived from CT. Significant differences were found between the three model types in HJC location, hip flexion–extension moment and inclination angle of the total contact force in the frontal plane. No model agreement was found between the three model types for the calculation of contact forces in terms of magnitude and orientations, and muscle moments. Therefore, we suggest that personalized models with individualized hip joint geometry and HJC location should be used for the quantification of hip loading. For biomechanical analyses aiming to understand modified hip joint loading, and planning hip surgery in patients with osteoarthritis, the amount of subject-specific detail, related to bone geometry and joint centre location in the musculoskeletal models used, needs to be considered.     

Prediction of the hip joint centre in adults, children, and patients with cerebral palsy based on magnetic resonance imaging

The location of the hip joint centre (HJC) is required for calculations of hip moments, the location and orientation of the femur, and muscle lengths and lever arms. In clinical gait analysis, the HJC is normally estimated using regression equations based on normative data obtained from adult populations. There is limited relevant anthropometric data available for children, despite the fact that clinical gait analysis is predominantly used for the assessment of children with cerebral palsy. In this study, pelvic MRI scans were taken of eight adults (ages 23-40), 14 healthy children (ages 5-13) and 10 children with spastic diplegic cerebral palsy (ages 6-13). Relevant anatomical landmarks were located in the scans, and the HJC location in pelvic coordinates was found by fitting a sphere to points identified on the femoral head. The predictions of three common regression equations for HJC location were compared to those found directly from MRI. Maximum absolute errors of 31 mm were found in adults, 26 mm in children, and 31 mm in the cerebral palsy group. Results from regression analysis and leave-one-out cross-validation techniques on the MRI data suggested that the best predictors of HJC location were: pelvic depth for the antero-posterior direction; pelvic width and leg length for the supero-inferior direction; and pelvic depth and pelvic width for the medio-lateral direction. For single-variable regression, the exclusion of leg length and pelvic depth from the latter two regression equations is proposed. Regression equations could be generalised across adults, children and the cerebral palsy group.     

Reduction of knee range of motion during continuous passive motion due to misaligned hip joint centre

When using continuous passive motion (CPM) devices, appropriate setting of the device and positioning of the patient are necessary to obtain maximum range of motion (ROM). In this study, the ROMs in both the knee joint and CPM device during CPM treatment were measured using a motion analysis system for three different CPM devices. Additionally, the trajectories of the angles at the knee for hip joint misalignments were evaluated using kinematic models of the three CPM devices. The results showed that discrepancies in ROM between the knee joints and the CPM device settings during CPM treatment were revealed regardless of the CPM device and that the effect of misalignment is dependent on the design of the CPM device. The present technology could be applied for the development of a better design configuration for the CPM device to reduce the discrepancy in ROM at the knee joint.     

In vivo estimation of the glenohumeral joint centre by functional methods: Accuracy and repeatability assessment

Several algorithms have been proposed for determining the centre of rotation of ball joints. These algorithms are used rather to locate the hip joint centre. Few studies have focused on the determination of the glenohumeral joint centre. However, no studies have assessed the accuracy and repeatability of functional methods for glenohumeral joint centre.This paper aims at evaluating the accuracy and the repeatability with which the glenohumeral joint rotation centre (GHRC) can be estimated in vivo by functional methods. The reference joint centre is the glenohumeral anatomical centre obtained by medical imaging. Five functional methods were tested: the algorithm of Gamage and Lasenby (2002), bias compensated (Halvorsen, 2003), symmetrical centre of rotation estimation (Ehrig et al., 2006), normalization method (Chang and Pollard, 2007), helical axis (Woltring et al., 1985). The glenohumeral anatomical centre (GHAC) was deduced from the fitting of the humeral head.Four subjects performed three cycles of three different movements (flexion/extension, abduction/adduction and circumduction). For each test, the location of the glenohumeral joint centre was estimated by the five methods. Analyses focused on the 3D location, on the repeatability of location and on the accuracy by computing the Euclidian distance between the estimated GHRC and the GHAC. For all the methods, the error repeatability was inferior to 8.25 mm. This study showed that there are significant differences between the five functional methods. The smallest distance between the estimated joint centre and the centre of the humeral head was obtained with the method of Gamage and Lasenby (2002).     

Accuracy of CT scans in identifying tumor tissue

The authors present their experience with stereotactic biopsy of brain tumors. Biopsies were obtained sequentially at different depths from the center of the tumor according to coordinates derived from computerized tomography (CT). Biopsies were obtained from 23 brain tumors: 17 gliomas, 5 metastases, and 1 lymphoma. In all a total of 137 biopsies were studied from both enhancing and nonenhancing areas. The tumor yield from these biopsies was 68 and 73%, respectively. It appears that tumor tissue may be obtained from both the enhancing periphery as well as the nonenhancing center of tumors.     

Detection and estimation of neural connectivity based on crosscorrelation analysis

Crosscorrelation analysis of simultaneously recorded activity of pairs of neurons is a common tool to infer functional neural connectivity. The adequacy of crosscorrelation procedures to detect and estimate neural connectivity has been investigated by means of computer simulations of small networks composed of fairly realistic modelneurons. If the mean interval of neural firings is much larger than the duration of postsynaptic potentials, which will be the case in many central brain areas excitatory connections are easier to detect than inhibitory ones. On the other hand, inhibitory connections are revealed better if the mean firing interval is much smaller than post-synaptic potential duration. In general the effects of external stimuli and the effects of neural connectivity do not add linearly. Furthermore, neurons may exhibit a certain degree of timelock to the stimulus. For these reasons the commonly applied shift predictor procedure to separate stimulus and neural effects appears to be of limited value. In case of parallel direct and indirect neural pathways between two neurons crosscorrelation analysis does not estimate the direct connection but instead an effective connectivity, which reflects the combined influences of the parallel pathways.Abbreviations ACH autocoincidence histogram - CCH crosscoincidence histogram - CPDF crossproduct density function - DCH difference crosscoincidence histogram - NACH nonsimultaneous autocoincidence histogram - NCCH nonsimultaneous crosscoincidence histogram - SCCH scaled crosscoincidence histogram - SDCH scaled difference crosscoincidence histogram     

A computer-based biomechanical analysis of the three-dimensional motion of cementless hip prostheses.

A computer-based mathematical technique was developed to measure and completely describe the migration and micromotion of a femoral hip prosthesis relative to the femur. This technique utilized the mechanics of rigid-body motion analysis and apparatus of seven linear displacement transducers to measure and describe the complete three-dimensional motion of the prosthesis during cyclic loading. Computer acquisition of the data and custom analysis software allowed one to calculate the magnitude and direction of the motion of any point of interest on the prostheses from information about the motion of two points on the device. The data were also used to replay the tests using a computer animation technique, which allowed a magnified view of the three-dimensional motion of the prosthesis. This paper describes the mathematical development of the rigid-body motion analysis, the experimental method and apparatus for data collection, the technique used to animate the motion, the sources of error and the effect of the assumptions (rigid bodies) on the results. Selected results of individual test runs of uncemented and cemented prostheses are presented to demonstrate the efficacy of the method. The combined effect of the vibration and electrical noise resulted in a resolution of the system of about 3-5 microns motion for each transducer. Deformation effects appear to contribute about 3-15 microns to the measurement error. This measurement and analysis technique is a very sensitive and powerful means of assessing the effects of different design parameters on the migration and micromotion of total joint prostheses and can be applied to any other case (knee, dental implant) where three-dimensional relative motion between two bodies is important.     

Reproducibility of geometrical acquisition of intra-thoracic organs of children on CT scans

This paper analyses geometry of intra-thoracic organs from computed tomography (CT) scans performed on 20 children aged from 4 months to 16 years. A set of two measurements on lungs and heart were performed by the same observer. A third set was performed by a second observer. Thus, the intra- and inter-observer relative deviation of measurements was analysed. Multiple regressions were used in order to study the relationship between the CT properties (scanner, voltage, dose, pixel size, slice increment) and the relative deviation of measurements. There is a very low systematic intra- and inter-observer bias in measurements except for the volume of the heart. None of the CT data properties has a significant influence on the relative deviation of measurement. In the present paper, the measurements and 3D reconstruction protocol described can be applied to characterise the growth of the intra-thoracic organs.     

MRI-based assessment of hip joint translations

To better understand movement limitations and, to some extent, the pathogenesis of osteoarthritis, it is important to quantitatively measure femoroacetabular translations to assess if any joint subluxation occurs. In this paper, we aim at measuring hip joint displacements from magnetic resonance images (MRI) based on a surface registration technique. Because this measurement is related to the location of the hip joint center (HJC), we investigate and compare different HJC estimation approaches based on patient-specific 3D bone models. We estimate the HJC based on a simulated circumduction while minimizing inter-articular distance changes. Measurements of femoroacetabular translations during low amplitude abductions ($80$ samples) and extreme flexions ($60$ samples) in female professional dancers, which is a population potentially exposed to femoroactebaluar impingements, do not show any significant subluxation.     

Effect of hip joint angle at seat-off on hip joint contact force during sit-to-stand movement: a computer simulation study

Background

Sit-to-stand movements are a necessary part of daily life, and excessive mechanical stress on the articular cartilage has been reported to encourage the progression of osteoarthritis. Although a change in hip joint angle at seat-off may affect hip joint contact force during a sit-to-stand movement, the effect is unclear. This study aimed to examine the effect of the hip joint angle at seat-off on the hip joint contact force during a sit-to-stand movement by using a computer simulation.

Methods

A musculoskeletal model was created for the computer simulation, and eight muscles were attached to each lower limb. Various sit-to-stand movements were generated using parameters (e.g., seat height and time from seat-off to standing posture) reported by previous studies. The hip joint contact force for each sit-to-stand movement was calculated. Furthermore, the effect of the hip joint angle at seat-off on the hip joint contact force during the sit-to-stand movement was examined. In this study, as the changes to the musculoskeletal model parameters affect the hip joint contact force, a sensitivity analysis was conducted.

Results and conclusions

The hip joint contact force during the sit-to-stand movement increased approximately linearly as the hip flexion angle at the seat-off increased. Moreover, the normal sit-to-stand movement and the sit-to-stand movement yielding a minimum hip joint contact force were approximately equivalent. The effect of the changes to the musculoskeletal model parameters on the main findings of this study was minimal. Thus, the main findings are robust and may help prevent the progression of hip osteoarthritis by decreasing mechanical stress, which will be explored in future studies.

     

Sensitivity analysis of a cemented hip stem to implant position and cement mantle thickness

Patient-specific finite element models of the implanted proximal femur can be built from pre-operative computed tomography scans and post-operative X-rays. However, estimating three-dimensional positioning from two-dimensional radiographs introduces uncertainty in the implant position. Further, accurately measuring the thin cement mantle and the degree of cement–bone interdigitation from imaging data is challenging. To quantify the effect of these uncertainties in stem position and cement thickness, a sensitivity study was performed. A design-of-experiment study was implemented, simulating both gait and stair ascent. Cement mantle stresses and bone–implant interface strains were monitored. The results show that small variations in alignment affect the implant biomechanics, especially around the most proximal and most distal ends of the stem. The results suggest that implant position is more influential than cement thickness. Rotation around the medial–lateral axis is the dominant factor in the proximal zones and stem translations are the dominant factors around the distal tip.     

A non-invasive acoustic and vibration analysis technique for evaluation of hip joint conditions

The performance evaluation of THA outcome is difficult and surgeons often use invasive methods to investigate effectiveness. A non-invasive acoustic and vibration analysis technique has recently been developed for more-in-depth evaluation of in vivo hip conditions.Gait kinematics, corresponding vibration and sound measurement of five THA subjects were analyzed post-operatively using video-fluoroscopy, sound and accelerometer measurements while walking on a treadmill. The sound sensor and a pair of tri-axial accelerometers, externally attached to the pelvic and femoral bone prominences, detected frequencies that are propagated through the femoral head and acetabular cup interactions. A data acquisition system was used to amplify the signal and filter out noise generated by undesired frequencies. In vivo kinematics and femoral head sliding quantified using video fluoroscopy were correlated to the sound and acceleration measurements.Distinct variations between the different subjects were identified. A correlation of sound and acceleration impulses with separation has been achieved. Although, in vivo sounds are quite variable in nature and all correlated well with the visual images.This is the first study to document and correlate visual and audible effects of THA under in-vivo conditions. This study has shown that the development of the acoustic and vibration technique provides a practical method and generates new possibilities for a better understanding of THA performance.     

Sensitivity of phylogeny estimation to taxonomic sampling

Recent studies have shown that addition or deletion of taxa from a data matrix can change the estimate of phylogeny. I used 29 data sets from the literature to examine the effect of taxon sampling on phylogeny estimation within data sets. I then used multiple regression to assess the effect of number of taxa, number of characters, homoplasy, strength of support, and tree symmetry on the sensitivity of data sets to taxonomic sampling. Sensitivity to sampling was measured by mapping characters from a matrix of culled taxa onto optimal trees for that reduced matrix and onto the pruned optimal tree for the entire matrix, then comparing the length of the reduced tree to the length of the pruned complete tree. Within-data-set patterns can be described by a second-order equation relating fraction of taxa sampled to sensitivity to sampling. Multiple regression analyses found number of taxa to be a significant predictor of sensitivity to sampling; retention index, number of informative characters, total support index, and tree symmetry were nonsignificant predictors. I derived a predictive regression equation relating fraction of taxa sampled and number of taxa potentially sampled to sensitivity to taxonomic sampling and calculated values for this equation within the bounds of the variables examined. The length difference between the complete tree and a subsampled tree was generally small (average difference of 0-2.9 steps), indicating that subsampling taxa is probably not an important problem for most phylogenetic analyses using up to 20 taxa.     

A three-dimensional investigation of temporomandibular joint loading

A mathematical model of the muscle groups applied to the human mandible is developed to study the forces developed on the condyles during maximum unilateral occlusion. The results show that the reaction forces are in approximately a 2:1 ratio with the balancing side condyle carrying the greater load. Furthermore, the direction in which these condylar reactions occur is presented.     

Statistical estimation of gene expression using multiple laser scans of microarrays

We propose a statistical model for estimating gene expression using data from multiple laser scans at different settings of hybridized microarrays. A functional regression model is used, based on a non-linear relationship with both additive and multiplicative error terms. The function is derived as the expected value of a pixel, given that values are censored at 65 535, the maximum detectable intensity for double precision scanning software. Maximum likelihood estimation based on a Cauchy distribution is used to fit the model, which is able to estimate gene expressions taking account of outliers and the systematic bias caused by signal censoring of highly expressed genes. We have applied the method to experimental data. Simulation studies suggest that the model can estimate the true gene expression with negligible bias. AVAILABILITY: FORTRAN 90 code for implementing the method can be obtained from the authors.     

Determining the location of hip joint centre: application of a conchoid's shape to the acetabular cartilage surface of magnetic resonance images

Preoperative planning, or intraoperative navigation of hip surgery, including joint-preserving procedures such as osteotomy or joint-replacing procedures such as total arthroplasty, needs to be performed with a high degree of accuracy to ensure a successful outcome. The ability to precisely localise the hip joint rotation centre may prove to be very useful in this context. The human hip joint has been shown to be a conchoid shape, and therefore the accurate location of the hip joint centre (HJC) cannot be computed simply as the centre of a sphere. This study describes a method for determining the HJC by applying a conchoid shape to the acetabular cartilage surface of magnetic resonance images, in order to increase the accuracy of the HJC location which had previously been calculated by a functional method using reconstructed three-dimensional surface bony models. By approximating a conchoid shape to the acetabulum, it was possible to compensate for HJC calculation errors.