Radiographic texture analysis of densitometer-generated calcaneus images differentiates postmenopausal women with and without fractures

Advanced 2D image processing technique to predict hip fracture risk in an older population based on single DXA scans

A new technique to enhance hip fracture risk prediction in older adults was presented and assessed. The new method dramatically improved prediction at high specificity levels using only a standard clinical diagnostic scan. This has the potential to be implemented in clinical practice to enhance patient fragility diagnosis.

INTRODUCTION

Diagnosis of osteoporosis is based on the measurement of bone mineral density (BMD) using dual-energy X-ray absorptiometry (DXA) scans. However, studies have shown this to be insufficient to accurately predict hip fractures. Therefore, complementary methods are needed to enhance hip fracture risk prediction to identify vulnerable patients.

METHODS

Hip DXA scans were obtained for 192 subjects from the Canadian Multicenter Osteoporosis Study (CaMos), 50 of whom had experienced a hip fracture within 5 years of the scan. 2D statistical shape and appearance modeling was performed to account for the effect of the femur's geometry and BMD distribution on hip fracture risk. Statistical shape modeling (SSM), and statistical appearance modeling (SAM) were also used separately to predict the fracture risk based solely on the femur's geometry and BMD distribution, respectively. Combined with BMD, age, and body mass index (BMI), logistic regression was performed to estimate the fracture risk over the 5-year period.

RESULTS

Using the new technique, hip fractures were correctly predicted in 78% of cases compared with 36% when using the T-score. The accuracy of the prediction was not greatly reduced when using SSM and SAM (78% and 74% correct, respectively). Various geometric and BMD distribution traits were identified in the fractured and non-fractured groups.

CONCLUSION

2D SSAM can dramatically improve hip fracture prediction at high specificity levels and estimate the year of the impending fracture using standard clinical images. This has the potential to be implemented in clinical practice to estimate hip fracture risk.

Ex Vivo Evaluation of Hip Fracture Risk by Proximal Femur Geometry and Bone Mineral Density in Elderly Chinese Women

Background

The incidence of hip fracture is steadily increasing. We aimed to establish a creative approach to precisely estimate the risk of hip fracture by exploring the relationship between hip fracture and bone mineral density (BMD)/femur geometry.

Material/Methods

Sixteen samples of cadaveric female proximal femora were randomly selected. Experiments were performed experimental measurement of the femoral neck BMD and geometric parameters (including neck length, neck diameter, head diameter, and neck-shaft angle). In addition, the experimental measurements contain the failure load, which represents the mechanical strength of the femoral neck, and we calculated the correlation coefficient among BMD, geometric parameters, and failure load.

Results

Significant correlations were discovered between femoral mechanical properties and femoral neck BMD (r=0.792, r²=0.628, P<0.001), trochanteric BMD (r=0.749, r²=0.560, P=0.001), and head diameter (r=0.706, r²=0.499, P=0.002). Multiple linear regression analyses indicated that the best predictor of hip fracture was the combination of femoral neck BMD, head diameter, and neck diameter (r²=0.844, P<0.001).

Conclusions

The results confirmed that, compared with BMD alone, the combination of BMD and geometric parameters of proximal femur is a better estimation of hip fracture. The geometry of the proximal femur played an important role in assessing the biomechanical strength of femur. This method greatly assists in predicting the risk of hip fracture in clinical trials and will assist studies on why the incidence of hip fracture varies among races.

Global sagittal alignment in elderly patients with osteoporosis and its relationship with severity of vertebral fracture and quality of life

We compared global sagittal alignment and quality of life in osteoporotic patients with and without vertebral compression fracture (VCF) and determined its relationship with VCF severity. The findings revealed osteoporotic patients with VCF showed decreased quality of life and worse global sagittal alignment, which was significantly associated with VCF severity.

INTRODUCTION

The aim of this study was to compare the global sagittal alignment and quality of life in elderly osteoporotic patients with and without vertebral compression fracture (VCF), and to investigate the relationship between global sagittal alignment and severity of VCF.

METHODS

A consecutive series of 72 female patients with osteoporosis aged over 60 years and 31 age-matched females without osteoporosis were prospectively enrolled. The patients were divided into VCF and non-VCF group. Patient's clinical demography, nature of VCF, and bone mineral density (BMD) were also recorded. Spinal deformity index was used to evaluate severity of VCF. EOS® biplanar imaging system was then used to evaluate global sagittal parameters: T1 pelvic angle (TPA) and global sagittal angle (GSA). In addition, quality of life was assessed with self-reported questionnaires: the Oswestry Disability Index (ODI) and Short-form 12 (SF-12).

RESULTS

Osteoporotic patients and controls were found to be significantly different in terms of TPA, GSA, and BMD. And in patients with VCF, they were found to have significantly higher TPA and GSA. TPA and GSA were significantly correlated with SF-12 and ODI. The number of VCF and SDI significantly correlated with global sagittal alignment. Using regression analysis, parameters significantly associated with abnormal global alignment were the number of VCF (OR = 1.13) and SDI (OR = 1.84).

CONCLUSION

Osteoporotic patients with VCF showed worse global sagittal alignment and decreased quality of life. The number and severity of VCF had a negative influence on global sagittal balance, which indicates that poorer sagittal global alignment may imply worse quality of life and more severe VCF.

Clinical Evaluation of Bone Strength and Fracture Risk

PURPOSE OF REVIEW

This paper seeks to evaluate and compare recent advances in the clinical assessment of the changes in bone mechanical properties that take place as a result of osteoporosis and other metabolic bone diseases and their treatments.

RECENT FINDINGS

In addition to the standard of DXA-based areal bone mineral density (aBMD), a variety of methods, including imaging-based structural measurements, finite element analysis (FEA)-based techniques, and alternate methods including ultrasound, bone biopsy, reference point indentation, and statistical shape and density modeling, have been developed which allow for reliable prediction of bone strength and fracture risk. These methods have also shown promise in the evaluation of treatment-induced changes in bone mechanical properties. Continued technological advances allowing for increasingly high-resolution imaging with low radiation dose, together with the expanding adoption of DXA-based predictions of bone structure and mechanics, as well as the increasing awareness of the importance of bone material properties in determining whole-bone mechanics, lead us to anticipate substantial future advances in this field.

Epidemiology of Vertebral Fractures

Vertebral fractures are one of the most common fractures associated with skeletal fragility and can cause as much morbidity as hip fractures. However, the epidemiology of vertebral fractures differs from that of osteoporotic fractures at other skeletal sites in important ways, largely because only one quarter to one-third of vertebral fractures are recognized clinically at the time of their occurrence and otherwise require lateral spine imaging to be recognized. This article first reviews the prevalence and incidence of clinical and radiographic vertebral fractures in populations across the globe and secular trends in the incidence of vertebral fracture over time. Next, associations of vertebral fractures with measures of bone mineral density and bone microarchitecture are reviewed followed by associations of vertebral fracture with various textural measures of trabecular bone, including trabecular bone score. Finally, the article reviews clinical risk factors for vertebral fracture and the association of vertebral fractures with morbidity, mortality, and other subsequent adverse health outcomes.

Utility of the trabecular bone score (TBS) in secondary osteoporosis

Altered bone micro-architecture is an important factor in accounting for fragility fractures. Until recently, it has not been possible to gain information about skeletal microstructure in a way that is clinically feasible. Bone biopsy is essentially a research tool. High-resolution peripheral Quantitative Computed Tomography, while non-invasive, is available only sparsely throughout the world. The trabecular bone score (TBS) is an imaging technology adapted directly from the Dual Energy X-Ray Absorptiometry (DXA) image of the lumbar spine. Thus, it is potentially readily and widely available. In recent years, a large number of studies have demonstrated that TBS is significantly associated with direct measurements of bone micro-architecture, predicts current and future fragility fractures in primary osteoporosis, and may be a useful adjunct to BMD for fracture detection and prediction. In this review, we summarize its potential utility in secondary causes of osteoporosis. In some situations, like glucocorticoid-induced osteoporosis and in diabetes mellitus, the TBS appears to out-perform DXA. It also has apparent value in numerous other disorders associated with diminished bone health, including primary hyperparathyroidism, androgen-deficiency, hormone-receptor positive breast cancer treatment, chronic kidney disease, hemochromatosis, and autoimmune disorders like rheumatoid arthritis. Further research is both needed and warranted to more clearly establish the role of TBS in these and other disorders that adversely affect bone.

Managing osteoporosis in ulcerative colitis: something new?

The authors revise the latest evidence in the literature regarding managing of osteoporosis in ulcerative colitis (UC), paying particular attention to the latest tendency of the research concerning the management of bone damage in the patient affected by UC. It is wise to assess vitamin D status in ulcerative colitis patients to recognize who is predisposed to low levels of vitamin D, whose deficiency has to be treated with oral or parenteral vitamin D supplementation. An adequate dietary calcium intake or supplementation and physical activity, if possible, should be guaranteed. Osteoporotic risk factors, such as smoking and excessive alcohol intake, must be avoided. Steroid has to be prescribed at the lowest possible dosage and for the shortest possible time. Moreover, conditions favoring falling have to been minimized, like carpets, low illumination, sedatives assumption, vitamin D deficiency. It is advisable to assess the fracture risk in all UC patient by the fracture assessment risk tool (FRAX(®) tool), that calculates the ten years risk of fracture for the population aged from 40 to 90 years in many countries of the world. A high risk value could indicate the necessity of treatment, whereas a low risk value suggests a follow-up only. An intermediate risk supports the decision to prescribe bone mineral density (BMD) assessment and a subsequent patient revaluation for treatment. Dual energy X-ray absorptiometry bone densitometry can be used not only for BMD measurement, but also to collect data about bone quality by the means of trabecular bone score and hip structural analysis assessment. These two indices could represent a method of interesting perspectives in evaluating bone status in patients affected by diseases like UC, which may present an impairment of bone quality as well as of bone quantity. In literature there is no strong evidence for instituting pharmacological therapy of bone impairment in UC patients for clinical indications other than those that are also applied to the patients with osteoporosis. Therefore, a reasonable advice is to consider pharmacological treatment for osteoporosis in those UC patients who already present fragility fractures, which bring a high risk of subsequent fractures. Therapy has also to be considered in patients with a high risk of fracture even if it did not yet happen, and particularly when they had long periods of corticosteroid therapy or cumulative high dosages. In patients without fragility fractures or steroid treatment, a medical decision about treatment could be guided by the FRAX tool to determine the intervention threshold. Among drugs for osteoporosis treatment, the bisphosphonates are the most studied ones, with the best and longest evidence of efficacy and safety. Despite this, several questions are still open, such as the duration of treatment, the necessity to discontinue it, the indication of therapy in young patients, particularly in those without previous fractures. Further, it has to be mentioned that a long-term bisphosphonates use in primary osteoporosis has been associated with an increased incidence of dramatic side-effects, even if uncommon, like osteonecrosis of the jaw and atypical sub-trochanteric and diaphyseal femoral fractures. UC is a long-lasting disease and the majority of patients is relatively young. In this scenario primary prevention of fragility fracture is the best cost-effective strategy. Vitamin D supplementation, adequate calcium intake, suitable physical activity (when possible), removing of risk factors for osteoporosis like smoking, and avoiding falling are the best medical acts.

Trabecular bone score: perspectives of an imaging technology coming of age

The trabecular bone score (TBS) is a new method to describe skeletal microarchitecture from the dual energy X-ray absorptiometry (DXA) image of the lumbar spine. While TBS is not a direct physical measurement of trabecular microarchitecture, it correlates with micro-computed tomography (µCT) measures of bone volume fraction, connectivity density, trabecular number, and trabecular separation, and with vertebral mechanical behavior in ex vivo studies. In human subjects, TBS has been shown to be associated with trabecular microarchitecture and bone strength by high resolution peripheral quantitative computed tomography (HRpQCT). Cross-sectional and prospective studies, involving a large number of subjects, have both shown that TBS is associated with vertebral, femoral neck, and other types of osteoporotic fractures in postmenopausal women. Data in men, while much less extensive, show similar findings. TBS is also associated with fragility fractures in subjects with secondary causes of osteoporosis, and preliminary data suggest that TBS might improve fracture prediction when incorporated in the fracture risk assessment system known as FRAX. In this article, we review recent advances that have helped to establish this new imaging technology.

Assessment of risk of femoral neck fracture with radiographic texture parameters: a retrospective study

PURPOSE

To investigate whether femoral neck fracture can be predicted retrospectively on the basis of clinical radiographs by using the combined analysis of bone geometry, textural analysis of trabecular bone, and bone mineral density (BMD).

MATERIALS AND METHODS

Formal ethics committee approval was obtained for the study, and all participants gave informed written consent. Pelvic radiographs and proximal femur BMD measurements were obtained in 53 women aged 79-82 years in 2006. By 2012, 10 of these patients had experienced a low-impact femoral neck fracture. A Laplacian-based semiautomatic custom algorithm was applied to the radiographs to calculate the texture parameters along the trabecular fibers in the lower neck area for all subjects. Intra- and interobserver reproducibility was calculated by using the root mean square average coefficient of variation to evaluate the robustness of the method.

RESULTS

The best predictors of hip fracture were entropy (P = .007; reproducibility coefficient of variation < 1%), the neck-shaft angle (NSA) (P = .017), and the BMD (P = .13). For prediction of fracture, the area under the receiver operating characteristic curve was 0.753 for entropy, 0.608 for femoral neck BMD, and 0.698 for NSA. The area increased to 0.816 when entropy and NSA were combined and to 0.902 when entropy, NSA, and BMD were combined.

CONCLUSION

Textural analysis of pelvic radiographs enables discrimination of patients at risk for femoral neck fracture, and our results show the potential of this conventional imaging method to yield better prediction than that achieved with dual-energy x-ray absorptiometry-based BMD. The combination of the entropy parameter with NSA and BMD can further enhance predictive accuracy.

Trabecular bone score: a noninvasive analytical method based upon the DXA image

The trabecular bone score (TBS) is a gray-level textural metric that can be extracted from the two-dimensional lumbar spine dual-energy X-ray absorptiometry (DXA) image. TBS is related to bone microarchitecture and provides skeletal information that is not captured from the standard bone mineral density (BMD) measurement. Based on experimental variograms of the projected DXA image, TBS has the potential to discern differences between DXA scans that show similar BMD measurements. An elevated TBS value correlates with better skeletal microstructure; a low TBS value correlates with weaker skeletal microstructure. Lumbar spine TBS has been evaluated in cross-sectional and longitudinal studies. The following conclusions are based upon publications reviewed in this article: 1) TBS gives lower values in postmenopausal women and in men with previous fragility fractures than their nonfractured counterparts; 2) TBS is complementary to data available by lumbar spine DXA measurements; 3) TBS results are lower in women who have sustained a fragility fracture but in whom DXA does not indicate osteoporosis or even osteopenia; 4) TBS predicts fracture risk as well as lumbar spine BMD measurements in postmenopausal women; 5) efficacious therapies for osteoporosis differ in the extent to which they influence the TBS; 6) TBS is associated with fracture risk in individuals with conditions related to reduced bone mass or bone quality. Based on these data, lumbar spine TBS holds promise as an emerging technology that could well become a valuable clinical tool in the diagnosis of osteoporosis and in fracture risk assessment.

Trabecular homogeneity index derived from plain radiograph to evaluate bone quality

Radiographic texture analysis has been developed lately to improve the assessment of bone architecture as a determinant of bone quality. We validate here an algorithm for the evaluation of trabecular homogeneity index (HI) in the proximal femur from hip radiographs, with a focus on the impact of the principal compressive system of the trabecular bone, and evaluate its correlation with femoral strength, bone mineral density (BMD), and volumetric trabecular structure parameters. A semiautomatic custom-made algorithm was applied to calculate the HI in the femoral neck and trochanteric areas from radiographs of 178 femoral bone specimens (mean age 79.3 ± 10.4 years). Corresponding neck region was selected in CT scans to calculate volumetric parameters of trabecular structure. The site-specific BMDs were assessed from dual-energy X-ray absorptiometry (DXA), and the femoral strength was experimentally tested in side-impact configuration. Regression analysis was performed between the HI and biomechanical femoral strength, BMD, and volumetric parameters. The correlation between HI and failure load was R(2)  = 0.50; this result was improved to R(2)  = 0.58 for cervical fractures alone. The discrimination of bones with high risk of fractures (load <3000 N) was similar for HI and BMD (AUC = 0.87). Regression analysis between the HIs versus site-specific BMDs yielded R(2)  = 0.66 in neck area, R(2)  = 0.60 in trochanteric area, and an overall of R(2)  = 0.66 for the total hip. Neck HI and BMD correlated significantly with volumetric structure parameters. We present here a method to assess HI that can explain 50% of an experimental failure load and determines bones with high fracture risk with similar accuracy as BMD. The HI also had good correlation with DXA and computed tomography-derived data.

Computational identification and quantification of trabecular microarchitecture classes by 3-D texture analysis-based clustering

High resolution peripheral quantitative computed tomography (HR-pQCT) permits the non-invasive assessment of cortical and trabecular bone density, geometry, and microarchitecture. Although researchers have developed various post-processing algorithms to quantify HR-pQCT image properties, few of these techniques capture image features beyond global structure-based metrics. While 3D-texture analysis is a key approach in computer vision, it has been utilized only infrequently in HR-pQCT research. Motivated by high isotropic spatial resolution and the information density provided by HR-pQCT scans, we have developed and evaluated a post-processing algorithm that quantifies microarchitecture characteristics via texture features in HR-pQCT scans. During a training phase in which clustering was applied to texture features extracted from each voxel of trabecular bone, three distinct clusters, or trabecular microarchitecture classes (TMACs) were identified. These TMACs represent trabecular bone regions with common texture characteristics. The TMACs were then used to automatically segment the voxels of new data into three regions corresponding to the trained cluster features. Regional trabecular bone texture was described by the histogram of relative trabecular bone volume covered by each cluster. We evaluated the intra-scanner and inter-scanner reproducibility by assessing the precision errors (PE), intra class correlation coefficients (ICC) and Dice coefficients (DC) of the method on 14 ultradistal radius samples scanned on two HR-pQCT systems. DC showed good reproducibility in intra-scanner set-up with a mean of 0.870±0.027 (no unit). Even in the inter-scanner set-up the ICC showed high reproducibility, ranging from 0.814 to 0.964. In a preliminary clinical test application, the TMAC histograms appear to be a good indicator, when differentiating between postmenopausal women with (n=18) and without (n=18) prevalent fragility fractures. In conclusion, we could demonstrate that 3D-texture analysis and feature clustering seems to be a promising new HR-pQCT post-processing tool with good reproducibility, even between two different scanners.

Prospective assessment of bone texture parameters at the hand in rheumatoid arthritis

OBJECTIVE

Fractal bone analysis (Hmean) is a texture parameter reflecting bone microarchitecture. The BMA device (D3A™ Medical Systems, Orléans, France) is a high-resolution X-ray device that allows assessment of bone texture analysis. We aimed to measure Hmean in rheumatoid arthritis patients at the second and third metacarpal bones, at baseline and after 1 year of follow-up, and to assess the relationship of Hmean and rheumatoid arthritis disease parameters.

METHODS

Patients with rheumatoid arthritis according to ACR criteria were included. They were assessed over 1 year, in the context of a prospective study conducted in Maastricht. For this substudy, activity of the disease was assessed by erythrocyte sedimentation rate, C-reactive protein and Disease Activity Score 28 performed at each visit. Radiographic bone damage was assessed using hand and feet radiographs at baseline and on a 1-year basis. The bone texture parameters were evaluated on the second and third metacarpal heads of the left hand using BMA device.

RESULTS

One hundred and sixty-five rheumatoid arthritis patients were included in this study. At baseline, Hmean was negatively correlated with age [r=-0.22 (P=0.013)] and erythrocyte sedimentation rate [r=-0.16 (P=0.039)]. No significant correlation was found between Hmean and Disease Activity Score, disease activity Visual Analog Scale, daily corticosteroid dose and C-reactive protein. There was a significant increase in Hmean of second and third metacarpal bones over 1 year (1.6% and 1.3%, P<0.01) except in patients with local second and third metacarpal bones erosion.

CONCLUSION

The bone texture parameter Hmean is influenced by age, inflammation and local erosions in rheumatoid arthritis.

Comparison of heterogeneity quantification algorithms for brain SPECT perfusion images

Background

Several algorithms from the literature were compared with the original random walk (RW) algorithm for brain perfusion heterogeneity quantification purposes. Algorithms are compared on a set of 210 brain single photon emission computed tomography (SPECT) simulations and 40 patient exams.

Methods

Five algorithms were tested on numerical phantoms. The numerical anthropomorphic Zubal head phantom was used to generate 42 (6 × 7) different brain SPECT simulations. Seven diffuse cortical heterogeneity levels were simulated with an adjustable Gaussian noise function and six focal perfusion defect levels with temporoparietal (TP) defects. The phantoms were successively projected and smoothed with Gaussian kernel with full width at half maximum (FWHM = 5 mm), and Poisson noise was added to the 64 projections. For each simulation, 5 Poisson noise realizations were performed yielding a total of 210 datasets. The SPECT images were reconstructed using filtered black projection (Hamming filter: α = 0.5).

The five algorithms or measures tested were the following: the coefficient of variation, the entropy and local entropy, fractal dimension (FD) (box counting and Fourier power spectrum methods), the gray-level co-occurrence matrix (GLCM), and the new RW.

The heterogeneity discrimination power was obtained with a linear regression for each algorithm. This regression line is a mean function of the measure of heterogeneity compared to the different diffuse heterogeneity and focal defect levels generated in the phantoms. A greater slope denotes a larger separation between the levels of diffuse heterogeneity.

The five algorithms were computed using 40 99mTc-ethyl-cysteinate-dimer (ECD) SPECT images of patients referred for memory impairment. Scans were blindly ranked by two physicians according to the level of heterogeneity, and a consensus was obtained. The rankings obtained by the algorithms were compared with the physicians' consensus ranking.

Results

The GLCM method (slope = 58.5), the fractal dimension (35.9), and the RW method (31.6) can differentiate the different levels of diffuse heterogeneity. The GLCM contrast parameter method is not influenced by a focal defect contrary to the FD and RW methods. A significant correlation was found between the RW method and the physicians' classification (r = 0.86; F = 137; p < 0.0001).

Conclusions

The GLCM method can quantify the different levels of diffuse heterogeneity in brain-simulated SPECT images without an influence from the focal cortical defects. However, GLCM classification was not correlated with the physicians' classification (Rho = −0.099). The RW method was significantly correlated with the physicians' heterogeneity perception but is influenced by the existence of a focal defect.

Bone texture analysis of human femurs using a new device (BMA™) improves failure load prediction

We measured bone texture parameters of excised human femurs with a new device (BMA™). We also measured bone mineral density by DXA and investigated the performance of these parameters in the prediction of failure load. Our results suggest that bone texture parameters improve failure load prediction when added to bone mineral density.

INTRODUCTION

Bone mineral density (BMD) is a strong determinant of bone strength. However, nearly half of the fractures occur in patients with BMD which does not reach the osteoporotic threshold. In order to assess fracture risk properly, other factors are important to be taken into account such as clinical risk factors as well as macro- and microarchitecture of bone. Bone microarchitecture is usually assessed by high-resolution QCT, but this cannot be applied in routine clinical settings due to irradiation, cost and availability concerns. Texture analysis of bone has shown to be correlated to bone strength.

METHODS

We used a new device to get digitized X-rays of 12 excised human femurs in order to measure bone texture parameters in three different regions of interest (ROIs). We investigated the performance of these parameters in the prediction of the failure load using biomechanical tests. Texture parameters measured were the fractal dimension (Hmean), the co-occurrence matrix, and the run length matrix. We also measured bone mineral density by DXA in the same ROIs as well as in standard DXA hip regions.

RESULTS

The Spearman correlation coefficient between BMD and texture parameters measured in the same ROIs ranged from -0.05 (nonsignificant (NS)) to 0.57 (p = 0.003). There was no correlation between Hmean and co-occurrence matrix nor Hmean and run length matrix in the same ROI (r = -0.04 to 0.52, NS). Co-occurrence matrix and run length matrix in the same ROI were highly correlated (r = 0.90 to 0.99, p < 0.0001). Univariate analysis with the failure load revealed significant correlation only with BMD results, not texture parameters. Multiple regression analysis showed that the best predictors of failure load were BMD, Hmean, and run length matrix at the femoral neck, as well as age and sex, with an adjusted r (2) = 0.88. Added to femoral neck BMD, Hmean and run length matrix at the femoral neck (without the effect of age and sex) improved failure load prediction (compared to femoral neck BMD alone) from adjusted r (2) = 0.67 to adjusted r (2) = 0.84.

CONCLUSION

Our results suggest that bone texture measurement improves failure load prediction when added to BMD.

Bone microarchitecture assessed by TBS predicts osteoporotic fractures independent of bone density: the Manitoba study

The measurement of BMD by dual-energy X-ray absorptiometry (DXA) is the "gold standard" for diagnosing osteoporosis but does not directly reflect deterioration in bone microarchitecture. The trabecular bone score (TBS), a novel gray-level texture measurement that can be extracted from DXA images, correlates with 3D parameters of bone microarchitecture. Our aim was to evaluate the ability of lumbar spine TBS to predict future clinical osteoporotic fractures. A total of 29,407 women 50 years of age or older at the time of baseline hip and spine DXA were identified from a database containing all clinical results for the Province of Manitoba, Canada. Health service records were assessed for the incidence of nontraumatic osteoporotic fracture codes subsequent to BMD testing (mean follow-up 4.7 years). Lumbar spine TBS was derived for each spine DXA examination blinded to clinical parameters and outcomes. Osteoporotic fractures were identified in 1668 (5.7%) women, including 439 (1.5%) spine and 293 (1.0%) hip fractures. Significantly lower spine TBS and BMD were identified in women with major osteoporotic, spine, and hip fractures (all p < 0.0001). Spine TBS and BMD predicted fractures equally well, and the combination was superior to either measurement alone (p < 0.001). Spine TBS predicts osteoporotic fractures and provides information that is independent of spine and hip BMD. Combining the TBS trabecular texture index with BMD incrementally improves fracture prediction in postmenopausal women.

Combination of radiograph-based trabecular and geometrical parameters can discriminate cervical hip fractures from controls in individuals with BMD in non-osteoporotic range

Majority of hip fractures occur in individuals with bone mineral density (BMD) in non-osteoporotic range. This suggests that factors other than BMD are associated with increased fracture risk in these individuals. The aim of this study was to investigate the combined ability of radiograph-based trabecular and geometrical parameters to discriminate cervical hip fractures from controls in individuals with non-osteoporotic BMD. A total of 39 postmenopausal females with non-pathologic cervical hip fracture were recruited to the study. Nineteen of the fracture patients (48.7%) had non-osteoporotic BMD and they constituted the fracture group. The control group consisted of 35 BMD-matched non-osteoporotic females. Several geometrical and trabecular parameters were extracted from plain pelvic radiographs, and their combined ability to discriminate fracture patients from controls was studied using a receiver operating characteristics (ROC) analysis. Significant differences in several radiograph-based geometrical and trabecular parameters were found between the fracture patients and controls, whereas no statistically significant difference in BMD was observed (p=0.92) between the groups. Area under the ROC curve was 0.993 (95% CI 0.977-1.008) for the combined multiple regression model, which included both trabecular and geometrical parameters as explanatory factors. Here, the sensitivity of 100% was achieved with the specificity of 94%. In a cross-validation of the model, 94.4% of the fracture patients, and 94.1% of the controls were classified correctly. The combination of radiograph-based trabecular and geometrical parameters was able to discriminate the cervical hip fracture cases from controls with similar BMD, showing that the method can provide additional information on bone structure and fracture risk beyond BMD.

Correlations between trabecular bone score, measured using anteroposterior dual-energy X-ray absorptiometry acquisition, and 3-dimensional parameters of bone microarchitecture: an experimental study on human cadaver vertebrae

Developing a novel technique for the efficient, noninvasive clinical evaluation of bone microarchitecture remains both crucial and challenging. The trabecular bone score (TBS) is a new gray-level texture measurement that is applicable to dual-energy X-ray absorptiometry (DXA) images. Significant correlations between TBS and standard 3-dimensional (3D) parameters of bone microarchitecture have been obtained using a numerical simulation approach. The main objective of this study was to empirically evaluate such correlations in anteroposterior spine DXA images. Thirty dried human cadaver vertebrae were evaluated. Micro-computed tomography acquisitions of the bone pieces were obtained at an isotropic resolution of 93μm. Standard parameters of bone microarchitecture were evaluated in a defined region within the vertebral body, excluding cortical bone. The bone pieces were measured on a Prodigy DXA system (GE Medical-Lunar, Madison, WI), using a custom-made positioning device and experimental setup. Significant correlations were detected between TBS and 3D parameters of bone microarchitecture, mostly independent of any correlation between TBS and bone mineral density (BMD). The greatest correlation was between TBS and connectivity density, with TBS explaining roughly 67.2% of the variance. Based on multivariate linear regression modeling, we have established a model to allow for the interpretation of the relationship between TBS and 3D bone microarchitecture parameters. This model indicates that TBS adds greater value and power of differentiation between samples with similar BMDs but different bone microarchitectures. It has been shown that it is possible to estimate bone microarchitecture status derived from DXA imaging using TBS.

Development and testing of texture discriminators for the analysis of trabecular bone in proximal femur radiographs

PURPOSE

Texture analysis of femur radiographs may serve as a potential low cost technique to predict osteoporotic fracture risk and has received considerable attention in the past years. A further application of this technique may be the measurement of the quality of specific bone compartments to provide useful information for treatment of bone fractures. Two challenges of texture analysis are the selection of the best suitable texture measure and reproducible placement of regions of interest (ROIs). The goal of this in vitro study was to automatically place ROIs in radiographs of proximal femur specimens and to calculate correlations between various different texture analysis methods and the femurs' anchorage strength.

METHODS

Radiographs were obtained from 14 femoral specimens and bone mineral density (BMD) was measured in the femoral neck. Biomechanical testing was performed to assess the anchorage strength in terms of failure load, breakaway torque, and number of cycles. Images were segmented using a framework that is based on the usage of level sets and statistical in-shape models. Five ROIs were automatically placed in the head, upper and lower neck, trochanteric, and shaft compartment in an atlas subject. All other subjects were registered rigidly, affinely, and nonlinearly, and the resulting transformation was used to map the five ROIs onto the individual femora.

RESULTS

In each ROI, texture features were extracted using gray level co-occurence matrices (GLCM), third-order GLCM, morphological gradients (MGs), Minkowski dimensions (MDs), Minkowski functionals (MFs), Gaussian Markov random fields, and scaling index method (SIM). Coefficients of determination for each texture feature with parameters of anchorage strength were computed. In a stepwise multiregression analysis, the most predictive parameters were identified in different models. Texture features were highly correlated with anchorage strength estimated by the failure load of up to R2=0.61 (MF and MG features, p<0.01) and were partially independent of BMD. The correlations were dependent on the choice of the ROI and the texture measure. The best predictive multiregression model for failure load R2adj=0.86 (p<0.001) included a set of recently developed texture methods (MF and SIM) but excluded bone mineral density and commonly used texture measures.

CONCLUSIONS

The results suggest that texture information contained in trabecular bone structure visualized on radiographs may predict whether an implant anchorage can be used and may determine the local bone quality from preoperative radiographs.

Radiographic texture analysis of densitometric calcaneal images: relationship to clinical characteristics and to bone fragility

Osteoporotic fractures are related not only to bone mineral density (BMD) but also to bone structure or microarchitecture, which is not assessed routinely with currently available methods. We have developed radiographic texture analysis (RTA) for calcaneal images from a peripheral densitometer as an easy, noninvasive method for assessing bone structure. We conducted a cross-sectional study of the relationship between RTA and prevalent vertebral fractures (n = 148) among 900 subjects (ages 19 to 99 years, 94 males) referred for bone densitometry as part of their routine medical care. RTA features were derived from Fourier-based image analysis of the radiographic texture pattern (including root mean square, first moment, and power spectral analyses). RTA features were associated with age, weight, gender, and race, as well as glucocorticoid use. When controlling for clinical risk factors and BMD (or a summary measure calculated using FRAX algorithms), RTA features were significantly different for subjects with and without prevalent vertebral fractures [adjusted odds ratio (OR) = 1.5 per 1 standard deviation (SD) decrease in RTA feature beta, 95% confidence interval (CI) 1.2-1.8, p = .001]. Gender and use of pharmacologic therapy for osteoporosis did not significantly affect this association, suggesting that RTA can be applied to a wide range of densitometry patients. We conclude that RTA obtained using a portable instrument has a potential as a noninvasive method to enhance identification of patients at increased risk of osteoporotic fractures.

Brain perfusion heterogeneity measurement based on Random Walk algorithm: choice and influence of inner parameters

A Random Walk (RW) algorithm was designed to quantify the level of diffuse heterogeneous perfusion in brain SPECT images in patients suffering from systemic brain disease or from drug-induced therapy. The goal of the present paper is to understand the behavior of the RW method on different kinds of images (extrinsic parameters) and also to understand how to choose the right parameters of the RW (intrinsic parameters) depending on the image characteristics (i.e. SPECT images). "Extrinsic parameters" are related to the image characteristics (level/size of defect and diffuse heterogeneity) and "intrinsic" parameters are related to the parameters of the method (number (N(rw)) and length of walk (L(rw)), temperature (T) and slowing parameter (S)). Two successive studies were conducted to test the influence of these parameters on the RW result. In the first study, calibrated checkerboard images are used to test the influence of "extrinsic parameters" (i.e. image characteristics) on the RW result (R-value). The R-value was tested as a function of (i) the size of black & white (B&W) squares simulating the size of a cortical defect, (ii) the intensity level gaps between the B&W squares simulating the intensity of the cortical defect and (iii) intensity (=variance) of noise, simulating the diffuse heterogeneity. The second study was constructed with simulated representative brain SPECT images, to test the "intrinsic" parameters. The R-value was tested regarding the influence of four parameters: S, T, N(rw) and L(rw). The third study is constructed so as to see if the classification by diffuse heterogeneity of real brain SPECT images is the same if it's made by senior clinicians or by RW algorithm.

RESULTS

Study 1: the RW was strongly influenced by all the characteristics of the images. Moreover, these characteristics interact with each other. The RW is influenced most by diffuse heterogeneity, then by intensity and finally by the size of a defect. Study 2: N(rw) and L(rw) values of 1000 give an optimal reproducibility of the measurement (mean standard deviation<0.1), a fast computation time (time<0.5s/image) and have a maximum difference in terms of R-value between the two extreme images corresponding to the range of the population studied. The best S and T values for SPECT images are 3 and 15, respectively. Study 3: A significant correlation was found between RW ranking and the physicians' consensus (rho=0.789; p<0.0001).

CONCLUSION

This study confirms that the RW method is able to measure the heterogeneity of brain SPECT images even in the presence of a large defect. However, the result of the method is strongly influenced by the "intrinsic" parameters, so the program should be calibrated for each different type of image.

Trabecular bone structure analysis in the osteoporotic spine using a clinical in vivo setup for 64-slice MDCT imaging: comparison to microCT imaging and microFE modeling

Assessment of trabecular microarchitecture may improve estimation of biomechanical strength, but visualization of trabecular bone structure in vivo is challenging. We tested the feasibility of assessing trabecular microarchitecture in the spine using multidetector CT (MDCT) on intact human cadavers in an experimental in vivo-like setup. BMD, bone structure (e.g., bone volume/total volume = BV/TV; trabecular thickness = Tb.Th; structure model index = SMI) and bone texture parameters were evaluated in 45 lumbar vertebral bodies using MDCT (mean in-plane pixel size, 274 microm(2); slice thickness, 500 microm). These measures were correlated with structure measures assessed with microCT at an isotropic spatial resolution of 16 microm and to microfinite element models (microFE) of apparent modulus and stiffness. MDCT-derived BMD and structure measures showed significant correlations to the density and structure obtained by microCT (BMD, R(2) = 0.86, p < 0.0001; BV/TV, R(2) = 0.64, p < 0.0001; Tb.Th, R(2) = 0.36, p < 0.01). When comparing microCT-derived measures with microFE models, the following correlations (p < 0.001) were found for apparent modulus and stiffness, respectively: BMD (R(2) = 0.58 and 0.66), BV/TV (R(2) = 0.44 and 0.58), and SMI (R(2) = 0.44 and 0.49). However, the overall highest correlation (p < 0.001) with microFE app. modulus (R(2) = 0.75) and stiffness (R(2) = 0.76) was achieved by the combination of QCT-derived BMD with the bone texture measure Minkowski Dimension. In summary, although still limited by its spatial resolution, trabecular bone structure assessment using MDCT is overall feasible. However, when comparing with microFE-derived bone properties, BMD is superior compared with single parameters for microarchitecture, and correlations further improve when combining with texture measures.

Advances in osteoporosis imaging

In the assessment of osteoporosis, the measurement of bone mineral density (BMD(a)) obtained from dual energy X-ray absorptiometry (DXA; g/cm(2)) is the most widely used parameter. However, bone strength and fracture risk are also influenced by parameters of bone quality such as micro-architecture and tissue properties. This article reviews the radiological techniques currently available for imaging and quantifying bone structure, as well as advanced techniques to image bone quality. With the recent developments in magnetic resonance (MR) techniques, including the availability of clinical 3T scanners, and advances in computed tomography (CT) technology (e.g. clinical Micro-CT), in-vivo imaging of the trabecular bone architecture is becoming more feasible. Several in-vitro studies have demonstrated that bone architecture, measured by MR or CT, was a BMD-independent determinant of bone strength. In-vivo studies showed that patients with, and without, osteoporotic fractures could better be separated with parameters of bone architecture than with BMD. Parameters of trabecular architecture were more sensitive to treatment effects than BMD. Besides the 3D tomographic techniques, projection radiography has been used in the peripheral skeleton as an additional tool to better predict fracture risk than BMD alone. The quantification of the trabecular architecture included parameters of scale, shape, anisotropy and connectivity. Finite element analyses required highest resolution, but best predicted the biomechanical properties of the bone. MR diffusion and perfusion imaging and MR spectroscopy may provide measures of bone quality beyond trabecular micro-architecture.

Local topological analysis of densitometer-generated scan images of the proximal femur for differentiation between patients with hip fracture and age-matched controls

SUMMARY

We evaluate densitometer-generated scan images of the proximal femur with respect to topological properties of bone mineral distribution patterns in selected regions of interest. In a population of 100 post-menopausal women, the method has a highly discriminative potential with a performance superior to standard densitometry. Results vary with anatomical location within the proximal femur.

INTRODUCTION

The objectives of the study were to evaluate densitometer-generated scan images of the proximal femur with respect to topological properties of bone mineral distribution patterns in selected regions of interest, to test the ability for differentiation between post-menopausal women hip fracture and controls, and to compare results with standard bone densitometry.

MATERIALS AND METHODS

We used dual-energy X-ray absorptiometry (DXA) to measure the femoral bone mineral density (BMD) of 100 post-menopausal women (73.4 +/- 12.2), 50 of whom had a recent hip fracture. Local bone mineral distribution in the scanner-generated images was analyzed in the standard DXA-regions of interest (ROIs; femoral neck, the shaft, the trochanteric area; and the total hip) using an optimized, local topological parameter MF2D. Performance of topological analysis and BMD was tested by receiver-operator characteristic and discriminant analysis.

RESULTS

Area under the curve (AUC) for correct differentiation between patients with and without fractures by BMD in the different ROIs ranged from 0.64 to 0.71; AUC of regional density-pattern analysis varied between 0.79 and 0.84. Using multivariate statistical models, between 71% and 84% of patients were correctly identified as fracture/non-fracture cases by regional topological analysis, whereas BMD reached levels from 58% to 68%.

CONCLUSION

Our analysis indicates that identification of patients with hip fracture by regional evaluation of density patterns varies with anatomical location within the proximal femur. In our study population, performance of the novel parameter was superior to densitometry.

Evaluation of the potential use of trabecular bone score to complement bone mineral density in the diagnosis of osteoporosis: a preliminary spine BMD-matched, case-control study

The trabecular bone score (TBS) is a new parameter that is determined from gray-level analysis of dual-energy X-ray absorptiometry (DXA) images. It relies on the mean thickness and volume fraction of trabecular bone microarchitecture. This was a preliminary case-control study to evaluate the potential diagnostic value of TBS as a complement to bone mineral density (BMD), by comparing postmenopausal women with and without fractures. The sample consisted of 45 women with osteoporotic fractures (5 hip fractures, 20 vertebral fractures, and 20 other types of fracture) and 155 women without a fracture. Stratification was performed, taking into account each type of fracture (except hip), and women with and without fractures were matched for age and spine BMD. BMD and TBS were measured at the total spine. TBS measured at the total spine revealed a significant difference between the fracture and age- and spine BMD-matched nonfracture group, when considering all types of fractures and vertebral fractures. In these cases, the diagnostic value of the combination of BMD and TBS likely will be higher compared with that of BMD alone. TBS, as evaluated from standard DXA scans directly, potentially complements BMD in the detection of osteoporotic fractures. Prospective studies are necessary to fully evaluate the potential role of TBS as a complementary risk factor for fracture.

Anniversary paper: History and status of CAD and quantitative image analysis: the role of Medical Physics and AAPM

The roles of physicists in medical imaging have expanded over the years, from the study of imaging systems (sources and detectors) and dose to the assessment of image quality and perception, the development of image processing techniques, and the development of image analysis methods to assist in detection and diagnosis. The latter is a natural extension of medical physicists' goals in developing imaging techniques to help physicians acquire diagnostic information and improve clinical decisions. Studies indicate that radiologists do not detect all abnormalities on images that are visible on retrospective review, and they do not always correctly characterize abnormalities that are found. Since the 1950s, the potential use of computers had been considered for analysis of radiographic abnormalities. In the mid-1980s, however, medical physicists and radiologists began major research efforts for computer-aided detection or computer-aided diagnosis (CAD), that is, using the computer output as an aid to radiologists-as opposed to a completely automatic computer interpretation-focusing initially on methods for the detection of lesions on chest radiographs and mammograms. Since then, extensive investigations of computerized image analysis for detection or diagnosis of abnormalities in a variety of 2D and 3D medical images have been conducted. The growth of CAD over the past 20 years has been tremendous-from the early days of time-consuming film digitization and CPU-intensive computations on a limited number of cases to its current status in which developed CAD approaches are evaluated rigorously on large clinically relevant databases. CAD research by medical physicists includes many aspects-collecting relevant normal and pathological cases; developing computer algorithms appropriate for the medical interpretation task including those for segmentation, feature extraction, and classifier design; developing methodology for assessing CAD performance; validating the algorithms using appropriate cases to measure performance and robustness; conducting observer studies with which to evaluate radiologists in the diagnostic task without and with the use of the computer aid; and ultimately assessing performance with a clinical trial. Medical physicists also have an important role in quantitative imaging, by validating the quantitative integrity of scanners and developing imaging techniques, and image analysis tools that extract quantitative data in a more accurate and automated fashion. As imaging systems become more complex and the need for better quantitative information from images grows, the future includes the combined research efforts from physicists working in CAD with those working on quantitative imaging systems to readily yield information on morphology, function, molecular structure, and more-from animal imaging research to clinical patient care. A historical review of CAD and a discussion of challenges for the future are presented here, along with the extension to quantitative image analysis.

Clinical interest of bone texture analysis in osteoporosis: a case control multicenter study

We demonstrate the clinical interest of bone texture analysis with a new high resolution X-ray device. We have found that the combination of BMD and texture parameter values provided a better assessment of the fracture risk than that obtainable solely by BMD measurement.

INTRODUCTION

Osteoporosis is characterized by BMD and trabecular bone microarchitecture. We have developed a new high-resolution X-ray device with direct digitization. The aim of this study was to demonstrate in a multicenter case control study the clinical interest of bone texture analysis with this new device.

METHODS

In this cross-sectional multicenter case-control population study in post-menopausal women, 159 osteoporotic fractures were compared with 219 control cases. Images were obtained on calcaneus with a direct digital X-ray device (BMA, D3A Medical Systems). Co-occurrence, run-length matrices and the fractal parameter Hmean were evaluated. BMD was measured at the lumbar spine (LS), femoral neck (FN) and total hip (TH) by DXA.

RESULTS

The three texture parameters were significantly lower in osteoporotic fracture cases than in control cases. These differences persisted after adjustment for TH BMD. Receiver operating characteristic curves were used to compare the discriminant capacity of texture parameters and BMD measurements for fracture. The highest areas under curve (AUC) were 0.721 for TH BMD and 0.706 for Hmean (AUC THBMD vs. AUC Hmean, p = NS). We determined the threshold between high and low Hmean parameter values and then the odds ratios (OR) of fracture for low Hmean, for BMD < or =2.5 SD in the T-score and for combinations of both parameters. The OR of fracture for low H was 2.72 (95% CI, 1.36-5.4). For a FN BMD < or = -2.5 SD, the OR of 4.78 (2.19-10.43) shifted to 14.06 (4.41-44.85) adding H.

CONCLUSIONS

These data confirmed the clinical interest of the combination of BMD and texture parameters to improve the assessment of the risk of fracture other that obtainable by the sole BMD measurement.

Laws' masks descriptors applied to bone texture analysis: an innovative and discriminant tool in osteoporosis

OBJECTIVE

The objective of this study was to explore Laws' masks analysis to describe structural variations of trabecular bone due to osteoporosis on high-resolution digital radiographs and to check its dependence on the spatial resolution. Laws' masks are well established as one of the best methods for texture analysis in image processing and are used in various applications, but not in bone tissue characterisation. This method is based on masks that aim to filter the images. From each mask, five classical statistical parameters can be calculated.

MATERIALS AND METHODS

The study was performed on 182 healthy postmenopausal women with no fractures and 114 age-matched women with fractures [26 hip fractures (HFs), 29 vertebrae fractures (VFs), 29 wrist fractures (WFs) and 30 other fractures (OFs)]. For all subjects radiographs were obtained of the calcaneus with a new high-resolution X-ray device with direct digitisation (BMA, D3A, France). The lumbar spine, femoral neck, and total hip bone mineral density (BMD) were assessed by dual-energy X-ray absorptiometry.

RESULTS

In terms of reproducibility, the best results were obtained with the TRE5E5 mask, especially for three parameters: "mean", "standard deviation" and "entropy" with, respectively, in vivo mid-term root mean square average coefficient of variation (RMSCV)%= 1.79, 4.24 and 2.05. The "mean" and "entropy" parameters had a better reproducibility but "standard deviation" showed a better discriminant power. Thus, for univariate analysis, the difference between subjects with fractures and controls was significant (P<10(-3)) and significant for each fracture group independently (P<10(-4) for HF, P=0.025 for VF and P< 10(-3) for OF). After multivariate analysis with adjustment for age and total hip BMD, the difference concerning the "standard deviation" parameter remained statistically significant between the control group and the HF and VF groups (P<5 x 10(-5), and P=0.04, respectively). No significant correlation between these Laws' masks parameters and BMD was obtained. In addition, this study showed the dependence of Laws' masks parameters on image resolution, which confirms the necessity to perform Laws' textural measurement on high-resolution images.

CONCLUSION

The reproducibility and discriminant power of the Laws' masks analysis has been demonstrated on bone images; thus, this method constitutes a promising routine technique for the determination of osteoporosis fracture risk from radiographs.

Reproducibility and sources of variability in radiographic texture analysis of densitometric calcaneal images

Radiographic texture analysis (RTA) is a computerized analysis of the spatial pattern of radiographic images used as a way of evaluating bone structure. We have shown that RTA performed on high-resolution heel images obtained using a portable densitometer differentiates subjects with and without osteoporotic fractures. In the present study, short-term precision of RTA was examined on densitometric heel images obtained from 33 subjects scanned 8 times each, with 3 observers placing a region of interest (ROI) 3 times on each image. The long-term precision was examined on images obtained from 10 subjects 3 times on each of 3 days separated by 1 week, with 2 observers placing an ROI on each image. The RTA features examined included the root mean square (RMS) variation, a measure of the contrast between the light and dark areas of the image, the first moment of the power spectrum, a measure of the spatial frequency of the trabecular pattern, and Minkowski fractal (MINK), a measure of roughness/smoothness of the trabecular pattern. The precision of the RTA features expressed as coefficient of variation ranged between the lowest of 0.5-0.7% for MINK and the highest of 14-16% for RMS. The short- and long-term precision was similar, and was not significantly influenced by repositioning and rescanning, or by ROI placement by the same or different observers. Significant sources of variability of RTA were the between-subject differences and differences between regions of the heel, but not differences due to repositioning, rescanning in the same position, or ROI placement by the same or different observers. We conclude that technical aspects of image acquisition and processing are adequate to allow further development of RTA of the densitometric images for clinical application as a method for noninvasive assessment of bone structure.

Experimental hip fracture load can be predicted from plain radiography by combined analysis of trabecular bone structure and bone geometry

Computerized analysis of the trabecular structure was used to test whether femur failure load can be estimated from radiographs. The study showed that combined analysis of trabecular bone structure and geometry predicts in vitro failure load with similar accuracy as DXA.

INTRODUCTION

Since conventional radiography is widely available with low imaging cost, it is of considerable interest to discover how well bone mechanical competence can be determined using this technology. We tested the hypothesis that the mechanical strength of the femur can be estimated by the combined analysis of the bone trabecular structure and geometry.

METHODS

The sample consisted of 62 cadaver femurs (34 females, 28 males). After radiography and DXA, femora were mechanically tested in side impact configuration. Fracture patterns were classified as being cervical or trochanteric. Computerized image analysis was applied to obtain structure-related trabecular parameters (trabecular bone area, Euler number, homogeneity index, and trabecular main orientation), and set of geometrical variables (neck-shaft angle, medial calcar and femoral shaft cortex thicknesses, and femoral neck axis length). Multiple linear regression analysis was performed to identify the variables that best explain variation in BMD and failure load between subjects.

RESULTS

In cervical fracture cases, trabecular bone area and femoral neck axis length explained 64% of the variability in failure loads, while femoral neck BMD also explained 64%. In trochanteric fracture cases, Euler number and femoral cortex thickness explained 66% of the variability in failure load, while trochanteric BMD explained 72%.

CONCLUSIONS

Structural parameters of trabecular bone and bone geometry predict in vitro failure loads of the proximal femur with similar accuracy as DXA, when using appropriate image analysis technology.

Temporal radiographic texture analysis in the detection of periprosthetic osteolysis

Periprosthetic osteolysis is one of the most serious long-term problems in total hip arthroplasty. It has been primarily attributed to the body's inflammatory response to submicron polyethylene particles worn from the hip implant, and it leads to bone loss and structural deterioration in the surrounding bone. It was previously demonstrated that radiographic texture analysis (RTA) has the ability to distinguish between osteolysis and normal cases at the time of clinical detection of the disease; however, that analysis did not take into account the changes in texture over time. The goal of this preliminary analysis, however, is to assess the ability of temporal radiographic texture analysis (tRTA) to distinguish between patients who develop osteolysis and normal cases. Two tRTA methods were used in the study: the RTA feature change from baseline at various follow-up intervals and the slope of the best-fit line to the RTA data series. These tRTA methods included Fourier-based and fractal-based features calculated from digitized images of 202 total hip replacement cases, including 70 that developed osteolysis. Results show that separation between the osteolysis and normal groups increased over time for the feature difference method, as the disease progressed, with area under the curve (AUC) values from receiver operating characteristic analysis of 0.65 to 0.72 at 15 years postsurgery. Separation for the slope method was also evident, with AUC values ranging from 0.65 to 0.76 for the task of distinguishing between osteolysis and normal cases. The results suggest that tRTA methods have the ability to measure changes in trabecular structure, and may be useful in the early detection of periprosthetic osteolysis.

Radiographic texture analysis in the characterization of trabecular patterns in periprosthetic osteolysis

RATIONALE AND OBJECTIVES

Periprosthetic osteolysis is a disease attributed to the body's reaction to fine polyethylene wear debris shed from total hip replacements. The purpose of this preliminary study was to investigate the ability of radiographic texture analysis (RTA) to characterize the trabecular texture patterns on pelvic images for osteolysis and normal total hip arthroplasty (THA) cases.

MATERIALS AND METHODS

Fourier-based and fractal-based texture features were calculated for a database of digitized radiographs from 202 THA cases, 70 of which developed osteolysis. The features were calculated from regions of interest selected at two time points: less than 1 month after surgery, and at the first clinical indication of osteolysis (or randomly selected follow-up time for normal cases). Receiver operating characteristic (ROC) analysis was used to compare feature performance at baseline and follow-up for osteolysis and normal cases.

RESULTS

Separation between the RTA features for osteolysis and normal cases was negligible at baseline and increased substantially for the follow-up images. The directional Fourier-based feature provided the best separation with an A(z) value from ROC analysis of 0.75 for the follow-up images, in the task of distinguishing between normal and osteolytic cases.

CONCLUSIONS

The results from this preliminary analysis indicate that qualitative changes in trabecular patterns from immediately after surgery to the eventual detection of osteolysis correspond to quantitative changes in RTA features. It therefore appears that RTA provides information that could potentially be useful to aid in the detection of this disease.