Is Micro-Computed Tomography Useful for Wear Assessment of Ceramic Femoral Heads? A Preliminary Evaluation of Volume Measurements

Abstract

Background

Wear associated with hip components represents the main clinical problem in these patients, and it is important to develop new techniques for more accurate measurements of that wear. Currently, the gravimetric method is the gold standard for assessing mass measurements in preclinical evaluations. However, this method does not give other information such as volumetric loss or surface deviation. This work aimed to develop and validate a new technique to quantify ceramic volume loss from in vitro experiments using micro-computed tomography (micro-CT).

Methods

An alumina (BIOLOX® forte) femoral head (Ø = 28 mm) was used. Mass and volume loss were approached by gravimetric method (using a four decimal place digital microbalance) and by using Skyscan 1176 microtomographic system, respectively.

Results

Standard error and coefficient of variance of both gravimetric and experimental groups demonstrated the reliability of the micro-CT analysis technique.

Conclusions

In conclusion, the findings of the present study suggest that this new protocol could be considered an important tool for wear assessment and that we have found a reliable metrological protocol for volumetric analysis of ceramic femoral head prostheses, demonstrating that the micro-CT technique can be an important tool for wear assessment.

Keywords

Dimensional metrology Gravimetric method Micro-CT Microtomography Volumetric wear

Introduction

Wear associated with hip components represents the main clinical problem of total hip arthroplasty (THA) treatment (1–2–3), and new techniques for more accurate measurements of wear are to be welcomed (4, 5). Currently, laboratory determination of wear rate of joint bearing constitutes an important aspect in preclinical validation of prostheses (6). Numerous measurement techniques are used to determine wear for both in vivo (i.e., radiostereometric analysis, radiography etc.) and in vitro (i.e., gravimetric method, coordinate measuring machine etc.) prosthetic hip joint components (4). In in vitro studies, wear can be measured quite accurately before, during and after the entire hip simulation. The gravimetric method, as suggested by Carmignato and coworkers (5), is the so-called gold standard practice adopted by orthopedic industry for wear evaluation of hip joint components. A microbalance is used for this procedure, and the mass loss of the specimens is measured following a procedure recommended by international guidelines (ISO 14242-2:2012). Currently the gravimetric method is the most commonly used measurement approach to measure wear in vitro (4–5–6–7).

Recently, volumetric methods such as the Coordinate Measuring Machine (CMM) and micro-computed tomography (micro-CT) techniques have been used as alternatives to overcome the limitations of the gravimetric method (5, 7–8–9). In particular, micro-CT is considered a promising technology for dimensional metrology (10, 11). In addition, geometrical approaches that use micro-CT would be the most direct ones to understand the effects of wear on contact geometry. In the case of a hip prosthesis, the effectiveness of micro-CT techniques to assess the volumetric mass loss during wear tests is not fully proven yet. Micro-CT is a nondestructive testing technique based on X-rays that gives quantitative 3-dimensional (3D) information for a whole object structure in a relatively short time (12–13–14). The most common micro-CT systems are equipped with a polyenergetic X-ray beam source (15). In these systems, low-energy photons are attenuated more rapidly than high-energy photons while they are penetrating dense materials (beam hardening effect), as reported by other authors (16–17–18–19–20–21). In addition to beam hardening, the accuracy of the images and the subsequent accuracy of micro-CT-based measurements, as reported by Kruth and coworkers (14), is influenced by the detector pixel size and by the quality of the reconstructed 2D images. In fact, both the choice of threshold value and the image-processing software techniques used to enhance sharpness or to cut image noise (e.g., speckles) directly affect the accuracy of dimensional measurements (22, 23).

The main goal of this study was to determine the appropriateness of micro-CT technique in assessing the volume of ceramic femoral heads in anticipation of a wear analysis. In particular, this study tried to determine the accuracy and precision of micro-CT volume measurements and to statistically compare the results with the gravimetric data.

Methods

An alumina (Biolox® Forte) femoral head (Ø = 28 mm) was used. For the gravimetric method, the mass of the sample was measured using a 4 decimal place digital microbalance (Denver Instruments, Bohemia, NY, USA) with a declared repeatability of 0.1 mg. The true density value, used to calculate the weight in a given volume, was measured by helium pycnometry (Quantachrome™ Ultrapycnometer; Quantachrome, Boynton Beach, FL, USA).

Afterwards, for the volumetric method, the alumina femoral head component was scanned 30 times with a Skyscan 1176 microtomographic system (Bruker Micro-CT; Bruker, Kontich, Belgium) using 2 different nominal resolutions (17.5- and 35-μm pixel size). The 2D detector resolution was set at 668 × 1,000 or 1,336 × 2,000 pixels depending on nominal resolution. In agreement with other protocols, source voltage and current were set to 90 kV and 278 μA, respectively, using a copper filter of 0.1 mm between the source and the sample to perform the scanning with only the hard spectrum of the beam. The X-ray beam and detector were rotated 360° with rotation steps of 0.4° and frame averaging of 3. The duration of each scan was nearly 60 minutes (Skyscan 1176 software; Bruker) to obtain nearly 1,500 images in 16-bit TIFF file format for each dataset.

The software NRecon (version 1.6.8.0; Bruker, Kontich, Belgium) was used to reconstruct these datasets to obtain the microtomographic sections using different specifications for beam hardening corrections – i.e., 25%, 40% and 50% (details in Fig. 1). For the resolution of 35 μm, the images were reconstructed with only the beam hardening of 40%, since this was found to be the most accurate percentage of beam hardening correction, from the analysis of datasets with a resolution of 17.5 μm. After reconstruction, 4 groups of 15 datasets were thus collected.

Fig. 1

Scheme of micro-computed tomography (micro-CT) analysis and group collection.

A new dedicated software (CTAn version 1.13.11.0; Bruker, Kontich, Belgium) was used to analyze reconstructed images. First, the density phase in the images was segmented considering a gray value threshold calculated using the method proposed by Otsu (24). This method automatically finds those thresholds that minimize interclass variance of the gray-scale histogram. The volume in cubic millimeters was directly calculated on these segmented images. Subsequently, to minimize the noise in the digital datasets, 2 levels of despeckle were applied: first all of the 2D objects except the largest one, as previously identified by the Otsu thresholding, were removed, and then the same operation was repeated in the 3D space. This procedure was also carried out on reconstructed datasets – i.e. on gray-level (nonbinarized) images after an unsharp mask filtering. This method enhances the edge details by adding the difference between the original and the blurred images to the original. With this image analysis process, 24 groups of 15 datasets each were collected (Fig. 1).

SPSS software (version 18.0; IBM Inc., Chicago, IL, USA) and OpenStat (open-source software by G.W. Miller, billmiller@gmail.com) were used to perform the statistical analysis. In particular, the statistical analyses were performed in 2 steps: during the first step, 1-way ANOVA with relevant post hoc tests (Dunnett and Bonferroni) (25, 26) was used to select a reduced number of groups to compare the 2 above-mentioned methods (gravimetric vs. micro-CT). The statistical level of significance was set at a p value <0.05. During the second phase, our analyses were performed to rule out statistical differences. Thus, 1-way ANOVA with relevant post hoc tests (Dunnett and Bonferroni), Student's t-test and Kruskal-Wallis analysis with the relevant post hoc test (pairwise) were used to find the best method equivalent to the gravimetric ones.

Results

To determine which micro-CT protocol, and consequently which group, was the most similar or equivalent to the gravimetric method in measuring volumes of alumina femoral head components, 15 repeated measurements of 24 different groups were analyzed. In the first phase of the analysis, descriptive statistics were generated (data were normally distributed), and the means with the relevant standard error (SE) and coefficient of variance (CV) of both gravimetric and experimental groups are presented in Table I.

TABLE I

Volume of femoral head component for gravimetric group and each of the 24 experimental micro-CT groups

Group	Volume ± SE (mm³)	Volume CV (%)
Gravimetric group	7,915.73 ± 0.00	0.00
1	7,964.82 ± 2.12	0.10
2	7,964.40 ± 2.12	0.10
3	7,964.79 ± 2.12	0.10
4	7,959.44 ± 2.25	0.11
5	7,957.24 ± 2.25	0.11
6	7,958.27 ± 2.25	0.11
7	7,925.78 ± 2.75	0.13
8	7,925.32 ± 2.75	0.13
9	7,925.74 ± 2.75	0.13
10	7,916.65 ± 2.95	0.14
11	7,913.59 ± 2.99	0.15
12	7,914.79 ± 2.98	0.15
13	7,889.99 ± 3.91	0.19
14	7,889.42 ± 3.91	0.19
15	7,889.93 ± 3.91	0.19
16	7,869.13 ± 4.17	0.21
17	7,865.06 ± 4.20	0.21
18	7,866.66 ± 4.19	0.21
19	7,954.83 ± 4.77	0.23
20	7,954.58 ± 4.77	0.23
21	7,954.82 ± 4.77	0.23
22	7,954.78 ± 4.76	0.23
23	7,954.06 ± 4.76	0.23
24	7,954.63 ± 4.76	0.23

Values are the means, standard error (SE) and coefficient of variance (CV) based on 15 repeated measurements derived from the same number of datasets acquired for each resolution.

Micro-CT = micro-computed tomography.

SE and CV values demonstrated the reliability of the micro-CT analysis technique. The variability in volume measurements followed an increasing trend as the group number incremented. The greatest variability was found for groups 19 to 24 (CV = 0.23%) corresponding to the micro-CT analysis protocols with the lower resolution of 35 µm. The lowest variability was found for volume measurements of groups 1 to 6 (range of CV 0.10%-0.11%) corresponding to the micro-CT analysis protocols with the lower beam hardening correction (25%).

The first phase of the statistical analysis using 1-way ANOVA with relevant post hoc tests (Dunnett and Bonferroni) and a p value <0.05 between all 24 groups and the control group resulted in the selection of 6 groups to be statistically reassessed: groups 7 to 12 (Fig. 2). In fact, these 6 groups were different without statistical significance, compared with the control group, whereas all of the other groups were significantly different, with a p value always close to zero. The 6 nonsignificantly different groups corresponded to the micro-CT analysis protocols with a beam hardening correction of 40%.

Fig. 2

Box plots for the nonsignificantly different experimental micro-computed tomography (micro-CT) groups, compared with the gravimetric group after the first statistical phase using 1-way ANOVA analysis with relevant post hoc tests (Bonferroni and Dunnett) and with a p value <0.05.

In the second phase of the statistical analysis, the p values calculated using 1-way ANOVA with relevant post hoc tests (Dunnett and Bonferroni), Student's t-test and Kruskal-Wallis analysis, with the relevant post hoc test (pairwise) between the 6 groups and the control group, were compared. The results of this second phase are shown in Table II. The greatest p value was found for group 10. This group corresponds to a 3D Otsu segmentation after an unsharp masking filtering of the image dataset. The dataset was obtained after a setting of 17.5 µm/pixel and a beam hardening correction of 40%. The micro-CT images of the most and least similar groups compared with the control group are shown in Figure 3.

TABLE II

The p values obtained using 1-way ANOVA with post hoc tests (Dunnett and Bonferroni), Student's t-test and Kruskal-Wallis analysis with relevant post hoc test (pairwise), between the 6 experimental micro-CT groups that were nonsignificantly different after the first phase of analysis, and the gravimetric group

Group number	Bonferroni	Dunnett	Student's t-test	Pairwise
7	0.1793	0.0416	0.0010	0.01
8	0.2518	0.0570	0.0016	0.08
9	0.1848	0.0428	0.0011	0.03
10	1.0000	0.9998	0.7576
11	1.0000	0.9835	0.4790	No values^*
12	1.0000	0.9998	0.7536

The group with the least difference from the gravimetric group is highlighted.

KW = Kruskal-Wallis; micro-CT = micro-computed tomography.

KW test resulted in a p value of 0.394, and a post hoc test was not completed.

Fig. 3

Micro-computed tomography (micro-CT) images of the most and least similar experimental micro-CT groups obtained in the statistical analysis, compared with the gravimetric group. Interclass variance of the gray-scale histogram is colored to highlight the pixels. Microtomographic sections along the 3 spatial planes, surface detail of micro-CT images in the axial plane and the segmented image: (A-C) settings of 17.5 µm/pixel, beam hardening correction of 40%, unsharp masking filter; (D-F) 17.5 µm/pixel, beam hardening correction of 40%, with no unsharp masking filter; (G-I) 17.5 µm/pixel, beam hardening correction of 25%, with no unsharp masking filter; (J-L) 35 µm/pixel, beam hardening correction of 40%, with no unsharp masking filter.

Discussion

Currently, prosthetic failures represent a clinical problem, and eliminating or reducing wear would play a crucial role in preventing them. In this work we asked whether a new micro-CT approach would result in agreement with the so-called gold standard gravimetric method, in measuring the volume of ceramic femoral heads. The gravimetric method is the standard procedure to measure loss of mass. Among the volumetric-based methods, micro-CT has some advantageous characteristics such as the relatively short scanning time for each specimen, the ability to evaluate specimens without touching their surfaces and the capability of analyzing the specimen structures even locally. However, the use of a micro-CT-based procedure to assess the wear of ceramic femur components has not been fully proven and validated yet. For this purpose, the volume of ceramic femoral heads was calculated with a commercial high-resolution micro-CT system, which is relatively frequently suited for routine laboratory or preclinical applications. Furthermore, state-of-the-art micro-CT analysis on prosthesis component wear considers only those made of polyethylene (8, 27, 28).

After 2 subsequent statistical analyses, we identified a new protocol in which volumetric mass loss was comparable between the micro-CT and the gravimetric method. The volume measurement of the best matching group differed from the measurement of the gravimetric group by less than 1 mm³ on average. In addition, the variability in volume measurements was 0.14% for CV, with an SE of 2.95 mm³.

Despite the statistical overlap, the relative SE was still higher than that for the gravimetric group. The aim of this preliminary study was to evaluate the possibility of using a nondestructive, imaging-based method. The appreciable advantages of using the micro-CT technique in wear assessments are the possibility of conducting longitudinal quantitative analyses and localizing the wear features on the hip replacement surface. The process of wear assessment will involve measurements based on the volume difference (before and after) as a result of in vitro wear tests. From this perspective, the difference of 1 mm³ observed between the chosen experimental method and the gravimetric one may not be relevant in determining the degree of wear. As regards the SE, the issue could be more problematic. Other studies have shown that the wear rate in the first phase of in vitro testing of alumina vs. alumina is between 0.1 and 0.2 mm³ for millions of cycles (29, 30).

Among the possible causes that affect the precision of micro-CT, we can identify (i) the possibly not perfect alignment of the specimens during the acquisition, (ii) the variability of the polychromatic X-ray beam and (iii) subsequent variability in beam hardening effect. Interestingly, using beam hardening correction algorithms with increasing intensity, whose effectiveness in improving the dimensional accuracy has not yet been verified (31), an increasing variability was detected.

This study has some limitations connected mainly with the small number of femoral heads analyzed. A study of many femoral prosthetic components would allow the elimination of systematic errors. This could result in a reduction in the measured uncertainty, and the samples could be analyzed varying the orientation in the 3 spatial planes during acquisition.

Further work will apply the protocol to the unworn and worn specimens of in vitro wear simulations to evaluate mass loss. Wear analysis conducted using micro-CT will be geometrically free. This means that no unworn reference surface geometry will be coaligned to measure the wear. The difference between the worn and the unworn volume of the specimen will just be calculated after only the image analysis operations. The advantage to this will lie in the fact that without a reference surface, the measurements will be influenced only by the micro-CT protocol adopted.

Conclusions

To date, there is still no literature available concerning international standards or guidelines to evaluate the uncertainty in CT measurements. The volume evaluation proposed in this study is repeatable and operator independent. The results could lay the foundations for international standardization of volumetric analysis by micro-CT, which is currently lacking due to the numerous influencing factors and their complex interactions in determining uncertain results.

In conclusion, we found a reliable metrological protocol for the volumetric analysis of ceramic femoral head prostheses, demonstrating that the micro-CT technique could be an important tool for wear assessment.

Footnotes

Acknowledgements

The authors wish to thank Dr. Giovanni Ridolfi from Centro Ceramico Bologna for femoral head density measurements.

Financial support: This work was partially supported by the Italian Program of Donation for Research “5 per Mille,” for the year 2011.

Conflict of interest: None of the authors has any financial interest related to this study to disclose.

References

Callaghan

Pedersen

Johnston

Brown

Clinical biomechanics of wear in total hip arthroplasty.

Iowa Orthop J 2003231–12

Affatato

Leardini

Jedenmalm

Ruggeri

Toni

Larger diameter bearings reduce wear in metal-on-metal hip implants.

Clin Orthop Relat Res 2007456456153–158

Sentürk

Perka

[Ceramic-on-ceramic bearings in total hip arthroplasty (THA)] [article in German].

Z Orthop Unfall 20151532198–202

Grillini

Affatato

How to measure wear following total hip arthroplasty.

Hip Int 2013233233–242

Carmignato

Spinelli

Affatato

Savio

Uncertainty evaluation of volumetric wear assessment from coordinate measurements of ceramic hip joint prostheses.

Wear 20112709-10 584–590

Affatato

Spinelli

Zavalloni

Mazzega-Fabbro

Viceconti

Tribology and total hip joint replacement: current concepts in mechanical simulation.

Med Eng Phys 200830101305–1317

Bills

Blunt

Jiang

Development of a technique for accurately determining clinical wear in explanted total hip replacements.

Wear 20072637-121133–1137

Bowden

Kurtz

Edidin

Validation of a micro-CT technique for measuring volumetric wear in retrieved acetabular liners.

J Biomed Mater Res B Appl Biomater 2005751205–209

Carmignato

Savio

Traceable volume measurements using coordinate measuring systems.

CIRP 2011601519–522

10.

Carmignato

Dreossi

Mancini

Marinello

Tromba

Savio

Testing of x-ray microtomography systems using a traceable geometrical standard.

Meas Sci Technol 2009 20 8 08421

11.

Tan

Kiekens

Kruth

Voet

Dewulf

Material dependent thresholding for dimensional X-ray computed tomography. International symposium on digital industrial radiology and computed tomography.

June 20-22, Berlin, Germany 2011. http://www.dir2011.com/Proceedings/Mo4-Computed-Tomography-Applications: 1-8.

12.

Parrilli

Pagani

Maltarello

et al.

Three-dimensional cellular distribution in polymeric scaffolds for bone regeneration: a microCT analysis compared to SEM, CLSM and DNA content.

J Microsc 2014255120–29

13.

Gandolfi

Parrilli

Fini

Prati

Dummer

3D micro-CT analysis of the interface voids associated with Thermafil root fillings used with AH Plus or a flowable MTA sealer.

Int Endod J 2013463253–263

14.

Kruth

Bartscher

Carmignato

Schmitt

De Chiffre

Weckenmann

Computed tomography for dimensional metrology.

CIRP Annals: Manufacturing Technology 201160821–842

15.

Stock

MicroComputed tomography: methodology and applications.

Boca Raton, FL:

CRC Press

2011.

16.

Meng

Zhang

Wang

Virtual experimentation of beam hardening effect in X-ray CT measurement of multiphase flow.

Powder Technol 20091941-2153–157

17.

Krumm

Kasperl

Franz

Reducing non-linear artifacts of multi-material objects in industrial 3D computed tomography.

NDT Int 2008414242–251

18.

Boas

Fleischmann

Evaluation of two iterative techniques for reducing metal artifacts in computed tomography.

Radiology 20112593894–902

19.

Brown

Zabic

Method for reducing windmill artifacts in multi-slice CT images.

Proc SPIE 201179611–5

20.

Stradiotti

Curti

Castellazzi

Zerbi

Metal-related artifacts in instrumented spine.

Techniques for reducing artifacts in CT and MRI: state of the art. Eur Spine J 200918Suppl 1 102–108

21.

Boas

Fleischmann

CT artifacts: causes and reduction techniques.

Imaging in Medicine 201242229–240

22.

Yagüe-Fabra

Ontiveros

Jiménez

et al.

A 3D edge detection technique for surface extraction in computed tomography for dimensional metrology applications.

CIRP Annals: Manufacturing Technology 201362531–534

23.

Whyms

Vorperian

Gentry

Schimek

Bersu

Chung

The effect of computed tomographic scanner parameters and 3-dimensional volume rendering techniques on the accuracy of linear, angular, and volumetric measurements of the mandible.

Oral Surgery, Oral Medicine, Oral Pathology and Oral Radiology 20131155682–691

24.

Otsu

A threshold selection method from gray-level histograms.

IEEE Trans Syst Man Cybern 19799162–66

25.

Armitage

Berry

Mathews

Statistical methods in medical research.

4^th ed. Oxford: Blackwell Science 2001.

26.

Bland

Altman

Statistical methods in medical research.

Stat Methods Med Res 19998161–179

27.

Teeter

Naudie

Milner

Holdsworth

Determination of reference geometry for polyethylene tibial insert wear analysis.

J Arthroplasty 2011263497–503

28.

Teeter

Naudie

McErlain

et al.

In vitro quantification of wear in tibial inserts using microcomputed tomography.

Clin Orthop Relat Res 20114691107–112

29.

Hannouche

Hamadouche

Nizard

Bizot

Meunier

Sedel

Ceramics in total hip replacement [review].

Clin Orthop Relat Res 200543062–71

30.

Clarke

Good

Williams

et al.

Ultra-low wear rates for rigid-on-rigid bearings in total hip replacements.

Proc Inst Mech Eng H 20002144331–347

31.

Kiekens

Tan

Kruth

Voet

Dewulf

Parameter dependent thresholding for dimensional X-ray computed tomography [poster].

International symposium on digital industrial radiology and computed tomography June 20-22, Berlin, Germany 2011: Poster 12.