Automated design of bone-preserving, insertable, and shape-matching patient-specific acetabular components.

Effective treatment of large acetabular defects remains among the most challenging aspects of revision total hip arthroplasty (THA), due to the deficiency of healthy bone stock and degradation of the support columns. Generic uncemented components, which are favored in primary THA, are often unsuitable in revision cases, where the bone-implant contact may be insufficient for fixation, without significant reaming of the limited residual bone. This study presents a computational design strategy for automatically generating patient-specific implants that simultaneously maximize the bone-implant contact area, and minimize bone reaming while ensuring insertability. These components can be manufactured using the same additive manufacturing methods as porous components and may reduce cost and operating-time, compared to existing patient-specific systems. This study compares the performance of implants generated via the proposed method to optimally fitted hemispherical implants, in terms of the achievable bone-implant contact surface, and the volume of reamed bone. Computer-simulated results based on the reconstruction of a set of 15 severe pelvic defects (Paprosky 2A-3B) suggest that the patient-specific components increase bone-implant contact by 63% (median: 63%; SD: 44%; 95% CI: 52.3%-74.0%; RMSD: 42%), and reduce the volume of reamed bone stock by 97% (median: 98%; SD: 4%; 95% CI: 95.9%-97.4%; RMSD: 3.7%).

Quantitative statistical shape model-based analysis of humeral head migration, Part 2: Shoulder osteoarthritis.

We wanted to investigate the quantitative characteristics of humeral head migration (HHM) in shoulder osteoarthritis (OA) and their possible associations with scapular morphology. We quantified CT-scan-based-HHM in 122 patients with a combination of automated 3D scapulohumeral migration (=HHM with respect to the scapula) and glenohumeral migration (=HHM with respect to the glenoid) measurements. We divided OA patients in Group 1 (without HHM), Group 2a (anterior HHM) and Group 2b (posterior HHM). We reconstructed and measured the prearthropathy scapular anatomy with a statistical shape model technique. HHM primarily occurs in the axial plane in shoulder OA. We found "not-perfect" correlation between subluxation distance AP and scapulohumeral migration values (rs = 0.8, p < 0.001). Group 2b patients had a more expressed prearthropathy glenoid retroversion (13° vs. 7°, p < 0.001) and posterior glenoid translation (4 mm vs. 6 mm, p = 0.003) in comparison to Group 1. Binary logistic regression analysis indicated prearthropathy glenoid version as a significant predictor of HHM (χ² = 27, p < 0.001). Multivariate regression analysis showed that the pathologic version could explain 56% of subluxation distance-AP variance and 75% of the scapulohumeral migration variance (all p < 0.001). Herewith, every degree increase in pathologic glenoid retroversion was associated with an increase of 1% subluxation distance-AP, and scapulohumeral migration. The occurrence of posterior HHM is associated with prearthropathy glenoid retroversion and more posterior glenoid translation. The reported regression values of HHM in the function of the pathologic glenoid version could form a basis toward a more patient-specific correction of HHM.

Advanced quantitative 3D imaging improves the reliability of the classification of acetabular defects.

INTRODUCTION: Classifying complex acetabular defects in revision total hip arthroplasty (THA) by means of conventional radiographs comes with significant limitations. Statistical shape modelling allows the virtual reconstruction of the native pelvic morphology, hereby enabling an analytic acetabular defect assessment. Our objective was to evaluate the effect of advanced imaging augmented with analytic representations of the defect on (1) intra- and inter-rater reliability, and (2) up- or downscaling of classification scores when evaluating acetabular defects in patients undergoing revision THA. MATERIALS AND METHODS: The acetabular defects of 50 patients undergoing revision THA were evaluated by three independent, fellowship-trained orthopaedic surgeons. Defects were classified according to the acetabular defect classification (ADC) using four different imaging-based representations, namely, standard radiographs, CT imaging, a virtual three-dimensional (3D) model and a quantitative analytic representation of the defect based on a statistical shape model reconstruction. Intra- and inter-rater reliabilities were quantified using Fleiss' and Cohen's kappa scores, respectively. Up- and downscaling of classification scores were compared for each of the imaging-based representations and differences were tested. RESULTS: Overall inter-rater agreement across all imaging-based representations for the classification was fair (κ 0.29 95% CI 0.28-0.30). Inter-rater agreement was lowest for radiographs (κ 0.21 95% CI 0.19-0.22) and increased for other representations with agreement being highest when using analytic defect models (κ 0.46 95% CI 0.43-0.48). Overall intra-rater agreement was moderate (κ 0.51 95% CI 0.42-0.60). Intra-rater agreement was lowest for radiographs (κ 0.40 95% CI 0.23-0.57), and highest for ratings including analytic defect models (κ 0.64:95% CI 0.46-0.82). Virtual 3D models with quantitative analytic defect representations upscaled acetabular defect scores in comparison to standard radiographs. CONCLUSIONS: Using 3D CT imaging with statistical shape models doubles the intra- and inter-rater reliability and results in upscaling of acetabular defect classification when compared to standard radiographs. This method of evaluating defects will aid in planning surgical reconstruction and stimulate the development of new classification systems based on advanced imaging techniques.

Quantitative SSM-based analysis of humeral head migration in rotator cuff tear arthropathy patients.

Rotator cuff tear arthropathy (RCTA) is characterized by massive rotator cuff tearing combined with humeral head migration (HHM). The aim of this study is to investigate the quantitative characteristics of this migration and its association with glenoid erosions and prearthropathy scapular anatomy. We quantified HHM and prearthropathy scapular anatomy of 64 RCTA patients with statistical shape modeling-based techniques. Glenoid erosion was classified according to Sirveaux et al. A cutoff value for confirming HHM was 5 mm based on a control group of 49 patients. Group 1 (RCTA without HHM) consisted of 21 patients, with a mean subluxation distance (SLD) of 3 mm. Group 2 (RCTA with HHM) consisted of 43 patients, with mean SLD of 9 mm, SLD in the anteroposterior plane of -1 mm (SD ± 4 mm), SLD in the superoinferior plane of 7 mm (SD ± 3 mm), and subluxation angle (SLA) of -5° (SD ± 40°). Analysis with Fisher's exact test showed a clear association between HHM and glenoid erosions (p = 0.002). Multivariate regression analysis of Group 2 showed that prearthropathy lateral acromial angle combined with critical shoulder angle (p = 0.004) explained 21% of the observed variability in SLD. The prearthropathy glenoid version explained 23% of the variability in SLA (p = 0.001). HHM in RCTA patients has a wide variation in both magnitude and direction leading to a distorted glenohumeral relationship in the coronal and axial plane. HHM is highly associated with the occurrence of glenoid erosions. There is a correlation between the prearthropathy scapular anatomy and the magnitude and direction of HHM.

The External Obturator Footprint Is a Usable, Accurate, and Reliable Landmark for Stem Depth in Direct Anterior THA.

BACKGROUND: Previous CT and cadaver studies have suggested that the external obturator footprint might be used as a landmark for stem depth in direct anterior THA. Instructions on where to template this structure with small variability in height have been developed but have not been tested in daily clinical practice. QUESTIONS/PURPOSES: In this study we sought to investigate the (1) usability, (2) accuracy, and (3) reliability of the external obturator footprint as a landmark for stem depth in direct anterior THA. METHODS: The distance between the superior border of the external obturator tendon and the shoulder of the stem was measured intraoperatively in all patients (n = 135) who underwent primary THA via a direct anterior approach performed by the senior author between November 2019 and October 2020. The landmark was considered useful when two of thre`e evaluators agreed that the intersection of the vertical line comprised of the lateral wall of the trochanteric fossa and the oblique line formed by the intertrochanteric crest was clearly visible on the preoperative planning radiograph, and when the landmark was furthermore identified with certainty during surgery. Accuracy was defined as the degree of agreement (categorical for thresholds of 2 and 5 mm, the latter representing the threshold for developing unphysiological gait parameters) between the intraoperative distance and radiographic distance as measured on intraoperative fluoroscopy images or postoperative radiographs, which were calibrated based on femoral head sizes in a software program commonly used for templating. Intrarater reliability was defined as the degree of agreement (categorical for thresholds of 1 mm, which we considered an acceptable measurement error) between the ratings of one observer, who measured the radiographic distance on two different occasions separated by a washout period of at least 2 weeks. Interrater reliability was defined as the degree of agreement (categorical for thresholds of 1 mm, which we considered an acceptable measurement error) between the ratings of three observers with varying levels of experience (a fellowship-trained hip surgeon, a hip surgery fellow, and a medical student). RESULTS: The landmark was considered useful in 77% (104 of 135) of patients who underwent direct anterior THA based on the observations that the trochanteric fossa was clearly visible on the planning radiograph in 117 patients and that the tendon was identified with certainty during surgery in 118 patients. There was good-to-excellent accuracy (intraclass correlation coefficient 0.75-087), and intrarater reliability (ICC 0.99) and interrater reliability (ICC 0.99) were both excellent. CONCLUSION: This clinical study showed that the external obturator footprint is a useful, accurate, and reliable landmark for stem depth in direct anterior THA. CLINICAL RELEVANCE: The external obturator landmark allows the surgeon to position the stem within a range of the templated depth that is beneath the threshold for the development of unphysiological gait parameters. Although strictly speaking it was found useful in 77% of patients in this study, we found that this percentage of usability can easily be improved to around 90% by providing the radiology lab technician with instructions to correct external rotation of the foot during the taking of the planning radiograph. Future studies could compare the established (in)equality in leg length in patients using the external obturator landmark with computer-assisted surgery.

Determination of predisposing scapular anatomy with a statistical shape model-Part II: shoulder osteoarthritis.

HYPOTHESIS AND BACKGROUND: Shoulder osteoarthritis can be divided into different glenoid types (A, B, C, and D) and subtypes. The aim of this study was to investigate if there is an association between the prearthropathy scapular anatomy, shoulder osteoarthritis, and different glenoid types and subtypes. METHODS: Using principal components analysis, a statistical shape model (SSM) of the scapula was constructed from a data set of 110 computed tomographic (CT) scans. These subjects formed the control group. Next, CT scan images of 117 patients with osteoarthritis were classified according to the modified Walch classification. A complete 3-dimensional (3D) scapular bone model was created for every patient, and using the SSM, a reconstruction of their prearthropathy scapular anatomy was performed. Automated 3D measurements were performed in both the patient and control group to obtain glenoid version and inclination, critical shoulder angle (CSA), posterior acromial slope (PAS), lateral acromion angle, scapular offset, and the rotational alignment of the coracoacromial complex. These parameters were compared between controls, patients with osteoarthritis, and glenoid types and subtypes. RESULTS: Mean version and inclination for the control group was 6° retroversion and 8° superior inclination (both SD 4°). The mean CSA, PAS, coracoid-posterior acromion angle, posterior acromion-scapular plane angle, and fulcrum axis ratio were 30° (SD 4°), 64° (SD 8°), 116° (SD 9°), 55° (SD 7°), and 46% (SD 4%), respectively. Patients with osteoarthritis had a significant lower CSA, posterior acromion-scapular plane angle, coracoid-posterior acromion angle, and fulcrum axis ratio (27°, 50°, 111°, and 44%, all P < .001). We found a significant difference between the control group and the respective glenoid types for the following parameters: mean CSA and coracoid-posterior acromion angle for A glenoids (27°, P = .001, and 111°, P = .007); mean version, CSA, PAS, coracoid-posterior acromion angle, posterior acromion-scapular plane angle, and fulcrum axis ratio for B glenoids (11°, 27°, 71°, 111°, 49°, and 43%, all P < .001); and mean version, CSA, and posterior acromion-scapular plane angle for D glenoids (2°, P = .002, 26°, P = .003, and 48°, P = .007). DISCUSSION: There seems to be an association between prearthropathy scapular anatomy and shoulder osteoarthritis. A small lateral extension and less posterior rotation of the acromion is associated with shoulder osteoarthritis and is present in almost all types and subtypes of glenoid morphology. Furthermore, B and D glenoids are associated with, respectively, a more and less pronounced prearthropathy glenoid retroversion.

The morphological variation of acetabular defects in revision total hip arthroplasty-A statistical shape modeling approach.

Classification and evaluation of acetabular defects remain challenging and are primarily based on qualitative classification methods. That is because quantitative techniques describing variations of acetabular defects and accompanying bone loss volume are not available. This study introduces a new method based on statistical shape models (SSMs) to quantitively describe acetabular defects. This method is then applied to 87 acetabular defects to objectively describe the variations in acetabular defects typically encountered during revision total hip arthroplasty. The absolute bone loss volume, relative bone loss volume, and relative bone loss surface area with respect to the SSM-based pre-diseased anatomy were used to quantify the acetabular bone defects in different segments of the acetabular surface. The absolute bone loss volume of the average defect shape was equal to 37.0 cm3 . The first three principal modes, accounting for 62% of the total shape variation, were found to represent variations in acetabular defect morphology. The first, second, and third principal modes described, respectively, the size of the bone defects, the difference between superomedially and superolaterally migrated defects, and the degree of involvement of the posterior or anterior column. The developed SSM and the introduced approach could be used to create automated and unbiased classification methods based on quantitative data. Moreover, the proposed model and the underlying data provide the basis for a quantitative design approach where the shape and size of new acetabular implants are determined according to clinical variation present in acetabular defects.

Determination of pre-arthropathy scapular anatomy with a statistical shape model: part I-rotator cuff tear arthropathy.

HYPOTHESIS AND BACKGROUND: Rotator cuff tear arthropathy (RCTA) is a pathology characterized by a massive rotator cuff tear combined with acromiohumeral and/or glenohumeral arthritis. The severity of RCTA can be staged according to the Hamada classification. Why RCTA develops in some patients is unknown. Furthermore, in RCTA patients, distinctly different articular damage patterns can develop on the glenoid side as categorized by the Sirveaux classification (glenoid erosion). The goal of this study was to determine whether an association exists between scapular anatomy and RCTA and different severity stages of RCTA, as well as the associated glenoid erosion types. METHODS: A statistical shape model of the scapula was constructed from a data set of 110 computed tomography scans using principal component analysis. Sixty-six patients with degenerative rotator cuff pathology formed the control group. The computed tomography scan images of 89 patients with RCTA were included and grouped according to the Hamada and Sirveaux classifications. A complete 3-dimensional scapular bone model was created, and statistical shape model reconstruction was performed. Next, automated 3-dimensional measurements of glenoid version and inclination, scapular offset, the critical shoulder angle (CSA), the posterior acromial slope (PAS), and the lateral acromial angle (LAA) were performed. All measurements were then compared between controls and RCTA patients. RESULTS: The control group had a median of 7° of retroversion (variance, 16°), 8° of superior inclination (variance, 19°), and 106 mm of scapular offset (variance, 58 mm). The median CSA, PAS, and LAA were 30° (variance, 14°), 65° (variance, 60°), and 90° (variance, 17°), respectively. In terms of inclination, version, scapular offset, and the PAS, we found no statistically significant differences between the RCTA and control groups. For RCTA patients, the median CSA and median LAA were 32° (P ≤ .01) and 86° (P ≤ .01), respectively. For all investigated parameters, we did not find any significant difference between the different stages of RCTA. Patients with type E3 erosion had a different pre-arthropathy anatomy with increased retroversion (12°, P = .006), an increased CSA (40°, P ≤ .001), and a reduced LAA (79°, P ≤ .001). DISCUSSION: Our results seem to indicate that a 4° more inferiorly tilted and 2° more laterally extended acromion is associated with RCTA. RCTA patients in whom type E3 erosion develops have a distinct pre-arthropathy scapular anatomy with a more laterally extended and more inferiorly tilted acromion and a more retroverted glenoid in comparison with RCTA patients with no erosion. The pre-arthropathy scapular anatomy does not seem to differ between patients with different stages of RCTA.

Accurate reconstructions of pelvic defects and discontinuities using statistical shape models.

Treatment of large acetabular defects and discontinuities remains challenging and relies on the accurate restoration of the native anatomy of the patient. This study introduces and validates a statistical shape model for the reconstruction of acetabular discontinuities with severe bone loss through a two-sided Markov Chain Monte Carlo reconstruction method. The performance of the reconstruction algorithm was evaluated using leave-one-out cross-validation in three defect types with varying severity as well as severe defects with discontinuities. The two-sided reconstruction method was compared to a one-sided methodology. Although, reconstruction errors increased with defect size and this increase was most pronounced for pelvic discontinuities, the two-sided reconstruction method was able to reconstruct the native anatomy with higher accuracy than the one-sided reconstruction method. These findings can improve the preoperative planning and custom implant design in patients with large pelvic defects, both with and without discontinuities.

Skull fracture prediction through subject-specific finite element modelling is highly sensitive to model parameters.

Reliable computer models are needed for a better understanding of the physical mechanisms of skull fracture in accidental hits, falls, bicycle - motor vehicle & car accidents and assaults. The performance and biofidelity of these models depend on the correct anatomical representation and material description of these structures. In literature, a strain energy criterion has been proposed to predict skull fractures. However, a broad range of values for this criterion has been reported. This study investigates if the impactor orientation, scalp thickness and material model of the skull could provide us with insight in the influencing factors of this criterion. 18 skull fracture experiments previously performed in our research group were reproduced in finite element simulations. Subject-specific skull geometries were derived from medical images and used to create high-quality finite element meshes. Based on local Hounsfield units, a subject-specific isotropic material model was assigned. The subject-specific models were able to predict fractures who matched visually with the corresponding experimental fracture patterns and provided detailed fracture patterns. The sensitivity study showed that small variations in impactor positioning as well as variations of the local geometry (frontal-temporal-occipital) strongly influenced the skull strain energy. Subject-specific modelling leads to a more accurate prediction of the force-displacement curve. The average error of the peak fracture force for all the 18 cases is 0.4190 for the subject-specific and 0.4538 for the homogeneous material model, for the displacement; 0.3368 versus 0.3844. But it should be carefully interpreted as small variations in the computational model significantly influence the outcome.