Assessment of bone density using ultrasonic backscatter

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.

Biomechanical structure-function relations for human trabecular bone - comparison of calcaneus, femoral neck, greater trochanter, proximal tibia, and vertebra

MicroCT-based finite element models were used to compute power law relations for uniaxial compressive yield stress versus bone volume fraction for 78 cores of human trabecular bone from five anatomic sites. The leading coefficient of the power law for calcaneus differed from those for most of the other sites (p < 0.05). However, after normalizing by site-specific mean values, neither the leading coefficient (p > 0.5) nor exponent (p > 0.5) differed among sites, suggesting that a given percentage deviation from mean bone volume fraction has the same mechanical consequence for all sites investigated. These findings help explain the success of calcaneal x-ray and ultrasound measurements for predicting hip fracture risk.

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 Backscatter Measurements of Human Cortical and Trabecular Bone Densities in a Head-Down Bed-Rest Study

This study aims to investigate the feasibility of quantitative ultrasonic backscatter in evaluating human cortical and trabecular bone densities in vivo based on a head-down-tilt bed rest study, with 36 participants tested through 90 d of bed rest and 180 d of recovery. Backscatter measurements were performed using an ultrasonic backscatter bone diagnostic instrument. Backscatter parameters were calculated with a dynamic signal-of-interest method, which was proposed to ensure the same ultrasonic interrogated volume in cortical and trabecular bones. The backscatter parameters exhibited significant correlations with site-matched bone densities provided by high-resolution peripheral quantitative computed tomography (0.33 < |R| < 0.72, p < 0.05). Some bone densities and backscatter parameters exhibited significant changes after the 90-d bed rest. The proposed method can be used to characterize bone densities, and the portable ultrasonic backscatter bone diagnostic device might be used to non-invasively reveal mean bone loss (across a group of people) after long-term bed rest and microgravity conditions of spaceflight missions.

Bayesian inference of human bone sample properties using ultrasonic reflected signals

The non-intrusiveness and low cost of ultrasonic interrogation is motivating the development of new means of detection of osteoporosis and other bone deficiencies. Bone is a porous media saturated with a viscous fluid and could thus be well characterized by the Biot model. The main purpose of this work is to present an in vitro methodology for the identification of the properties and structural parameters of the bone, adopting a statistical Bayesian inference technique using ultrasonic reflected signals at normal incidence. It is, in this respect, a companion paper to a previous work [J. Acoust. Soc. Am. 146, 3 (2019), pp. 1629-1640], where ultrasonic transmitted signals were considered. This approach allows the retrieval of some important parameters that characterize the bone structure and associated uncertainties. The method was applied to seven samples of bone extracted from femoral heads, immersed in water, and exposed to ultrasonic signals with a center frequency of ≈500 kHz. For all seven samples, signals at different sites were acquired to check the method robustness. The porosity, pore mean size and standard deviation, and the porous frame bulk modulus were all successfully identified using only ultrasonic reflected signals.

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.

Ultrasonic Backscatter Difference Measurement of Bone Health in Preterm and Term Newborns

Metabolic bone disease of prematurity remains a significant problem for preterm infants. Quantitative ultrasound (QUS) has potential as a non-invasive tool for assessing bone health of newborns. The aim of this study was to assess bone health in preterm and term newborns using ultrasonic backscatter difference measurement. This study analyzed a total of 493 neonates, including 239 full-term infants (gestational age [GA] >37 wk), 201 preterm I infants (GA: 32-37 wk) and 53 extreme preterm II infants (GA <32 wk). Ultrasonic backscatter measurements were performed on the calcaneus of infants at birth, and the normalized mean of the backscatter difference spectrum (nMBD) was calculated as an ultrasonic index of neonatal bone status. Simple and multiple linear regressions were performed to determine the association of ultrasonic nMBD with GA, anthropometric characteristics and biochemical markers. Statistically significant differences in GA, anthropometric characteristics (birth weight, birth length [BL], birth head circumference and body mass index [BMI]) and biochemical markers (alkaline phosphatase, serum calcium and serum phosphate) were observed among preterm and term infants. The nMBD for term infants (median = 3.72 dB/μs, interquartile range [IR] = 1.95 dB/μs) was significantly higher than that for preterm I infants (median = 1.95 dB/μs, IR = 3.12 dB/μs), which was, in turn, significantly higher than that for preterm II infants (median = 0.19 dB/μs, IR = 3.50 dB/μs). The nMBD yielded moderate correlations (ρ = 0.57-0.62, p < 0.001) with GA and anthropometric characteristics and weak correlations (|ρ| = 0.08-0.21, p < 0.001 or not significant) with biochemical markers. Multivariate regressions revealed that only BL (p = 0.002) and BMI (p = 0.032) yielded significantly independent contributions to the nMBD measurement, and combinations of BL and BMI could explain up to 42% of the variation of nMBD in newborn infants. This study found that ultrasonic backscatter difference measurement might be helpful in bone health evaluation in preterm and term newborns. The utility of ultrasonic backscatter measurement in diagnosis of metabolic bone disease in infants should be investigated further.

Artificial neural network to estimate micro-architectural properties of cortical bone using ultrasonic attenuation: A 2-D numerical study

The goal of this study is to estimate micro-architectural parameters of cortical porosity such as pore diameter (φ), pore density (ρ) and porosity (ν) of cortical bone from ultrasound frequency dependent attenuation using an artificial neural network (ANN). First, heterogeneous structures with controlled pore diameters and pore densities (mono-disperse) were generated, to mimic simplified structure of cortical bone. Then, more realistic structures were obtained from high resolution CT scans of human cortical bone. 2-D finite-difference time-domain simulations were conducted to calculate the frequency-dependent attenuation in the 1-8 MHz range. An ANN was then trained with the ultrasonic attenuation at different frequencies as the input feature vectors while the output was set as the micro-architectural parameters (pore diameter, pore density and porosity). The ANN is composed of three fully connected dense layers with 24, 12 and 6 neurons, connected to the output layer. The dataset was trained over 6000 epochs with a batch size of 16. The trained ANN exhibits the ability to predict the micro-architectural parameters with high accuracy and low losses. ANN approaches could potentially be used as a tool to help inform physics-based modelling of ultrasound propagation in complex media such as cortical bone. This will lead to the solution of inverse-problems to retrieve bone micro-architectural parameters from ultrasound measurements for the non-invasive diagnosis and monitoring osteoporosis.

Monitoring the Setting of Calcium Sulfate Bone-Graft Substitute Using Ultrasonic Backscattering

We report a method to monitor the setting process of bone-graft substitutes (calcium sulfate) using ultrasonic backscattering techniques. Analyzing the backscattered fields using a pulse-echo technique, we show that it is possible to dynamically describe the acoustic properties of the material which are linked to its setting state. Several experiments were performed to control the setting process of calcium sulfate using a 3.5-MHz transducer. The variation of the apparent integrated backscatter (AIB) with time during the setting process is analyzed and compared with measurements of the speed of sound (SOS) and temperature of the sample. The correlation of SOS and AIB allows us to clearly identify two different states of the samples, liquid and solid, in addition to the transition period. Results show that using backscattering analysis, the setting state of the material can be estimated with a threshold of 15 dB. This ultrasonic technique is indeed the first step to develop real-time monitoring systems for time-varying complex media as those present in bone regeneration for dental implantology applications.

Improving broadband ultrasound attenuation assessment in cancellous bone by mitigating the influence of cortical bone: Phantom and in-vitro study

The purpose of this work was to present a new approach that allows the influence of cortical bone on noninvasive measurement of broadband ultrasound attenuation (BUA) to be corrected. The method, implemented here at 1 MHz makes use of backscattered signal and once refined and clinically confirmed, it would offer an alternative to ionizing radiation based methods, such as DEXA (Dual-energy X-ray absorptiometry), quantitative computed tomography (QCT), radiographic absorptiometry (RA) or single X-ray absorptiometry (SXA), which are clinically approved for assessment of progress of osteoporosis. In addition, as the method employs reflected waves, it might substantially enhance the applicability of BUA - from being suitable to peripheral bones only it would extend this applicability to include such embedded bones as hip and femoral neck. The proposed approach allows the cortical layer parameters used for correction and the corrected value and parameter of the cancellous bone (BUA) to be determined simultaneously from the single (pulse-echo) bone backscattered wave; to the best of the authors' knowledge such approach was not previously reported. The validity of the method was tested using acoustic data obtained from a custom-designed bone-mimicking phantom and a calf femur. The relative error of the attenuation coefficient assessment was determined to be 3.9% and 4.7% for the bone phantom and calf bone specimens, respectively. When the cortical shell influence was not taken into account the corresponding errors were considerably higher 8.3% (artificial bone) and 9.2% (calf femur). As indicated above, once clinically proven, the use of this BUA measurement technique in reflection mode would augment diagnostic power of the attending physician by permitting to include bones, which are not accessible for transmission mode evaluation, e.g. hip, spine, humerus and femoral neck.

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.

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%).

Measurement of the Human Calcaneus In Vivo Using Ultrasonic Backscatter Spectral Centroid Shift

OBJECTIVES

The aim of this study was to determine the relationship between the backscattered spectral centroid shift and the bone mineral density (BMD) in vivo and investigate the feasibility of using the backscattered spectral centroid shift to characterize the cancellous bone status.

METHODS

Ultrasonic backscatter measurements were performed in vivo on 1216 participants at the right calcaneus using an ultrasonic backscattered bone diagnostic system, and the backscattered spectral centroid shift was calculated at central frequencies of 3.5 and 5.0 MHz. The BMD values were measured at the sites of the lumbar spine and left hip by dual energy x-ray absorptiometry.

RESULTS

The study population included 592 male and 624 female participants aged 20 to 89 years. The correlations between the backscattered spectral centroid shift in the calcaneus and the spine and hip BMD were found to be statistically significant in both the male and female groups (P < .0001). Linear regression showed that the spectral centroid shift at 3.5 MHz had negative correlations with the spine BMD (R = -0.65 for male participants; R = -0.67 for female participants) and hip BMD (R = -0.64 for male participants; R = -0.64 for female participants). The spectral centroid shift at 5.0 MHz was also found to be closely related to the spine BMD (R = -0.68 for male participants; R = -0.68 for female participants) and hip BMD (R = -0.66 for male participants; R = -0.64 for female participants).

CONCLUSIONS

The moderate correlations observed between the spectral centroid shift and the spine and hip BMD demonstrate that the ultrasonic backscattered spectral centroid shift may be a useful measurement for assessment of the cancellous bone status.

Backscatter-difference Measurements of Cancellous Bone Using an Ultrasonic Imaging System

Backscatter-difference measurements may be used to detect changes in bone caused by osteoporosis. The backscatter-difference technique measures the power difference between two portions of an ultrasonic backscatter signal. The goal of this study is to evaluate the feasibility of using an ultrasonic imaging system to perform backscatter-difference measurements of bone. Ultrasonic images and backscatter signals were acquired from 24 specimens of human cancellous bone. The signals were analyzed in the frequency domain to determine the normalized mean backscatter-difference (nMBD) and in the time domain to determine the normalized backscatter amplitude ratio (nBAR). The images were analyzed to determine the normalized pixel value difference (nPVD), which measures the difference in average pixel brightness between regions of interest placed at two different depths in the image. All three parameters were found to increase with bone mineral density. The signal-based parameters, nMBD and nBAR, correlated well with bone mineral density, yielding linear correlation coefficients that ranged from 0.74 to 0.87. The image based parameter, nPVD, performed somewhat less well, yielding correlation coefficients that ranged from 0.42 to 0.81. These results suggest that ultrasonic imaging systems may be used to perform backscatter-difference measurements for the purpose of ultrasonic bone assessment.

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.

Numerical Analysis of Ultrasound Backscattered Waves in Cancellous Bone Using a Finite-Difference Time-Domain Method: Isolation of the Backscattered Waves From Various Ranges of Bone Depths

Using a finite-difference time-domain method, ultrasound backscattered waves inside cancellous bone were numerically analyzed to investigate the backscatter mechanism. Two bone models with different thicknesses were modeled with artificial absorbing layers positioned at the back surfaces of the model, and an ultrasound pulse wave was transmitted toward the front surface. By calculating the difference between the simulated waveforms obtained using the two bone models, the backscattered waves from a limited range of depths in cancellous bone could be isolated. The results showed that the fast and slow longitudinal waves, which have previously been observed only in the ultrasound waveform transmitted through the bone, could be distinguished in the backscattered waveform from a deeper bone depth when transmitting the ultrasound wave parallel to the main orientation of the trabecular network. The amplitudes of the fast and slow backscattered waves were more closely correlated with the bone porosity [R2 = 0.84 and 0.66 (p < 0.001), respectively] than the amplitude of the whole (nonisolated) backscattered waves [R2 = 0.48 (p < 0.001)]. In conclusion, the nonisolated backscattered waves could be regarded as the superposition of the fast and slow waves reflected from various bone depths, returning at different times.

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.

Feasibility of bone assessment with ultrasonic backscatter signals in neonates

The objective of this study was to assess the value of ultrasonic backscatter signals and the backscatter coefficient (BSC) in the analysis of bone status in neonates and to analyze the relationships between the BSC and gestational age, birth weight, length, head circumference and gender. A total of 122 neonates participated in the study, including 83 premature infants and 39 full-term infants. Their BSCs were measured by ultrasound after birth. The results revealed a significant correlation between the BSC and gestational age (R = 0.47, p < 0.001), birth weight (R = 0.47, p < 0.0001) and length at birth (R = 0.43, p < 0.001) at a frequency of 5.0 MHz. This study suggests that the use of ultrasonic backscattering and the BSC is feasible for assessment of the bone status of neonates.

Ultrasonic scanner for in vivo measurement of cancellous bone properties from backscattered data

A dedicated ultrasonic scanner for acquiring RF echoes backscattered from the trabecular bone was developed. The design of device is based on the goal of minimizing of custom electronics and computations executed solely on the main computer processor and the graphics card. The electronic encoder-digitizer module executing all of the transmission and reception functions is based on a single low-cost field programmable gate array (FPGA). The scanner is equipped with a mechanical sector-scan probe with a concave transducer with 50 mm focal length, center frequency of 1.5 MHz and 60% bandwidth at -6 dB. The example of femoral neck bone examination shows that the scanner can provide ultrasonic data from deeply located bones with the ultrasound penetrating the trabecular bone up to a depth of 20 mm. It is also shown that the RF echo data acquired with the scanner allow for the estimation of attenuation coefficient and frequency dependence of backscattering coefficient of trabecular bone. The values of the calculated parameters are in the range of corresponding in vitro data from the literature but their variation is relatively high.

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.

Modeling and analysis of multiple scattering of acoustic waves in complex media: application to the trabecular bone

The integral equations that describe scattering in the media with step-rise changing parameters have been numerically solved for the trabecular bone model. The model consists of several hundred discrete randomly distributed elements. The spectral distribution of scattering coefficients in subsequent orders of scattering has been presented. Calculations were carried on for the ultrasonic frequency ranging from 0.5 to 3 MHz. Evaluation of the contribution of the first, second, and higher scattering orders to total scattering of the ultrasounds in trabecular bone was done. Contrary to the approaches that use the μCT images of trabecular structure to modeling of the ultrasonic wave propagation condition, the 3D numerical model consisting of cylindrical elements mimicking the spatial matrix of trabeculae, was applied. The scattering, due to interconnections between thick trabeculae, usually neglected in trabecular bone models, has been included in calculations when the structure backscatter was evaluated. Influence of the absorption in subsequent orders of scattering is also addressed. Results show that up to 1.5 MHz, the influence of higher scattering orders on the total scattered field characteristic can be neglected while for the higher frequencies, the relatively high amplitude interference peaks in higher scattering orders clearly occur.

Statistics of the envelope of ultrasonic backscatter from human trabecular bone

The paper describes the investigations intended to compare the results of experimental measurements of backscattering properties of the trabecular bone with the results of computer simulations. Ultrasonic RF echoes were collected using two bone scanners operating at 0.58 and 1.3 MHz. The simulations of the backscattered RF echoes were performed using the scattering model of the trabecular bone that consisted of cylindrical and spherical elements uniformly distributed in water-like medium. For each measured or simulated RF backscatter the statistical properties of the signal envelope were determined. Experimental results suggest deviations of the backscattering properties from the Rayleigh distribution. The results of simulation suggest that deviation from Rayleigh distribution depends on the variation of trabeculae diameters and the number of thin trabeculae. Experimentally determined deviations corresponded well to the deviations calculated from simulated echoes assuming trabeculae thickness variation equaled to the earlier published histomorphometric study results.

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.

Semi-empirical bone model for determination of trabecular structure properties from backscattered ultrasound

A novel semi-empirical scattering model of trabecular bone facilitating its characterization and allowing optimization of the interrogating pulse-echo transducer performance was developed. The model accounts for spatial density distribution of the trabeculae and includes measurement conditions such as pressure-time waveform of the probing ultrasound wave, the emitted field structure, and the transfer function and limited bandwidth of the acoustic source operating in pulse-echo mode. These measurement conditions are of importance as they modify the scattered echoes, which in turn are linked to the micro-architecture of the bone. The bone was modeled by a random distribution of long and thin cylindrical scatterers having randomly varying diameters and mechanical properties, and oriented perpendicularly to the ultrasound beam axis. To mimic clinically encountered conditions the relevant empirical data obtained at 1 MHz were input to the model. The data included pulse-echo source pressure field distribution in the focal zone and the above mentioned transfer function. With these data the model allowed frequency dependent backscattering coefficient of the simulated bone structure and its statistical properties to be determined. The results obtained indicated that the computer simulation is of particular relevance in studying scattering properties of the cancellous bone and holds promise as a tool to determine the relationship between the physical dimensions and shape of the scatterers and for monitoring of osteoporosis. The results of simulations also indicated that the new bone model proposed is well suited to mimic clinically relevant conditions. In contrast to the existing bone models, which usually assume scatterers to be randomly distributed as infinitely long identical cylinders with a cross-section much smaller than the probing ultrasound wave, the new model includes two populations of scatterers having different physical dimensions and also allows the mechanical properties of the scatterers to be varied.

Analysis of frequency dependence of ultrasonic backscatter coefficient in cancellous bone

The ultrasonic scattering mechanism in cancellous bone is investigated theoretically and a model describing the frequency dependence of ultrasonic scattering from cancellous bone is presented. The ultrasonic backscatter coefficient (BSC) of bovine tibiae, human calcanei in vitro and in vivo, were measured and discussed. The data of BSC were also fitted by polynomial. The results demonstrate that BSC is a nonlinear function of frequency and increases with frequency. A good agreement was obtained between BSC values from theory and experiment. Also, the high correlation coefficient between BSC and bone mineral density was obtained, r=0.85+/-0.07 (mean+/-SD) (n=15, p<0.001). Based on the values of BSC, the status of cancellous bone and the degree of osteoporotic fracture risk may be assessed.

Ultrasonic attenuation in parallel-nylon-wire cancellous-bone-mimicking phantoms

Attenuation coefficients between 1.5 and 3.5 MHz were measured on four parallel-nylon-wire arrays (simulating cancellous bone) with four different wire diameters (150, 200, 250, and 300 microm). Interwire spacing was 800 microm for all four parallel-nylon-wire arrays. The measured frequency dependencies of attenuation were consistent with theoretical predications based on Faran's theory, which considers the component of attenuation due to scattering of longitudinal waves.

Mechanisms for attenuation in cancellous-bone-mimicking phantoms

Broadband ultrasound attenuation (BUA) in cancellous bone is useful for prediction of osteoporotic fracture risk, but its causes are not well understood. To investigate attenuation mechanisms, 9 cancellous-bone-mimicking phantoms containing nylon filaments (simulating bone trabeculae) embedded within soft-tissue-mimicking fluid (simulating marrow) were interrogated. The measurements of frequency-dependent attenuation coefficient had 3 separable components: 1) a linear (with frequency) component attributable to absorption in the soft-tissue-mimicking fluid, 2) a quasilinear (with frequency) component, which may include absorption in and longitudinal-shear mode conversion by the nylon filaments, and 3) a nonlinear (with frequency) component, which may be attributable to longitudinal-longitudinal scattering by the nylon filaments. The slope of total linear (with frequency) attenuation coefficient (sum of components #1 and #2) versus frequency was found to increase linearly with volume fraction, consistent with reported measurements on cancellous bone. Backscatter coefficient measurements in the 9 phantoms supported the claim that the nonlinear (with frequency) component of attenuation coefficient (component #3) was closely associated with longitudinal-longitudinal scattering. This work represents the first experimental separation of these 3 components of attenuation in cancellous bone-mimicking phantoms.

Ultrasonic scattering from cancellous bone: a review

This paper reviews theory, measurements, and computer simulations of scattering from cancellous bone reported by many laboratories. Three theoretical models (binary mixture, Faran cylinder, and weak scattering) for scattering from cancellous bone have demonstrated some consistency with measurements of backscatter. Backscatter is moderately correlated with bone mineral density in human calcaneus in vitro (r(2) = 0.66 - 0.68). Backscatter varies approximately as frequency cubed and trabecular thickness cubed in human calcaneus and femur in vitro. Backscatter from human calcaneus and bovine tibia exhibits substantial anisotropy. So far, backscatter has demonstrated only modest clinical utility. Computer simulation models have helped to elucidate mechanisms underlying scattering from cancellous bones.

Ultrasonic characterization of human cancellous bone in vitro using three different apparent backscatter parameters in the frequency range 0.6-15.0 mhz

Ultrasonic techniques based on measurements of apparent backscatter may provide a useful means for diagnosing bone diseases such as osteoporosis. The term "apparent" means that the backscattered signals are not compensated for the frequency-dependent effects of attenuation and diffraction. We performed in vitro apparent backscatter measurements on 23 specimens of human cancellous bone prepared from the left and right femoral heads of seven donors. A mechanical scanning system was used to obtain backscattered signals from each specimen at several sites. Scans were performed using five different ultrasonic transducers with center frequencies of 1, 2.25, 5, 7.5, and 10 MHz. The -6 dB bandwidths of these transducers covered a frequency range of 0.6-15.0 MHz. The backscattered signals were analyzed to determine three ultrasonic parameters: apparent integrated backscatter (AIB), frequency slope of apparent backscatter (FSAB), and time slope of apparent backscatter (TSAB). Linear regression analysis was used to examine the correlation of these ultrasonic parameters with five measured physical characteristics of the specimens: mass density, X-ray bone mineral density, Young's modulus, yield strength, and ultimate strength. A total of 75 such correlations were examined (3 ultrasonic parameters x 5 specimen characteristics x 5 transducers). Good correlations were observed for AIB using the 5 MHz (r = 0.70 - 0.89) and 7.5 MHz (r = 0.75-0.93) transducers; for FSAB using the 2.25 MHz (r = 0.70 - 0.88), 5 MHz (r = 0.79 - 0.94), and 7.5 MHz (r = 0.80 - 0.92) transducers; and for TSAB using the 5 MHz (r = 0.68 - 0.89), 7.5 MHz (r = 0.75 - 0.89), and 10 MHz (r = 0.75 - 0.92) transducers.

Instrumentation for in vivo ultrasonic characterization of bone strength

Although it has been more than 20 years since the first recorded use of a quantitative ultrasound (QUS) technology to predict bone fragility, the field has not yet reached its maturity. QUS has the potential to predict fracture risk in several clinical circumstances and has the advantages of being nonionizing, inexpensive, portable, highly acceptable to patients, and repeatable. However, the wide dissemination of QUS in clinical practice is still limited and suffering from the absence of clinical consensus on how to integrate QUS technologies in bone densitometry armamentarium. Several critical issues need to be addressed to develop the role of QUS within rheumatology. These include issues of technologies adapted to measure the central skeleton, data acquisition, and signal processing procedures to reveal bone properties beyond bone mineral quantity and elucidation of the complex interaction between ultrasound and bone structure. This article reviews the state-of-the art in technological developments applied to assess bone strength in vivo. We describe generic measurement and signal processing methods implemented in clinical ultrasound devices, the devices and their practical use, and performance measures. The article also points out the present limitations, especially those related to the absence of standardization, and the lack of comprehensive theoretical models. We conclude with suggestions of future lines and trends in technology challenges and research areas such as new acquisition modes, advanced signal processing techniques, and modelization.

Spatial variation of acoustic properties is related with mechanical properties of trabecular bone

In clinical applications, ultrasound parameters are measured as an average value over a region of interest (ROI) or as a value at a single measurement point. Due to natural adaptation to loading conditions, trabecular bone is structurally, compositionally and mechanically heterogeneous and anisotropic. Thus, spatial variation of ultrasound parameters within ROI may contain valuable information on the mechanical integrity of trabecular bone. However, this issue has not been thoroughly investigated. In the present study, we aimed at investigating the significance of the spatial variation of ultrasound parameters for the prediction of mechanical properties of human trabecular bone. For this aim, parametric maps of apparent integrated backscattering (AIB), integrated reflection coefficient (IRC), speed of sound (SOS), average attenuation (AA) and normalized broadband ultrasound attenuation (nBUA) were calculated for femoral and tibial bone cylinders (n = 19-20). Further, the effect of time window length on the AIB, variation of AIB within ROI and association between AIB and bone mechanical properties were characterized. Based on linear correlation analysis, spatial variation of AIB, assessed as standard deviation of measurements within ROI, was a strong predictor of bone ultimate strength (r = -0.82, n = 19, p < 0.01). Further, the time window length affected absolute values of AIB and strength of correlation between AIB and bone ultimate strength. Interestingly, linear combination of mean IRC and spatial variation of AIB within ROI was the strongest predictor of bone ultimate strength (r = 0.92, n = 19, p < 0.01). In conclusion, our findings suggest that the measurement of two-dimensional parametric maps of ultrasound parameters could yield information on bone status not extractable from single point measurements. This highlights the potential of parametric imaging in osteoporosis diagnostics.

Acoustic properties of trabecular bone--relationships to tissue composition

In osteoporosis, changes in tissue composition and structure reduce bone strength and expose it to fractures. The current primary diagnostic technique, i.e., dual energy X-ray absorptiometry, measures areal bone mineral density (BMD) but provides no direct information on trabecular structure or organic composition. Although still poorly characterized, ultrasound techniques may bring about information on bone composition and structure. In this study, relationships of 2.25-MHz ultrasound speed, attenuation, reflection and backscattering with composition of human trabecular bone (n=26) were characterized experimentally, as well as by using numerical analyses. We also determined composition of the trabecular sample (fat and water content, bone volume fraction) and that of the calcified matrix (mineral, proteoglycan and collagen content of trabeculae). In experimental analyses, bone volume fraction and mineral content of the calcified matrix were the only determinants of BMD. Further, bone volume fraction served as the strongest determinant of ultrasound parameters (r=0.51-0.87). In numerical simulations, density and mechanical properties of the calcified matrix systematically affected ultrasound speed, attenuation, reflection and backscattering. However, partial correlation coefficients revealed only low associations(|r|<or=0.4) between the composition of calcified matrix and ultrasound parameters in experimental measurements. To conclude, the content and structure of calcified matrix, rather than its composition, affect more significantly acoustic properties of healthy trabecular bone.

Characterization of dense bovine cancellous bone tissue microstructure by ultrasonic backscattering using weak scattering models

A weak scattering model was proposed for the ultrasonic frequency-dependent backscatter in dense bovine cancellous bone, using two autocorrelation functions to describe the medium: one with discrete homogeneities (spherical distribution of equal spheres) and another, which considers tissue as an inhomogeneous continuum (densely populated medium). The inverse problem to estimate trabecular thickness of bone tissue has been addressed. A combination of the two autocorrelation functions was required to closely approximate the backscatter from bovine bone with various microarchitecture, given that the shape of trabeculae ranges from a rodlike to a platelike shape. Because of the large variation in trabecular thickness, both at an intraspecimen and an interspecimen level, thickness distributions for individual trabeculae for each bone specimen were obtained, and dominant trabecular sizes were determined. Comparison of backscatter measurements to theoretical predictions indicated that there were more than one dominant trabecular sizes that scatter sound for most specimens. Linear regression, performed between dominant trabecular thickness and estimated correlation length, showed significant linear correlation (R(2)=0.81). Attenuation due to scattering by a continuous distribution of scatterers was predicted to be linear over a frequency range from 0.3 to 0.9 MHz, suggesting a possibility that scattering may be a significant source of attenuation.

Wavelet-based signal processing of in vitro ultrasonic measurements at the proximal femur

To estimate osteoporotic fracture risk, several techniques for quantitative ultrasound (QUS) measurements at peripheral sites have been developed. As these techniques are limited in the prediction of fracture risk of the central skeleton, such as the hip, we are developing a QUS device for direct measurements at the femur. In doing so, we noticed the necessity to improve the conventional signal processing because it failed in a considerable number of measurements due to multipath transmission. Two sets of excised human femurs (n = 6 + 34) were scanned in transmission mode. Instead of using the conventional methods, the radio-frequency signals were processed with the continuous wavelet transform to detect their time-of-flights for the calculation of speed-of-sound (SOS) in bone. The SOS-values were averaged over a region similar to the total hip region of dual X-ray absorptiometry (DXA) measurements and compared with bone mineral density (BMD) measured with DXA. Testing six standard wavelets, this algorithm failed for only 0% to 6% of scan in test set 1 compared with 29% when using conventional algorithms. For test set 2, it failed for 2% to 12% compared with approximately 40%. SOS and BMD correlated significantly in both test sets (test set 1: r2 = 0.87 to 0.92, p < 0.007; test set 2: r2 = 0.68 to 0.79, p < 0.0001). The correlations are comparable with correlations recently reported. However, the number of evaluable signals could be substantially increased, which improves the perspectives of the in vivo measurements.

Predict ultrasonic backscatter coefficient in cancellous bone by theory and experiment

The scattering mechanism of ultrasound in cancellous bone is investigated theoretically. The relationship of backscatter coefficient (BSC) in cancellous bone with frequency is analyzed in theoretical and experimental. The results of theory and experiment for cancellous bone of bovine tibiae, human calcaneus in vitro and in vivo showed that BSC is a non-linear function of frequency, increasing with frequency. In general, all curve of BSC can be divided into three sections with different slope. The slopes of the first and third section have a large value, and the slope of second section is flat. A good agreement was obtained in the averaged BSC of experiment and cellular model. Those results suggest that the backscatter signal and BSC have a particularly important action in assessment of cancellous bone status and diagnosis of osteoporosis.

Influence of the precision of spectral backscatter measurements on the estimation of scatterers size in cancellous bone

The goal of this study is to propose a model for the ultrasonic frequency-dependent backscatter coefficient in femoral cancellous bone. This model has been developed with success to predict backscatter in human calcaneal bone [Jenson, Ultr. Med. Biol. 2003]. A weak scattering model is used and the backscatter coefficient is expressed in terms of a Gaussian autocorrelation function of the medium. The backscatter coefficient is computed and comparison is made with experimental data for 37 specimens and for frequency ranging from 0.4 to 1.2 MHz. An excellent agreement between experimental data and predictions is found for both the magnitude and the frequency-dependence of the backscatter coefficient. Then, a nonlinear regression is performed for each specimen, and the mean trabecular thickness is estimated. Experimental data and theoretical predictions are averaged over the 37 specimens. We also find a close agreement between theoretical predictions obtained using the Gaussian autocorrelation function (scatterer size=134+/-15 microm) and the mean trabecular thickness (Tb.Th=132+/-12 microm) derived from the analysis of bone 3-D micro-architecture using high-resolution micro-tomography. However, the correlation between individual experimental and estimated Tb.Th values is moderate (R(2)=0.44). The performance of the estimator are limited mainly by two factors: interference noise due to random positioning of the scatterers and attenuation. We show that the fundamental limitation of our estimator due to the speckle noise is around 5 microm for trabecular thickness estimation. This limitation is lower than the observed biological variability which is around 30 microm and should not be a limiting factor for individual prediction. A second limitation is the tremendous attenuation encountered in highly scattering media such as cancellous bone, which results in highly damped backscatter signals. The compensation for attenuation is difficult to perform, and it may be a critical point that limits the precision of the estimator.

Ultrasonic characterization of cancellous bone using apparent integrated backscatter

Apparent integrated backscatter (AIB) is a measure of the frequency-averaged (integrated) backscattered power contained in some portion of a backscattered ultrasonic signal. AIB has been used extensively to study soft tissues, but its usefulness as a tissue characterization technique for cancellous bone has not been demonstrated. To address this, we performed measurements on 17 specimens of cancellous bone over two different frequency ranges using a 1 MHz and 5 MHz broadband ultrasonic transducer. Specimens were obtained from bovine tibiae and prepared in the shape of cubes (15 mm side length) with faces oriented along transverse (anterior, posterior, medial and lateral) and longitudinal (superior and inferior) principal anatomic directions. A mechanical scanning system was used to acquire multiple backscatter signals from each direction for each cube. AIB demonstrated highly significant linear correlations with bone mineral density (BMD) for both the transverse (R2 = 0.817) and longitudinal (R2 = 0.488) directions using the 5 MHz transducer. In contrast, the correlations with density were much weaker for the 1 MHz transducer (R2 = 0.007 transverse, R2 = 0.228 longitudinal). In all cases where a significant correlation was observed, AIB was found to decrease with increasing BMD.

Ultrasonic characterization of human trabecular bone microstructure

New quantitative ultrasound (QUS) techniques involving ultrasound backscattering have been introduced for the assessment of bone quality. QUS parameters are affected by the transducer characteristics, e.g. frequency range, wave and pulse length. Although frequency-dependent backscattering has been studied extensively, understanding of the ultrasound scattering phenomenon in trabecular bone is still limited. In the present study, the relationships between QUS parameters and the microstructure of human trabecular bone were investigated experimentally and by using numerical simulations. Speed of sound (SOS), normalized broadband ultrasound attenuation (nBUA), average attenuation, integrated reflection coefficient (IRC) and broadband ultrasound backscatter (BUB) were measured for 26 human trabecular bone cylinders. Subsequently, a high-resolution microCT system was used to determine the microstructural parameters. Moreover, based on the sample-specific microCT data, a numerical model for ultrasound propagation was developed for the simulation of experimental measurements. Experimentally, significant relationships between the QUS parameters and microstructural parameters were demonstrated. The relationships were dependent on the frequency, and the strongest association (r = 0.88) between SOS and structural parameters was observed at a centre frequency of 5 MHz. nBUA, average attenuation, IRC and BUB showed somewhat lower linear correlations with the structural properties at a centre frequency of 5 MHz, as compared to those determined at lower frequencies. Multiple regression analyses revealed that the variation of acoustic parameters could best be explained by parameters reflecting the amount of mineralized tissue. A principal component analysis demonstrated that the strongest determinants of BUB and IRC were related to the trabecular structure. However, other structural characteristics contributed significantly to the prediction of the acoustic parameters as well. The two-dimensional numerical model introduced in the present study demonstrated good agreement with the experimental measurements. However, further studies with the simulation model are warranted to systematically investigate the relation between the structural parameters and ultrasound scattering.

Estimation of trabecular thickness using ultrasonic backcatter

We present a method to estimate trabecular thickness (Tb.Th) in trabecular bones from ultrasound backscatter measurements. The estimation scheme is based on a nonlinear adjustment of predictions from a model to experimental data. The model assumes weak scattering from bone, where scattering is assumed to arise from the elastic solid trabeculae. The fluctuations of acoustical properties between bone tissue and the saturating fluid are assumed to be random and are described by the 3-D spatial autocorrelation function of the medium. In this paper, a Gaussian autocorrelation function is used. The inversion procedure is applied to a set of data measured on 33 femoral bone specimens. Results show that the model can predict both the magnitude and the frequency-dependence of the backscatter coefficient (root mean square error RMSE = 1 dB). The estimated trabecular thickness values are compared to the true trabecular thickness measured on high resolution microcomputed tomography 3-D reconstruction of bones microarchitecture. A close agreement is obtained on average over the group of specimens between predictions and the reference values: true Tb.Th is 132 +/- 12 microm and estimated Tb.Th is 134 +/- 15 microm. However, a moderate correlation between actual and estimated Tb.Th values is found (R2 = 0.44, p<10(-4), RMSE = 8.7 microm) suggesting a modest predictability at the individual level. Sources for the variability of the estimator are studied. Using synthetic rf signals, we demonstrate that the fundamental limitation of the estimator due to speckle noise is approximately 5 microm. Taking into account the measurement errors, the total uncertainty on Tb.Th estimates is of the order of 7 microm. The influence of the attenuation compensation function used to derive the backscatter coefficient is studied. In particular, we demonstrate the necessity of compensating for the effect of the gating time window. The results are discussed with respect to their meaningful clinical value. The requirements to be fulfilled by the performance of the technique change with regard to the question being posed. Two different strategies are examined: 1. characterize trabecular thickness without consideration of bone quantity (or bone mineral density) and 2. estimate trabecular thickness after adjustment for BMD. Considering the first strategy, a comparison between the precision of our estimator and the biological variability leads us to the conclusion that our estimator should only permit to distinguish between micro-architectures characterized by extreme values of trabecular thickness (i.e., very thin or very thick trabecular thickness). In this respect, it would be interesting to test whether the estimator is able to discriminate between rod-like (thin) and plate-like (thick) structures that are known to influence differently bone strength. The second strategy is more demanding in terms of technique performance and our estimator is not able yet to catch small differences in Tb.Th values expected after adjustment to bone density. Progress in the field will require a significant reduction in speckle noise and measurement errors and/or the development of other and more efficient microstructural estimators.

Kramers-Kronig analysis of attenuation and dispersion in trabecular bone

A restricted-bandwidth form of the Kramers-Kronig dispersion relations is applied to in vitro measurements of ultrasonic attenuation and dispersion properties of trabecular bone specimens from bovine tibia. The Kramers-Kronig analysis utilizes only experimentally measured properties and avoids extrapolation of ultrasonic properties beyond the known bandwidth. Compensation for the portions of the Kramers-Kronig integrals over the unknown bandwidth is partially achieved by the method of subtractions, where a subtraction frequency acts as an adjustable parameter. Good agreement is found between experimentally measured and Kramers-Kronig reconstructed dispersions. The restricted-bandwidth approach improves upon other forms of the Kramers-Kronig relations and may provide further insight into how ultrasound interacts with trabecular bone.

Optimal prediction of bone mineral density with ultrasonic measurements in excised human femur

Bone mineral density (BMD) measured with dual energy X-ray absorptiometry (DXA) techniques is the current gold standard for osteoporotic fracture risk prediction. Quantitative ultrasound (QUS) techniques in transmission measurements are, however, increasingly recognized as an alternative approach. It is feasible to select different QUS methods, one type being optimized to assess microarchitectural properties of bone structure and another to assess BMD. Broadband ultrasonic attenuation (BUA) and ultrasonic velocity (UV) measured on the proximal human femur have been shown to be both significantly correlated with BMD. However, a great diversity of algorithms has been reported to measure the time-of-flight used to derive UV values. The purpose of this study was to determine which procedure results in the optimal BMD prediction at the proximal femur from ultrasound measurements. Thirty-eight excised human femurs were measured in transmission with a pair of focused 0.5-MHz central frequency transducers. Two-dimensional scans were performed and radiofrequency (RF) signals were recorded digitally at each scan position. BUA was estimated and eight different signal processing techniques were performed to estimate UV. For each signal-processing technique UV was compared to BMD. We show that the best prediction of BMD was obtained with signal-processing techniques taking into account only the first part of the transmitted signal (r2BMD-SOS = 0.86). Moreover, we show that a linear multiple regression using both BUA and speed of sound (SOS) and applied to site-matched regions of interest improved the accuracy of BMD predictions (r2BMD-SOS/BUA = 0.95). Our results demonstrate that selecting specific signal-processing methods for QUS variables allows optimal assessment of BMD. Correlation is sufficiently high that this specific QUS method can be considered as a good surrogate of BMD.

Prediction of density and mechanical properties of human trabecular bone in vitro by using ultrasound transmission and backscattering measurements at 0.2-6.7 MHz frequency range

The ultrasound (US) backscattering method has been introduced as an alternative for the through-transmission measurement of sound attenuation and speed in diagnosis of osteoporosis. Both attenuation and backscattering depend strongly on the US frequency. In this study, 20 human trabecular bone samples were measured in transmission and pulse-echo geometry in vitro. The aim of the study was to find the most sensitive frequency range for the quantitative ultrasound (QUS) analyses. Normalized broadband US attenuation (nBUA), speed of sound (SOS), broadband US backscatter (BUB) and integrated reflection coefficient (IRC) were determined for each sample. The samples were spatially scanned with five pairs of US transducers covering a frequency range of 0.2-6.7 MHz. Furthermore, mechanical properties and density of the same samples were determined. At all frequencies, SOS, BUB and IRC showed statistically significant linear correlations with the mechanical properties or density of human trabecular bone (0.51 < r < 0.82, 0.54 < r < 0.81 and 0.70 < r < 0.85, respectively). In contrast to SOS, IRC and BUB, nBUA showed statistically significant correlations with mechanical parameters or density at the centre frequency of 1 MHz only. Our results suggest that frequencies up to 5 MHz can be useful in QUS analyses for the prediction of bone mechanical properties and density. Since the use of higher frequencies provides better axial and spatial resolution, improved structural analyses may be possible. While extensive attenuation of high frequencies in trabecular bone limits the clinically feasible frequency range, selection of optimal frequency range for in vivo QUS application should be carefully considered.

Ability of ultrasound backscattering to predict mechanical properties of bovine trabecular bone

Ultrasound (US) backscatter measurements have been proposed for the quantitative evaluation of bone quality. In this study, we explored the ability of broadband US backscatter (BUB) and integrated reflection coefficient (IRC) to predict density and mechanical properties of trabecular bone, as compared to normalized broadband US attenuation (nBUA) and speed of sound (SOS). These acoustic parameters were measured in 41 in vitro samples of bovine trabecular bone and correlated with a number of mechanical parameters and with volumetric bone mineral density (BMDvol). BUB correlated statistically significantly with the volumetric bone mineral density (r = 0.61, p < 0.01), Young's modulus (r = 0.40, p < 0.01) and ultimate strength (r = 0.40, p < 0.01). IRC was even more strongly correlated with BMD(vol) (r = 0.92, p < 0.01) and most of the mechanical parameters (0.81 < r < 0.85). Strong correlations were also found between mechanical parameters and SOS (0.87 < r < 0.90). No significant correlation was found between attenuation (nBUA) and either BMD(vol) or mechanical parameters. Reproducibilities (standardized CV%) of BUB (3.5%) and IRC (1.5%) were comparable to those of nBUA (2.3%) and SOS (0.5%). To conclude, BUB and IRC are promising parameters for the evaluation of density and mechanical properties of trabecular bone. Advantageously, BUB and IRC can be determined with a single transducer, hypothetically enabling measurements at many clinically relevant fracture sites.

Singular spectrum analysis applied to backscattered ultrasound signals from in vitro human cancellous bone specimens

Mean scatterer spacing (MSS) holds particular promise for the detection of changes in quasiperiodic tissue microstructures such as may occur during development of disease in the liver, spleen, or bones. Many techniques that may be applied for MSS estimation (temporal and spectral autocorrelation, power spectrum and cepstrum, higher order statistics, and quadratic transformation) characterize signals that contain a mixture of periodic and nonperiodic contributions. In contrast, singular spectrum analysis (SSA), a method usually applied in nonlinear dynamics, first identifies components of signals corresponding to periodic structures and, second, identifies dominant periodicity. Thus, SSA may better separate periodic structures from nonperiodic structures and noise. Using an ultrasound echo simulation model, we previously demonstrated SSA's potential to identify MSS of structures in quasiperiodic scattering media. The current work aims to observe the behavior of MSS estimation by SSA using ultrasound measurements in phantom materials (two parallel, nylon-line phantoms and four foam phantoms of different densities). The SSA was able to estimate not only the nylon-line distances but also nylon-line thickness. The method also was sensitive to the average pore-size differences of the four sponges. The algorithms then were applied to characterize human cancellous bone microarchitectures. Using 1-MHz center-frequency, radio-frequency ultrasound signals, MSS was measured in 24 in vitro bone samples and ranged from 1.0 to 1.7 mm. The SSA MSS estimates correlate significantly to MSS measured independently from synchrotron microtomography, r2 = 0.68. Thus, application of SSA to backscattered ultrasound signals seems to be useful for providing information linked to tissue microarchitecture that is not evident from clinical images.