A computational tool for automatic selection of total knee replacement implant size using X-ray images

Purpose: The aim of this study was to outline a fully automatic tool capable of reliably predicting the most suitable total knee replacement implant sizes for patients, using bi-planar X-ray images. By eliminating the need for manual templating or guiding software tools via the adoption of convolutional neural networks, time and resource requirements for pre-operative assessment and surgery could be reduced, the risk of human error minimized, and patients could see improved outcomes.

Methods: The tool utilizes a machine learning-based 2D—3D pipeline to generate accurate predictions of subjects’ distal femur and proximal tibia bones from X-ray images. It then virtually fits different implant models and sizes to the 3D predictions, calculates the implant to bone root-mean-squared error and maximum over/under hang for each, and advises the best option for the patient. The tool was tested on 78, predominantly White subjects (45 female/33 male), using generic femur component and tibia plate designs scaled to sizes obtained for five commercially available products. The predictions were then compared to the ground truth best options, determined using subjects’ MRI data.

Results: The tool achieved average femur component size prediction accuracies across the five implant models of 77.95% in terms of global fit (root-mean-squared error), and 71.79% for minimizing over/underhang. These increased to 99.74% and 99.49% with ±1 size permitted. For tibia plates, the average prediction accuracies were 80.51% and 72.82% respectively. These increased to 99.74% and 98.98% for ±1 size. Better prediction accuracies were obtained for implant models with fewer size options, however such models more frequently resulted in a poor fit.

Conclusion: A fully automatic tool was developed and found to enable higher prediction accuracies than generally reported for manual templating techniques, as well as similar computational methods.

1 Introduction

Prior to total knee replacement (TKR) procedures, surgeons often manually evaluate the size and morphology of patient femur and tibia bones using X-rays taken during preoperative assessment (Tanzer and Makhdom, 2016). This information is used to estimate the most appropriately sized TKR implant components that will be required to achieve both a good global fit, and minimize local areas of over/underhang (OUH), typically from between 5 – 8 sizes (Hitt et al., 2003; Wernecke et al., 2012). Poor implant size selection can result in increased rates of complications, revisions, and patient pain post-surgery (Culler et al., 2017; Schroeder and Martin, 2019; Buller et al., 2018). The presence of regions of OUH ≥3 mm in particular is highlighted as being clinically significant and directly linked with cases of increased soft-tissue irritation, bleeding, osteolysis, laxity in flexion, subsidence, and instability (Schroeder and Martin, 2019; Dai et al., 2014; Shao et al., 2020). Despite the importance of achieving a good fit, the accuracy of manual X-ray assessment and templating, for both femur components and tibia plates, can be poor. Hernández-Vaquero et al. (2019) summarized the reported accuracies from 10 different studies and found mean selection accuracies of 59.2% for femur components and 60.7% for tibia plates. With ±1 size permitted, mean scores of 97.4% and 96.4% were recorded respectively. Consequently, surgeons often opt to implant a different size than estimated during preoperative assessment (Sheth et al., 2017). This could lead to a higher chance of human error, as well as longer procedures.

To attempt to improve the accuracy and reliability of pre-operative size selection, other techniques besides manual X-ray assessment and templating have been explored. Using computerized tomography (CT) imaging was investigated by Vaishya et al. (2018) and Kobayashi et al. (2012). The former reported accuracies of 66% for femur components and 72% for tibia plates, whilst the latter reported an accuracy of just 59% for both. Using demographic factors such as patient height and shoe size as a means of predicting implant size was investigated by Trainor et al. (2018). They reported accuracies of 56/58% for femur components and 56/63% for tibia plates respectively. Sershon et al. (2017) used a combination of demographic variables including height, weight, and sex to build a multivariate linear regression model that achieved ±1 size prediction accuracies of 71–92% and 81–97% for a range of femur component and tibia plate models. Optimum size prediction accuracies were not reported. Limited computational tools capable of predicting TKR implant sizes from X-rays were identified. Zheng et al. (2018) developed an X-ray based tool for pre-operative knee prosthesis planning named “3X.” The authors reported size selection accuracies of 78% for femur components and 70% for tibia plates, based on a study featuring 23 subjects. Massé and Ghate (2021) more recently developed a similar tool named “X-Atlas”, primarily aimed at creating patient-matched cutting guides from X-ray images. The study involved 45 subjects and utilized the Zimmer Biomet “Persona” knee implant. Size selection accuracies of 53.3% for femur components and 57.8% for tibia plates were reported. Both these tools feature “semi-automatic” workflows that require users to manually identify landmarks on the inputted X-ray images and guide “live-wire algorithms” to extract the required bone contours.

In recent years, prediction and classification tools have been developed in other biomedical imaging applications using machine learning techniques to enable automatic workflows. U-Net convolutional neural networks (CNNs) in particular are growing in popularity as highly accurate classification and segmentation models can be trained with a relatively low number of reference images (Ronneberger et al., 2015). Example applications that have adopted this technology include in the diagnosis of COVID-19 (Mohammed et al., 2021) (Narin et al., 2021) and for detecting tuberculosis (Liu et al., 2018). In these studies, CNNs were trained to recognize features and patterns in X-ray images and extract the necessary information for diagnosis, all without requiring any guidance from the user. Furthermore, Cernazanu-Glavan and Holban (2013) demonstrated the benefits of usings CNNs for X-ray-based bone segmentation over alternative methods such as artificial neural networks, principal component analysis and fuzzy clustering. The authors highlighted the superior accuracy of CNNs and that all other methods investigated required some level of human intervention. In the case of computer assisted TKR implant size selection, the tools outlined in previous studies all featured workflows requiring trained users’ feedback (Zheng et al., 2018; Massé and Ghate, 2021). Consequently, the authors do not believe that these “semi-automatic” solutions adequately address the issues preventing widespread adoption of computer assisted sizing tools in TKR. The aim of this study was therefore to develop a proof-of-concept, “fully” automatic alternative, facilitated by employing CNNs in the workflow. Such a tool could lower the possibility of human error, whilst reducing the burden on hospitals relating to training, completing manual X-ray templating, and driving software packages. Necessitating no more imaging requirements over conventional manual X-ray based methods would help avoid additional costs, radiation exposure, and resource/time. Lastly, an automatic, X-ray based tool could also be of benefit in less developed countries where implant variety and access to 3D imaging would likely be more limited. The workflow of the developed size prediction tool, from inputted X-rays to size selection, is outlined. Its performance for the dimensions of five manufacturer models is then analyzed via various fit metrics and compared to the manual and semi-automatic solutions detailed above.

2 Materials and methods

2.1 Datasets

Two datasets were utilized to train and test the size prediction tool. Anterior-posterior (AP) and lateral X-ray images, as well as high resolution (3T) magnetic resonance imaging (MRI) scans, of knee joints were obtained from the Osteoarthritis Initiative (OAI) (Nevitt et al., 2006). This facilitated data to train image segmentation models within the tool and provided the necessary pairing of X-rays with 3D ground truth data to test the accuracy of the full workflow. 3D Slicer was used to segment the OAI MRI data and generate 3D bone models (Fedorov et al., 2012). Scaling of the X-ray data retrieved from the OAI dataset was completed to address geometric magnification effects and ensure consistency with the MRI data in lieu of calibration artifacts being included. CT scans with pre-segmented 3D femur and tibia models were sourced from the Korea Institute of Science and Technology Information (KISTI) (Lee and Lee, 2010). This data was used for training aspects of the size selection tool as detailed in Section 2.3. Both datasets consisted of male and female subjects with a broad range of ages. The KISTI dataset consisted purely of Asian Korean subjects, while the OAI contained White, Black and Asian Americans. No subjects used in training or testing had above mild levels of arthritis (above grade 2 out of 4) according to Kellgren and Lawrence (1956).

2.2 Test subjects

Due to the retrospective nature of the study, it was not possible to control the X-ray imaging process. Subjects’ X-rays were instead visually screened to ensure appropriate anatomical alignment and clarity in both AP and lateral projections. Figure 1 illustrates correctly captured X-rays with the scan taken straight on with both the medial and lateral condyles pointing forward in the AP projection, and side on with the posterior surfaces of the femur condyles aligned in the lateral view. 78 subjects were found to be suitable and were taken from the OAI dataset to test the size prediction tool. The test subjects ranged from 46 to 79 years old, and their demographics are summarized in Table 1. KISTI subjects were used purely for training the tool as the dataset only contained CT data without accompanying X-ray images.

FIGURE 1

FIGURE 1. Top level workflow of the automatic 2D—3D reconstruction pipeline (top row–femur, bottom row–tibia): (A) CNN segmentation of input X-rays, (B) extracted bone contours, (C) contours aligned in 3D space, (D) PDMs used to transform contours into sparse point clouds which are then fitted to base shapes of SSMs, (E) 3D model predictions.

TABLE 1

TABLE 1. Test subject demographical information.

2.3 2D—3D pipeline

Like in Zheng et al. (2018) and Massé and Ghate (2021), the size prediction tool initially utilizes inputted AP and lateral X-ray images to generate 3D estimations of patients’ femur and tibia bones before component size predictions can be made. To achieve this automatically, the 2D—3D pipeline developed in Burge et al. (2022a), Burge et al. (2022b) was utilized. The key aspects of the workflow, built using Python 3, are summarized below and illustrated in Figure 1.

The first step of the 2D—3D process utilizes CNNs to isolate the femur and tibia bones from the surrounding soft tissue in the inputted X-ray images, as shown in Figure 1A. Four U-Net CNN image segmentation models (one for each X-ray projection of each bone) are used in the tool. The models were built using TensorFlow and a U-Net architecture, based on Ronneberger et al. (2015), was adopted (illustrated in Figure 2). To better capture the bone profiles occupying a large proportion of the X-ray images, an additional filter resolution level (16 x 16 pixels) was incorporated. The number of filter channels used in each layer and learning rate were adjusted to achieve the best validation results. Batch size was set at 10 with 100 epochs and model weights were saved after each iteration. Once training had finished, the weights with the lowest validation loss were selected for each model to minimize overfitting. 176 X-ray image/mask pairs were developed from the OAI dataset to train each model. 20 of these (∼10%) were reserved for validation (separate to those used to test the full tool). It is noted that to build well generalized CNN models, capable of reliably segmenting any inputted X-ray, thousands of varied training images would likely be required (Liu et al., 2018). Due to the proof-of-concept nature of this study, and lack of suitable training data, it was deemed acceptable to mitigate this by adjusting the contrast/brightness settings of the inputted subject X-ray images to improve compatibility with the CNN models when testing the size prediction tool.

FIGURE 2

FIGURE 2. Architecture of U-Net CNN X-ray segmentation models from inputted X-ray image (left) to predicted mask (right). Each box corresponds to a multi-channel feature map with those shaded grey indicating where channels have been copied from previous layers. The number of channels is denoted on top of each box with the size provided at the lower left (256² = 256 × 256 pixels). Arrows denote different operations such as convolutions.

After segmentation, a Cany edge detector is applied to extract the four bone contours (Figure 1B) which are then aligned in 3D space as shown in Figure 1C. Point depth models (PDMs) are applied to estimate the third dimension of the points along each contour and create a sparse 3D point cloud of reference coordinates (Figure 1D). The femur and tibia reference point clouds are independently fitted to two statistical shape models (SSMs) (one for each bone) via an iterative point method. The SSMs then morph base femur and tibia shape models to fit the reference points as closely as possible. For the PDMs and SSMs, 3D bone models (20 and 100 respectively) were used from the KISTI dataset for training.

The outputs from the 2D—3D pipeline are predictions of subjects’ distal femur and proximal tibia anatomies in the form of 3D surface mesh models (Figure 1E). These can subsequently be used by the tool to approximate the most appropriately sized implant components.

2.4 Implant designs and sizes

Generic base models of TKR femur components and tibia plates were designed based on widely used commercial products (Figure 3A). The base models were scaled (Figure 3B) to the sizes listed for five manufacturers’ models as reported in size charts within their respective surgical technique manuals (Tables 2 and 3). This method was used as access to the official geometry of all the manufacturer models and sizes utilized within the study was not possible. Nor was acquiring physical samples for reverse engineering. Utilizing generic implant shapes also allowed for the test models to be easily controlled/edited for the purposes of the study.

FIGURE 3

FIGURE 3. Method for evaluating fit of implant model on subject ground truth anatomy (top row–femur, bottom row–tibia): (A) Generic implant component designs, (B) edges of scaled components, (C) scaled components fitted to subject 3D model predictions, (D) scaled components fitted to subject ground truth anatomy. Dashed box in (A) shows detail not included in fit analysis. Red squares in (C) and (D) show the location of the maximum OUH.

TABLE 2

TABLE 2. Femur component sizing chart, detailing the manufacturer, model, size identifier, and ML/AP dimensions (obtained from respective surgical technique manuals).

TABLE 3

TABLE 3. Tibia plate sizing chart, detailing the manufacturer, model, size identifier, and ML/AP dimensions (obtained from respective surgical technique manuals).

The height of the 3D femur components was kept proportional by scaling the Z dimension in line with the AP and ML dimensions. The design of the base femur component was subtly adjusted for each model to reflect differences in transepicondylar (ML) width between the designs and ensure consistency with the dimensioning used in each of the surgical manuals. The design of the base tibia plate was kept constant and only the 2D profile was used like in prior studies (Shao et al., 2020; Clary et al., 2014). A 2D analysis was sufficient for the tibia plate because the components only interface with the resected bone on a singular 2D face.

2.5 Size prediction and accuracy calculation

To predict which femur component and tibia plate sizes are the most suited for a subject, the tool systematically fits each implant model and size to the 3D estimates of the subject’s anatomy using an iterative closest point method (Figure 3C). Two fit metrics were calculated after each size fitting including the global root-mean-squared error (RMSE), and local maximum OUH. The RMSE calculation was performed between the surface of the positioned component and the subject anatomy as described by Eq. 1:

$\mathit{R}\mathit{M}\mathit{S}\mathit{e}\mathit{r}\mathit{r}\mathit{o}\mathit{r}=\sqrt{\frac{{\sum }_{\mathit{i}=1}^{\mathit{N}}{\left({\mathit{x}}_{\mathit{i}}-{\widehat{\mathit{x}}}_{\mathit{i}}\right)}^2}{\mathit{N}}}(1)$

where N is the number of points and $x_i-{\widehat{x}}_i$ is the Euclidian distance between each point of the component surface and the bone. The maximum OUH was reported as the Hausdorff distance (h) anywhere between the edges of the component (C) to the edges of the bone (B), as described by Eq. 2 and shown in Figure 3. The distance (d) between each point along the component edges (c) and the bone (b) was calculated as the Euclidian distance.

$\mathit{h}\left(\mathit{C},\mathit{B}\right)=\mathit{max}{}_{\mathbf{c}\in \mathbf{C}}\left\{{\mathit{min}}_{\mathbf{b}\in \mathbf{B}}\left\{\mathit{d}\left(\mathit{c},\mathit{b}\right)\right\}\right\}(2)$

For the femur, the fit analysis was completed in 3D assuming that the bone would be resected from the lateral view to the same standard dimensions as the mating faces on the femur component being fitted. For the tibia, the RMSE and maximum OUH were calculated using the 2D profile of the tibia plate and a cross-section profile taken at the intended resection point on the tibia bone. In practice, the necessary resection depth is determined by the surgeon, however, for the purposes of this study, the resection plane was located 2 mm below the height of the widest point of the tibia medial condyle, parallel with the surface of the medial plateau. Rotation of the tibia plate was limited to about its central axis and no flexion or extension was permitted. This approach consistently facilitated the largest possible surface area for stability, minimized bone loss (Schnurr et al., 2011), and enabled a continuous, flat resection plane to be created through the bone. Once RMSE and maximum OUH were recorded for each component size, the options for each model that achieved the lowest errors were outputted as the tool predictions.

2.6 Results analysis

After predictions were made, all sizes for each implant model were positioned on 3D ground truth models via an iterative closest point method and the resultant RMSE and maximum OUH values were calculated (Figure 3D). Ground truth models were created for each test subject by manually segmenting their MRI data as described previously. The “ground truth best” size for the individual (in terms of each fit metric, implant component, and model) was recorded as the option which resulted in the lowest calculated error when fitted to the subject’s 3D ground truth model. If the tool predicted the same as the ground truth best, the prediction was deemed to be correct. Spearman’s correlation coefficients were calculated and used to evaluate the impact on performance due to continuous variables such as subject age.

3 Results

3.1 2D—3D model generation

The accuracy of the 2D—3D process was evaluated by Burge et al. (2022b) by comparing the tool’s bone surface estimations with subjects’ 3D ground truth models. A mean RMSE of 1.09 mm (SD 0.18 mm) for the femur and 0.98 mm (SD 0.15 mm) for the tibia across the 78 subjects were reported.

3.2 Size selection accuracy

The accuracy of the femur component size selection for each model is reported in Table 4 and in Table 5 for tibia plates. Tables 4 and 5 demonstrate that high levels of accuracy for the size selection tool across the five implant model sizes for both implant components were obtained in terms of RMSE and maximum OUH. Minimal difference in prediction accuracy was observed across the five tibia plate models for males and females. For femur components however, males achieved 87.27% for both RMSE and maximum OUH, whilst accuracies of 71.11% and 60.44% were recorded for females respectively. Calculating Spearman’s correlation coefficients showed no strong correlations between prediction accuracies and subject age for both femur components (RMSE = −0.03, maximum OUH = 0.02) and tibia plates (RMSE = 0.21, maximum OUH = 0.25) were present.

TABLE 4

TABLE 4. Femur component size prediction accuracies for each manufacturer model sizings. Size prediction accuracy is shown in terms of RMSE and maximum OUH, as well as allowing for ± one model size.

TABLE 5

TABLE 5. Tibia plate size prediction accuracies for each manufacturer model sizings. Size prediction accuracy is shown in terms of RMSE and maximum OUH, as well as allowing for ± one model size.

3.3 Implant fit

For the sizes deemed to be the best by the size prediction tool for each model and each of the test subjects, the resulting fit on the ground truth models were recorded for both metrics. These results are shown in Table 6 for femur components and Table 7 for tibia plates, alongside the best possible results. The tables show that the results obtained using the predictions made by the tool for both component types were on average only marginally inferior to the best possible outcomes in terms of both RMSE and maximum OUH.

TABLE 6

TABLE 6. Femur component mean RMSE and mean maximum OUH results for ground truth (GT) best sizes and predictions, split by manufacturer model sizings.

TABLE 7

TABLE 7. Tibia plate mean RMSE and mean maximum OUH results for ground truth (GT) best sizes and predictions, split by manufacturer model sizings.

Figures 4, 5 show the proportion of subjects achieving maximum OUH <3 mm for the various femur component and tibia plate model sizes. The figures demonstrate that the level of subjects attaining below the threshold was on average 10% less for femur components compared to tibia plates and the levels varied significantly between the different manufacturer model sizes. The proportion of subjects seeing maximum OUH of ≥3 mm was on average 12% higher than the ground truth best when the tool’s size predictions were used.

FIGURE 4

FIGURE 4. Boxplot illustrating proportion of test subjects achieving maximum OUH < 3 mm for various femur component model sizings. Both ground truth (GT) best possible size and tool prediction results shown.

FIGURE 5

FIGURE 5. Boxplot illustrating proportion of test subjects achieving maximum OUH < 3 mm for various tibia plate model sizings. Both ground truth (GT) best possible size and tool prediction results shown.

4 Discussion

The most important outcome from this study is that a fully automatic, X-ray based TKR size prediction tool was developed and shown to be reliable for a range of implant model sizes and test subjects. For both femur components and tibia plates, the tool more accurately predicted the ground truth best size in terms of RMSE compared to maximum OUH. Nevertheless, the accuracies for both metrics were consistently high with the ground truth best size predicted on average 78% of the time (in terms of RMSE) and 72% of the time (in terms of maximum OUH) for femur components. For tibia plates the prediction accuracies were 81% in terms of RMSE and 73% in terms of maximum OUH. These increased to 99–100% for ±1 size for both metrics and component types. The prediction accuracy of the tool was not found to be sensitive to subject age, nor was it sensitive to sex for predicting tibia plate size. However, 16–27% higher accuracies were achieved for male subject femur component size selection. This was likely due to the large White American male dimensions in the test population– often requiring the upper limits of the femur component size ranges or beyond. Females on the other hand utilized a broader range of smaller sizes which made predicting the correct option more challenging. Using a more balanced group of subjects, featuring ethnicities typically requiring smaller implant sizes such as Asian Chinese (Li et al., 2019), would likely reduce this effect and result in similar prediction accuracies between sexes.

Another key finding was the proportion of subjects obtaining clinically significant levels of OUH, even when the best possible sizes were used, was above 30% for all implant models evaluated. This aligns with results reported by previous studies (Wernecke et al., 2012; Mahoney and Kinsey, 2010) where the performance of non-customised implants, particularly with a limited number of size options, was shown to be poor. In this study, the model with the fewest number of sizes (the Zimmer Biomet NexGen), performed the worst in terms of the resulting mean RMSE, maximum OUH, and proportion of subjects achieving maximum OUH <3mm, whilst models with more sizes generally performed better. Nevertheless, for the tool’s size prediction accuracies, the opposite was found. The Zimmer Biomet NexGen (5/6 sizes) obtained the best results, whilst the Stryker Scorpio (9 sizes) achieved the lowest accuracy for femur components (both in terms of RMSE). Comparing the results for the Smith & Nephew Legion and the Maxx Orthopedics Freedom tibia plate sizes showed considerable differences in prediction accuracies and resulting fits, despite an equal number of sizes. It is again emphasized that these results do not reflect the true performance of the named manufacturer models as generic base models were used in the analysis. However, the importance of selecting the optimum sizing dimensions, not just enough options, is highlighted. For models with many similar sizes, it was often observed that when the tool failed to predict the best possible size, the optimal component dimensions would have been close to the center of two options. As a result, the resulting RMSE and/or maximum OUH for the predicted size were only marginally worse. To better align with the approach taken by surgeons in said scenarios, the tool could be configured to always downsize the component when minimal difference in fit is predicted between two sizes (Dai et al., 2014).

Comparing the results to those available in the literature for manual planning/templating, the size prediction accuracy of the tool was shown to be considerably better than the average reported across the 10 studies summarized by Hernández-Vaquero et al. (2019). This was also true when compared to studies that utilized CT data for manual templating (Vaishya et al., 2018; Kobayashi et al., 2012). Therefore, due to the high accuracies obtained without the need for user feedback or training, it is anticipated that hospitals could realize considerable time and resource savings and reduce the possibility of human error by adopting the tool. In terms of other computational approaches, comparing the tool to results detailed previously for prediction models based on subject demographics (Trainor et al., 2018; Sershon et al., 2017) showed significantly improved performance was possible by using the approach outlined in this study. The tool achieved a better average tibia plate size prediction accuracy and matched the level reported for femur components (in terms of RMSE) when compared to the solution published by Zheng et al. (2018). When compared to the results reported by Massé and Ghate (2021), the tool was on average >20% more accurate in selecting the best size for both component types. This study has therefore shown that fully automating the size prediction tool did not compromise accuracy, but similar or better results were in fact obtained. Moreover, the tools developed by Zheng et al. (2018) and Massé and Ghate (2021) were only tested on 23 and 45 subjects respectively, compared to 78 in this work. The authors also only used single implant models and one global fit metric in the cited studies. By exploring a range of implant model sizes and employing both global and clinical fit metrics to evaluate performance, this study has provided a more comprehensive assessment of the developed tool’s robustness and applicable use with various products.

The performance of the tool could be further improved by controlling the alignment and quality of the inputted X-ray images to ensure compatibility with the 2D—3D pipeline (Burge et al., 2022b). Due to the datasets available, the testing of the tool was limited to predominantly White Americans without severe arthritis degeneration, whilst the modules of the tool were principally trained using Asian Korean data. Differences in performance between sexes and with age were evaluated in the study, however, with the data available it was not possible to analyze potential biases across ethnicities and/or the impact of arthritis severity on performance. For the development of future implant size prediction tools, it is encouraged that larger datasets, featuring hundreds of subjects with various ethnicities, ages, sexes and Kellgren and Lawrence grades, all imaged with consistent alignments/scan settings, be obtained. This would help better generalize the CNN models and improve robustness across the full range of potential patient demographics. Moreover, to further reduce overfitting, regularization techniques such as drop-out layers could be incorporated into the CNN model architecture. It should be noted that the method for determining the best fitting implant designs varied between the studies referenced. The best size in most of the referenced studies was determined by a medical professional during or after surgery. This study however was completed without clinical work and used computational fit metrics to determine performance (like in Zheng et al. (2018)). Future work could seek to test the size prediction tool against size choices made during real surgical procedures to provide a more realistic means for comparison. Furthermore, generic implant designs, scaled to published manufacturer model size charts, were used instead of official geometries. Going forwards, it is hoped manufacturers will adapt the framework of the tool for use with commercially available products.

5 Conclusion

In this study a computational TKR size prediction tool was developed that uses two X-ray images to assist clinicians with selecting the best implant sizes for patients. The tool achieved selection accuracies superior to those reported for manual templating, as well as when compared to other computational alternatives. By removing the need for manual templating or guiding semi-automatic software tools, the tool could minimize the time and resource required for TKR preoperative assessment, increase surgeons’ confidence in the outcomes, and minimize the possibility of size changes during surgery, whilst still achieving reliable size prediction results. This could also help reduce surgery time, minimize the level of component inventory required, lessen surgeon accountability for determining the right size implants, and ultimately improve outcomes for patients. Finally, the tool could be used to assess the need for customised solutions, or for a bespoke selection of implants depending on patient morphology.

Data availability statement

The datasets presented in this article are not readily available because access to the data used in the study is restricted and controlled by the OAI and KISTI. Requests to access the datasets should be directed to https://nda.nih.gov/oai/and https://www.kisti.re.kr/eng/.

Author contributions

TB: study design, algorithm development, data generation, results processing, manuscript writing. GJ and CJ: study design, advice on current templating and surgical procedures, reviewing manuscript/results. JJ and CM: project supervision, study design, reviewing manuscript/results.

Funding

The authors would like to thank GSK for funding and supporting the research.

Acknowledgments

Data and/or research tools used in the preparation of this manuscript were obtained and analyzed from the controlled access datasets distributed from the Osteoarthritis Initiative (OAI), a data repository housed within the NIMH Data Archive (NDA). OAI is a collaborative informatics system created by the National Institute of Mental Health and the National Institute of Arthritis, Musculoskeletal and Skin Diseases (NIAMS) to provide a worldwide resource to quicken the pace of biomarker identification, scientific investigation and OA drug development. Dataset identifier(s): (10.15154/1524779). The study also used data from the Digital Korean provided by the Korea Institute of Science and Technology Information (KISTI).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

Buller, L. T., Menken, L., and Rodriguez, J. A. (2018). The custom total knee replacement: A bespoke solution. Seminars Arthroplasty 29 (3), 209–213. doi:10.1053/j.sart.2019.01.006

CrossRef Full Text | Google Scholar

Burge, T. A., Jeffers, J. R. T., and Myant, C. W. (2022a). Development of an automated mass-customization pipeline for knee replacement surgery using biplanar X-rays J. Mech. Des. 144 (2), 1–11. doi:10.1115/1.4052192

CrossRef Full Text | Google Scholar

Burge, T. A., Jeffers, J. R. T., and Myant, C. W. (2022b). Performance and sensitivity analysis of an automated X-ray based total knee replacement mass-customization pipeline. J. Med. Devices., 1–12. doi:10.1115/1.4055000

CrossRef Full Text | Google Scholar

Cernazanu-Glavan, C., and Holban, S. (2013). Segmentation of bone structure in X-ray images using convolutional neural network. Adv. Electr. Comput. Eng. 13 (1), 87–94. doi:10.4316/AECE.2013.01015

CrossRef Full Text | Google Scholar

Clary, C., Aram, L., Deffenbaugh, D., and Heldreth, M. (2014). Tibial base design and patient morphology affecting tibial coverage and rotational alignment after total knee arthroplasty. Knee Surg. Sports Traumatol. Arthrosc. 22 (12), 3012–3018. doi:10.1007/s00167-014-3402-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Culler, S. D., Martin, G. M., and Swearingen, A. (2017). Comparison of adverse events rates and hospital cost between customized individually made implants and standard off-the-shelf implants for total knee arthroplasty. Arthroplasty Today 3 (4), 257–263. doi:10.1016/j.artd.2017.05.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Dai, Y., Scuderi, G. R., Penninger, C., Bischoff, J. E., and Rosenberg, A. (2014). Increased shape and size offerings of femoral components improve fit during total knee arthroplasty. Knee Surg. Sports Traumatol. Arthrosc. 22, 2931–2940. doi:10.1007/s00167-014-3163-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Fedorov, A., Beichel, R., Kalpathy-Cramer, J., Finet, J., Fillion-Robin, J. C., Pujol, S., et al. (2012). 3D slicer as an image computing platform for the quantitative imaging network. Magn. Reson. Imaging 30 (9), 1323–1341. doi:10.1016/j.mri.2012.05.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Hernández-Vaquero, D., Noriega-Fernandez, A., Roncero-Gonzalez, S., Perez-Coto, I., Sierra-Pereira, A. A., and Sandoval-Garcia, M. A. (2019). Agreement in component size between preoperative measurement, navigation and final implant in total knee replacement. J. Orthop. Transl. 18, 84–91. doi:10.1016/j.jot.2018.10.006

CrossRef Full Text | Google Scholar

Hitt, K., Shurman, J. R., Greene, K., McGarthy, J., Moskal, J., Hoeman, T., et al. (2003). Anthropometric measurements of the human knee: Correlation to the sizing of current knee arthroplasty systems. J. Bone Jt. Surgery-American 85, 115–122. doi:10.2106/00004623-200300004-00015

PubMed Abstract | CrossRef Full Text | Google Scholar

Kellgren, J. H., and Lawrence, J. S. (1956) ‘Radiological assessment of osteo-arthrosis’, Ann. Rheumatic Dis., 16(4), pp. 494–502. doi:10.1136/ard.16.4.494

PubMed Abstract | CrossRef Full Text | Google Scholar

Kobayashi, A., Ishii, Y., Takeda, M., Noguchi, H., Higuchi, H., and Toyabe, S. (2012). Comparison of analog 2D and digital 3D preoperative templating for predicting implant size in total knee arthroplasty. Comput. Aided Surg. 17 (2), 96–101. doi:10.3109/10929088.2011.651488

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, S. H., and Lee, S. B. (2010). Production and usage of Korean human information in KISTI. J. Korea Contents Assoc. 10 (5), 416–421. doi:10.5392/jkca.2010.10.5.416

CrossRef Full Text | Google Scholar

Li, K., Saffarini, M., Valluy, J., Desseroit, M. C., Morvan, Y., Telmon, N., et al. (2019). Sexual and ethnic polymorphism render prosthetic overhang and under-coverage inevitable using off-the shelf TKA implants. Knee Surg. Sports Traumatol. Arthrosc. 27 (7), 2130–2139. doi:10.1007/s00167-019-05410-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, C., Cao, Y., Alcantara, M., Liu, B., Brunette, M., Peinado, J., et al. (2018). “TX-CNN: Detecting tuberculosis in chest X-ray images using convolutional neural network,” in Proceedings - International Conference on Image Processing, Beijing, China, 17-20 September 2017. IEEE, 2314–2318. doi:10.1109/ICIP.2017.8296695

CrossRef Full Text | Google Scholar

Mahoney, O. M., and Kinsey, T. (2010). Overhang of the femoral component in total knee arthroplasty: Risk factors and clinical consequences. J. Bone Jt. Surgery-American 92 (5), 1115–1121. doi:10.2106/JBJS.H.00434

PubMed Abstract | CrossRef Full Text | Google Scholar

Massé, V., and Ghate, R. S. (2021). Using standard X-ray images to create 3D digital bone models and patient-matched guides for aiding implant positioning and sizing in total knee arthroplasty. Comput. Assist. Surg. 26 (1), 31–40. doi:10.1080/24699322.2021.1894239

CrossRef Full Text | Google Scholar

Mohammed, M. A., Abdulkareem, K. H., Garcia-Zapirain, B., Mostafa, S. A., Maashi, M. S., Al-Waisy, A. S., et al. (2021). A comprehensive investigation of machine learning feature extraction and classification methods for automated diagnosis of COVID-19 based on X-ray images. Comput. Mat. Contin. 66 (3), 3289–3310. doi:10.32604/cmc.2021.012874

CrossRef Full Text | Google Scholar

Narin, A., Kaya, C., and Pamuk, Z. (2021). Automatic detection of coronavirus disease (COVID-19) using X-ray images and deep convolutional neural networks. Pattern Anal. Appl. 24 (3), 1207–1220. doi:10.1007/s10044-021-00984-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Nevitt, M., Felson, D., and Lester, G. (2006). The osteoarthritis initiative: Protocol for the cohort study. Osteoarthr. Initiative 16, 1–74.

Google Scholar

Ronneberger, O., Fischer, P., and Brox, T. (2015). U-Net: Convolutional networks for biomedical image segmentation. Med. Image Comput. Computer-Assisted Intervention – MICCAI 2015 9351, 234–241. doi:10.1007/978-3-319-24574-4

CrossRef Full Text | Google Scholar

Schnurr, C., Csécsei, G., Nessler, J., Eysel, P., and König, D. P. (2011). How much tibial resection is required in total knee arthroplasty? Int. Orthop. 35 (7), 989–994. doi:10.1007/s00264-010-1025-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Schroeder, L., and Martin, G. (2019). In vivo tibial fit and rotational analysis of a customized, patient-specific TKA versus off-the-shelf TKA. J. Knee Surg. 32 (6), 499–505. doi:10.1055/s-0038-1653966

PubMed Abstract | CrossRef Full Text | Google Scholar

Sershon, R. A., Courtney, P. M., Rosenthal, B. D., Sporer, S. M., and Levine, B. R. (2017). Can demographic variables accurately predict component sizing in primary total knee arthroplasty? J. Arthroplasty 32 (10), 3004–3008. doi:10.1016/j.arth.2017.05.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Shao, L., Wu, X. D., Wang, T., Liu, X. K., Xu, W., Huang, W., et al. (2020). Approximating the maximum tibial coverage in total knee arthroplasty does not necessarily result in implant malrotation. Sci. Rep. 10, 10529–9. doi:10.1038/s41598-020-67613-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Sheth, N. P., Husain, A., and Nelson, C. L. (2017). Surgical techniques for total knee arthroplasty: Measured resection. J. Am. Acad. Orthop. Surg. 25 (7), 499–508. doi:10.5435/JAAOS-D-14-00320

PubMed Abstract | CrossRef Full Text | Google Scholar

Tanzer, M., and Makhdom, A. M. (2016). Preoperative planning in primary total knee arthroplasty. J. Am. Acad. Orthop. Surg. 24 (4), 220–230. doi:10.5435/JAAOS-D-14-00332

PubMed Abstract | CrossRef Full Text | Google Scholar

Trainor, S., Collins, J., Mulvey, H., and Fitz, W. (2018). Total knee replacement sizing: Shoe size is a better predictor for implant size than body height. Arch. Bone Jt. Surg. 6 (2), 100–104. doi:10.22038/abjs.2017.22499.1603

PubMed Abstract | CrossRef Full Text | Google Scholar

Vaishya, R., Vijay, V., Krishnan, M., and Agarwal, A. K. (2018). Fallacies of CT based component size prediction in total knee arthroplasty – are patient specific instruments the answer? J. Clin. Orthop. Trauma 9 (1), 34–39. doi:10.1016/j.jcot.2017.11.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Wernecke, G. C., Harris, I. A., Houang, M. T. W., Seeto, B. G., Chen, D. B., and MacDessi, S. J. (2012). Comparison of tibial bone coverage of 6 knee prostheses: A magnetic resonance imaging study with controlled rotation. J. Orthop. Surg. Hong. Kong. 20 (2), 143–147. doi:10.1177/230949901202000201

PubMed Abstract | CrossRef Full Text | Google Scholar

Zheng, G., Hommel, H., Akcoltekin, A., Thelen, B., Stifter, J., and Peersman, G. (2018). A novel technology for 3D knee prosthesis planning and treatment evaluation using 2D X-ray radiographs: A clinical evaluation. Int. J. Comput. Assist. Radiol. Surg. 13 (8), 1151–1158. doi:10.1007/s11548-018-1789-4

PubMed Abstract | CrossRef Full Text | Google Scholar