Anisotropy of Trabecular Bone from Ultra-Distal Radius Digital X-Ray Imaging: Effects on Bone Mineral Density and Age ()
1. Introduction
Osteoporosis is a bone disease characterized by low bone mineral density (BMD) and micro-architectural deterioration of bone tissue [1] . The BMD measurement by using dual energy X-ray absorptiometry (DXA) is defined by the World Health Organization (WHO) Working Group as the gold standard and is clinically used to diagnose osteoporosis [2] . BMD is clearly one of the major determinants of bone strength and fracture risk [3] [4] [5] , but the assessment of fracture risk by BMD sometime may still lack enough sensitivity [6] . It is clear that other factors in addition to BMD may account for bone strength and fracture risk. They could include bone micro-architecture, bone geometry, and other extra-skeletal factors [6] . The trabecular structure with specific reference to bone micro-architecture, known to affect bone fragility, may not be completely accounted for by standard DXA measurements, but alternative approaches, such as Trabecular Bone Score (TBS), have been proposed in order to overcome this limitation [6] [7] [8] . TBS originally was defined as the slope at the origin of the log-log representation of 2D projected image evaluated from least square regression line calculation, and is a novel texture parameter that evaluated pixel gray-level local variations in a DXA lumbar spine or human femur image and may directly related to bone micro-architecture and fragility fracture risk [6] , and might improve the prediction accuracy for major osteoporotic fractures in elderly people [8] . Conceptually, a dense trabecular network, associated with strong mechanical bone strength, produces a projection image with many fine gray-level texture variations of small amplitude and therefore a steep slope. In contrast, a loose trabecular network, associated with weak mechanical bone strength, produces a projection image with fewer fine gray-level texture variations of large amplitude and therefore a shallow slop. Furthermore, mechanical and structural properties of trabecular bone are anisotropic rather than isotropic in general [9] [10] , which has been observed through variation of the mechanical behavior of trabecular bone following specific testing directions. The anisotropy of trabecular bone corresponds to these variations and may be considered to predict the fracture risk. In our previous paper [11] , we only studied the initial Slope of variogram of digital X-Ray images along the direction of perpendicular to the forearm. In the current study, first, the values of trabecular bone score along the direction of perpendicular to the forearm, TBSx, and along the direction of parallel to the forearm, TBSy, were calculated respectively. The statistics of TBSx and TBSy were calculated and analyzed. The ratio of TBSy/TBSx was used to study the trabecular bone anisotropy. If the ratio deviates from one, trabecular bone will demonstrate greater anisotropy; otherwise, if the ratio is closer to one, it will demonstrate greater isotropy. This anisotropy changed with subjects’ BMD and age will be reported.
2. Materials and Methods
Study Subjects
In order to study the properties of TBSx and TBSy of trabecular bone located at ultra-distal radius, both radiographic and BMD data from a group of 98 male and 128 female adults, who were recruited and participated in the project of BMD measurements and osteoporosis assessments with radiographer at two hospitals in China, were studied retrospectively [11] . The participants were originally screened by detailed questionnaire, disease history, and physical examination. Their ages ranged from 23 to 87 years old, with an average age of 62 years. Participants were excluded if they had a history of forearm trauma and cancer.
A typical digital X-ray image of the ultra-distal radius, as shown in Figure 1, was obtained by using a commercial available digital radiographic system. These digital X-ray images could qualitatively show the image texture difference between normal and abnormal bones due to deteriorate microarchitecture of bone. As a typical setup for X-ray imaging of extremities, the acquisition parameters for these digital X-ray images were set as 55 kVp and 10 mAs (100 mA and 100 msec) with a specific SID (source-to-image distance) of 100 cm, and a X-ray focal spot of 0.6 mm. No any added X-ray filtration was used. According to the detector specifications, the digital image detector has 3 k × 3 k pixels matrix with a pixel size of 0.139 mm, and the image signal has a gray level depth of 12 bits. In order to reduce Heel effect of X-ray imaging, the region of interest (ROI) was put in the center of X-ray beam (around the original point of X-ray image).
Calculation of the Trabecular Bone Scores TBSx and TBSy
In our study, a digital 2D X-ray image signal,
, was used for both TBSx and TBSy calculation, where
were the spatial coordinates. First, the experimental variogram Vx(k) vs k is calculated by averaging the squared differences of the image signal values at several pairs of points with a specified separation, k, along the direction perpendicular to the forearm, as shown in Figure 2(a):
(1)
Figure 1. Representative image depicting the ultra-distal radius of subjects obtained using a digital X-ray imaging system.
where M × N is the total number of data pairs
and
within the region of interest (ROI). ROI starts at the bifurcation of the radius and ulna, and the coverage range is 1/12 length of radius, and is similar to the ROI of ultra-distal radius defined for GE Prodigy, and has an enclosed area within ultra-distal radius [12] , as shown in Figure 2(b). Similarly, the experimental variogram Vy(k) vs k is calculated by averaging the squared differences of image signal values at several pairs of points with a specified separation, k, along the direction of parallel to the forearm, as shown in Figure 2(a):
. (2)
Then, the values of trabecular bone scores, TBSx and TBSy, were calculated from the slopes at the origin of log-log representation of Vx(k) and Vy(k), evaluated from the slope of least-square regression line (by using a line of best fit) calculation over the four points of k = 1, 2, 3, and 4, or log10(k) = 0, 0.301, 0.477, and 0.602, as shown in Figure 3.
BMD Measurements with DXA
The measurement of BMD at ultra-distal radius was performed for all participants using a GE Lunar Prodigy DXA scanner (GE Healthcare, Madison, WI) running enCORE software version 13.31, which is designed to make DXA even more powerful [12] . During the study, the DXA scanner was calibrated with a GE Lunar aluminum spine phantom to avoid any affect from significant drift and/or shift on our study. For BMD measurements, each subject was asked to sit beside the scanning table in a chair without any arms or wheels. The subject was asked to keep the forearm immobile during scanning and the forearm was scanned parallel to the long axis of the table. The scanning range was from the center of wrist to the proximal forearm according to the operation manual provided by the manufacturer.
Figure 2. ROI used for our analysis starts from the reference point of the bifurcation of the radius and ulna, and the coverage range is 1/12 length of radius as shown in (a), and is similar to the ROI of ultra-distal radius area defined for GE Prodigy as shown in (b) [11] . Here x-axis is in the direction of perpendicular to the forearm, and y-axis is in the direction of parallel to the forearm.
Figure 3. The calculation of TBSx and TBSy of and as the slopes at the origin of the experimental variogram, where k is the separation (in pixels) between a pair of imaging position.
Statistical Analysis
Pearson correlation coefficients for TBSx and TBSy of trabecular bone located at ultra-distal radius were calculated. Person correlation coefficient is given by [13] .
(3)
where
is the Pearson correlation coefficient,
is the TBSx, and
is the TBSy for the subject i.
is the mean of
,
is the mean of
, and
is the total number of subjects. Then, both linear regressions and their R-squared values are also calculated. Here R-squared value shows how well the data fit the regression model.
3. Results
Participant Characteristics
Among 228 subjects, 43% were men (98/228) with an average age 57.7 ± 19.5 years, and 57% were women (130/228) with an average age 65.3 ± 13.3 years. The BMD of ultra-distal radius was 0.522 ± 0.088 g/cm2 for men, and 0.350 ± 0.103 g/cm2 for women. The average TBSx of ultra-distal radius was 1.081 ± 0.096 cm−1 in both men and women. The average TBSy of ultra-distal radius was 1.184 ± 0.071 cm−1 in men, and 1.163 ± 0.069 cm−1 in both men and women. The Pearson correlation coefficient between TBSx and TBSy for the trabecular bone located at ultra-distal radius is 0.72 (p < 0.001), as shown in Figure 4. The ratio of TBSy/TBSx has a slight positive correlation with the subject’s BMD, as shown in Figure 5, and a slight negative correlation with the subject’s age, as shown in Figure 6, and will approach to 1, when subject’s age increases, as shown in Figure 6. These
Figure 4. The correlation coefficient between TBSx and TBSy of trabecular bone located at ultra-distal radius, and linear regression lines with R-squared value = 0.52.
Figure 5. The ratio of TBSy/TBSx as a function of subjects’ BMD.
results suggest that healthier trabecular bone located at the ultra-distal radius will have greater anisotropy.
4. Discussion
Osteoporosis is defined as a systemic skeletal disease characterized by low bone mass density and micro-architectural degradation of bone tissue, with an increase in bone fragility to fracture. Current study shows that a decrease in BMD could often be accompanied by bone micro-architectural deterioration. Pothuaud
Figure 6. The ratio of TBSy/TBSx as a function of subjects’ age.
et al. introduced a parameter called trabecular bone score (TBS) which is a texture index that evaluates pixel gray-level variations of 2D projection images of the trabecular bone from human lumbar spines acquired by DXA scanners [6] , providing an indirect index of trabecular architecture. TBS is defined as the initial slope of the log-log transform of variogram, which describes the rate of changes of variance with respect to a lag distance. Significant correlations were observed between TBS and microarchitecture parameters (e.g. bone volume fraction and trabecular number) of trabecular bone from human femurs and lumbar spines [6] . A dense trabecular structure produces a 2D image with a large number of pixel value variations of small amplitude, and consequently a steep slope at the origin of the variogram and a high TBS value. Conversely, a 2D projection of deteriorated bone architecture produces an image with a low number of pixel value variations of high amplitude, and therefore a mild slope at the origin of the variogram and a low TBS. It was found that TBS correlated positively with the BMD and decreased with age [6] .
More studies found that trabecular bone microarchitecture showed the spatial organization and morphology of the trabecular network. Anisotropy of the trabecular bone may also be used to predict fracture risk [14] . Different methods have been used to study characterized the structure anisotropy on trabecular radiographs. A semi-quantitative index was applied to femoral neck radiographs [15] . This index is based on the presence of several arches of trabeculae in the femoral neck. Some arches preferentially disappear with age and osteoporosis, and counting of these arch systems can be used to determine that semi-quantitative index. Fractal analysis on bone radiographs was proposed [16] . The measured parameters provide a global assessment of the image irregularities as an indicator of roughness. The Maximum Likelihood Estimator is applied first on a specific direction. Then, the process is repeated in 36 directions in steps of 10 degrees. The orientation results can be represented by a polar diagram, which itself can fit an ellipse; the shape of the polar diagram can be characterized to determine textural anisotropy. Fast Fourier transforms (FFT) were applied to regions of interest in trabecular bone radiographs, with all FFT spectra involving horizontal and vertical components corresponding to longitudinal and transversal trabeculae respectively [10] .
In current paper, the directional trabecular bone scores, TBSx and TBSy, along the directions perpendicular and parallel to the forearm were calculated from X-ray images of trabeculae located at the ultra-distal radius. The ratio of TBSx and TBSy was used to evaluate the anisotropy of the trabecular bone.
The results of current studies show that anisotropy of TBS has the potential to be used quantitatively to assess age-related and/or pathological changes in the microstructures of ultra-distal radius, which may influence bone fragility. If the anisotropy of TBS proposed in this study was applied, it can be directly evaluation quantitatively on the projection X-ray images obtained in routine clinical practice. In additional to quantifying the microstructures of ultra-distal radius as shown in current study, the anisotropy of TBS can also be used to quantitatively assess the microstructures of other trabecular bones. The overall purpose of such a study is to understand whether radiographically determined trabecular bone micro-architecture can improve clinicians’ ability to estimate patients’ risk of osteoporotic fracture, beyond that obtained using bone mass measurements alone.
5. Limitations
This study has several limitations. First, the irregular shape of the trabecular bone located at ultra-distal radius may produce uneven bone image due to thickness variations along the projection direction in clinical settings. Second, all measurements of BMD reported in this study were obtained from a GE Lunar Prodigy DXA scanner, so our findings can not necessarily be extrapolated to scanners from other manufacturers due to the ultra-distal radius sites of different DXA scanners. The ROI may vary slightly from scanner to scanner. Third, standardized settings for X-ray image acquisition procedures need to be developed because the image signal received by the X-ray detector depends on the energy and intensity of X-rays. This should include the development of automatic compensation techniques related to digital detector’s characteristics, as well as methods to generate good reproducibility. In the current study, we only analyzed images of trabecular bone and didn’t process any images of cortical bone.
6. Conclusion
The results of this study demonstrate that healthier trabecular bone located at the ultra-distal radius will have greater anisotropy. Imaging analysis of the anisotropy of trabecular bone texture may potentially play an alternative role in the study of osteoporosis and bone health.