Abstract
The implant design of the talar component for total ankle replacement (TAR) should match the surface morphology of the talus so that the replaced ankle can restore the natural motion of the tibiotalar joint and may reduce postoperative complications. The purpose of this study was to introduce a new 3D fitting method (the two-sphere fitting method of the talar trochlea with three fitting resection planes) to approximate the shape of the upper part of the talus for the Chinese population. 90 models of the tali from CT images of healthy volunteers were used in this study. Geometrical fitting and morphological measurements were performed for the surface morphology of the upper part of the talus. The accuracy of the two-sphere fitting method of the talar trochlea was assessed by a comparison of previously reported data. Parameters of the fitting geometries with different sizes were recorded and compared. Results showed that compared with previously reported one-sphere, cylinder, and bitruncated cone fitting methods, the two-sphere fitting method presented the smallest maximum distance difference, indicating that talar trochlea can be approximated well as two spheres. The radius of the medial fitting sphere was mm which was significantly smaller than the radius of the lateral fitting sphere of mm. After grouping all data by the average radius of fitting spheres, the result showed that different sizes of the upper part of the talus presented significantly different parameters except the orientation of the lateral cutting plane, indicating that the orientation of the lateral cutting plane may keep consistent for all upper part of the talus and have no relationship with the size. The linear regression analyses demonstrated a weak correlation () between the majority of parameters and the average radius of the fitting spheres. Therefore, different sizes of the upper part of the talus presented unique morphological features, and the design of different sizes of talar components for TAR should consider the size-specific characteristics of the talus. The parameters measured in this study provided a further understanding of the talus and can guide the design of different sizes of the talar components of the TAR implant.
1. Background
The kinematics of the tibiotalar joint during the gait is highly influenced by the surface morphology of ankle bones [1]. Thus, the implant design for total ankle replacement (TAR) should be based on the surface morphology of the tibia and talus, so that the replaced ankle can restore the natural motion of the tibiotalar joint and may reduce postoperative complications [2–5]. The implantation of the talar component of the TAR implant requires the resection of the proximal part of the talus. Therefore, investigation of the surface morphology of the upper part of the talus, especially the talar trochlea, is essential to the design of the talar component of the TAR implant. However, the design of the talar component varies among existing implant systems, and no consensus has been established [2, 6].
Shape approximation is one of the useful methods to study bone morphology and design ankle implants. Although several studies have measured the morphology of the talar trochlea [7–12], only a few studies [13–16] approximated the shape of the talar trochlea from the perspective of TAR implant design. Note that the majority of studies analyzed the shape of the talus in CT images considering that the thickness of the articular cartilage layer of the talar trochlea was shown in the previous literature [17] to be small and vary slightly, averaging mm in males and in females. Therefore, not including the articular cartilage layer in measuring models derived from CT images would produce limited errors.
Pioneered work by Inman [18] measured the shape of the talar trochlea of cadaver specimens and found out that the talar trochlea can be approximated as a truncated cone with its axis oriented from medial-distal end to lateral-proximal end. Although concerns have been raised, the concept of a truncated cone with one fixed axis has been used in the design of several early implant systems, such as the STAR implant system (Figure 1(a)) and the Hintegra implant system (Figure 1(b)). Implant systems with such designs have a natural instability along the medial/lateral direction, which requires additional central ridge or edge rims to constrain the medial/lateral motion of the bearing. However, unknown biomechanical risks might be brought by these structures for stabilization. To further approximate the shape of the talar trochlea and avoid additional structures for stabilization, some researchers considered the talar trochlea as a saddle-shaped truncated cone [8, 19]. Few implant systems, such as INBONE Total Ankle System (INBONE Technologies, U.S.) (Figure 1(c)), adopted this design. But the saddle-shaped talar component of INBONE Total Ankle System was later revised to INBONE II Total Ankle System (INBONE Technologies, U.S.) (Figure 1(d)) due to the high implant failure rate [2]. It may be explained by the fact that the articular surface of the latter implant system with the shape of two intersected spheres has a deeper sulcus in the center of the articular surface, thus providing higher stability of the implant system.

(a)

(b)

(c)

(d)
Fitting the medial and lateral surfaces of the talar trochlea with two spheres seems to be a good way to approximate the surface of the talar trochlea. However, in previous literature, this two-sphere fitting method has only been used to define the rotational axis of the tibiotalar joint [20, 21]. To define a reliable 3D modeling method for quantifying the shape of the upper part of the talus of the Chinese population, the present study approximated the shape of the talar trochlea using the two-sphere fitting method and the medial and lateral surfaces using the plane fitting method. The key parameters of these fitting geometries, including the radius of the spheres, the position of sphere origins, the inclination angle, and the position of the medial and lateral surfaces, were recorded. The accuracy of the two-sphere fitting method was assessed by a comparison of previously reported data. Also, the size difference of fitting geometries was further investigated to guide the design of the implant components with different sizes. We hypothesized that the parameters of the fitting geometries have a strong relationship with the size of the upper part of the talus.
2. Methods
2.1. Subjects
The protocol of this study was approved by the Ethics Committee of Huashan Hospital, Fudan University. The CT scan data of 90 intact ankles from previous studies were used to approximate the shape of the talar trochlea (46 males, 44 females, years, 19-37 years of age range, cm of height, and kg of body weight, 48 left ankles and 42 right ankles) [22, 23]. No a priori power analysis was performed since there were no data available to perform such analysis. CT examination was performed on a CT scanner (Brilliance iCT, Philips, Cleveland, U.S.) with 120 kV of voltage, 250 mA of current, 0.625 mm of slice thickness, and pixels of matrix. The assessments of all CT images were performed by senior foot and ankle surgeons (X.W. and X.M.), and no signs of previous trauma, severe deformity, or degenerative changes in the ankle, such as ankle arthritis, osteochondral lesions of the talus, and talar cyst, were observed in the CT images of all ankles.
2.2. 3D Model Reconstruction
The Digital Imaging and Communications in Medicine (DICOM) file of each ankle was imported into Mimics (Materialise NV, Belgium) for the segmentation of the talus with a threshold of more than 203 Hounsfield units (HU) [24]. The 3D model of each talus was imported into 3-Matic Medical (Materialise NV, Belgium) for shape approximation and surface morphology measurement of the talar trochlea.
2.3. Shape Approximation of the Talar Trochlea
The two-sphere fitting method was described in previous studies [20, 21] to determine the rotation axis of the talus. We tried to apply this method to approximate the surface shape of the talar trochlea. This method used spheres to least-squares fit the medial or lateral part of the trochlear surface of the talus. The medial or lateral part of the trochlear surface of the talus was defined and manually selected as the facet surface between the central trochlea groove and the medial or lateral rim based on a curvature analysis of each talus. The curvature analysis built in the 3-Matic Medical (Materialise NV, Belgium) projected a color map of the surface curvature on the talar model which helped us identify the edge of each selection. The radius of the lateral sphere , the radius of the medial sphere , and the distance between sphere origins were recorded. To define the size of each talar trochlea, we used the average radius which was the mean of and of each talus. The model reconstruction process and sphere fitting method are illustrated in Figure 2.

(a)

(b)

(c)

(d)

(e)

(f)

(g)
2.4. Coordinate System Definition
The method of a coordinate system definition was based on the subject-specific bone geometry of the talus and similar to that previously described by Yamaguchi et al. [25]. For each talus, the medial-lateral axis (-axis) was defined as the line connecting the sphere origins (Point and ). The intersection of two spheres created an intersection circle. The origin of the coordinate system was defined as the center of the intersection circle, and the sagittal plane was defined as the intersection plane of the two spheres. The most anterior and posterior points of the medial and lateral edges of the talar trochlea were defined, and the midpoints of the most anterior or posterior points of both edges of the talar trochlea were calculated (Points and ). Then, a horizontal plane (plane ) perpendicular to the sagittal plane which contained these two midpoints was created, and the superior-inferior axis (-axis) was defined as the line perpendicular to the horizontal plane passing through the origin. The height () of the plane was recorded. The anterior-posterior axis (-axis) was a cross product of the superior-inferior and medial-lateral axes. Therefore, the coronal plane was formed by the intersection of the -axis and -axis, while the sagittal plane was formed by the intersection of the -axis and -axis. The Cartesian coordinate system was established in 3-Matics Medical (Materialise NV, Leuven, Belgium).
2.5. Shape Approximation of the Medial and Lateral Articular Surfaces
The medial or lateral articular surfaces of the upper part of the talus were defined and manually selected as the facet surface between the medial or lateral rim of the talus and the horizontal plane. The implantation of the talar component for total ankle replacement only requires the resection of the upper part of the talus. So, each talus was resected by the previously defined horizontal plane, and the medial and lateral articular surfaces can be fitted by two fitting planes (planes and ). To define these two fitting planes, the position of intersection points (Point ) among the coronal plane, horizontal plane, and was calculated (), while the intersection point (Point ) among the coronal plane, horizontal plane, and was also calculated (). The angles between the fitting plane to the horizontal or sagittal plane were also recorded (, , , and ). The establishment of the coordinate system and geometry parameters of the talar trochlea were measured in 3-Matics (Materialise NV, Leuven, Belgium) and presented in Figure 3.

(a)

(b)

(c)

(d)
2.6. Evaluating the Accuracy of the 3D Fitting Method and Statistical Analysis
To evaluate the accuracy of the current 3D fitting method, sensitivity tests were performed in Supplementary File 1 to assess the stability of both the two-sphere fitting method to the talar trochlea and the plane fitting method to the medial and lateral articular surfaces. Also, the Euclidean distances between vertices on the talar trochlear surface and the geometries of two spheres were calculated and compared with the reported value of the single sphere, cylinder, or bitruncated cone fitting method by Huang et al. [16].
All data were presented as mean values with standard deviation (SD). To investigate the size difference of each parameter, we chose to divide all data into 4 groups based on the mean and standard deviation of the average radius (). All key parameters (: radius of the lateral sphere; : radius of the lateral sphere; : the distance between sphere origins; : the height of the plane ; : the position of the medial fitting plane; : the position of the lateral fitting plane;: the angle between the horizontal plane and the medial fitting plane; : the angle between the horizontal plane and the lateral fitting plane; : the angle between the sagittal plane and the medial fitting plane; : the angle between the sagittal plane and the lateral fitting plane) both overall or within each size group were checked for normal distribution and homogeneity of variance. A paired -test was used to determine the significance of the difference between the radius of the medial and lateral spheres or between two poses. One-way repeated measures Analysis of Variance (ANOVA) was conducted to compare the distances or angles among each group. Statistical significance level was set at 0.05. The linear regressions between the average radii and other parameters of the fitting models (, , , , , , , and ) were performed, and the correlation coefficients () were computed. Data were processed in MATLAB (2018b, MathWorks Inc., Natick, MA).
3. Results
The sensitivity analyses in Supplementary File 1 showed that changing the area of the selected trochlear surface of the talus (see Figure S1) made a small impact on the radius of the spheres and the position of the sphere origin (maximum difference of the radii of the medial and lateral fitting spheres was both less than 5%). Also, changing the area of the selected medial articular surfaces of the talus (Figure S2) made a small impact on the orientation of the medial fitting plane (the maximum difference of the inclination angle was less than 5%). However, due to the curved shape of the lateral articular surface, the inclination angle of the lateral fitting plane varied as the selection area changed. Therefore, the selection of the lateral articular surface should be carefully selected above the horizontal plane.
Compared with fitting methods of one sphere, cylinder, and bitruncated cone, the mean, standard deviation, and maximum of the distances between the talar trochlea and two-sphere model are given in Table 1 [16]. The two-sphere fitting method presented a relatively small mean and the smallest standard deviation and maximum (, ,, respectively).
The mean and SD of all parameters grouped by the mean and SD of the average radius are listed in Table 2. The overall average radius was . Thus, the four groups were average mm, average and <21.00 mm, average and <23.09 mm, and average mm. All values of the geometry parameters of the talar trochlea (, , ,, , , , , , ) are presented in Supplementary File 2, and they were normally distributed in both overall and within each size group category. Also, the variance of each parameter was not significantly different. The radius of the medial fitting sphere was mm which was significantly smaller than the radius of the lateral fitting sphere of mm (). Statistically significant differences among four groups were found in , , ,, , , , , , and . The mean of , , ,, , , and increased with the average radius, while the mean of decreased with the average radius.
The result of linear regressions between the average radii and other parameters of the fitting models is illustrated in Figure 4, and the correlation coefficients for ,, , , , , , and were 0.1376, 0.6534, 0.4074, 0.1706, 0.1116, 0.0871, 0.0864, and 0.0294, respectively.

4. Discussion
In this study, we proposed a method for shape approximation of talar trochlea with two spheres resected by three fitting planes and investigated the size difference of the parameters for this fitting method (Figure 3). Calculation of the distances between the surface of the talar trochlea and the two-sphere model showed that the current fitting method was an accurate shape approximate method to the talar trochlea. The sensitivity test also showed that fitting the talar trochlea with two spheres was a stable method with a small selection bias.
Since the mean and standard deviation of the distance difference for all fitting methods were below the slicing distance of the CT scan, the maximums of the distances between vertices on the talar trochlea and geometrical models were primarily used to assess the accuracy of the fitting method. Compared with other fitting methods, the two-sphere fitting method presented the smallest maximum (), suggesting that the similarity between the talar trochlea and the two-sphere fitting method was more prominent and this method seems to be a reasonable choice for shape approximation of the talar trochlea.
Several studies have previously reported the radii of the talar trochlea, and their results are listed in Table 3 [7–9, 15, 22, 26, 27]. The disagreements on the measurement of the medial and lateral curvature radii of the talar trochlea seem to be well explained by the truncated cone model presented by Inman [18]. The cross-sections of a cone intersected by two planes can be two ellipses. Therefore, the medial and lateral radii of curvature of the talar trochlea may not be equal. Similarly, the anterior and posterior radii of curvature of one side of the talar trochlea may not be the same. Modified from the cone fitting method, the two-sphere fitting method shared the rotational axis of the cone but provided a closer shape approximation at the central groove and medial and lateral rims. The average and standard deviation of the radii of the talar trochlea measured in this study were in the range of previously reported values [7–9, 15, 22, 26, 27].
However, the radii measured in this study showed a high variety, indicating high anatomical variability of the size of talar trochlea in the Chinese population. The ranges of the medial and lateral radii were 16.39 to 26.27 degrees and 16.46 to 25.18 degrees, respectively. Although statistical analysis showed that the radius of the lateral fitting sphere was found to be significantly larger than that of the medial fitting sphere, 27 out of 90 tali (30%) had a larger medial sphere radius . Therefore, the design of the talar component for total ankle replacement should prepare extra configuration for the need of special cases.
Only one of the existing total ankle implant systems (Cadence implant system, Integra Lifesciences, Cordova, TN, U.S.) had different medial and lateral sphere radii in the design of talar component for TAR with 8 degrees of conical angle, and the lateral sphere is larger than the medial side. This value was derived from the mean angle between the axis and the horizontal plane of the ankle measured by Inman [18]. The current study also found that the lateral sphere presented a statistically significantly larger mean radius than the medial sphere. The mean conical angle of the current study can be calculated from the value of mean lateral radius (), mean medial radius (), and the mean distance between sphere origins () as , which was 2.29 degrees. The difference between these may be due to the population difference.
After grouping all data by the average radius of fitting spheres, the result showed that different sizes of the upper part of the talus presented significantly different parameters except the orientation of the lateral cutting plane, indicating that the orientation of the lateral cutting plane may keep consistent for all upper part of the talus and have no relationship with the size. The result of the linear regression showed that the majority of parameters showed a weak relationship (less than 0.5) to the average radius except for the height of with a value of 0.6534. was determined by the medial and lateral edges of the talar trochlea which were strongly related to the medial and lateral radii. Small also rejected our conjecture that parameters of the fitting geometries have a strong relationship with the size of the upper part of the talus. Therefore, the design of the different sizes of talar components for TAR should consider the size-specific characteristics of the talus instead of simply scaling the shape of one general model.
There were some limitations of the current study. First, the sample size of this study is limited. Although our sample size was larger than the majority of previous studies [7–9, 15, 22, 26, 27], future studies should recruit a large number of subjects to increase the power of the study and expand the participants to older population or patients with damaged talar trochlea to improve the clinical relevance. Second, the two-sphere fitting method presented here for shape approximation only studied the surface of the talus. The mated surfaces of the tibia and fibula within the tibiotalar joint and the anatomy of the tibiofibular syndesmosis should also be considered in future studies. The kinematics of the tibiotalar should also be considered. Future studies should evaluate whether the axis by connecting the fitting sphere origins was the rotational axis of the tibiotalar joint during the gait. Third, the two-sphere fitting method proposed in the current study has not been proved as the best method for talar shape approximation with the most accuracy. Future studies should compare the two-sphere fitting method with more new geometric fitting methods to better approximate the morphology of the talus. Last, the current fitting geometry of two spheres with three resection planes is only the first step of implant design for TAR. Before a finished product, additional manufacturing processes will be required, such as fillets will be added to the edges to obtain smooth transitions from surfaces to surfaces to avoid shape edges of the fitting geometry. Still, the results of this study could be a useful reference for implant component design for TAR.
5. Conclusion
The two-sphere fitting method of the talar trochlear was a stable and accurate shape approximation method with the smallest difference of maximum distance. The medial fitting sphere has a significantly smaller radius than the lateral fitting sphere. The presented fitting method with two spheres and three resection planes can be used as a standardized definition for shape approximation of the upper part of the talus. Also, different sizes of the upper part of the talus presented unique morphological features, and the design of different sizes of talar components for TAR should consider the size-specific characteristics of the talus. The parameters measured in this study provided a further understanding of the talus and can guide the design of different sizes of the talar components of the TAR implant.
Data Availability
The datasets used and analyzed during the current study are available from the corresponding author on reasonable request.
Conflicts of Interest
The authors received no financial support for the research, authorship, and/or publication of this article.
Authors’ Contributions
Jian Yu and Dahang Zhao contributed equally.
Acknowledgments
This project was supported by National Natural Science Foundation of China (Grant No. 82072388), Clinical Research Plan of SHDC (Grant No. SHDC2020CR3071B), and Natural Science Foundation of Shanghai (Grant No. 19ZR1407400).
Supplementary Materials
The present study included two supplementary files, which were all referred to in the manuscript. (Supplementary Materials)