This retrospective study used de-identified data from archived, clinically indicated CT examinations. The study was conducted in accordance with the Declaration of Helsinki and was approved by the Ethics Committee of Beijing Chaoyang Hospital, Capital Medical University (approval no. 2021-Ke-441). The requirement for informed consent was formally waived by the ethics committee. The tools used in this protocol are listed in the Table of Materials.
1. Patient selection
This retrospective study included 100 unique patients who underwent clinically indicated CT examinations at the study hospital between 2016 and 2023. One shoulder was included per patient, yielding 100 shoulder CT datasets and 100 corresponding scapular models. When both shoulders fulfilled all eligibility criteria, the right shoulder was selected. The cohort comprised 50 male and 50 female patients aged 60–87 years (mean, 69.97 ± 7.91 years).
Patients were eligible if they were 60–90 years old and an archived clinical CT examination included the complete scapula. Selection was not based on an orthopedic shoulder complaint. Exclusion criteria were radiographically evident osteoarthritis, defined by at least moderate joint-space narrowing, osteophyte formation, or subchondral sclerosis; rheumatoid arthritis; previous shoulder infection; previous humeral head fracture or dislocation; or congenital shoulder deformity. These criteria were chosen to isolate relatively preserved glenoid geometry rather than to represent the pathological glenoids typically treated with RTSA.
2. CT acquisition and 3D reconstruction
Patients were positioned according to the routine shoulder CT protocol, and the entire scapula was included in the scan field. Images were acquired with a 64-row spiral CT scanner. The series was reconstructed at a slice thickness of 1.25 mm and an interval of 0.625 mm.
The axial series was reviewed to confirm inclusion of the glenoid, scapular spine, medial border, superior and inferior angles, and the complete scapular body. A dataset was excluded if any of these structures was truncated or affected by motion artifact that prevented segmentation. The selected image series was exported from the picture archiving and communication system in Digital Imaging and Communications in Medicine (DICOM) format. Native orientation and voxel spacing were retained during export. The DICOM series was imported into three-dimensional medical image-processing software. An initial bone mask was created using threshold-based segmentation. Because CT attenuation and image quality varied among examinations, the threshold was adjusted case by case until the cortical outline and cancellous portion of the entire scapula were included without incorporating adjacent soft tissue.
Region-growing segmentation was applied from a seed point within the scapula to retain the connected scapular mask. The mask was inspected in the axial, coronal, and sagittal planes. When the humerus, clavicle, ribs, or residual soft tissue remained connected, the mask was edited to separate these structures slice by slice, and internal cancellous-bone cavities were filled. The complete image stack was reviewed, and discontinuities or leakage were corrected manually, with particular attention to the glenoid rim, articular surface, scapular neck, scapular spine, and medial border. A triangular surface mesh was then generated.
Smoothing was applied only to small stair-step irregularities, and mesh reduction was used to decrease excessive mesh density while preserving the glenoid rim, articular surface, scapular spine, and medial border. The original reconstruction workflow did not retain fixed numerical smoothing or reduction settings; therefore, these tools were applied conservatively, and the anatomical landmarks were confirmed in all three orthogonal planes before the model was exported in stereolithography (STL) format. The STL mesh was imported into three-dimensional surface-processing software. Mesh defects identified by the software were repaired. Local relaxation and spike-removal functions were applied to correct residual sharp projections or irregular triangles, and mesh repair was repeated until no defect that prevented surface construction remained.
During exact surfacing, automatic surface generation was selected, an organic geometry type was used, the target number of surface patches was set to 500, and adaptive fitting tolerance was applied. The resulting surface was inspected for continuity and preservation of the glenoid and scapular-spine landmarks, and the solid model was exported as a STEP AP203 file. The STEP AP203 file was imported into three-dimensional computer-aided design software. The model was confirmed to be a closed solid body, the original scale in millimeters was retained, and the model was saved for trajectory construction and virtual peg simulation.
3. Definition of the coordinate system and four placement strategies
In the three-dimensional computer-aided design software, the glenoid was displayed in an en face lateral view. The superior and inferior glenoid poles were marked, and line AB was created between them.
The anterior and posterior glenoid rims were marked at the widest anteroposterior portion of the articular surface. Line CD was created between these points and constrained perpendicular to AB. Lines AB and CD were used to define the glenoid reference plane. The intersection of AB and CD was created and projected onto the glenoid articular surface to define the glenoid center O. The projected point was verified in the en face and orthogonal views.
The standard centerline (SCL) was created as a 3D reference line passing through O and perpendicular to the glenoid reference plane, directed medially into the scapula. The longitudinal axis of the scapular spine toward its medial insertion at the medial scapular border was identified. The alternative centerline (ACL) was created as a 3D reference line passing through O and directed medially along this scapular-spine axis.
On line AB, a point 12 mm superior to the inferior glenoid rim was identified and projected onto the articular surface to create the second entry point. Through this point, one trajectory parallel to the SCL and another parallel to the ACL were created. At each entry point, a local coordinate system was created. In the glenoid plane, +x was defined as posterior and +y as superior, and +z was defined along the initial peg trajectory, directed medially into the scapula.
A straight cylindrical central peg with a diameter of 7.5 mm was created and assembled with the scapular model (Figure 1). The four strategies (Figure 2) were defined as follows: method a, O/SCL; method b, O/ACL; method c, 12 mm superior entry/SCL-parallel; and method d, 12 mm superior entry/ACL-parallel.

Figure 1: Three-dimensional reconstructed scapular model with a simulated central peg. A cylindrical 7.5 mm-diameter central peg is shown within the reconstructed scapula. The red segment represents the portion of the peg contained within the scapular bone; the gray segment shows the peg axis outside the entry surface. Please click here to view a larger version of this figure.

Figure 2: Four simulated central peg placement strategies. (A) Method a: Glenoid-center entry along the standard centerline (SCL). (B) Method b: Glenoid-center entry along the alternative centerline (ACL). (C) Method c: Entry 12 mm superior to the inferior glenoid rim along a trajectory parallel to the SCL. (D) Method d: Entry 12 mm superior to the inferior glenoid rim along a trajectory parallel to the ACL. The red cylinder represents the 7.5 mm-diameter central peg, and the gray extension indicates its trajectory. Please click here to view a larger version of this figure.
4. Measurement of maximum contained peg depth and cortical penetration
For each scapular model and placement method, an assembly constraint was used to place the center of the peg end face at the selected entry point. The peg axis was aligned with the SCL, ACL, or the corresponding parallel trajectory.
The entry point and trajectory were confirmed in the lateral en face, axial, and superior or inferior views. A section view was used when necessary to inspect the peg within the scapular vault. The peg was translated medially along its axis while the entry-point and trajectory constraints were maintained. Interference analysis was used to inspect the peg–bone relationship after each adjustment. Translation was stopped at the greatest depth at which the complete cylindrical peg remained within or tangent to the external cortical boundary and no portion projected beyond the cortex.
The distance from the glenoid entry point to the limiting distal end of the peg was recorded as the maximum contained peg depth in millimeters. The value was recorded before the assembly configuration was changed. The assembly was repeated with nominal peg lengths of 15, 20, and 25 mm. A configuration was classified as contained when the entire peg was inside or tangent to the cortical boundary and as cortical penetration when any part of the peg extended beyond the external cortex.
5. Angular-deviation simulations
The entry point, initial trajectory, and local coordinate system established for the selected method were retained (Figure 3). The scapula was displayed from the lateral side of the glenoid while viewed medially along the +z trajectory, and the same working view was used for Figure 3 and Figure 4.
With the entry point fixed, a deviated reference trajectory was created at a polar angle α of 5° from the initial +z trajectory. The construction was repeated for α = 10° and 15°. The entry point remained unchanged when the polar angle was altered. The azimuth θ was defined in the local x–y plane and increased clockwise in the lateral-to-medial working view. The directional convention was defined as θ = 0° superiorly (+y), 90° posteriorly (+x), 180° inferiorly (−y), and 270° anteriorly (−x).
For each polar angle α, the deviation direction was rotated clockwise about the initial +z trajectory from θ = 0° to 345° in 15° increments, generating 24 azimuthal directions. The viewing direction remained lateral-to-medial throughout the rotation. At every α–θ combination, the 7.5 mm peg was aligned with the deviated trajectory, and the maximum contained peg depth was determined using the procedure described above. The 24-direction sequence was repeated for all three polar deviations, all four placement methods, and all 100 scapular models. The same clockwise azimuth convention was applied to the directional results and Figure 4.

Figure 3: Local coordinate system and working view for angular-deviation analysis. The scapula is displayed from the lateral side of the glenoid while looking medially along the initial +z peg trajectory. In the plane parallel to the glenoid surface, +x points posteriorly and +y superiorly; −x and −y indicate anterior and inferior. Azimuth increases clockwise in this lateral-to-medial view: 0° is superior (+y), 90° posterior (+x), 180° inferior (−y), and 270° anterior (−x). Figures 3 and 4 use the same viewing direction and azimuth convention. Please click here to view a larger version of this figure.

Figure 4: Directional maximum contained peg depth under simulated angular deviations. Panels A–D correspond to methods a–d. The plots use the same lateral-to-medial view as Figure 3, and azimuth increases clockwise from 0° superior through 90° posterior, 180° inferior, and 270° anterior in 15° increments. The radial axis is the mean maximum contained peg depth (mm) across 100 models. Solid yellow, blue, and green curves denote imposed polar deviations of 5°, 10°, and 15°, respectively; deviation describes trajectory offset and is not a statistical error estimate. Points are directional means. No SD, confidence interval, or error bars are displayed. Please click here to view a larger version of this figure.
6. Repeated measurements and resolution of disagreements
Two observers independently repeated all continuous depth measurements in two sessions. The final continuous value for each model and configuration was the mean of the four accepted readings (two observers × two sessions). If the review identified an incorrect landmark, trajectory, or cortical boundary, a senior radiologist confirmed the anatomical setup, and the affected observer repeated the measurement. The corrected reading replaced the invalid reading before calculation of the four-reading mean. The senior radiologist did not supply a fifth depth value.
Each observer/session also classified cortical penetration for each nominal peg length. When the four binary classifications were discordant, the senior radiologist re-examined the 3D model and assigned the final binary status used in the analysis. Inter-rater reliability for continuous depth was evaluated from each observer’s two-session mean using a two-way random-effects, absolute-agreement intraclass correlation coefficient (ICC). Intra-rater reliability was evaluated by comparing sessions using a two-way mixed-effects, absolute-agreement ICC.
Figure 5 summarizes patient selection, CT acquisition, 3D reconstruction, coordinate definition, the four placement strategies, nominal peg-length simulation, angular-deviation testing, outcome determination, and paired statistical analysis.

Figure 5: Overall 3D CT-based simulation workflow. The schematic summarizes selection of 100 unique patients (one shoulder and one scapular model per patient), CT acquisition, DICOM-based 3D reconstruction, coordinate and strategy definition, simulation of 15, 20, and 25 mm peg lengths, 5°, 10°, and 15° angular deviations at 15° azimuthal intervals, determination of maximum contained peg depth and adjudicated cortical penetration, and paired statistical analysis. Please click here to view a larger version of this figure.
7. Data analysis
For each nominal peg length and placement method, the adjudicated cortical-penetration outcome was summarized as n (%). Maximum contained peg depth was summarized as mean ± standard deviation (SD) and range. For angular-deviation analyses, the mean maximum contained peg depth across the 24 azimuths was calculated for each method and polar deviation, and these directional means were used to construct polar plots.
Data were managed, and descriptive summaries were calculated using spreadsheet software. Inferential analyses were performed using statistical analysis software. Because all four placement methods were evaluated in every scapular model, initial maximum contained peg depths were compared using one-factor repeated-measures analysis of variance (ANOVA). After a significant omnibus result, paired post hoc comparisons were conducted, and the family-wise error rate across the six method pairs was controlled with the Holm procedure. At each nominal peg length, the four paired binary penetration outcomes were compared using Cochran’s Q test. Pairwise exact McNemar tests with Holm adjustment were performed only if the omnibus test was significant. Two-sided tests were used with α = 0.05.