Sources of error in bone mineral density estimates from quantitative CT

Age-dependent changes in CT vertebral attenuation values in opportunistic screening for osteoporosis: a nationwide multi-center study

OBJECTIVES

To examine how vertebral attenuation changes with aging, and to establish age-adjusted CT attenuation value cutoffs for diagnosing osteoporosis.

MATERIALS AND METHODS

This multi-center retrospective study included 11,246 patients (mean age ± standard deviation, 50 ± 13 years; 7139 men) who underwent CT and dual-energy X-ray absorptiometry (DXA) in six health-screening centers between 2022 and 2023. Using deep-learning-based software, attenuation values of L1 vertebral bodies were measured. Segmented linear regression in women and simple linear regression in men were used to assess how attenuation values change with aging. A multivariable linear regression analysis was performed to determine whether age is associated with CT attenuation values independently of the DXA T-score. Age-adjusted cutoffs targeting either 90% sensitivity or 90% specificity were derived using quantile regression. Performance of both age-adjusted and age-unadjusted cutoffs was measured, where the target sensitivity or specificity was considered achieved if a 95% confidence interval encompassed 90%.

RESULTS

While attenuation values declined consistently with age in men, they declined abruptly in women aged > 42 years. Such decline occurred independently of the DXA T-score (p < 0.001). Age adjustment seemed critical for age ≥ 65 years, where the age-adjusted cutoffs achieved the target (sensitivity of 91.5% (86.3-95.2%) when targeting 90% sensitivity and specificity of 90.0% (88.3-91.6%) when targeting 90% specificity), but age-unadjusted cutoffs did not (95.5% (91.2-98.0%) and 73.8% (71.4-76.1%), respectively).

CONCLUSION

Age-adjusted cutoffs provided a more reliable diagnosis of osteoporosis than age-unadjusted cutoffs since vertebral attenuation values decrease with age, regardless of DXA T-scores.

KEY POINTS

Question How does vertebral CT attenuation change with age? Findings Independent of dual-energy X-ray absorptiometry T-score, vertebral attenuation values on CT declined at a constant rate in men and abruptly in women over 42 years of age. Clinical relevance Age adjustments are needed in opportunistic osteoporosis screening, especially among the elderly.

Regions of interest in opportunistic computed tomography-based screening for osteoporosis: impact on short-term in vivo precision

OBJECTIVE

To determine an optimal region of interest (ROI) for opportunistic screening of osteoporosis in terms of short-term in vivo diagnostic precision.

MATERIALS AND METHODS

We included patients who underwent two CT scans and one dual-energy X-ray absorptiometry scan within a month in 2022. Deep-learning software automatically measured the attenuation in L1 using 54 ROIs (three slice thicknesses × six shapes × three intravertebral levels). To identify factors associated with a lower attenuation difference between the two CT scans, mixed-effect model analysis was performed with ROI-level (slice thickness, shape, intravertebral levels) and patient-level (age, sex, patient diameter, change in CT machine) factors. The root-mean-square standard deviation (RMSSD) and area under the receiver-operating-characteristic curve (AUROC) were calculated.

RESULTS

In total, 73 consecutive patients (mean age ± standard deviation, 69 ± 9 years, 38 women) were included. A lower attenuation difference was observed in ROIs in images with slice thicknesses of 1 and 3 mm than that in images with a slice thickness of 5 mm (p < .001), in large elliptical ROIs (p = .007 or < .001, respectively), and in mid- or cranial-level ROIs than that in caudal-level ROIs (p < .001). No patient-level factors were significantly associated with the attenuation difference. Large, elliptical ROIs placed at the mid-level of L1 on images with 1- or 3-mm slice thicknesses yielded RMSSDs of 12.4-12.5 HU and AUROCs of 0.90.

CONCLUSION

The largest possible regions of interest drawn in the mid-level trabecular portion of the L1 vertebra on thin-slice images may yield improvements in the precision of opportunistic screening for osteoporosis via CT.

Misinterpretations about CT numbers, material decomposition, and elemental quantification

BACKGROUND

Quantitative CT imaging, particularly iodine and calcium quantification, is an important CT-based biomarker.

PURPOSE

This study quantifies sources of errors in quantitative CT imaging in both single-energy and spectral CT.

MATERIALS AND METHODS

This work examines the theoretical relationship between CT numbers, linear attenuation coefficient, and material quantification. We derive four understandings: (1) CT numbers are not proportional with element mass in vivo, (2) CT numbers are proportional with element mass only when contained in a voxel of pure water, (3) iodine-water material decomposition is never accurate in vivo, and (4) for error-free material decomposition a voxel must only consist of the basis decomposition vectors. Misinterpretation-based errors are calculated using the National Institute of Standards and Technology (NIST) XCOM database for: tissue chemical compositions, clinical concentrations of hydroxyapatite (HAP), and iodine. Quantification errors are also demonstrated experimentally using phantoms.

RESULTS

In single-energy CT, misinterpretation-induced errors for HAP density in adipose, muscle, lung, soft tissue, and blood ranged from 0-132%, i.e., a mass error of 0-749 mg/cm3. In spectral CT, errors with iodine in the same tissues resulted in a range of < 0.1-33% error, resulting in a mass error of < 0.1-1.2 mg/mL.

CONCLUSION

Our work demonstrates material quantification is fundamentally limited when measured in vivo due to measurement conditions differing from assumed and the errors are at or above detection limits for bone mineral density (BMD) and spectral iodine quantification. To define CT-derived biomarkers, the errors we demonstrate should either be avoided or built into uncertainty bounds.

CLINICAL RELEVANCE STATEMENT

Improving error bounds in quantitative CT biomarkers, specifically in iodine and BMD quantification, could lead to improvements in clinical care aspects based on quantitative CT.

KEY POINTS

CT numbers are only proportional with element mass only when contained in a voxel of pure water, therefore iodine-water material decomposition is never accurate in vivo. Misinterpretation-induced errors ranged from 0-132% for HAP density and < 0.1-33% in spectral CT with iodine. For error-free material decomposition, a voxel must only consist of the basis decomposition vectors.

Asynchronous calibration of a CT scanner for bone mineral density estimation: sources of error and correction

The estimation of BMD with CT scans requires a calibration method, usually based on a phantom. In asynchronous calibration, the phantom is scanned separately from the patient. A standardized acquisition protocol must be used to avoid variations between patient and phantom. However, variations can still be induced, for example, by temporal fluctuations or patient characteristics. Based on the further use of 739 forensic and 111 clinical CT scans, this study uses the proximal femur BMD value ("total hip") to assess asynchronous calibration accuracy, using in-scan calibration as ground truth. It identifies the parameters affecting the calibration accuracy and quantifies their impact. For time interval and table height, the impact was measured by calibrating the CT scan twice (once using the phantom scan with closest acquisition parameters and once using a phantom scan with standard values) and comparing the calibration accuracy. For other parameters such as body weight, the impact was measured by computing a linear regression between parameter values and calibration accuracy. Finally, this study proposes correction methods to reduce the effect of these parameters and improve the calibration accuracy. The BMD error of the asynchronous calibration, using the phantom scan with the closest acquisition parameters, was -1.2 ± 1.7% for the forensic and - 1.6 ± 3.5% for the clinical dataset. Among the parameters studied, time interval and body weight were identified as the main sources of error for asynchronous calibration, followed by table height and reconstruction kernel. Based on these results, a correction method was proposed to improve the calibration accuracy.

Comparison of Hounsfield Units within the Humeral Trochlea and Medial Coronoid Process in a Population of Labrador X Golden Retriever Guide Dogs and Border Collies

OBJECTIVE

 The aim of this study was to determine if variations in Hounsfield units (HU) are present within the elbow between asymptomatic dogs of two breeds of dogs susceptible to elbow dysplasia.

STUDY DESIGN

 Guide Dogs and Border Collies that presented for routine computed tomography (CT) screening for elbow dysplasia prior to breeding were evaluated. All dogs had no documented history of lameness. Dogs diagnosed with CT as being free of elbow dysplasia were included. The CT images were randomized and assessed by three blinded observers. A standardised approach to CT image reconstruction to create consistent image planes was used. Hounsfield units were measured within a standardised region of interest (ROI) at the humeral trochlea and medial coronoid process. The minimum, mean and maximum HU within each ROI was recorded.

RESULTS

 Eighty-six elbows were included in the study with 32 Guide Dogs, and 11 Border Collies. Guide Dogs had significantly higher minimum (99.75 HU, 95% confidence interval [CI]: 15.02-184.48, p = 0.022), mean (115.09 HU, 95% CI: 80.53-149.64, p < 0.01) and maximum (74.00 HU, 95% CI: 44.58-103.42, p < 0.01) difference in HU within the medial coronoid process ROI, and significantly higher mean (146.49 HU, 95% CI: 100.12-192.87, p < 0.01) and maximum (147.77 HU, 95% CI: 102.57-192.97, p < 0.01) difference in HU within the humeral trochlea ROI.

CONCLUSION

 In this dataset breed variations in elbow HU were present between asymptomatic Guide Dogs and Border Collies. This needs to be considered in breeding screening programmes to avoid over-interpretation of elbow sclerosis, in the absence of elbow pathology.

Phantomless calibration of CT scans for hip fracture risk prediction in silico: Comparison with phantom-based calibration

Finite element models built from quantitative computed tomography images rely on element-wise mapping of material properties starting from Hounsfield Units (HU), which can be converted into mineral densities upon calibration. While calibration is preferably carried out by scanning a phantom with known-density components, conducting phantom-based calibration may not always be possible. In such cases, a phantomless procedure, where the scanned subject's tissues are used as a phantom, is an interesting alternative. The aim of this study was to compare a phantom-based and a phantomless calibration method on 41 postmenopausal women. The proposed phantomless calibration utilized air, adipose, and muscle tissues, with reference equivalent mineral density values of -797, -95, and 38 mg/cm3, extracted from a previously performed phantom-based calibration. A 9-slice volume of interest (VOI) centred between the femoral head and knee rotation centres was chosen. Reference HU values for air, adipose, and muscle tissues were extracted by identifying HU distribution peaks within the VOI, and patient-specific calibration was performed using linear regression. Comparison of FE models calibrated with the two methods showed average relative differences of 1.99% for Young's modulus1.30% for tensile and 1.34% for compressive principal strains. Excellent correlations (R2 > 0.99) were identified for superficial maximum tensile and minimum compressive strains. Maximum normalised root mean square relative error (RMSRE) values settled at 4.02% for Young's modulus, 2.99% for tensile, and 3.22% for compressive principal strains, respectively. The good agreement found between the two methods supports the adoption of the proposed methodology when phantomless calibration is needed.

Coronal plane in opportunistic screening of osteoporosis using computed tomography: comparison with axial and sagittal planes

OBJECTIVE

To compare the coronal plane with axial and sagittal planes in opportunistic screening of osteoporosis using computed tomography (CT).

MATERIALS AND METHODS

A total of 100 patients aged ≥ 50 years who underwent both lumbar spine CT and dual-energy X-ray absorptiometry within 3 months were included. Osteoporosis was diagnosed based on dual-energy X-ray absorptiometry results. The CT number was measured at the center of the vertebral body in coronal, axial, and sagittal planes. To compare the coronal plane with axial and sagittal planes in diagnosing osteoporosis, the areas under the receiver operating characteristic curve (AUC) were compared and intraclass correlation coefficient (ICC) was calculated. The optimal cutoff values were calculated using Youden's index.

RESULTS

The AUC of the coronal plane (0.80; 95% confidence interval [CI], 0.71-0.89) was not significantly different from that of the axial plane (0.78; 95% CI, 0.68-0.87; P = 0.39) and that of the sagittal plane (0.78; 95% CI, 0.69-0.87; P = 0.68). Excellent concordance rates were observed between coronal and axial planes with ICC of 0.95 (95% CI, 0.92-0.96) and between coronal and sagittal planes with ICC of 0.93 (95% CI, 0.85-0.96). The optimal cutoff values for the coronal, axial, and sagittal planes were 110, 112, and 112 HU, respectively.

CONCLUSION

The coronal plane does not significantly differ from axial and sagittal planes in opportunistic screening of osteoporosis. Thus, the coronal plane as well as axial and sagittal planes can be used interchangeably in measuring bone mineral density using CT.

Differential effect of atorvastatin and pravastatin on thoracic spine attenuation: A sub-analysis of a randomized clinical trial

BACKGROUND

Statins reduce cardiovascular events and may improve bone mineral density.

METHODS

We conducted a sub-analysis of a randomized clinical trial that investigated the differential effect of moderate vs intensive low-density lipoprotein cholesterol (LDL-C) lowering therapies on coronary artery calcium (CAC) scores, and used the acquired images to assess the change in radiological attenuation of selected thoracic vertebrae. Baseline and 12-month unenhanced chest CT scans were performed in 420 hyperlipidemic, postmenopausal women randomized to atorvastatin (ATV) 80 mg/day or pravastatin (PRV) 40 mg/day in the Beyond Endorsed Lipid Lowering with Electron Beam Tomography Scanning (BELLES) trial. Bone attenuation was measured in three contiguous thoracic vertebrae at baseline and 12 months.

RESULTS

There were no differences in baseline demographic and clinical characteristics between treatment arms. The median percent lowering (interquartile range) in LDL-C was significantly greater with ATV than PRV [-53 (-69 to 20)% vs -28 (-55 to 74)%, p < 0.001], although the CAC score change was similar [12 (-63 to 208)% vs 13 (-75 to 358)%; p = 0.44]. At follow-up, the median bone attenuation loss was significantly greater with PRV than with ATV [-2.6 (-27 to 11)% vs 0 (-11 to 25)%; p < 0.001]. The attenuation loss in the PRV group was comparable to that of a historical untreated general population sample. In the entire cohort, the changes in LDL-C and total cholesterol were inversely correlated with bone attenuation change (p < 0.01). In adjusted multivariable linear regression analyses, race and percent change in LDL-C were independent predictors of bone attenuation change. Age, body mass index, history of smoking, diabetes mellitus, hypertension, peripheral vascular disease, or hormone replacement therapy did not affect percent change in BMD.

CONCLUSIONS

These findings support the hypothesis that there is an interaction between bone and cardiometabolic health and that intensive lipid lowering has a beneficial effect on bone health.

Bone marrow metabolism is affected by body weight and response to exercise training varies according to anatomical location

AIM

High body weight is a protective factor against osteoporosis, but obesity also suppresses bone metabolism and whole-body insulin sensitivity. However, the impact of body weight and regular training on bone marrow (BM) glucose metabolism is unclear. We studied the effects of regular exercise training on bone and BM metabolism in monozygotic twin pairs discordant for body weight.

METHODS

We recruited 12 monozygotic twin pairs (mean ± SD age 40.4 ± 4.5 years; body mass index 32.9 ± 7.6, mean difference between co-twins 7.6 kg/m2 ; eight female pairs). Ten pairs completed the 6-month long training intervention. We measured lumbar vertebral and femoral BM insulin-stimulated glucose uptake (GU) using 18 F-FDG positron emission tomography, lumbar spine bone mineral density and bone turnover markers.

RESULTS

At baseline, heavier co-twins had higher lumbar vertebral BM GU (p < .001) and lower bone turnover markers (all p < .01) compared with leaner co-twins but there was no significant difference in femoral BM GU, or bone mineral density. Training improved whole-body insulin sensitivity, aerobic capacity (both p < .05) and femoral BM GU (p = .008). The training response in lumbar vertebral BM GU was different between the groups (time × group, p = .02), as GU tended to decrease in heavier co-twins (p = .06) while there was no change in leaner co-twins.

CONCLUSIONS

In this study, regular exercise training increases femoral BM GU regardless of weight and genetics. Interestingly, lumbar vertebral BM GU is higher in participants with higher body weight, and training counteracts this effect in heavier co-twins even without reduction in weight. These data suggest that BM metabolism is altered by physical activity.

The correlation between osteoporotic vertebrae fracture risk and bone mineral density measured by quantitative computed tomography and dual energy X-ray absorptiometry: a systematic review and meta-analysis

This paper presents a comparison of quantitative computed tomography (QCT) and dual-energy X-ray absorptiometry (DXA) in osteoporosis with vertebral fracture and osteoporosis without fracture. It has been proved that the volumetric bone mineral density (vBMD) measured by QCT exhibits a stronger correlation with fracture risk than areal bone mineral density (aBMD) measured by DXA.

PURPOSE

This study aims to systematically evaluate the ability of QCT and DXA to distinguish between osteoporosis with vertebral fracture and osteoporosis without fracture according to vBMD and aBMD.

METHODS

We conducted a primary literature search of the online databases up to 3 July, 2022, in both English and Chinese publications, combining synonyms for "QCT", "DXA" and "osteoporosis". The Newcastle-Ottawa scale (NOS) was employed to evaluate the quality of the selected articles. vBMD obtained through QCT and aBMD obtained through DXA were extracted, and were analyzed by Review Manager 5.4 and RStudio.

RESULTS

Six studies with 610 individuals aged 45 to 90, of which 179 had vertebral fractures, were included in the final analysis. The weighted mean difference (WMD) between osteoporosis with vertebral fracture and osteoporosis without fracture for vBMD was - 27.08 (95% CI  - 31.24 to  - 22.92), while for aBMD was - 0.05 (95% CI - 0.08 to - 0.03).

CONCLUSIONS

Both vBMD detected by QCT and aBMD detected by DXA could discriminate fracture status in the spine, and vBMD performed a stronger correlation with fracture risk.

TRIAL REGISTRATION

PROSPERO 2022 CRD42022349185.

Hidden blood loss and its risk factors in percutaneous vertebroplasty surgery for osteoporotic vertebral compression fractures

OBJECTIVES

Percutaneous vertebroplasty (PVP) is a percutaneous interventional procedure for osteoporotic vertebral compression fractures (OVCFs). However, hidden blood loss (HBL) during the surgery is easily disregarded. This study aimed to evaluate HBL and its possible risk factors in the patients following PVP for OVCFs.

METHODS

Patients with OVCFs who underwent PVP surgery between January 2019 and November 2022 at our hospital were retrospectively analyzed. Patients' demographics, laboratory data, and imaging and clinical date were also collected. Preoperative and postoperative hematocrit were recorded, the hidden blood loss was calculated according to Sehat formula, and the risk factors were analyzed by multivariate linear regression analysis.

RESULTS

One hundred and fifty-five patients (26 males and 129 females) were retrospectively enrolled in this study. 85.2% of patients had one segment vertebral fracture and the mean surgical time was 30.5 ± 11.0 min. No intraspinal cement leakage occurred. The mean HBL was 204.0 ± 89.6 ml. Multivariate linear regression analysis revealed that HBL was positively associated with number of fracture segments (P = 0.001), degree of vertebral height restoration (P = 0.001), surgical time (P = 0.000), number of puncture (P = 0.002), and cement leakage (P = 0.038).

CONCLUSIONS

Multiple vertebral fractures, higher degree of vertebral height restoration, longer surgical time, more number of puncture, and cement leakage are independent risk factors for HBL. Therefore, HBL should not be neglected in the patients with OVCFs undergoing PVP surgery, especially in those with poor preoperative physical condition and presence of anemia.

Feasibility analysis of bone density evaluation with Hounsfield unit value after fibula flap reconstruction of jaw defect

BACKGROUND

Implant-supported dentures have become an essential means of restoring occlusal function after jaw reconstruction. Bone mineral density (BMD) may influence the success rate of implant denture restorations. This study aimed to explore whether the Hounsfield unit (HU) value can be used to monitor the changing trend of fibular BMD after jaw reconstruction.

RESULTS

A total of 54 patients who underwent maxillar/mandibular reconstruction with a fibula flap were included in this study. There was a significant correlation between the HU value and BMD at 1 week, 3 months, and 6 months after surgery, and both were significantly correlated with follow-up time. The difference between each pair of absorption rates (DAR) was less than 10% in 66.7% and 75.9% of patients at 3 and 6 months; however, the DAR was more than 20% in 12% and 13.8% of patients at 3 and 6 months, respectively.

CONCLUSIONS

There is a significant correlation between HU value and BMD. The HU value can be used to roughly reflect the fibular BMD changing trend in a group of patients as opposed to an individual, and the HU value is not equivalent to BMD.

TRIAL REGISTRATION

ChiCTR, ChiCTR2300069661, retrospectively registered on 22 March 2023. Retrospectively registered, https://www.chictr.org.cn/showproj.html?proj=188953 .

Inflammation, bone loss and 2-year bone formation at the same vertebra in axial spondyloarthritis: a multilevel MRI and low-dose CT analysis

OBJECTIVE

To investigate whether in radiographic axial spondyloarthritis (r-axSpA) inflammation is associated with lower trabecular bone density (TBD), and subsequently, if a lower TBD increases the likelihood of 2-year bone formation at the same vertebra.

METHODS

Whole spine (C3-L5) data from patients included in the multicentre 2-year Sensitive Imaging in Ankylosing Spondylitis cohort was used. Two readers measured baseline TBD by Hounsfield units (HU) on low-dose CT (ldCT). Baseline MRI bone marrow oedema (BME) status scores and ldCT syndesmophyte formation and/or growth change-from-baseline scores were assessed by three and two readers, respectively. Average of readers' continuous measurements or readers' agreement in binary scores generated within the same vertebra (1-present in ≥1 quadrant/0-absent in all quadrants) were used. Multilevel generalised estimating equations models were used, the unit of analysis being the vertebra.

RESULTS

In 50 patients with r-axSpA, TBD HU decreased from cranial to caudal vertebrae. Baseline MRI-BME was present in 300/985 (30%) and syndesmophytes in 588/910 (65%) vertebrae, both most prevalent at thoracolumbar region. Syndesmophyte formation or growth was observed in 18% of at-risk vertebrae (124/691). A significant confounder-adjusted association was found between inflammation and lower TBD (regression coefficient=-51; 95% CI-63 to -39). TBD was not associated with 2-year syndesmophyte formation or growth (adjusted OR 1.00; 95% CI 0.99 to 1.00).

CONCLUSION

In r-axSpA, while vertebral inflammation was associated with lower vertebral TBD, lower vertebral TBD itself did not increase the risk for new bone formation at the same vertebra. In preventing syndesmophyte progression, targeting local inflammation seems more important than targeting vertebral trabecular bone loss.

Asynchronous calibration of quantitative computed tomography bone mineral density assessment for opportunistic osteoporosis screening: phantom-based validation and parameter influence evaluation

Asynchronous calibration could allow opportunistic screening based on routine CT for early osteoporosis detection. In this phantom study, a bone mineral density (BMD) calibration phantom and multi-energy CT (MECT) phantom were imaged on eight different CT scanners with multiple tube voltages (80-150 kV_p) and image reconstruction settings (e.g. soft/hard kernel). Reference values for asynchronous BMD estimation were calculated from the BMD-phantom and validated with six calcium composite inserts of the MECT-phantom with known ground truth. Relative errors/changes in estimated BMD were calculated and investigated for influence of tube voltage, CT scanner and reconstruction setting. Reference values for 282 acquisitions were determined, resulting in an average relative error between calculated BMD and ground truth of - 9.2% ± 14.0% with a strong correlation (R² = 0.99; p < 0.0001). Tube voltage and CT scanner had a significant effect on calculated BMD (p < 0.0001), with relative differences in BMD of 3.8% ± 28.2% when adapting reference values for tube voltage, - 5.6% ± 9.2% for CT scanner and 0.2% ± 0.2% for reconstruction setting, respectively. Differences in BMD were small when using reference values from a different CT scanner of the same model (0.0% ± 1.4%). Asynchronous phantom-based calibration is feasible for opportunistic BMD assessment based on CT images with reference values adapted for tube voltage and CT scanner model.

Low-dose CT hounsfield units: a reliable methodology for assessing vertebral bone density in radiographic axial spondyloarthritis

Objective

Studying vertebral bone loss in radiographic axial spondyloarthritis (r-axSpA) has been challenging due to ectopic bone formation. We cross-sectionally analysed low-dose CT (ldCT) trabecular bone density Hounsfield units (HU) measurements and calculated inter-reader reliability at the vertebral level in patients with r-axSpA.

Methods

LdCT scans of 50 patients with r-axSpA from the sensitive imaging in ankylosing spondylitis study, a multicentre 2-year prospective cohort were included. Trabecular bone HU taken from a region of interest at the centre of each vertebra (C3-L5) were independently assessed by two trained readers. HU mean (SD), and range were provided at the vertebral level, for each reader and centre separately. Inter-reader reliability and agreement were assessed using intraclass correlation coefficients (ICC; single measurements, absolute agreement, two-way mixed effects models); smallest detectable difference and Bland-Altman plots.

Results

Overall, 1100 vertebrae were assessed by each reader. HU values decreased from cranial to caudal vertebrae. For readers 1 and 2 respectively, the highest mean (SD) HU value was obtained at C3 (354(106) and 355(108)), and the lowest at L3 (153(65) and 150 (65)). Inter-reader reliability was excellent (ICC(2,1):0.89 to 1.00). SDD varied from 4 to 8. For most vertebrae, reader 1 scored somewhat higher than reader 2 (mean difference of scores ranging from −0.6 to 2.9 HU). Bland-Altman plots showed homoscedasticity.

Conclusion

LdCT measurement of HU is a feasible method to assess vertebral bone density in r-axSpA with excellent inter-reader reliability from C3 to L5. These results warrant further validation and longitudinal assessment of reliability.