Frequency dependence of apparent ultrasonic backscatter from human cancellous bone

Backscatter measurement of cancellous bone using the ultrasound transit time spectroscopy

Recently, ultrasound transit time spectroscopy (UTTS) was proposed as a promising method for bone quantitative ultrasound measurement. Studies have showed that UTTS could estimate the bone volume fraction and other trabecular bone structure in ultrasonic through-transmission measurements. The goal of this study was to explore the feasibility of UTTS to be adapted in ultrasonic backscatter measurement and further evaluate the performance of backscattered ultrasound transit time spectrum (BS-UTTS) in the measurement of cancellous bone density and structure. First, taking ultrasonic attenuation into account, the concept of BS-UTTS was verified on ultrasonic backscatter signals simulated from a set of scatterers with different positions and intensities. Then, in vitro backscatter measurements were performed on 26 bovine cancellous bone specimens. After a logarithmic compression of the BS-UTTS, a linear fitting of the log-compressed BS-UTTS versus ultrasonic propagated distance was performed and the slope and intercept of the fitted line for BS-UTTS were determined. The associations between BS-UTTS parameters and cancellous bone features were analyzed using simple linear regression. The results showed that the BS-UTTS could make an accurate deconvolution of the backscatter signal and predict the position and intensity of the simulated scatterers eliminating phase interference, even the simulated backscatter signal was with a relatively low signal-to-noise ratio. With varied positions and intensities of the scatterers, the slope of the fitted line for the log-compressed BS-UTTS versus ultrasonic propagated distance (i.e., slope of BS-UTTS for short) yield a high agreement (r2 = 99.84%-99.96%) with ultrasonic attenuation in simulated backscatter signal. Compared with the high-density cancellous bone, the low-density specimen showed more abundant backscatter impulse response in the BS-UTTS. The slope of BS-UTTS yield a significant correlation with bone mineral density (r = 0.87; p < 0.001), BV/TV (r = 0.87; p < 0.001), and cancellous bone microstructures (r up to 0.87; p < 0.05). The intercept of BS-UTTS was also significantly correlated with bone densities (r = -0.87; p < 0.001) and trabecular structures (|r|=0.43-0.80; p < 0.05). However, the slope of the BS-UTTS underestimated attenuation when measurements were performed experimentally. In addition, a significant non-linear relationship was observed between the measured attenuation and the attenuation estimated by the slope of the BS-UTTS. This study demonstrated that the UTTS method could be adapted to ultrasonic backscatter measurement of cancellous bone. The derived slope and intercept of BS-UTTS could be used in the measurement of bone density and microstructure. The backscattered ultrasound transit time spectroscopy might have potential in the diagnosis of osteoporosis in the clinic.

The Protocol of Ultrasonic Backscatter Measurements of Musculoskeletal Properties

This study aims to introduce the protocol for ultrasonic backscatter measurements of musculoskeletal properties based on a novel ultrasonic backscatter bone diagnostic (UBBD) instrument. Dual-energy X-ray absorptiometry (DXA) can be adopted to measure bone mineral density (BMD) in the hip, spine, legs and the whole body. The muscle and fat mass in the legs and the whole body can be also calculated by DXA body composition analysis. Based on the proposed protocol for backscatter measurements by UBBD, ultrasonic backscatter signals can be measured in vivo, deriving three backscatter parameters [apparent integral backscatter (AIB), backscatter signal peak amplitude (BSPA) and the corresponding arrival time (BSPT)]. AIB may provide important diagnostic information about bone properties. BSPA and BSPT may be important indicators of muscle and fat properties. The standardized backscatter measurement protocol of the UBBD instrument may have the potential to evaluate musculoskeletal characteristics, providing help for promoting the application of the backscatter technique in the clinical diagnosis of musculoskeletal disorders (MSDs), such as osteoporosis and muscular atrophy.

In Vivo Comparison of Backscatter Techniques for Ultrasonic Bone Assessment at the Femoral Neck

Ultrasonic techniques are being developed to detect changes in cancellous bone caused by osteoporosis. The goal of this study was to test the relative in vivo performance of eight backscatter parameters developed over the last several years for ultrasonic bone assessment: apparent integrated backscatter (AIB), frequency slope of apparent backscatter (FSAB), frequency intercept of apparent backscatter (FIAB), normalized mean of the backscatter difference (nMBD), normalized slope of the backscatter difference (nSBD), normalized intercept of the backscatter difference (nIBD), normalized backscatter amplitude ratio (nBAR) and backscatter amplitude decay constant (BADC). Backscatter measurements were performed on the left and right femoral necks of 80 adult volunteers (age = 25 ± 11 y) using an imaging system equipped with a convex array transducer. For comparison, additional ultrasonic measurements were performed at the left and right heel using a commercially available heel-bone ultrasonometer that measured the stiffness index. Six of the eight backscatter parameters (all but nSBD and nIBD) exhibited similar and highly significant (p < 0.000001) left-right correlations (0.51 ≤ R ≤ 0.68), indicating sensitivity to naturally occurring variations in bone tissue. Left-right correlations for the stiffness index measured at the heel (R = 0.75) were not significantly better than those produced by AIB, FSAB and FIAB. The short-term precisions of AIB, nMBD, nBAR and BADC (7.8%-11.7%) were comparable to that of the stiffness index measured with the heel-bone ultrasonometer (7.5%).

Ultrasonic Bone Assessment: Ability of Apparent Backscatter Techniques to Detect Changes in the Microstructure of Human Cancellous Bone

Ultrasonic backscatter techniques may offer a useful approach for detecting changes in bone caused by osteoporosis. The goal of this study was to investigate how bone mineral density (BMD) and the microstructure of human cancellous bone affect three ultrasonic backscatter parameters that have been identified as potentially useful for ultrasonic bone assessment purposes: the apparent integrated backscatter (AIB), the frequency slope of apparent backscatter (FSAB), and the frequency intercept of apparent backscatter (FIAB). Ultrasonic measurements were performed with a 3.5-MHz broadband transducer on 54 specimens of human cancellous bone prepared from the proximal femur. Microstructural parameters and BMD were measured using X-ray microcomputed tomography (micro-CT). Relationships between AIB, FSAB, FIAB, and the micro-CT parameters were investigated using univariate and multivariate statistical analysis techniques. Moderate-to-strong univariate correlations were observed between the backscatter parameters and microstructure and BMD in many cases. The partial correlation analysis indicated that the backscatter parameters are dependent on microstructure independently of BMD in some cases. Multiple stepwise linear regression analysis used to generate multivariate models found that microstructure was a significant predictor of the backscatter parameters in most cases.

Ultrasonic Bone Assessment Using the Backscatter Amplitude Decay Constant

Ultrasonic backscatter techniques are being developed to detect changes in bone caused by osteoporosis. The present study introduces a new technique that measures the exponential decay in the amplitude of the backscatter signal quantified by a parameter called the backscatter amplitude decay constant (BADC). Measurements were performed on 54 specimens of cancellous bone from 14 human femurs using a 3.5-MHz transducer. Six methods were tested to determine BADC. The recommended method measures the time slope of the natural log of the rectified signal. Measured values of BADC ranged from approximately 0.1 μs^-1 to 0.6 μs^-1. Moderate to strong correlations (Spearman's ρ >0.7) were found between BADC and the density and microstructural characteristics of the specimens determined using X-ray microcomputed tomography. The results of this study suggest that BADC may be able to detect changes in the density and microstructure of cancellous bone caused by osteoporosis and other diseases.

Relationships of the ultrasonic backscatter measurements with the bone mineral density and the microarchitectural parameters in bovine trabecular bone in vitro

Relationships of the backscatter coefficient (BC), the apparent integrated backscatter (AIB), and the integrated reflection coefficient (IRC) with the bone mineral density (BMD) and the microarchitectural parameters were investigated in 28 bovine femoral trabecular bone samples. The BC was highly correlated with the BMD and the microarchitectural parameters (R = -0.66 to 0.71). In contrast, the AIB and the IRC exhibited high correlations with the BMD and the bone volume fraction (R = -0.68 to 0.77) and relatively lower correlations with the remaining microarchitectural parameters (R = -0.62 to 0.60). The multiple regression models yielded the adjusted squared correlation coefficients of 0.54-0.76.

A Combined Ultrasonic Backscatter Parameter for Bone Status Evaluation in Neonates

Metabolic bone disease (MBD) is one of the major complications of prematurity. Ultrasonic backscatter technique has the potential to be a portable and noninvasive method for early diagnosis of MBD. This study firstly applied CAS to neonates, which was defined as a linear combination of the apparent integrated backscatter coefficient (AIB) and spectral centroid shift (SCS). The objective was to evaluate the feasibility of ultrasonic backscatter technique for assessing neonatal bone health using AIB, SCS, and CAS. Ultrasonic backscatter measurements at 3.5 MHz, 5.0 MHz, and 7.5 MHz were performed on a total of 505 newborns within 48 hours after birth. The values of backscatter parameters were calculated and compared among gestational age groups. Correlations between backscatter parameters, gestational age, anthropometric indices, and biochemical markers were analyzed. The optimal predicting models for CAS were determined. The results showed term infants had lower SCS and higher AIB and CAS than preterm infants. Gestational age and anthropometric indices were negatively correlated with SCS (|r| = 0.45 – 0.57, P < 0.001), and positively correlated with AIB (|r| = 0.36 – 0.60, P < 0.001) and CAS (|r| = 0.56 – 0.69, P < 0.001). Biochemical markers yielded weak or nonsignificant correlations with backscatter parameters. CAS had relatively stronger correlations with the neonatal variables than AIB and SCS. At 3.5 MHz and 5.0 MHz, only gestational age (P < 0.001) independently contributed to the measurements of CAS, and could explain up to 40.5% – 44.3% of CAS variation. At 7.5 MHz, the combination of gestational age (P < 0.001), head circumference (P = 0.002), and serum calcium (P = 0.037) explained up to 40.3% of CAS variation. This study suggested ultrasonic backscatter technique was feasible to evaluate neonatal bone status. CAS was a promising parameter to provide more information about bone health than AIB or SCS alone.

Mechanisms of Interaction of Ultrasound With Cancellous Bone: A Review

Ultrasound is now a clinically accepted modality in the management of osteoporosis. The most common commercial clinical devices assess fracture risk from measurements of attenuation and sound speed in cancellous bone. This review discusses fundamental mechanisms underlying the interaction between ultrasound and cancellous bone. Because of its two-phase structure (mineralized trabecular network embedded in soft tissue-marrow), its anisotropy, and its inhomogeneity, cancellous bone is more difficult to characterize than most soft tissues. Experimental data for the dependencies of attenuation, sound speed, dispersion, and scattering on ultrasound frequency, bone mineral density, composition, microstructure, and mechanical properties are presented. The relative roles of absorption, scattering, and phase cancellation in determining attenuation measurements in vitro and in vivo are delineated. Common speed of sound metrics, which entail measurements of transit times of pulse leading edges (to avoid multipath interference), are greatly influenced by attenuation, dispersion, and system properties, including center frequency and bandwidth. However, a theoretical model has been shown to be effective for correction for these confounding factors in vitro and in vivo. Theoretical and phantom models are presented to elucidate why cancellous bone exhibits negative dispersion, unlike soft tissue, which exhibits positive dispersion. Signal processing methods are presented for separating "fast" and "slow" waves (predicted by poroelasticity theory and supported in cancellous bone) even when the two waves overlap in time and frequency domains. Models to explain dependencies of scattering on frequency and mean trabecular thickness are presented and compared with measurements. Anisotropy, the effect of the fluid filler medium (marrow in vivo or water in vitro), phantoms, computational modeling of ultrasound propagation, acoustic microscopy, and nonlinear properties in cancellous bone are also discussed.

The Ability of Ultrasonic Backscatter Parametric Imaging to Characterize Bovine Trabecular Bone

The ultrasonic backscatter technique holds the promise of characterizing bone density and microstructure. This paper conducts ultrasonic backscatter parametric imaging based on measurements of apparent integrated backscatter (AIB), spectral centroid shift (SCS), frequency slope of apparent backscatter (FSAB), and frequency intercept of apparent backscatter (FIAB) for representing trabecular bone mass and microstructure. We scanned 33 bovine trabecular bone samples using a 7.5 MHz focused transducer in a 20 mm × 20 mm region of interest (ROI) with a step interval of 0.05 mm. Images based on the ultrasonic backscatter parameters (i.e., AIB, SCS, FSAB, and FIAB) were constructed to compare with photographic images of the specimens as well as two-dimensional (2D) μ-CT images from approximately the same depth and location of the specimen. Similar structures and trabecular alignments can be observed among these images. Statistical analyses demonstrated that the means and standard deviations of the ultrasonic backscatter parameters exhibited significant correlations with bone density (|R| = 0.45-0.78, p < 0.01) and bone microstructure (|R| = 0.44-0.87, p < 0.001). Some bovine trabecular bone microstructure parameters were independently associated with the ultrasonic backscatter parameters (ΔR² = 4.18%-44.45%, p < 0.05) after adjustment for bone apparent density (BAD). The results show that ultrasonic backscatter parametric imaging can provide a direct view of the trabecular microstructure and can reflect information about the density and microstructure of trabecular bone.

Relationships of Ultrasonic Backscatter With Bone Densities and Microstructure in Bovine Cancellous Bone

This study was designed to investigate the associations among ultrasonic backscatter, bone densities, and microstructure in bovine cancellous bone. Ultrasonic backscatter measurements were performed on 33 bovine cancellous bone specimens with a 2.25-MHz transducer. Ultrasonic apparent backscatter parameters ("apparent" means not compensating for ultrasonic attenuation and diffraction) were calculated with optimal signals of interest. The results showed that ultrasonic backscatter was significantly related to bone densities and microstructure ( R² = 0.17 -0.88 and ). After adjusting the correlations by bone mineral density (BMD), the bone apparent density (BAD) and some trabecular structural features still contributed significantly to the adjusted correlations, with moderate additional variance explained ( ∆R² = 9.7 % at best). Multiple linear regressions revealed that both BAD and trabecular structure contributed significantly and independently to the prediction of ultrasound backscatter (adjusted R² = 0.75 -0.89 and ), explaining an additional 14% of the variance at most, compared with that of BMD measurements alone. The results proved that ultrasonic backscatter was primarily determined by BAD, not BMD, but the combination of bone structure and densities could achieve encouragingly better performances (89% of the variance explained at best) in predicting backscatter properties. This study demonstrated that ultrasonic apparent backscatter might provide additional density and structural features unrelated to current BMD measurement. Therefore, we suggest that ultrasonic backscatter measurement could play a more important role in cancellous bone evaluation.

Ultrasonic backscatter difference measurements of cancellous bone from the human femur: Relation to bone mineral density and microstructure

Ultrasonic backscatter techniques are being developed to detect changes in cancellous bone caused by osteoporosis. One technique, called the backscatter difference technique, measures the power difference between two portions of a backscatter signal. The goal of the present study is to investigate how bone mineral density (BMD) and the microstructure of human cancellous bone influence four backscatter difference parameters: the normalized mean of the backscatter difference (nMBD) spectrum, the normalized slope of the backscatter difference spectrum, the normalized intercept of the backscatter difference spectrum, and the normalized backscatter amplitude ratio (nBAR). Ultrasonic measurements were performed with a 3.5 MHz broadband transducer on 54 specimens of human cancellous bone from the proximal femur. Volumetric BMD and the microstructural characteristics of the specimens were measured using x-ray micro-computed tomography. Of the four ultrasonic parameters studied, nMBD and nBAR demonstrated the strongest univariate correlations with density and microstructure. Multivariate analyses indicated that nMBD and nBAR depended on trabecular separation and possibly other microstructural characteristics of the specimens independently of BMD. These findings suggest that nMBD and nBAR may be sensitive to changes in the density and microstructure of bone caused by osteoporosis.

Characterization of a polymer, open-cell rigid foam that simulates the ultrasonic properties of cancellous bone

Materials that simulate the ultrasonic properties of tissues are used widely for clinical and research purposes. However, relatively few materials are known to simulate the ultrasonic properties of cancellous bone. The goal of the present study was to investigate the suitability of using a polymer, open-cell rigid foam (OCRF) produced by Sawbones^®. Measurements were performed on OCRF specimens with four different densities. Ultrasonic speed of sound and normalized broadband ultrasonic attenuation were measured with a 0.5 MHz transducer. Three backscatter parameters were measured with a 5 MHz transducer: apparent integrated backscatter, frequency slope of apparent backscatter, and normalized mean of the backscatter difference. X-ray micro-computed tomography was used to measure the microstructural characteristics of the OCRF specimens. The trabecular thickness and relative bone volume of the OCRF specimens were similar to those of human cancellous bone, but the trabecular separation was greater. In most cases, the ultrasonic properties of the OCRF specimens were similar to values reported in the literature for cancellous bone, including dependence on density. In addition, the OCRF specimens exhibited an ultrasonic anisotropy similar to that reported for cancellous bone.

Relationships among ultrasonic and mechanical properties of cancellous bone in human calcaneus in vitro

Clinical bone sonometers applied at the calcaneus measure broadband ultrasound attenuation and speed of sound. However, the relation of ultrasound measurements to bone strength is not well-characterized. Addressing this issue, we assessed the extent to which ultrasonic measurements convey in vitro mechanical properties in 25 human calcaneal cancellous bone specimens (approximately 2×4×2cm). Normalized broadband ultrasound attenuation, speed of sound, and broadband ultrasound backscatter were measured with 500kHz transducers. To assess mechanical properties, non-linear finite element analysis, based on micro-computed tomography images (34-micron cubic voxel), was used to estimate apparent elastic modulus, overall specimen stiffness, and apparent yield stress, with models typically having approximately 25-30 million elements. We found that ultrasound parameters were correlated with mechanical properties with R=0.70-0.82 (p<0.001). Multiple regression analysis indicated that ultrasound measurements provide additional information regarding mechanical properties beyond that provided by bone quantity alone (p≤0.05). Adding ultrasound variables to linear regression models based on bone quantity improved adjusted squared correlation coefficients from 0.65 to 0.77 (stiffness), 0.76 to 0.81 (apparent modulus), and 0.67 to 0.73 (yield stress). These results indicate that ultrasound can provide complementary (to bone quantity) information regarding mechanical behavior of cancellous bone.

Effect of gate choice on backscatter difference measurements of cancellous bone

A variety of ultrasonic techniques have been developed to detect changes in bone caused by osteoporosis. One approach, called the backscatter difference technique, analyzes the power difference between two different portions of a backscatter signal. Analysis gates with a certain delay τ_d, width τ_w, and separation τ_s are used to define portions of the backscatter signal for analysis. The goal of the present study was to investigate how different choices of τ_d, τ_w, and τ_s affect four backscatter difference parameters: the normalized mean of the backscatter difference (nMBD), the normalized slope of the backscatter difference (nSBD), the normalized intercept of the backscatter difference (nIBD), and the normalized backscatter amplitude ratio (nBAR). Backscatter measurements were performed on 54 cube shaped specimens of human cancellous bone. nMBD, nSBD, nIBD, and nBAR were determined for 34 different combinations of τ_d, τ_w, and τ_s for each specimen. nMBD and nBAR demonstrated the strongest correlations with apparent bone density (0.48 ≤ R_s ≤ 0.90). Generally, the correlations were found to improve as τ_w + τ_s was increased and as τ_d was decreased. Among the four backscatter difference parameters, the measured values of nMBD were least sensitive to gate choice (<16%).

Numerical evaluation of the backward propagating acoustic field in healing long bones

The propagation of ultrasound in healing long bones induces complex scattering phenomena due to the interaction of an ultrasonic wave with the composite nature of callus and osseous tissues. This work presents numerical simulations of ultrasonic propagation in healing long bones using the boundary element method aiming to provide insight into the complex scattering mechanisms and better comprehend the state of bone regeneration. Numerical models of healing long bones are established based on scanning acoustic microscopy images from successive postoperative weeks considering the effect of the nonhomogeneous callus structure. More specifically, the scattering amplitude and the acoustic pressure variation are calculated in the backward direction to investigate their potential to serve as quantitative and qualitative indicators for the monitoring of the bone healing process. The role of the excitation frequency is also examined considering frequencies in the range 0.2-1 MHz. The results indicate that the scattering amplitude decreases at later stages of healing compared to earlier stages of healing. Also, the acoustic pressure could provide supplementary qualitative information on the interaction of the scattered energy with bone and callus.

Ultrasonic backscatter from cancellous bone: the apparent backscatter transfer function

Ultrasonic backscatter techniques are being developed to detect changes in cancellous bone caused by osteoporosis. Many techniques are based on measurements of the apparent backscatter transfer function (ABTF), which represents the backscattered power from bone corrected for the frequency response of the measurement system. The ABTF is determined from a portion of the backscatter signal selected by an analysis gate of width τw delayed by an amount τd from the start of the signal. The goal of this study was to characterize the ABTF for a wide range of gate delays (1 μs ≤ τd ≤ 6 μs) and gate widths (1 μs ≤ τw ≤ 6 μs). Measurements were performed on 29 specimens of human cancellous bone in the frequency range 1.5 to 6.0 MHz using a broadband 5-MHz transducer. The ABTF was found to be an approximately linear function of frequency for most choices of τd and τw. Changes in τd and τw caused the frequency-averaged ABTF [quantified by apparent integrated backscatter (AIB)] and the frequency dependence of the ABTF [quantified by frequency slope of apparent backscatter (FSAB)] to change by as much as 24.6 dB and 6.7 dB/MHz, respectively. τd strongly influenced the measured values of AIB and FSAB and the correlation of AIB with bone density (-0.95 ≤ R ≤ +0.68). The correlation of FSAB with bone density was influenced less strongly by τd (-0.97 ≤ R ≤ -0.87). τw had a weaker influence than τd on the measured values of AIB and FSAB and the correlation of these parameters with bone density.

Effect of intervening tissues on ultrasonic backscatter measurements of bone: An in vitro study

Ultrasonic backscatter techniques are being developed to diagnose osteoporosis. Tissues that lie between the transducer and the ultrasonically interrogated region of bone may produce errors in backscatter measurements. The goal of this study is to investigate the effects of intervening tissues on ultrasonic backscatter measurements of bone. Measurements were performed on 24 cube shaped specimens of human cancellous bone using a 5 MHz transducer. Measurements were repeated after adding a 1 mm thick plate of cortical bone to simulate the bone cortex and a 3 cm thick phantom to simulate soft tissue at the hip. Signals were analyzed to determine three apparent backscatter parameters (apparent integrated backscatter, frequency slope of apparent backscatter, and frequency intercept of apparent backscatter) and three backscatter difference parameters [normalized mean backscatter difference (nMBD), normalized slope of the backscatter difference, and normalized intercept of the backscatter difference]. The apparent backscatter parameters were impacted significantly by the presence of intervening tissues. In contrast, the backscatter difference parameters were not affected by intervening tissues. However, only one backscatter difference parameter, nMBD, demonstrated a strong correlation with bone mineral density. Thus, among the six parameters tested, nMBD may be the best choice for in vivo backscatter measurements of bone when intervening tissues are present.

Signal of Interest Selection Standard for Ultrasonic Backscatter in Cancellous Bone Evaluation

The aim of this study was to examine the effect of the backscattered signal of interest (SOI) on ultrasonic cancellous bone evaluation. In vitro backscatter measurements were performed using 16 bovine cancellous bone specimens and six different transducers with central frequencies of 0.5, 1, 2.25, 3.5, 5 and 10 MHz. The SOI for signal analysis was selected by a rectangular window. The delay (T1) and duration (T2) of the time window were varied, and the apparent integrated backscatter (AIB) and its correlation to bone volume fraction (BV/TV) were calculated. The results indicate that in addition to affecting the measured value of AIB, the SOI influences the observed correlation between AIB and BV/TV. Strong positive correlations were observed for short T1 (0.5 MHz: ≤6 μs, 1 MHz: ≤3 μs, 2.25 and 3.5 MHz: ≤2 μs, 5 and 10 MHz: ≤1 μs). However, strong negative correlations were observed when T1 was long (0.5 MHz: >9 μs, 1 MHz: >7 μs, 2.25 and 3.5 MHz: >3 μs, 5 and 10 MHz: >2 μs). The T2 value, especially low values (≤3 μs), also influenced the correlation coefficients. Positive correlations were more commonly observed at lower frequencies (i.e., 0.5-1 MHz), whereas negative correlations were more common at higher frequencies (i.e., 2.25-10 MHz). An explicit standard for in vitro SOI selection and cancellous bone assessment was proposed for a broad frequency range (0.5-10 MHz). Current conflicting findings are explained, and constructive suggestions for ultrasonic backscatter cancellous bone evaluation are provided.

Photoacoustic and ultrasound imaging of cancellous bone tissue

We used ultrasound (US) and photoacoustic (PA) imaging modalities to characterize cattle trabecular bones. The PA signals were generated with an 805-nm continuous wave laser used for optimally deep optical penetration depth. The detector for both modalities was a 2.25-MHz US transducer with a lateral resolution of ~1 mm at its focal point. Using a lateral pixel size much larger than the size of the trabeculae, raster scanning generated PA images related to the averaged values of the optical and thermoelastic properties, as well as density measurements in the focal volume. US backscatter yielded images related to mechanical properties and density in the focal volume. The depth of interest was selected by time-gating the signals for both modalities. The raster scanned PA and US images were compared with microcomputed tomography (μCT) images averaged over the same volume to generate similar spatial resolution as US and PA. The comparison revealed correlations between PA and US modalities with the mineral volume fraction of the bone tissue. Various features and properties of these modalities such as detectable depth, resolution, and sensitivity are discussed.

The application of backscattered ultrasound and photoacoustic signals for assessment of bone collagen and mineral contents

BACKGROUND

This study examines the backscattered ultrasound (US) and back-propagating photoacoustic (PA) signals from trabecular bones, and their variations with reduction in bone minerals and collagen content. While the collagen status is directly related to the strength of the bone, diagnosis of its condition using US remains a challenge.

METHODS

For both PA and US methods, coded-excitation signals and matched filtering were utilized to provide high sensitivity of the detected signal. The optical source was a 805-nm CW laser and signals were detected employing a 2.2-MHz ultrasonic transducer. Bone decalcification and decollagenization were induced with mild ethylenediaminetetraacetic acid (EDTA) and sodium hypochlorite solutions, respectively.

RESULTS

The PA and US signals were measured on cattle bones, and apparent integrated backscatter/back-propagating (AIB) parameters were compared before and after demineralization and decollagenization.

CONCLUSIONS

The results show that both PA and US are sensitive to mineral changes. In addition, PA is also sensitive to changes in the collagen content of the bone, but US is not significantly sensitive to these changes.

A novel ultrasound methodology for estimating spine mineral density

We investigated the possible clinical feasibility and accuracy of an innovative ultrasound (US) method for diagnosis of osteoporosis of the spine. A total of 342 female patients (aged 51-60 y) underwent spinal dual X-ray absorptiometry and abdominal echographic scanning of the lumbar spine. Recruited patients were subdivided into a reference database used for US spectral model construction and a study population for repeatability and accuracy evaluation. US images and radiofrequency signals were analyzed via a new fully automatic algorithm that performed a series of spectral and statistical analyses, providing a novel diagnostic parameter called the osteoporosis score (O.S.). If dual X-ray absorptiometry is assumed to be the gold standard reference, the accuracy of O.S.-based diagnoses was 91.1%, with k = 0.859 (p < 0.0001). Significant correlations were also found between O.S.-estimated bone mineral densities and corresponding dual X-ray absorptiometry values, with r(2) values up to 0.73 and a root mean square error of 6.3%-9.3%. The results obtained suggest that the proposed method has the potential for future routine application in US-based diagnosis of osteoporosis.

Ultrasound backscatter measurements of intact human proximal femurs--relationships of ultrasound parameters with tissue structure and mineral density

Ultrasound reflection and backscatter parameters are related to the mechanical and structural properties of bone in vitro. However, the potential of ultrasound reflection and backscatter measurements has not been tested with intact human proximal femurs ex vivo. We hypothesize that ultrasound backscatter can be measured from intact femurs and that the measured backscattered signal is associated with cadaver age, bone mineral density (BMD) and trabecular bone microstructure. In this study, human femoral bones of 16 male cadavers (47.0±16.1 years, range: 21-77 years) were investigated using pulse-echo ultrasound measurements at the femoral neck in the antero-posterior direction and at the trochanter major in the anteroposterior and lateromedial directions. Recently introduced ultrasound backscatter parameters, independent of cortical thickness, e.g., time slope of apparent integrated backscatter (TSAB) and mean of the backscatter difference technique (MBD) were obtained and compared with the structural properties of trabecular bone samples, extracted from the locations of ultrasound measurements. Moreover, more conventional backscatter parameters, e.g., apparent integrated backscatter (AIB) and frequency slope of apparent integrated backscatter (FSAB) were analyzed. Bone mineral density of the intact femurs was evaluated using dual energy X-ray absorptiometry (DXA). AIB and MDB measured from the femoral neck correlated significantly (p<0.01) with the neck BMD (R2=0.44 and 0.45), cadaver age (R2=0.61 and 0.41) and several structural parameters, e.g., bone volume fraction (R2=0.33 and 0.39, p<0.05 and p<0.01), respectively. To conclude, ultrasound backscatter parameters, measured from intact proximal femurs, are significantly related (p<0.05) to structural properties and mineral density of trabecular bone.

Time-domain separation of interfering waves in cancellous bone using bandlimited deconvolution: simulation and phantom study

In through-transmission interrogation of cancellous bone, two longitudinal pulses ("fast" and "slow" waves) may be generated. Fast and slow wave properties convey information about material and micro-architectural characteristics of bone. However, these properties can be difficult to assess when fast and slow wave pulses overlap in time and frequency domains. In this paper, two methods are applied to decompose signals into fast and slow waves: bandlimited deconvolution and modified least-squares Prony's method with curve-fitting (MLSP + CF). The methods were tested in plastic and Zerdine(®) samples that provided fast and slow wave velocities commensurate with velocities for cancellous bone. Phase velocity estimates were accurate to within 6 m/s (0.4%) (slow wave with both methods and fast wave with MLSP + CF) and 26 m/s (1.2%) (fast wave with bandlimited deconvolution). Midband signal loss estimates were accurate to within 0.2 dB (1.7%) (fast wave with both methods), and 1.0 dB (3.7%) (slow wave with both methods). Similar accuracies were found for simulations based on fast and slow wave parameter values published for cancellous bone. These methods provide sufficient accuracy and precision for many applications in cancellous bone such that experimental error is likely to be a greater limiting factor than estimation error.

A backscatter difference technique for ultrasonic bone assessment

Ultrasonic backscatter techniques may offer a useful approach for detecting changes in cancellous bone caused by osteoporosis and other diseases. The goal of this study was to investigate the utility of a backscatter difference technique for ultrasonic bone assessment. Measurements were performed on 22 cube-shaped specimens of human cancellous bone using four broadband transducers with center frequencies 2.25, 5, 7.5, and 10 MHz. The backscatter difference spectrum D(f) was obtained by subtracting power spectra (in dB) from two different portions of the same backscatter signal. D(f) was found to be a monotonically increasing, quasi-linear function of frequency when averaged over multiple measurement sites on multiple specimens. The frequency slope of D(f) demonstrated weak to moderate correlations with specimen density (R = 0.21-0.80). The frequency averaged mean of D(f) demonstrated moderate to good correlations with density (R = 0.70-0.95). These results suggest that parameters based on the frequency averaged mean of the backscatter difference spectrum may be useful for bone assessment purposes.

Multi-site bone ultrasound measurements in elderly women with and without previous hip fractures

About 75% of patients suffering from osteoporosis are not diagnosed. This study describes a multi-site bone ultrasound method for osteoporosis diagnostics. In comparison with axial dual energy X-ray absorptiometry (DXA), the ultrasound method showed good diagnostic performance and could discriminate fracture subjects among elderly females.

INTRODUCTION

Axial DXA, the gold standard diagnostic method for osteoporosis, predicts fractures only moderately. At present, no reliable diagnostic methods are available at the primary health care level. Here, a multi-site ultrasound method is proposed for osteoporosis diagnostics.

METHODS

Thirty elderly women were examined using the ultrasound backscatter measurements in proximal femur, proximal radius, proximal and distal tibia in vivo. First, we predicted the areal bone mineral density (BMD) at femoral neck by ultrasound measurements in tibia combined with specific subject characteristics (density index, DI) and, second, we tested the ability of ultrasound backscatter measurements at proximal femur to discriminate between individuals with previously fractured hips from those without fractures. Areal BMD was determined by axial DXA.

RESULTS

Combined ultrasound parameters, cortical thickness at distal and proximal tibia, with age and weight of the subject, provided a significant estimate of BMD(neck) (r = 0.86, p < 0.001, n = 30). When inserted into FRAX (World Health Organization fracture risk assessment tool), the DI indicated the same treatment proposal as the BMD(neck) with 86% sensitivity and 100% specificity. The receiver operating characteristic analyses, with a combination of ultrasound parameters and patient characteristics, discriminated fracture subjects from the controls similarly as the model combining BMD(neck) and patient characteristics.

CONCLUSIONS

For the first time, ultrasound backscatter measurements of proximal femur were conducted in vivo. The results indicate that ultrasound parameters, combined with patient characteristics, may provide a means for osteoporosis diagnostics.

Relationships of quantitative ultrasound parameters with cancellous bone microstructure in human calcaneus in vitro

Ultrasound parameters (attenuation, phase velocity, and backscatter), bone mineral density (BMD), and microarchitectural features were measured on 29 human cancellous calcaneus samples in vitro. Regression analysis was performed to predict ultrasound parameters from BMD and microarchitectural features. The best univariate predictors of the ultrasound parameters were the indexes of bone quantity: BMD and bone volume fraction (BV/TV). The most predictive univariate models for attenuation, phase velocity, and backscatter coefficient yielded adjusted squared correlation coefficients of 0.69-0.73. Multiple regression models yielded adjusted correlation coefficients of 0.74-0.83. Therefore attenuation, phase velocity, and backscatter are primarily determined by bone quantity, but multiple regression models based on bone quantity plus microarchitectural features achieve slightly better predictive performance than models based on bone quantity alone.

Effect of the cortex on ultrasonic backscatter measurements of cancellous bone

Ultrasonic backscatter techniques offer a promising new approach for detecting changes in bone caused by osteoporosis. However, several challenges impede clinical implementation of backscatter techniques. This study examines how the dense outer surface of bone (the cortex) affects backscatter measurements of interior regions of porous (cancellous) bone tissue. Fifty-two specimens of bone were prepared from 13 human femoral heads so that the same region of cancellous bone could be ultrasonically interrogated through the cortex or along directions that avoided the cortex. Backscatter signals were analyzed over a frequency range of 0.8-3.0 MHz to determine two ultrasonic parameters: apparent integrated backscatter (AIB) and frequency slope of apparent backscatter (FSAB). The term 'apparent' means that the parameters are sensitive to the frequency-dependent effects of diffraction and attenuation. Significant (p < 0.001) changes in AIB and FSAB indicated that measurements through the cortex decreased the apparent backscattered power and increased the frequency dependence of the power. However, the cortex did not affect the correlation of AIB and FSAB with the x-ray bone mineral density of the specimens. This suggests that results from many previous in vitro backscatter studies of specimens of purely cancellous bone may be extrapolated with greater confidence to in vivo conditions.