Pre-operative 3D glenoid planning improves component placement in terms of version, inclination, offset and orientation. Version and inclination measurements require the position of the inferior angle. As a consequence, current planning tools require a 3D model of the full scapula to accurately determine the glenoid parameters. Statistical shape models (SSMs) can be used to reconstruct the missing anatomy of bones. Therefore, the objective of this study is to develop and validate an SSM for the reconstruction of the inferior scapula, hereby reducing the irradiation exposure for patients. The training dataset for the statistical shape consisted of 110 CT images from patients without observable scapulae pathologies as judged by an experienced shoulder surgeon. 3D scapulae models were constructed from the segmented images. An open-source non-rigid B-spline-based registration algorithm was used to obtain point-to-point correspondences between the models. A statistical shape model was then constructed from the dataset using principal component analysis. Leave-one-out cross-validation was performed to evaluate the accuracy of the predicted glenoid parameters from virtual partial scans. Five types of virtual partial scans were created on each of the training set models, where an increasing amount of scapular body was removed to mimic a partial CT scan. The statistical shape model was reconstructed using the leave-one-out method, so the corresponding training set model is no longer incorporated in the shape model. Reconstruction was performed using a Monte Carlo Markov chain algorithm, random walk proposals included both shape and pose parameters, the closest fitting proposal was selected for the virtual reconstruction. Automatic 3D measurements were performed on both the training and reconstructed 3D models, including glenoid version, inclination, glenoid centre point position and glenoid offset. In terms of inclination and version we found a mean absolute difference between the complete model and the different virtual partial scan models of 0.5° (SD 0.4°). The maximum difference between models was 3° for inclination and 2° for version. For offset and centre point position the mean absolute difference was 0 mm with an absolute maximum of 1 mm. The magnitude of the mean and maximum differences for all anatomic measurements between the partial scan and complete models is smaller than the current surgical accuracy. Considering these findings, we believe a SSM based reconstruction technique can be used to accurately reconstruct the glenoid parameters from partial CT scans.
A large number of total hip arthroplasties (THA) are performed each year, of which 60 % use cementless femoral fixation. This means that the implant is press-fitted in the bone by hammer blows. The initial fixation is one of the most important factors for a long lasting fixation [Gheduzzi 2007]. It is not easy to obtain the point of optimal initial fixation, because excessively press-fitting the implant by the hammer blows can cause peak stresses resulting in femoral fracture. In order to reduce these peak stresses during reaming, IMT Integral Medizintechnik (Luzern, Switzerland) designed the Woodpecker, a pneumatic reaming device using a vibrating tool. This study explores the feasibility of using this Woodpecker for implant insertion and detection of optimal fixation by analyzing the vibrational response of the implant and Woodpecker. The press-fit of the implant is quantified by measuring the strain in the cortical bone surrounding the implant. An in vitro study is presented. Two replica femur models (Sawbones Europe AB, Malmo Sweden) were used in this study. One of the femur models was instrumented with three rectangular strain gauge rosettes (Micro-Measurements, Raleigh, USA). The rosettes were placed medially, posteriorly and anteriorly on the proximal femur. Five paired implant insertions were performed on both bone models, alternating between standard hammer blow insertions and using the Woodpecker. The vibrational response was measured during the insertion process, at the implant and Woodpecker side using two shock accelerometers (PCB Piezotronics, Depew, NY, USA). The endpoint of insertion was defined as the point when the static strain stopped increasing. Significant trends were observed in the bandpower feature that was calculated from the vibrational spectrum at the implant side during the Woodpecker insertion. The bandpower is defined as the percentage power of the spectrum in the band 0–1000 Hz. Peak stress values calculated from the strain measurement during the insertion showed to be significantly (p < 0.05) lower at two locations using the Woodpecker compared to the hammer blows at the same level of static strain. However, the final static strain at the endpoint of insertion was approximately a factor two lower using the Woodpecker compared to the hammer. A decreasing trend was observed in the bandpower feature, followed by a stagnation. This point of stagnation was correlated with the stagnation of the periprosthetic stress in the bone measured by the strain gages. The behavior of this bandpower feature shows the possibility of using vibrational measurements during insertion to assess the endpoint of insertion. However it needs to be taken into account that it was not possible to reach the same level of static strain using the Woodpecker as with the hammer insertion. This could mean that either extra hammer blows or a more powerful pneumatic device could be needed for proper implant insertion.