Database of segmentations and surface models of bones of the entire lower body created from cadaver CT scans.

The range of applications of digital surface models of the bones in science and industry is wide. Three-dimensional reconstructions of bones are used in biomechanics, biomedical engineering, medical image processing, orthopedics, traumatology, radiology, patient education, anatomy, anthropometry, forensic anthropology, ergonomics, usability and human factors engineering, or accident and injury analysis and prevention. No open access database or repository of skeletal surface models of the full lower extremities exists. Therefore, the objective of this publication was to provide access to consistent complete bone models of the pelvis and lower limbs of multiple subjects, including biometric data. Segmentations and surface models of the bones of the lower extremities of more than twenty subjects were created from open access postmortem whole-body computed tomography scans. The database provides a broad range of applications by giving access to the data of the complete process chain, from the raw medical imaging data through the segmentations to the surface models.

Implications of the uncertainty of postoperative functional parameters for the preoperative planning of total hip arthroplasty.

The functional parameters pelvic tilt (PT) and hip joint force (HJF) are required to calculate patient-specific target zones based on the range of motion (ROM) and implant loading for preoperative planning of total hip arthroplasty (THA). Both functional parameters may change after THA. The preoperative prediction of the postoperative PT and HJF is associated with a specific amount of uncertainty. The prediction uncertainty has to be considered in the preoperative planning process to avoid a suboptimal implantation. So far, very little attention has been paid to the necessary reduction of patient-specific target zones by the prediction uncertainties of postoperative functional parameters. Prediction models for the postoperative PT in standing position and for the HJF during one-leg stance as a surrogate for the peak force phase during level walking were used to quantify the reduction of the ROM- and load-based target zones of 196 Japanese THA patients. The prediction uncertainty was about 14° for the postoperative standing PT and ranged from 17% body weight to 37% body weight for the components of the HJF. On average, the prosthetic ROM-based target zone had to be significantly reduced by 43% and the load-based target zone by 39%. This led to a median reduction of the combined prosthetic ROM- and load-based target zone of 96%. The study sharpens the awareness for the substantial reduction of ROM- and load-based target zones by prediction uncertainties of the postoperative PT and HJF and highlights the importance of further research to improve prediction models for both functional parameters.

Effect of the underlying cadaver data and patient-specific adaptation of the femur and pelvis on the prediction of the hip joint force estimated using static models.

The prediction of the hip joint force (HJF) is a fundamental factor for the prevention of edge loading in total hip arthroplasty. Naturally, the loading of the liner of the acetabular component depends on the HJF acting on the artificial joint. In contrast to dynamic musculoskeletal models, static models for HJF prediction do not require motion analysis of the patient. However, patient-specific adaptability and validity of static models have to be scrutinized. In this study, a modular framework for HJF prediction using static models is introduced to compare the results of different cadaver templates that are the basis of most static and dynamic models, and different scaling laws for the patient-specific adaptation with in vivo HJF of ten patients for one-leg stance and level walking. The results revealed the significant effect of the underlying cadaver template used for the prediction of the HJF (p < 0.01). A higher degree of patient-specific scaling of the cadaver template often did not significantly reduce the prediction error. Three static models with the lowest prediction errors were compared to results of dynamic models from literature. The prediction error of the peak HJF of the static models (median absolute errors below 15% body weight in magnitude and below 5° in direction) was similar in magnitude and even smaller in direction compared to dynamic models. The necessary reduction of a load-based target zone for the prevention of edge loading due to the uncertainty of the HJF prediction has to be considered in the preoperative planning. The framework for HJF prediction is openly accessible at https://github.com/RWTHmediTEC/HipJointForceModel.

The Patient-Specific Combined Target Zone for Morpho-Functional Planning of Total Hip Arthroplasty.

Background Relevant criteria for total hip arthroplasty (THA) planning have been introduced in the literature which include the hip range of motion, bony coverage, anterior cup overhang, leg length discrepancy, edge loading risk, and wear. The optimal implant design and alignment depends on the patient's anatomy and patient-specific functional parameters such as the pelvic tilt. The approaches proposed in literature often consider one or more criteria for THA planning. but to the best of our knowledge none of them follow an integrated approach including all criteria for the definition of a patient-specific combined target zone (PSCTZ). Questions/purposes (1) How can we calculate suitable THA implant and implantation parameters for a specific patient considering all relevant criteria? (2) Are the resulting target zones in the range of conventional safe zones? (3) Do patients who fulfil these combined criteria have a better outcome score? Methods A method is presented that calculates individual target zones based on the morphology, range of motion and load acting on the hip joint and merges them into the PSCTZ. In a retrospective analysis of 198 THA patients, it was calculated whether the patients were inside or outside the Lewinnek safe zone, Dorr combined anteversion range and PSCTZ. The postoperative Harris Hip Scores (HHS) between insiders and outsiders were compared. Results 11 patients were inside the PSCTZ. Patients inside and outside the PSCTZ showed no significant difference in the HHS. However, a significant higher HHS was observed for the insiders of two of the three sub-target zones incorporated in the PSCTZ. By combining the sub-target zones in the PSCTZ, all PSCTZ insiders except one had an HHS higher than 90. Conclusions The results might suggest that, for a prosthesis implanted in the PSCTZ a low outcome score of the patient is less likely than using the conventional safe zones by Lewinnek and Dorr. For future studies, a larger cohort of patients inside the PSCTZ is needed which can only be achieved if the cases are planned prospectively with the method introduced in this paper. Clinical Relevance The method presented in this paper could help the surgeon combining multiple different criteria during THA planning and find the suitable implant design and alignment for a specific patient.

A robust method for automatic identification of femoral landmarks, axes, planes and bone coordinate systems using surface models.

The identification of femoral landmarks is a common procedure in multiple academic fields. Femoral bone coordinate systems are used particularly in orthopedics and biomechanics, and are defined by landmarks, axes and planes. A fully automatic detection overcomes the drawbacks of a labor-intensive manual identification. In this paper, a new automatic atlas- and a priori knowledge-based approach that processes femoral surface models, called the A&A method, was evaluated. The A&A method is divided in two stages. Firstly, a single atlas-based registration maps landmarks and areas from a template surface to the subject. In the second stage, landmarks, axes and planes that are used to construct several femoral bone coordinate systems are refined using a priori knowledge. Three common femoral coordinate systems are defined by the landmarks detected. The A&A method proved to be very robust against a variation of the spatial alignment of the surface models. The results of the A&A method and a manual identification were compared. No significant rotational differences existed for the bone coordinate system recommended by the International Society of Biomechanics. Minor significant differences of maximally 0.5° were observed for the two other coordinate systems. This might be clinically irrelevant, depending on the context of use and should, therefore, be evaluated by the potential user regarding the specific application. The entire source code of the A&A method and the data used in the study is open source and can be accessed via https://github.com/RWTHmediTEC/FemoralCoordinateSystem .

Preoperative factors improving the prediction of the postoperative sagittal orientation of the pelvis in standing position after total hip arthroplasty.

The aims of this study were to investigate if the sagittal orientation of the pelvis (SOP) in the standing position changes after total hip arthroplasty (THA) and evaluate what preoperative factors may improve the prediction of the postoperative standing SOP in the context of a patient-specific functional cup orientation. 196 primary THA patients from Japan were retrospectively selected for this study. Computed tomography imaging of the pelvis, EOS imaging of the lower body and lateral radiographs of the lumbar spine in the standing position were taken preoperatively. Common biometrics and preoperative Harris Hip Score were recorded. The EOS imaging in the standing position was repeated three months following THA. A 3D/2.5D registration process was used to determine the standing SOP. Thirty-three preoperative biometric, morphological and functional parameters were measured. Important preoperative parameters were identified that significantly improve the prediction of the postoperative standing SOP by using multiple linear LASSO regression. On average, the SOP changed significantly (p < 0.001) between the preoperative and postoperative standing position three months after THA by 3° ± 4° in the posterior direction. The age, standing lumbar lordosis angle (LLA) and preoperative supine and standing SOP significantly (p < 0.001) improve the prediction of the postoperative standing SOP. The linear regression model for the prediction of the postoperative standing SOP is significantly (p < 0.001) improved by adding the parameters preoperative standing SOP and LLA, in addition to the preoperative supine SOP, reducing the root mean square error derived from a leave-one-out cross-validation by more than 1°. The mean standing SOP in Japanese patients changes already three months after THA in comparison to the preoperative value. The preoperative factors age, LLA, supine and standing SOP can significantly improve the prediction of the postoperative standing SOP and should be considered within the preoperative planning process of a patient-specific functional cup orientation.

A robust method for automatic identification of landmarks on surface models of the pelvis.

The recognition of bony landmarks of the pelvis is a required operation in patient-specific orthopedics, subject-specific biomechanics or morphometrics. A fully automatic detection is preferable to a subjective and time-consuming manual identification. In this paper, a new approach, called the iterative tangential plane (ITP) method, for fully automatic identification of landmarks on surface models of the pelvis is introduced. The method includes the landmarks to construct the two most established anatomical reference frames of the pelvis: the anterior pelvic plane (APP) coordinate system and superior inferior spine plane (SISP) coordinate system. The ITP method proved to be robust against the initial alignment of the pelvis in space. A comparison to a manual identification was performed that showed minor but significant (p < 0.05) median differences below 3 mm for the position of the landmarks and below 1° for the orientation of the APP coordinate system. Whether these differences are acceptable, has to be evaluated for each specific use case. There were no significant differences for the orientation of the SISP coordinate system recommended by the International Society of Biomechanics.

Patient-specific musculoskeletal modeling of the hip joint for preoperative planning of total hip arthroplasty: A validation study based on in vivo measurements.

Validation of musculoskeletal models for application in preoperative planning is still a challenging task. Ideally, the simulation results of a patient-specific musculoskeletal model are compared to corresponding in vivo measurements. Currently, the only possibility to measure in vivo joint forces is to implant an instrumented prosthesis in patients undergoing a total joint replacement. In this study, a musculoskeletal model of the AnyBody Modeling System was adapted patient-specifically and validated against the in vivo hip joint force measurements of ten subjects performing one-leg stance and level walking. The impact of four model parameters was evaluated; hip joint width, muscle strength, muscle recruitment, and type of muscle model. The smallest difference between simulated and in vivo hip joint force was achieved by using the hip joint width measured in computed tomography images, a muscle strength of 90 N/cm2, a third order polynomial muscle recruitment, and a simple muscle model. This parameter combination reached mean deviations between simulation and in vivo measurement during the peak force phase of 12% ± 14% in magnitude and 11° ± 5° in orientation for one-leg stance and 8% ± 6% in magnitude and 10° ± 5° in orientation for level walking.