Use of internal references for assessing CT density measurements of the pelvis as replacement for use of an external phantom

Internal calibration for opportunistic computed tomography muscle density analysis

INTRODUCTION

Muscle weakness can lead to reduced physical function and quality of life. Computed tomography (CT) can be used to assess muscle health through measures of muscle cross-sectional area and density loss associated with fat infiltration. However, there are limited opportunities to measure muscle density in clinically acquired CT scans because a density calibration phantom, allowing for the conversion of CT Hounsfield units into density, is typically not included within the field-of-view. For bone density analysis, internal density calibration methods use regions of interest within the scan field-of-view to derive the relationship between Hounsfield units and bone density, but these methods have yet to be adapted for muscle density analysis. The objective of this study was to design and validate a CT internal calibration method for muscle density analysis.

METHODOLOGY

We CT scanned 10 bovine muscle samples using two scan protocols and five scan positions within the scanner bore. The scans were calibrated using internal calibration and a reference phantom. We tested combinations of internal calibration regions of interest (e.g., air, blood, bone, muscle, adipose).

RESULTS

We found that the internal calibration method using two regions of interest, air and adipose or blood, yielded accurate muscle density values (< 1% error) when compared with the reference phantom. The muscle density values derived from the internal and reference phantom calibration methods were highly correlated (R2 > 0.99). The coefficient of variation for muscle density across two scan protocols and five scan positions was significantly lower for internal calibration (mean = 0.33%) than for Hounsfield units (mean = 6.52%). There was no difference between coefficient of variation for the internal calibration and reference phantom methods.

CONCLUSIONS

We have developed an internal calibration method to produce accurate and reliable muscle density measures from opportunistic computed tomography images without the need for calibration phantoms.

The effect of variations in CT scan protocol on femoral finite element failure load assessment using phantomless calibration

Recently, it was shown that fracture risk assessment in patients with femoral bone metastases using Finite Element (FE) modeling can be performed using a calibration phantom or air-fat-muscle calibration and that non-patient-specific calibration was less favorable. The purpose of this study was to investigate if phantomless calibration can be used instead of phantom calibration when different CT protocols are used. Differences in effect of CT protocols on Hounsfield units (HU), calculated bone mineral density (BMD) and FE failure loads between phantom and two methods of phantomless calibrations were studied. Five human cadaver lower limbs were scanned atop a calibration phantom according to a standard scanning protocol and seven additional commonly deviating protocols including current, peak kilovoltage (kVp), slice thickness, rotation time, field of view, reconstruction kernel, and reconstruction algorithm. The HUs of the scans were calibrated to BMD (in mg/cm³) using the calibration phantom as well as using air-fat-muscle and non-patient-specific calibration, resulting in three models for each scan. FE models were created, and failure loads were calculated by simulating an axial load on the femur. HU, calculated BMD and failure load of all protocols were compared between the three calibration methods. The different protocols showed little variation in HU, BMD and failure load. However, compared to phantom calibration, changing the kVp resulted in a relatively large decrease of approximately 10% in mean HU and BMD of the trabecular and cortical region of interest (ROI), resulting in a 13.8% and 13.4% lower failure load when air-fat-muscle and non-patient-specific calibrations were used, respectively. In conclusion, while we observed significant correlations between air-fat-muscle calibration and phantom calibration as well as between non-patient-specific calibration and phantom calibration, our sample size was too small to prove that either of these calibration approaches was superior. Further studies are necessary to test whether air-fat-muscle or non-patient-specific calibration could replace phantom calibration in case of different scanning protocols.

Monetite-based composite cranial implants demonstrate long-term clinical volumetric balance by concomitant bone formation and degradation

The use of calcium phosphates (CaPs) as synthetic bone substitutes should ideally result in a volumetric balance with concomitant bone formation and degradation. Clinical data on such properties is nevertheless lacking, especially for monetite-based CaPs. However, a monetite-based composite implant has recently shown promising cranial reconstructions, with both CaP degradation and bone formation. In this study, the volumetric change at the implant site was quantified longitudinally by clinical computed tomography (CT). The retrospective CT datasets had been acquired postoperatively (n = 10), in 1-year (n = 9) and 3-year (n = 5) follow-ups. In the 1-year follow-up, the total volumetric change at the implant site was -8 ± 8%. A volumetric increase (bone formation) was found in the implant-bone interface, and a volumetric decrease was observed in the central region (CaP degradation). In the subjects with 2- or 3-year follow-ups, the rate of volumetric decrease slowed down or plateaued. The reported degradation rate is lower than previous clinical studies on monetite, likely due to the presence of pyrophosphate in the monetite-based CaP-formulation. A 31-months retrieval specimen analysis demonstrated that parts of the CaP had been remodeled into bone. The CaP phase composition remained stable, with 6% transformation into hydroxyapatite. In conclusion, this study demonstrates successful bone-bonding between the CaP-material and the recipient bone, as well as a long-term volumetric balance in cranial defects repaired with the monetite-based composite implant, which motivates further clinical use. The developed methods could be used in future studies for correlating spatiotemporal information regarding bone regeneration and CaP degradation to e.g. patient demographics. STATEMENT OF SIGNIFICANCE: In bone defect reconstructions, the use of calcium phosphate (CaP) bioceramics ideally results in a volumetric balance between bone formation and CaP degradation. Clinical data on the volumetric balance is nevertheless lacking, especially for monetite-based CaPs. Here, this concept is investigated for a composite cranial implant. The implant volumes were quantified from clinical CT-data: postoperatively, one year and three years postoperatively. In total, -8 ± 8% (n = 9) volumetric change was observed after one year. But the change plateaued, with only 2% additional decrease at the 3-year follow-up (n = 5), indicating a lower CaP degradation rate. Osseointegration was seen at the bone-implant interface, with a 9 ± 7% volumetric change after one year. This study presented the first quantitative spatiotemporal CT analysis of monetite-based CaPs.

Evaluation of patient tissue selection methods for deriving equivalent density calibration for femoral bone quantitative CT analyses

Osteoporosis affects an increasing number of people every year and patient specific finite element analysis of the femur has been proposed to identify patients that could benefit from preventative treatment. The aim of this study was to demonstrate, verify, and validate an objective process for selecting tissues for use as the basis of phantomless calibration to enable patient specific finite element analysis derived hip fracture risk prediction. Retrospective reanalysis of patient computed tomography (CT) scans has the potential to yield insights into more accurate prediction of osteoporotic fracture. Bone mineral density (BMD) specific calibration scans are not typically captured during routine clinical practice. Tissue-based BMD calibration can therefore empower the retrospective study of patient CT scans captured during routine clinical practice. Together the method for selecting tissues as the basis for phantomless calibration coupled with the post-processing steps for deriving a calibration equation using the selected tissues provide an estimation of quantitative equivalent density results derived using calibration phantoms. Patient tissues from a retrospective cohort of 211 patients were evaluated. The best phantomless calibration resulted in a femoral strength (FS) [N] bias of 0.069 ± 0.07% over FS derived from inline calibration and a BMD [kg/cm³] bias of 0.038 ± 0.037% over BMD derived from inline calibration. The phantomless calibration slope for the best method presented was within the range of patient specific calibration curves available for comparison and demonstrated a small bias of 0.028 ± 0.054 HU/(mg/cm³), assuming the Mindways Model 3 BMD inline calibration phantom as the gold standard. The presented method of estimating a calibration equation from tissues showed promise for CT-based femoral fracture analyses of retrospective cohorts without readily available calibration data.

A Registration Method for Three-Dimensional Analysis of Bone Mineral Density in the Proximal Tibia

Although alterations in bone mineral density (BMD) at the proximal tibia have been suggested to play a role in various musculoskeletal conditions, their pathophysiological implications and their value as markers for diagnosis remain unclear. Improving our understanding of proximal tibial BMD requires novel tools for three-dimensional (3D) analysis of BMD distribution. Three-dimensional imaging is possible with computed tomography (CT), but computational anatomy algorithms are missing to standardize the quantification of 3D proximal tibial BMD, preventing distribution analyses. The objectives of this study were to develop and assess a registration method, suitable with routine knee CT scans, to allow the standardized quantification of 3D BMD distribution in the proximal tibia. Second, as an example of application, the study aimed to characterize the distribution of BMD below the tibial cartilages in healthy knees. A method was proposed to register both the surface (vertices) and the content (voxels) of proximal tibias. The method combines rigid transformations to account for differences in bone size and position in the scanner's field of view and to address inconsistencies in the portion of the tibial shaft included in routine CT scan, with a nonrigid transformation locally matching the proximal tibias. The method proved to be highly reproducible and provided a comprehensive description of the relationship between bone depth and BMD. Specifically it reported significantly higher BMD in the first 6 mm of bone than deeper in the proximal tibia. In conclusion, the proposed method offers promising possibilities to analyze BMD and other properties of the tibia in 3D.

Use of Bone Health Evaluation in Orthopedic Surgery: 2019 ISCD Official Position

This position development conference (PDC) Task Force examined the assessment of bone status in orthopedic surgery patients. Key questions included which orthopedic surgery patients should be evaluated for poor bone health prior to surgery and which subsets of patients are at high risk for poor bone health and adverse outcomes. Second, the reliability and validity of using bone densitometry techniques and measurement of specific geometries around the hip and knee before and after arthroplasty was determined. Finally, the use of computed tomography (CT) attenuation coefficients (Hounsfield units) to estimate bone quality at anatomic locations where orthopedic surgery is performed including femur, tibia, shoulder, wrist, and ankle were reviewed. The literature review identified 665 articles of which 198 met inclusion exclusion criteria and were selected based on reporting of methodology, reliability, or validity results. We recommend that the orthopedic surgeon be aware of established ISCD guidelines for determining who should have additional screening for osteoporosis. Patients with inflammatory arthritis, chronic corticosteroid use, chronic renal disease, and those with history of fracture after age 50 are at high risk of osteoporosis and adverse events from surgery and should have dual energy X-ray absorptiometry (DXA) screening before surgery. In addition to standard DXA, bone mineral density (BMD) measurement along the femur and proximal tibia is reliable and valid around implants and can provide valuable information regarding bone remodeling and identification of loosening. Attention to positioning, selection of regions of interest, and use of special techniques and software is required. Plain radiographs and CT provide simple, reliable methods to classify the shape of the proximal femur and to predict osteoporosis; these include the Dorr Classification, Cortical Index, and critical thickness. Correlation of these indices to central BMD is moderate to good. Many patients undergoing orthopedic surgery have had preoperative CT which can be utilized to assess regional quality of bone. The simplest method available on most picture archiving and communications systems is to simply measure a regions of interest and determine the mean Hounsfield units. This method has excellent reliability throughout the skeleton and has moderate correlation to DXA based on BMD. The prediction of outcome and correlation to mechanical strength of fixation of a screw or implant is unknown.

Calibration with or without phantom for fracture risk prediction in cancer patients with femoral bone metastases using CT-based finite element models

The objective of this study was to develop a new calibration method that enables calibration of Hounsfield units (HU) to bone mineral densities (BMD) without the use of a calibration phantom for fracture risk prediction of femurs with metastases using CT-based finite element (FE) models. Fifty-seven advanced cancer patients (67 femurs with bone metastases) were CT scanned atop a separate calibration phantom using a standardized protocol. Non-linear isotropic FE models were constructed based on the phantom calibration and on two phantomless calibration methods: the “air-fat-muscle” and “non-patient-specific” calibration. For air-fat-muscle calibration, peaks for air, fat and muscle tissue were extracted from a histogram of the HU in a standardized region of interest including the patient’s right leg and surrounding air. These CT peaks were linearly fitted to reference “BMD” values of the corresponding tissues to obtain a calibration function. For non-patient-specific calibration, an average phantom calibration function was used for all patients. FE failure loads were compared between phantom and phantomless calibrations. There were no differences in failure loads between phantom and air-fat-muscle calibration (p = 0.8), whereas there was a significant difference between phantom and non-patient-specific calibration (p<0.001). Although this study was not designed to investigate this, in four patients who were scanned using an aberrant reconstruction kernel, the effect of the different kernel seemed to be smaller for the air-fat-muscle calibration compared to the non-patient-specific calibration. With the air-fat-muscle calibration, clinical implementation of the FE model as tool for fracture risk assessment will be easier from a practical and financial viewpoint, since FE models can be made using everyday clinical CT scans without the need of concurrent scanning of calibration phantoms.

Quantitative regional and sub-regional analysis of femoral and tibial subchondral bone mineral density (sBMD) using computed tomography (CT): comparison of non-osteoarthritic (OA) and severe OA knees

OBJECTIVE

This study aimed to compare subchondral bone mineral density (sBMD) between non-radiographic osteoarthritic (OA) and medial femorotibial OA knees, using computed tomography (CT).

DESIGN

CT exams from 16 non-radiographic OA (KL grade < 2) and 16 severe medial OA (KL grade ≥ 3) knees (average age of 61.7 ± 3 and 62.2 ± 5 years old respectively, 50% male in each group), were retrospectively analyzed. CT exams were segmented and 3D maps of sBMD based on the CT number in the most superficial 3 mm of femoral and tibial subchondral bone were computed. Average sBMD and medial-to-lateral sBMD ratios were calculated for total load-bearing regions and for sub-regions of interest in the femur and tibia.

RESULTS

The analysis of total load-bearing regions did not reveal any significant difference between groups, except for the lateral tibia, where OA knees had lower sBMD. Sub-regional analysis unveiled differences with some sub-regions of the femur and tibia presenting significantly lower (in the lateral compartment) or higher (in the medial compartment) sBMD in OA knees compared to non-OA knees. The M/L sBMD ratios were significantly higher for OA knees compared to non-OA knees for all regions and sub-regions, except for the internal sub-regions.

CONCLUSIONS

sBMD locally differs between non-OA and OA knees, in agreement with prior knowledge on biomechanics. CT proved to be a valuable tool for 3D analysis of femoral and tibial sBMD, which can be used in future studies to describe the chronology of sBMD alterations and improve our understanding of the role of subchondral bone in knee OA.

CT-based quantification of bone stock in large head metal-on-metal unilateral total hip replacements

PURPOSE

To explore ipsilateral and contralateral acetabular roof bone stock density in unilateral large head MoM THA whether there is a significant lower acetabular bone stock in the hip with a metal-on-metal (MoM) total hip replacement compared to the contralateral side. Second part of this study is to examine if there are any associates with regard to potential bone stock density difference.

MATERIALS & METHODS

A database of 317 patients with unilateral metal-on-metal (MoM) total hip replacements was set up retrospectively for this study. On computed tomography scans, conducted after a relative short in situ time period averaging 2.8 years, regions-of-interests were drawn in the trabecular bone of the acetabulum to measure average Hounsfield Units (HU). HU differences were calculated and tested by Wilcoxon signed-rank test. Univariate analysis was conducted to examine associates of potential bone loss.

RESULTS

In a population of 317 patients (156 male, 161 female) with an average age of 61.9 ± 7.8, the median HU on the side of the MoM replacement was 123.3 (7.6-375.4). On the contralateral side, median HU was 144.7 (-0.4 to 332.8). The median HU difference was 21.4 after a mean post-operative in situ time of 2.8 years. The Wilcoxon signed-rank test proved a significant difference (p<0.001). Univariate analyses show that the in situ time of the MoM THA has a significant correlation with the bone density difference.

CONCLUSION

Results show a significant lower bone density at the acetabular roof at the side of the prosthesis compared with the contralateral side after short in situ time of the MoM THA in patients with unilateral MoM total hip replacements. In our patient population, the in situ time showed a significant association with the acetabular bone density difference. As acetabular roof bone stock measurements are feasible and show temporal decline this could become an important parameter to be used in orthopedic decision making for revision surgery.