A Digital Model to Simulate Effects of Bone Architecture Variations on Texture at Spatial Resolutions of CT, HR-pQCT, and μCT Scanners

Commonality in the microarchitecture of trabecular bone: A preliminary study

Understanding the relationship between the microstructure and mechanical function of trabecular bone is critical for prediction and prevention of bone fragility fractures. However, a detailed understanding of the structural design of trabecular microarchitecture is still missing. This study hypothesized that there exists a commonality in the underlying probabilistic distributions of microstructural features of trabecular bones, whereas the microstructural differences among individuals are primarily describe by a set of scalar parameters. To test the hypothesis, twenty-three trabecular bone specimens were obtained from two anatomic locations (i.e., femoral neck and vertebral body) and a diverse group of seventeen donors of different age and sex. The number, size, spatial location, and orientation of individual plates and rods in the trabecular bone specimens were determined via volumetric decomposition of 3D μCT images using the Individual Trabecula Segmentation (ITS) technique. Then, m/n bootstrap Kolmogorov-Smirnov tests were performed to compare the normalized distributions of size, orientation, and spatial arrangement of trabecular plates and rods in the specimens. The results showed that 100% of the twenty-three normalized distributions of each microstructural feature were statistically equivalent irrespective of individual differences among the bone specimens, except the distributions of rod spatial arrangement (<100%). On the other hand, nonparametric Mann-Whitney U tests showed that a set of scalar parameters (i.e., the number, average size, and average nearest neighbor distance of trabecular plates and rods) were statistically different among the individual specimens (p<0.05). Due to the commonality of the underlying distributions, the individual differences in the trabecular microstructure among the specimens seemed to be reflected primarily by changes in the scalar parameters. The above results strongly support the hypothesis of this study and may shed more light on understanding the natural design of trabecular bone microstructures.

The robustness of T2 value as a trabecular structural index at multiple spatial resolutions of 7 Tesla MRI

PURPOSE

To evaluate the robustness of MR transverse relaxation times of trabecular bone from spin-echo and gradient-echo acquisitions at multiple spatial resolutions of 7 T.

METHODS

The effects of MRI resolutions to T₂ and T2* of trabecular bone were numerically evaluated by Monte Carlo simulations. T₂ , T2*, and trabecular structural indices from multislice multi-echo and UTE acquisitions were measured in defatted human distal femoral condyles on a 7 T scanner. Reference structural indices were extracted from high-resolution microcomputed tomography images. For bovine knee trabecular samples with intact bone marrow, T₂ and T2* were measured by degrading spatial resolutions on a 7 T system.

RESULTS

In the defatted trabecular experiment, both T₂ and T2* values showed strong ( |r| > 0.80) correlations with trabecular spacing and number, at a high spatial resolution of 125 µm³ . The correlations for MR image-segmentation-derived structural indices were significantly degraded ( |r| < 0.50) at spatial resolutions of 250 and 500 µm³ . The correlations for T2* rapidly dropped ( |r| < 0.50) at a spatial resolution of 500 µm³ , whereas those for T₂ remained consistently high ( |r| > 0.85). In the bovine trabecular experiments with intact marrow, low-resolution (approximately 1 mm³ , 2 minutes) T₂ values did not shorten ( |r| > 0.95 with respect to approximately 0.4 mm³ , 11 minutes) and maintained consistent correlations ( |r| > 0.70) with respect to trabecular spacing (turbo spin echo, 22.5 minutes).

CONCLUSION

T₂ measurements of trabeculae at 7 T are robust with degrading spatial resolution and may be preferable in assessing trabecular spacing index with reduced scan time, when high-resolution 3D micro-MRI is difficult to obtain.

Advanced Knee Structure Analysis (AKSA): a comparison of bone mineral density and trabecular texture measurements using computed tomography and high-resolution peripheral quantitative computed tomography of human knee cadavers

Background

A change of loading conditions in the knee causes changes in the subchondral bone and may be a cause of osteoarthritis (OA). However, quantification of trabecular architecture in vivo is difficult due to the limiting spatial resolution of the imaging equipment; one approach is the use of texture parameters. In previous studies, we have used digital models to simulate changes of subchondral bone architecture under OA progression. One major result was that, using computed tomography (CT) images, subchondral bone mineral density (BMD) in combination with anisotropy and global homogeneity could characterize this progression.

The primary goal of this study was a comparison of BMD, entropy, anisotropy, variogram slope, and local and global inhomogeneity measurements between high-resolution peripheral quantitative CT (HR-pQCT) and CT using human cadaveric knees. The secondary goal was the verification of the spatial resolution dependence of texture parameters observed in the earlier simulations, two important prerequisites for the interpretation of in vivo measurements in OA patients.

Method

The applicability of texture analysis to characterize bone architecture in clinical CT examinations was investigated and compared to results obtained from HR-pQCT. Fifty-seven human knee cadavers (OA status unknown) were examined with both imaging modalities. Three-dimensional (3D) segmentation and registration processes, together with automatic positioning of 3D analysis volumes of interest (VOIs), ensured the measurement of BMD and texture parameters at the same anatomical locations in CT and HR-pQCT datasets.

Results

According to the calculation of dice ratios (>0.978), the accuracy of VOI locations between methods was excellent. Entropy, anisotropy, and global inhomogeneity showed significant and high linear correlation between both methods (0.68 < R² < 1.00). The resolution dependence of these parameters simulated earlier was confirmed by the in vitro measurements.

Conclusion

The high correlation of HR-pQCT- and CT-based measurements of entropy, global inhomogeneity, and anisotropy suggests interchangeability between devices regarding the quantification of texture. The agreement of the experimentally determined resolution dependence of global inhomogeneity and anisotropy with earlier simulations is an important milestone towards their use to quantify subchondral bone structure. However, an in vivo study is still required to establish their clinical relevance.

Characterization of knee osteoarthritis-related changes in trabecular bone using texture parameters at various levels of spatial resolution-a simulation study

Articular cartilage and subchondral bone are the key tissues in osteoarthritis (OA). The role of the cancellous bone increasingly attracts attention in OA research. Because of its fast adaptation to changes in the loading distribution across joints, its quantification is expected to improve the diagnosis and monitoring of OA. In this study, we simulated OA progression-related changes of trabecular structure in a series of digital bone models and then characterized the potential of texture parameters and bone mineral density (BMD) as surrogate measures to quantify trabecular bone structure. Five texture parameters were studied: entropy, global and local inhomogeneity, anisotropy and variogram slope. Their dependence on OA relevant structural changes was investigated for three spatial resolutions typically used in micro computed tomography (CT; 10 μm), high-resolution peripheral quantitative CT (HR-pQCT) (90 μm) and clinical whole-body CT equipment (250 μm). At all resolutions, OA-related changes in trabecular bone architecture can be quantified using a specific (resolution dependent) combination of three texture parameters. BMD alone is inadequate for this purpose but if available reduces the required texture parameter combination to anisotropy and global inhomogeneity. The results are summarized in a comprehensive analysis guide for the detection of structural changes in OA knees. In conclusion, texture parameters can be used to characterize trabecular bone architecture even at spatial resolutions below the dimensions of a single trabecula and are essential for a detailed classification of relevant OA changes that cannot be achieved with a measurement of BMD alone.