Trabecular homogeneity index derived from plain radiograph to evaluate bone quality

Quantitative Skeletal Imaging and Image-Based Modeling in Pediatric Orthopaedics

PURPOSE OF REVIEW

Musculoskeletal imaging serves a critical role in clinical care and orthopaedic research. Image-based modeling is also gaining traction as a useful tool in understanding skeletal morphology and mechanics. However, there are fewer studies on advanced imaging and modeling in pediatric populations. The purpose of this review is to provide an overview of recent literature on skeletal imaging modalities and modeling techniques with a special emphasis on current and future uses in pediatric research and clinical care.

RECENT FINDINGS

While many principles of imaging and 3D modeling are relevant across the lifespan, there are special considerations for pediatric musculoskeletal imaging and fewer studies of 3D skeletal modeling in pediatric populations. Improved understanding of bone morphology and growth during childhood in healthy and pathologic patients may provide new insight into the pathophysiology of pediatric-onset skeletal diseases and the biomechanics of bone development. Clinical translation of 3D modeling tools developed in orthopaedic research is limited by the requirement for manual image segmentation and the resources needed for segmentation, modeling, and analysis. This paper highlights the current and future uses of common musculoskeletal imaging modalities and 3D modeling techniques in pediatric orthopaedic clinical care and research.

Methods for bone quality assessment in human bone tissue: a systematic review

Background

For biomechanical investigations on bone or bone implants, bone quality represents an important potential bias. Several techniques for assessing bone quality have been described in the literature. This study aims to systematically summarize the methods currently available for assessing bone quality in human bone tissue, and to discuss the advantages and limitations of these techniques.

Methods

A systematic review of the literature was carried out by searching the PubMed and Web of Science databases from January 2000 to April 2021. References will be screened and evaluated for eligibility by two independent reviewers as per PRISMA (Preferred Reporting Items for Systematic Reviews and Meta-Analyses) guidelines. Studies must apply to bone quality assessment with imaging techniques, mechanical testing modalities, and compositional characterization. The terms used for the systematic search were: “(bone quality”. Ti,ab.) AND “(human bone specimens)”.

Results

The systematic review identified 502 relevant articles in total. Sixty-eight articles met the inclusion criteria. Among them, forty-seven articles investigated several imaging modalities, including radiography, dual-energy X-ray absorptiometry (DEXA), CT-based techniques, and MRI-based methods. Nineteen articles dealt with mechanical testing approaches, including traditional testing modalities and novel indentation techniques. Nine articles reported the correlation between bone quality and compositional characterization, such as degree of bone mineralization (DBM) and organic composition. A total of 2898 human cadaveric bone specimens were included.

Conclusions

Advanced techniques are playing an increasingly important role due to their multiple advantages, focusing on the assessment of bone morphology and microarchitecture. Non-invasive imaging modalities and mechanical testing techniques, as well as the assessment of bone composition, need to complement each other to provide comprehensive and ideal information on the bone quality of human bone specimens.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13018-022-03041-4.

Discrimination of Low-Energy Acetabular Fractures from Controls Using Computed Tomography-Based Bone Characteristics

The incidence of low-energy acetabular fractures has increased. However, the structural factors for these fractures remain unclear. The objective of this study was to extract trabecular bone architecture and proximal femur geometry (PFG) measures from clinical computed tomography (CT) images to (1) identify possible structural risk factors of acetabular fractures, and (2) to discriminate fracture cases from controls using machine learning methods. CT images of 107 acetabular fracture subjects (25 females, 82 males) and 107 age-gender matched controls were examined. Three volumes of interest, one at the acetabulum and two at the femoral head, were extracted to calculate bone volume fraction (BV/TV), gray-level co-occurrence matrix and histogram of the gray values (GV). The PFG was defined by neck shaft angle and femoral neck axis length. Relationships between the variables were assessed by statistical mean comparisons and correlation analyses. Bayesian logistic regression and Elastic net machine learning models were implemented for classification. We found lower BV/TV at the femoral head (0.51 vs. 0.55, p = 0.012) and lower mean GV at both the acetabulum (98.81 vs. 115.33, p < 0.001) and femoral head (150.63 vs. 163.47, p = 0.005) of fracture subjects when compared to their matched controls. The trabeculae within the femoral heads of the acetabular fracture sides differed in structure, density and texture from the corresponding control sides of the fracture subjects. Moreover, the PFG and trabecular architectural variables, alone and in combination, were able to discriminate fracture cases from controls (area under the receiver operating characteristics curve 0.70 to 0.79). In conclusion, lower density in the acetabulum and femoral head with abnormal trabecular structure and texture at the femoral head, appear to be risk factors for low-energy acetabular 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.

Does image normalization and intensity resolution impact texture classification?

Image texture is a very important component in many types of images, including medical images. Medical images are often corrupted by noise and affected by artifacts. Some of the texture-based features that should describe the structure of the tissue under examination may also reflect, for example, the uneven sensitivity of the scanner within the tissue region. This in turn may lead to an inappropriate description of the tissue or incorrect classification. To limit these phenomena, the analyzed regions of interest are normalized. In texture analysis methods, image intensity normalization is usually followed by a reduction in the number of levels coding the intensity. The aim of this work was to analyze the impact of different image normalization methods and the number of intensity levels on texture classification, taking into account noise and artifacts related to uneven background brightness distribution. Analyses were performed on four sets of images: modified Brodatz textures, kidney images obtained by means of dynamic contrast-enhanced magnetic resonance imaging, shoulder images acquired as T2-weighted magnetic resonance images and CT heart and thorax images. The results will be of use for choosing a particular method of image normalization, based on the types of noise and distortion present in the images.

Enhancing hip fracture risk prediction by statistical modeling and texture analysis on DXA images

Each year in the US more than 300,000 older adults suffer from hip fractures. While protective measures exist, identification of those at greatest risk by DXA scanning has proved inadequate. This study proposed a new technique to enhance hip fracture risk prediction by accounting for many contributing factors to the strength of the proximal femur. Twenty-two isolated cadaveric femurs were DXA scanned, 16 of which had been mechanically tested to failure. A function consisting of the calculated modes from the statistical shape and appearance modeling (to consider the shape and BMD distribution), homogeneity index (representing trabecular quality), BMD, age and sex of the donor was created in a training set and used to predict the fracture load in a test group. To classify patients as "high risk" or "low risk", fracture load thresholds were investigated. Hip fracture load estimation was significantly enhanced using the new technique in comparison to using t-score or BMD alone (average R² of 0.68, 0.32, and 0.50, respectively) (P < 0.05). Using a fracture cut-off of 3400 N correctly predicted risk in 94% of specimens, a substantial improvement over t-score classification (38%). Ultimately, by identifying patients at high risk more accurately, devastating hip fractures can be prevented through applying protective measures.

Structural risk factors for low-energy acetabular fractures

In this study, we aimed to clarify proximal femur and acetabular structural risk factors associated with low-energy acetabular fractures in the elderly using three-dimensional (3D) computed tomography (CT). Pelvic bones and femurs were segmented and modeled in 3D from abdominopelvic CT images of 121 acetabular fracture patients (mean age 72 ± 12 years, range 50-98 years, 31 females and 90 males) and 121 age-gender matched controls with no fracture. A set of geometric parameters of the proximal femur and the acetabulum was measured. An independent-samples t-test or a Mann-Whitney U test was used for statistical analyses. The fractured side was used for proximal femur geometry, while the contralateral side was used for acetabular geometry. The neck shaft angle (NSA) was significantly smaller (mean 122.1° [95% CI 121.1°-123.2°] vs. 124.6° [123.6°-125.6°], p = 0.001) and the femoral neck axis length (FNALb) was significantly longer (78.1 mm [77.0-79.2 mm] vs. 76.0 mm [74.8-77.2 mm], p = 0.026) in the fracture group than in the controls when genders were combined. The NSA was significantly smaller both for females (120.2° [117.8°-122.6°] vs. 124.7° [122.5°-127.0°], p = 0.007) and for males (122.7° [121.5°-123.8°] vs. 124.6° [123.4°-125.7°], p = 0.006) in the fracture group. However, only males showed a significantly longer FNALb (80.0 mm [78.9-81.1 mm] vs. 77.8 mm [76.6-79.0 mm], p = 0.025). No statistically significant associations of acetabular geometry with fractures were found. However, the mean values of the acetabular angle of Sharp (34°), the lateral center-edge angle (40°), the anterior center-edge angle (62°), and the posterior center-edge angle (105°) indicated possible over-coverage. In conclusion, our findings suggest that proximal femur geometry is associated with low-energy acetabular fractures. Especially elderly subjects with an NSA smaller than normal have an increased risk of acetabular fractures.

Technologies for assessment of bone reflecting bone strength and bone mineral density in elderly women: an update

Reduced bone mineral density is a strong risk factor for fracture. The WHO's definition of osteoporosis is based on bone mineral density measurements assessed by dual x-ray absorptiometry. Several on other techniques than dual x-ray absorptiometry have been developed for quantitative assessment of bone, for example, quantitative ultrasound and digital x-ray radiogrammetry. Some of these techniques may also capture other bone properties than bone mass that contribute to bone strength, for example, bone porosity and microarchitecture. In this article we give an update on technologies which are available for evaluation primarily of bone mass and bone density, but also describe methods which currently are validated or are under development for quantitative assessment of other bone properties.

Quantification of differences in bone texture from plain radiographs in knees with and without osteoarthritis

OBJECTIVE

To quantify differences in bone texture between subjects with different stages of knee osteoarthritis (OA) and age- and gender-matched controls from plain radiographs using advanced image analysis methods.

DESIGN

Altogether 203 knees were imaged using constant X-ray parameters and graded according to Kellgren-Lawrence (KL) grading scale (KL0: n = 110, KL1: n = 28, KL2: n = 27, KL3: n = 31, KL4: n = 7). Bone density-related and structure-related parameters were calculated from medial and lateral tibial subchondral bone plate and trabecular bone and from femur. Density-related parameters were derived from grayscale values and structure-related parameters from Laplacian- and local binary patterns (LBP)-based images.

RESULTS

Reproducibilities of structure-related parameters were better than bone density-related parameters. Bone density-related parameters were significantly (P < 0.05) higher in KL2-4 groups than in control group (KL0) in medial tibial subchondral bone plate and trabecular bone. LBP-based structure parameters differed significantly between KL0 and KL2-4 groups in medial subchondral bone plate, between KL0 and KL1-4 groups in medial and lateral trabecular bone, and between KL0 and KL1-4/KL2-4 in medial and lateral femur. Laplacian-based parameters differed significantly between KL0 and KL2-4 groups in medial side regions-of-interest (ROIs).

CONCLUSIONS

Our results indicate that the changes in bone texture in knee OA can be quantitatively evaluated from plain radiographs using advanced image analysis. Based on the results, increased bone density can be directly estimated if the X-ray imaging conditions are constant between patients. However, structural analysis of bone was more reproducible than direct evaluation of grayscale values, and is therefore better suited for quantitative analysis when imaging conditions are variable.

Measuring bone quality

Osteoporosis is defined as a reduction in bone mass and impairment of bone quality that lead to bone fragility and fracture risk. Bone quality includes a hierarchy of properties from macroscopic to nanoscale level. Several techniques have been developed in an attempt to measure these non-density properties. Densitometry, high-resolution images (radiography, CT scan), and MRI can measure the geometry and microarchitecture of bone. Tissue mineralization and composition can be assessed by use of microradiography, Fourier-transform infrared spectroscopy, or Raman microspectroscopy. Finite-element analysis is an image-based method that enables calculation of bone strength. More recently, microindentation has enabled direct estimation of bone material strength, measured in the cortical bone of the tibia. Most of these techniques are of limited use to clinics, although finite-element analysis and microindentation have high potential for clinical use and can enable more comprehensive and accurate evaluation of bone fragility and fracture susceptibility.

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.