Frequency dependence of ultrasonic backscatter from human trabecular bone: theory and experiment

Quantitative ultrasound techniques for assessing thermal ablation: Measurement of the backscatter coefficient from ex vivo human liver

BACKGROUND

Understanding the changes occurring in biological tissue during thermal ablation is at the heart of many current challenges in both therapy and medical imaging research.

PURPOSE

The objective of this work is to quantitatively interpret the scattering response of human liver samples, before and after thermal ablation. We report acoustic measurements performed involving n = 21 human liver samples. Thermal ablation is achieved at temperatures between 45 and 80°C and quantification of the irreversible changes in acoustic attenuation and Backscattering Coefficient (BSC) is reported, with a particular attention to the latter.

METHODS

Both attenuation coefficient and BSCs were measured in the frequency range from 10 to 52 MHz. Scans were performed before heating and after cooling down. Attenuation coefficients were calculated using spectral difference method and BSC estimated using the reference phantom method.

RESULTS

Strong increases of attenuation coefficients and BSCs with heating temperature were observed. Quantitative ultrasonic parameters obtained with the polydisperse structure factor model (poly-SFM)are compared to histological observations and seen to be close to hepatocyte mean diameter (HMD).

CONCLUSIONS

The results presented in this study provide a description of the impact of thermal ablation in human liver tissue on acoustic attenuation and the BSC. For the first time, quantitative agreement between the Effective Scatterer Diameter (ESD) estimated from BSC and HMD was shown, highlighting the important role of cellular network in the scattering response of the medium. This core result is an important step toward the determination of the nature of scattering sources in biological tissues.

Extracting Quantitative Ultrasonic Parameters from the Backscatter Coefficient

The ultrasonic backscatter coefficient (BSC) is a fundamental quantitative ultrasound (QUS) parameter that contains rich information about the underlying tissue. Deriving parameters from the BSC is essential for fully utilizing the information contained in BSC for tissue characterization. In this chapter, we review two primary approaches for extracting parameters from the BSC versus frequency curve: the model-based approach and the model-free approach, focusing on the model-based approach, where a scattering model is fit to the observed BSC to yield model parameters. For this approach, we will attempt to unite commonly used models under a coherent theoretical framework. We will focus on the underlying assumptions and conditions for various BSC models. Computer code is provided to facilitate the use of some of the models. The strengths and weaknesses of various models are also discussed.

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.

Estimation of Cortical Bone Microstructure From Ultrasound Backscatter

Multichannel pulse-echo ultrasound using linear arrays and single-channel data acquisition systems opens new perspectives for the evaluation of cortical bone. In combination with spectral backscatter analysis, it can provide quantitative information about cortical microstructural properties. We present a numerical study, based on the finite-difference time-domain method, to estimate the backscatter cross section of randomly distributed circular pores in a bone matrix. A model that predicts the backscatter coefficient using arbitrary pore diameter distributions was derived. In an ex vivo study on 19 human tibia bones (six males, 13 females, 83.7 ± 8.4 years), multidirectional ultrasound backscatter measurements were performed using an ultrasound scanner equipped with a 6-MHz 128-element linear array with sweep motor control. A normalized depth-dependent spectral analysis was performed to derive backscatter and attenuation coefficients. Site-matched reference values of tissue acoustic impedance Z , cortical thickness (Ct.Th), pore density (Ct.Po.Dn), porosity (Ct.Po), and characteristic parameters of the pore diameter (Ct.Po.Dm) distribution were obtained from 100-MHz scanning-acoustic microscopy images. Proximal femur areal bone mineral density (aBMD), stiffness S , and ultimate force Fu from the same donors were available from a previous study. All pore structure and material properties could be predicted using linear combinations of backscatter parameters with a median to high accuracy (0.28 ≤ adjusted R² ≤ 0.59). The combination of cortical thickness and backscatter parameter provided similar or better prediction accuracies than aBMD. For the first time, a method for the noninvasive assessment of the pore diameter distribution in cortical bone by ultrasound is proposed. The combined assessment of cortical thickness, sound velocity, and pore size distribution in a mobile, nonionizing measurement system could have a major impact on preventing osteoporotic fractures.

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.

Safety Assessment of a Wearable Low-Intensity Pulsed Ultrasound Device for Relieving Mental Illness Symptoms

Depression is expected to be one of the significant global medical burdens. Ultrasound therapy, with much-encouraging evidence, has been demonstrated to have a beneficial effect on alleviating mental illness symptoms (neuropsychiatric conditions). However, the study of estimating the risk of using transcranial therapeutic ultrasound has barely been investigated. In this experiment, we develop a wearable head-mounted LIPUS device and assessed the possible tissue damage when applying the LIPUS stimulation on the brain. Our computer simulation and in-vitro experiment results demonstrate that the low-intensity ultrasound (LIPUS) device can safely deliver small doses of low-intensity pulsed ultrasound through the skull into the brain without significant thermal injury. The preliminary results indicate that this modality has the potential for the transcranial treatment of neurological diseases in safety.

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.

Evaluation of Structural Anisotropy in a Porous Titanium Medium Mimicking Trabecular Bone Structure Using Mode-Converted Ultrasonic Scattering

The mode-converted (longitudinal to transverse, L-T) ultrasonic scattering method was utilized to characterize the structural anisotropy of a phantom mimicking the structural properties of trabecular bone. The sample was fabricated using metal additive manufacturing from high-resolution computed tomography (CT) images of a sample of trabecular horse bone with strong anisotropy. Two focused transducers were used to perform the L-T ultrasonic measurements. A normal incidence transducer was used to transmit longitudinal ultrasonic waves into the sample, while the scattered transverse signals were received by an oblique incidence transducer. At multiple locations on the sample, four L-T measurements were performed by collecting ultrasonic scattering from four directions. The amplitude of the root mean square (rms) of the collected ultrasonic scattering signals was calculated for each L-T measurement. The ratios of rms amplitudes for L-T measurements in different directions were calculated to characterize the anisotropy of sample. The results show that the amplitude of L-T converted scattering is highly dependent on the direction of microstructural anisotropy. A strong anisotropy of the microstructure was observed, which coincides with simulation results previously published on the same structure as well as with the anisotropy estimated from the CT images. These results suggest the potential of mode-converted ultrasonic scattering methods to assess the anisotropy of materials with porous, complex structures, including trabecular bone.

Patient-Specific 3-Dimensional Model for High-Intensity Focused Ultrasound Treatment Through the Rib Cage: A Preliminary Study

OBJECTIVES

The purpose of this study was to develop a patient-specific 3-dimensional model for high-intensity focused ultrasound (HIFU) treatment through the rib cage using patient data.

METHODS

Experimental testing to derive parameters used in defining the amount of energy and alteration needed in treatment protocols for upper abdominal disorders under the rib cage was performed. Reconstructed rib cage models based on patient data, tissue-mimicking material phantoms, and magnetic resonance imaging-guided HIFU using a multielement phased array transducer were used in the experiments. Changes in the focal temperature, acoustic power, and acoustic pressure distribution were investigated with and without the presence of the rib cage model. An ExAblate system (InSightec Ltd, Tirat Carmel, Israel) was used to sonicate phantoms by varying the target phantom or rib cage model location.

RESULTS

The effect of the rib cage on the acoustic pressure distribution and acoustic power was closely related to the anatomic structures of the ribs. Thermometry revealed that heating at the focus could be controlled by changing either the power or duration of HIFU application to improve the focal temperature change. The focal temperature change was found to be related to the distance between the rib cage model and focus and the shadow area on the transducer elements covered by the rib cage model in the beam path.

CONCLUSIONS

Experimental results suggest that the rib cage model is a valuable and useful tool that can provide realistic human anatomic structures and properties for evaluating the effects of the rib cage on ultrasound propagation.

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.

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.

Influence of cancellous bone microstructure on ultrasonic attenuation: a theoretical prediction

Background

Quantitative ultrasound has been used for the assessment of cancellous bone status. The attenuation mechanisms of cancellous bone, however, have not been well understood, because the microstructure of cancellous bone is significantly inhomogeneous and the interaction between ultrasound and the microstructure of cancellous bone is complex. In this study, a theoretical approach was applied to investigate the influence of the microstructure of cancellous bone on ultrasonic attenuation.

Results

The scattering from a trabecular cylinder was significantly angle dependent. The dependencies of the ultrasonic attenuation on frequency, scatterer size, and porosity were explored from the theoretical calculation. Prediction results showed that the ultrasonic attenuation increased with the increase of frequency and decreased linearly with the increase in porosity, and the broadband ultrasound attenuation decreased with the increase in porosity. All these predicted trends were consistent with published experimental data. In addition, our model successfully explained the principle of broadband ultrasound attenuation measurement (i.e., the attenuation over the frequency range 0.3–0.65 MHz was approximately linearly proportional to frequency) by considering the contributions of scattering and absorption to attenuation.

Conclusion

The proposed theoretical model may be a potentially valuable tool for understanding the interaction of ultrasound with cancellous bone.

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.

Correlations of the frequency dependence of the ultrasonic backscatter coefficient with the bone volume fraction and the trabecular thickness in bovine trabecular bone: Application of the binary mixture model

The ultrasonic backscatter coefficient and the exponent n (frequency dependence of the backscatter coefficient) were measured in 24 bovine femoral trabecular bone samples. The binary mixture model for ultrasonic scattering from trabecular bone was applied to predict the variations of the ultrasound parameters with the bone volume fraction (BV/TV) and the trabecular thickness (Tb.Th) in trabecular bone. The backscatter coefficient exhibited significant, positive correlations with the BV/TV (R = 0.82) and the Tb.Th (R = 0.79). In contrast, the exponent n was found to be significantly, negatively correlated with the BV/TV (R = -0.77) and the Tb.Th, (R = -0.71).

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.

Velocity dispersion and backscatter in marrow-filled and water-filled trabecular bone samples in vitro

The phase velocity and the backscatter coefficient were measured in 28 bovine femoral trabecular bone samples filled with marrow and water in vitro from 0.2 to 0.6 MHz. The phase velocities decreased approximately linearly with increasing frequency and the average dispersion rate of -34 ms^-1MHz^-1 in the marrow-filled samples was higher than that of -42 ms^-1MHz^-1 in the water-filled samples. The backscatter coefficients exhibited nonlinear, monotonically increasing dependences on the frequency and the average value of the exponent n = 2.92 (frequency dependence) in the marrow-filled samples was higher than the value of n = 2.79 in the water-filled samples.

Acoustic properties of porous microlattices from effective medium to scattering dominated regimes

Microlattices are architected materials that allow for an unprecedented control of mechanical properties (e.g., stiffness, density, and Poisson's coefficient). In contrast to their quasi-static mechanical properties, the acoustic properties of microlattices remain largely unexplored. This paper analyzes the acoustic response of periodic millimeter-sized microlattices immersed in water using experiments and numerical simulations. Microlattices are fabricated using high-precision stereolithographic three-dimensional printing in a large variety of porosities and lattice topologies. This paper shows that the acoustic propagation undergoes a frequency dependent transition from a classic poroelastic behaviour that can be described by Biot's theory to a regime that is dominated by scattering effects. Biot's acoustic parameters are derived from direct simulations of the microstructure using coupled fluid and solid finite elements. The wave speeds predicted with Biot's theory agree well with the experimental measures. Within the scattering regime, the signals show a strong attenuation and dispersion, which is characterized by a cut-off frequency. The strong dispersion results in a frequency dependent group velocity. A simplified model of an elastic cylindrical scatterer allows predicting the signal attenuation and dispersion observed experimentally. The results in this paper pave the way for the creation of microlattice materials for the control of ultrasonic waves across a wide range of frequencies.

Variability in Ultrasound Backscatter Induced by Trabecular Microstructure Deterioration in Cancellous Bone

To determine the relationship between the ultrasonic backscatter parameters and trabecular microstructural variations in cancellous bone, three erosion procedures were performed to simulate various changes in the cancellous bone microstructure. The finite difference time domain (FDTD) method was used to simulate the backscatter signal in cancellous bone. Ultrasonic backscatter properties were derived as functions of the porosity when the ultrasound incident directions were perpendicular and parallel to the major trabeculae direction (MTD), respectively. The variability in the apparent backscatter coefficient (ABC) and apparent integrated backscatter (AIB) due to the trabecular microstructure was revealed. Significant negative correlations between the backscatter parameters (ABC and AIB) and the porosity of the cancellous bone were observed. The simulations showed that the ABC and AIB were influenced by the direction of the trabecular microstructural variations. The linear regressions between the ultrasonic backscatter parameters (ABC and AIB) and the porosity showed significantly different slopes for three erosion procedures when they are ultrasonically perpendicular (for ABC, −1.22 dB, −0.98 dB, and −0.46 dB; for AIB, −0.74 dB, −0.69 dB, and −0.25 dB) and parallel (for ABC, −1.87 dB, −0.69 dB, and −0.51 dB; for AIB, −0.9 dB, −0.5 dB, and −0.34 dB) to the MTD. This paper investigated the relationship between ultrasonic backscatter and cancellous bone microstructure deterioration and indicated that the ultrasonic backscatter could be affected by cancellous bone microstructure deterioration direction.

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.

Monitoring Microwave Ablation of Ex Vivo Bovine Liver Using Ultrasonic Attenuation Imaging

Thermal ablation of soft tissue changes the tissue microstructure and, consequently, induces changes in its acoustic properties. Although B-mode ultrasound provides high-resolution and high-frame-rate images of ablative therapeutic procedures, it is not particularly effective at delineating boundaries of ablated regions because of poor contrast in echogenicity between ablated and surrounding normal tissue. Quantitative ultrasound techniques can provide quantitative estimates of acoustic properties, such as backscatter and attenuation coefficients, and differentiate ablated and unablated regions more effectively, with the potential for monitoring minimally invasive thermal therapies. In this study, a previously introduced attenuation estimation method was used to create quantitative attenuation coefficient maps for 11 microwave ablation procedures performed on refrigerated ex vivo bovine liver. The attenuation images correlate well with the pathologic images of the ablated region. The mean attenuation coefficient for regions of interest drawn inside and outside the ablated zones were 0.9 (±0.2) and 0.45 (±0.15) dB/cm/MHz, respectively. These estimates agree with reported values in the literature and establish the usefulness of non-invasive attenuation imaging for monitoring therapeutic procedures in the liver.

Lower Bound on Estimation Variance of the Ultrasonic Attenuation Coefficient Using the Spectral-Difference Reference-phantom Method

Ultrasonic attenuation is one of the primary parameters of interest in Quantitative Ultrasound (QUS). Non-invasive monitoring of tissue attenuation can provide valuable diagnostic and prognostic information to the physician. The Reference Phantom Method (RPM) was introduced as a way of mitigating some of the system-related effects and biases to facilitate clinical QUS applications. In this paper, under the assumption of diffuse scattering, a probabilistic model of the backscattered signal spectrum is used to derive a theoretical lower bound on the estimation variance of the attenuation coefficient using the Spectral-Difference RPM. The theoretical lower bound is compared to simulated and experimental attenuation estimation statistics in tissue-mimicking (TM) phantoms. Estimation standard deviation (STD) of the sample attenuation in a region of interest (ROI) of the TM phantom is measured for various combinations of processing parameters, including Radio-Frequency (RF) data block length (i.e., window length) from 3 to 17 mm, RF data block width from 10 to 100 A-lines, and number of RF data blocks per attenuation estimation ROI from 3 to 10. In addition to the Spectral-Difference RPM, local attenuation estimation for simulated and experimental data sets was also performed using a modified implementation of the Spectral Fit Method (SFM). Estimation statistics of the SFM are compared to theoretical variance predictions from the literature.1 Measured STD curves are observed to lie above the theoretical lower bound curves, thus experimentally verifying the validity of the derived bounds. This theoretical framework benefits tissue characterization efforts by isolating processing parameter ranges that could provide required precision levels in estimation of the ultrasonic attenuation coefficient using Spectral Difference methods.

Optimum Diffraction-Corrected Frequency-Shift Estimator of the Ultrasonic Attenuation Coefficient

The ultrasonic attenuation coefficient is an important parameter that has been studied extensively in Quantitative Ultrasound and Tissue Characterization. There are various methods described in the literature that estimate this parameter by measuring either spectral difference (i.e., decay) or spectral shift of the backscattered echo signal. Under ideal conditions, i.e., in the absence of abrupt changes in tissue backscattering, Spectral Difference methods can produce estimates with high accuracy and precision. On the other hand, diffraction-corrected Spectral Shift methods (e.g., the Hybrid method) are better suited for application in practical settings using clinical ultrasound scanners. However, current Spectral Shift methods use inefficient frequency shift estimators that ultimately degrade the quality of attenuation coefficient estimates. In this paper, a probabilistic model of the backscattered radiofrequency (RF) echo is used to derive the Cramér-Rao lower bound (CRLB) on estimation variance of the spectral centroid. Next, an efficient correlation-based shift estimator is presented that achieves the CRLB. Used in conjunction with a well-characterized reference phantom to correct for diffraction and other system-related effects, this estimator greatly improves the accuracy and precision of Spectral- Shift attenuation estimation. A theoretical analysis of this method is provided, and its performance is quantitatively compared with that of the Hybrid method using simulated and experimental phantom studies. A minimum of 3-fold reduction in the standard deviation of attenuation coefficient estimates is observed using the new method.

Measuring mass density and ultrasonic wave velocity: A wavelet-based method applied in ultrasonic reflection mode

When assessing ultrasonic measurements of material parameters, the signal processing is an important part of the inverse problem. Measurements of thickness, ultrasonic wave velocity and mass density are required for such assessments. This study investigates the feasibility and the robustness of a wavelet-based processing (WBP) method based on a Jaffard-Meyer algorithm for calculating these parameters simultaneously and independently, using one single ultrasonic signal in the reflection mode. The appropriate transmitted incident wave, correlated with the mathematical properties of the wavelet decomposition, was determined using a adapted identification procedure to build a mathematically equivalent model for the electro-acoustic system. The method was tested on three groups of samples (polyurethane resin, bone and wood) using one 1-MHz transducer. For thickness and velocity measurements, the WBP method gave a relative error lower than 1.5%. The relative errors in the mass density measurements ranged between 0.70% and 2.59%. Despite discrepancies between manufactured and biological samples, the results obtained on the three groups of samples using the WBP method in the reflection mode were remarkably consistent, indicating that it is a reliable and efficient means of simultaneously assessing the thickness and the velocity of the ultrasonic wave propagating in the medium, and the apparent mass density of material.

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.

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.

Nonlinear attenuation and dispersion in human calcaneus in vitro: statistical validation and relationships to microarchitecture

Through-transmission measurements were performed on 30 human calcaneus samples in vitro. Nonlinear attenuation and dispersion measurements were investigated by estimating 95% confidence intervals of coefficients of polynomial expansions of log magnitude and phase of transmission coefficients. Bone mineral density (BMD) was measured with dual x-ray absorptiometry. Microarchitecture was measured with microcomputed tomography. Statistically significant nonlinear attenuation and nonzero dispersion were confirmed for a clinical bandwidth of 300-750 kHz in 40%-43% of bone samples. The mean linear coefficient for attenuation was 10.3 dB/cm MHz [95% confidence interval (CI): 9.0-11.6 dB/cm MHz]. The mean quadratic coefficient for attenuation was 1.6 dB/cm MHz(2) (95% CI: 0.4-2.8 dB/cm MHz(2)). Nonlinear attenuation provided little information regarding BMD or microarchitecture. The quadratic coefficient for phase (which is related to dispersion) showed moderate correlations with BMD (r = -0.65; 95% CI: -0.82 to -0.36), bone surface-to-volume ratio (r = 0.47; 95% CI: 0.12-0.72) and trabecular thickness (r = -0.40; 95% CI: -0.67 to -0.03). Dispersion was proportional to bone volume fraction raised to an exponent of 2.1 ± 0.2, which is similar to the value for parallel nylon-wire phantoms (2.4 ± 0.2) and supports a multiple-scattering model for dispersion.

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.

Numerical investigation of ultrasound reflection and backscatter measurements in cancellous bone on various receiving areas

In this study, new ultrasound reflection and backscatter measurements in cancellous bone using a membrane-type hydrophone are proposed. A membrane hydrophone made of a piezoelectric polymer film mounted on an annular frame allows an incident ultrasound wave to pass through its aperture because it has no backing material. Therefore, in measurements using the membrane hydrophone, the receiving area could be located independently from the transmitting area. In addition, the size and shape of the receiving area, which corresponded to those of the electrode deposited on the piezoelectric film, could be arranged in various ways. To investigate the validity of the proposed measurements, before bench-top experiments, the reflected and backscattered waves from cancellous bone were numerically simulated using a finite-difference time-domain method. The reflection and backscatter parameters were measured on various receiving areas, and their correlation coefficients with the structural parameters in the cancellous bone were derived. The simulated results suggested that appropriate receiving areas for the reflection and backscatter measurements could exist and that the proposed measurements could be more effective for evaluating bone properties than conventional measurements.

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.

Dependences of quantitative ultrasound parameters on frequency and porosity in water-saturated nickel foams

The frequency-dependent phase velocity, attenuation coefficient, and backscatter coefficient were measured from 0.8 to 1.2 MHz in 24 water-saturated nickel foams as trabecular-bone-mimicking phantoms. The power law fits to the measurements showed that the phase velocity, the attenuation coefficient, and the backscatter coefficient were proportional to the frequency with exponents n of 0.95, 1.29, and 3.18, respectively. A significant linear correlation was found between the phase velocity at 1.0 MHz and the porosity. In contrast, the best regressions for the normalized broadband ultrasound attenuation and the backscatter coefficient at 1.0 MHz were obtained with the polynomial fits of second order.

Numerical investigation of ultrasonic attenuation through 2D trabecular bone structures reconstructed from CT scans and random realizations

In this paper, we compare ultrasound interrogations of actual CT-scanned images of trabecular bone with artificial randomly constructed bone. Even though it is known that actual bone does not have randomly distributed trabeculae, we find that the ultrasound attenuations are close enough to cast doubt on any microstructural information, such as trabeculae width and distance between trabeculae, being gleaned from such experiments. More precisely, we perform numerical simulations of ultrasound interrogation on cancellous bone to investigate the phenomenon of ultrasound attenuation as a function of excitation frequency and bone porosity. The theoretical model is based on acoustic propagation equations for a composite fluid-solid material and is solved by a staggered-grid finite-difference scheme in the time domain. Numerical experiments are performed on two-dimensional bone samples reconstructed from CT-scanned images of real human calcaneus and from random distributions of fluid-solid particles generated via the turning bands method. A detailed comparison is performed on various parameters such as the attenuation rate and speed of sound through the bone samples as well as the normalized broadband ultrasound attenuation coefficient. Comparing results from these two types of bone samples allows us to assess the role of bone microstructure in ultrasound attenuation. It is found that the random model provides suitable bone samples for ultrasound interrogation in the transverse direction of the trabecular network.

Frequency specific ultrasound attenuation is sensitive to trabecular bone structure

This study investigated the efficacy of frequency modulated ultrasound attenuation in the assessment of the trabecular structural properties. Four frequency modulated signals were created to represent four frequency bands centered at 500 kHz, 900 kHz, 1.3 MHz and 1.7 MHz with the bandwidth of 400 kHz. Five 1-cm trabecular cubes were harvested from fresh bovine distal femur. The cubes underwent four steps of demineralization process to expand the sample size to 25 with the greater variations of the structural properties for the better correlation study. Pearson correlation study was performed between the ultrasound attenuation in four frequency bands and the trabecular structural properties. The results showed that correlations of frequency modulated ultrasound attenuation to the trabecular structural properties are dependent on frequency bands. The attenuation in proximal-distal orientation had the highest correlation to BV/TV (R(2) = 0.73, p < 0.001) and trabecular thickness (R(2) = 0.50, p < 0.001) at the frequency band centered at 1.7 MHz. It was equivalent in the four frequency bands in correlation to the trabecular number (average R(2) = 0.80, p < 0.001) and to the trabecular separation (average R(2) = 0.83, p < 0.001). The attenuation in anterio-posterial orientation had the highest correlation to BV/TV (R(2) = 0.80, p < 0.001) and trabecular thickness (R(2) = 0.71, p < 0.001) at the frequency band centered at 1.3 MHz. The attenuation in the first frequency band was the most sensitive to the trabecular number (R(2) = 0.71, p < 0.001) and trabecular separation (R(2) = 0.80, p < 0.001). No significant correlation was observed for the attenuation in medial-lateral orientation across the four frequency bands.

Comparison of phase velocity in trabecular bone mimicking-phantoms by time domain numerical (EFIT) and analytical multiple scattering approaches

The corrected Waterman-Truell model and the Elastodynamic Finite Integration Technique were used to analyze the ultrasonic wave dispersion in trabecular bones mimicking phantoms. A simple two-phase model of the trabecular bone is assumed; the trabeculae structure and the bone marrow. The phase velocity for frequencies within the range from 400kHz to 800kHz were computed for different scatterer arrays varying their dimensions and number. The theoretical and numerical results were compared to experimental published data, obtained from a mimicking phantom composed by a periodic array of nylon shreds (trabeculae array) immersed in a water tank. Our results showed an excellent consistency when compared to experimental data. The negative dispersions of -8.48m/s/MHz and -9.16m/s/MHz were computed by the multiple scattering method and the numerical approach, respectively, where the latter is closer to the experimental dispersion of -12.09m/s/MHz. Similar result has been reported in the literature, where the dispersion predicted by the Generalized Self-Consistent Method [J. Acoust. Soc. Am. 124 (2008) 4047] is -9.96m/s/MHz.

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.

Investigation of the influence of reflection on the attenuation of cancellous bone

The model proposed in this paper is based on the fact that the reflection might have a significant contribution to the attenuation of the acoustic waves propagating through the cancellous bone. The numerical implementation of the mentioned effect is realized by the development of a new representative volume element that includes an infinitesimally thin 'transient' layer on the contact surface of the bone and the marrow. This layer serves to model the amplitude transformation of the incident wave by the transition through media with different acoustic impedances and to take into account the energy loss due to the reflection. The proposed representative volume element together with the multiscale finite element is used to simulate the wave propagation and to evaluate the attenuation coefficient for samples with different effective densities in the dependence of the applied excitation frequency. The obtained numerical values show a very good agreement with the experimental results. Moreover, the model enables the determination of the upper and the lower bound for the attenuation coefficient.

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.

Frequency-dependent attenuation and backscatter coefficients in bovine trabecular bone from 0.2 to 1.2 MHz

The frequency-dependent attenuation and backscatter coefficients were measured in 25 bovine femoral trabecular bone samples from 0.2 to 1.2 MHz. When the average attenuation coefficient was fitted to a nonlinear power law α(f)=α(0)+α(1)f(n), the exponent n was found to be 1.65. In contrast, the average backscatter coefficient was fitted to a power law η(f)=η(1)f(n) and the exponent n was measured as 3.25. The apparent bone density was significantly correlated with the parameter α(1) (0.2-0.7 MHz: r = 0.852, 0.6-1.2 MHz: r = 0.832) as well as the backscatter coefficient (0.5 MHz: r = 0.751, 1.0 MHz: r = 0.808).

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.

Development and characterization of a tissue-mimicking material for high-intensity focused ultrasound

A tissue-mimicking material (TMM) for the acoustic and thermal characterization of high-intensity focused ultrasound (HIFU) devices has been developed. The material is a high-temperature hydrogel matrix (gellan gum) combined with different sizes of aluminum oxide particles and other chemicals. The ultrasonic properties (attenuation coefficient, speed of sound, acoustical impedance, and the thermal conductivity and diffusivity) were characterized as a function of temperature from 20 to 70°C. The backscatter coefficient and nonlinearity parameter B/A were measured at room temperature. Importantly, the attenuation coefficient has essentially linear frequency dependence, as is the case for most mammalian tissues at 37°C. The mean value is 0.64f(0.95) dB·cm(-1) at 20°C, based on measurements from 2 to 8 MHz. Most of the other relevant physical parameters are also close to the reported values, although backscatter signals are low compared with typical human soft tissues. Repeatable and consistent temperature elevations of 40°C were produced under 20-s HIFU exposures in the TMM. This TMM is appropriate for developing standardized dosimetry techniques, validating numerical models, and determining the safety and efficacy of HIFU devices.

The 25th Anniversary of BUA for the Assessment of Osteoporosis: Time for a New Paradigm?

The measurement of broadband ultrasonic attenuation (BUA) in cancellous bone at the calcaneus for the assessment of osteoporosis was first described within this journal 25 years ago. It was recognized in 2006 by Universities UK as being one of the ‘100 discoveries and developments in UK Universities that have changed the world’ over the past 50 years. In 2008, the UK's Department of Health also recognized BUA assessment of osteoporosis in a publication highlighting 11 projects that have contributed to ‘60 years of NHS research benefiting patients’. The BUA technique has been extensively clinically validated and is utilized worldwide, with at least seven commercial systems currently providing calcaneal BUA measurement. However, there is still no fundamental understanding of the dependence of BUA upon the material and structural properties of cancellous bone. This review aims to provide an ‘engineering in medicine’ perspective and proposes a new paradigm based upon phase cancellation due to variation in propagation transit time across the receive transducer face to explain the non-linear relationship between BUA and bone volume fraction in cancellous bone.

Padé approximants and their application to scattering from fluid media

In this work, a numerical method for modeling the scattered acoustic pressure from fluid occlusions is described. The method is based on the asymptotic series expansion of the pressure expressed in terms of sound speed contrast between the host medium and entrained fluid occlusions. Padé approximants are used to extend the applicability of the result for larger values of sound speed contrast. For scattering from a circular cylinder, an improvement in convergence between the exact and numerical solutions is demonstrated. In the case of scattering from an inhomogeneous medium, a numerical solution with reduced order of Padé approximants is presented.

Frequency dependence of average phase shift from human calcaneus in vitro

If dispersion in a medium is weak and approximately linear with frequency (over the experimental band of frequencies), then it can be shown that the constant term in a polynomial representation of phase shift as a function of frequency can produce errors in measurements of phase-velocity differences in through-transmission, substitution experiments. A method for suppressing the effects of the constant phase shift in the context of the single-wave-model was tested on measurements from 30 cancellous human calcaneus samples in vitro. Without adjustment for constant phase shifts, the estimated phase velocity at 500 kHz was 1516+/-6 m/s (mean+/-standard error), and the estimated dispersion was -24+/-4 m/s MHz (mean+/-standard error). With adjustment for constant phase shifts, the estimated mean velocity decreased by 4-9 m/s, and the estimated magnitude of mean dispersion decreased by 50%-100%. The average correlation coefficient between the measured attenuation coefficient and frequency was 0.997+/-0.0026 (mean+/-standard deviation), suggesting that the signal for each sample was dominated by one wave. A single-wave, linearly dispersive model conformed to measured complex transfer functions from the 30 cancellous-bone samples with an average root-mean-square error of 1.9%+/-1.0%.

Influence of the filling fluid on frequency-dependent velocity and attenuation in cancellous bones between 0.35 and 2.5 MHz

The paper is focused on experiments on human cancellous bones filled with different fluids with the goal of evaluating their contribution to velocity dispersion, absorption, and scattering mechanisms. The specimens were measured first filled with marrow and subsequently, after marrow removal, with water and alcohol. No significant influence of the fluids was evidenced on the attenuation coefficient. Given the absence of impact of viscosity of the saturating fluid, the authors hypothesized that the source of attenuation is associated with viscoelastic absorption in the solid trabeculae and with scattering. Alteration of scattering obtained by changing the acoustic impedance mismatch between the fluid (alcohol vs water) and the trabeculae was reflected neither in the attenuation nor in its slope. This led the authors to suggest that longitudinal-to-shear scattering together with absorption in the solid phase are candidates as main sources for the attenuation. The differences in velocity values indicate that the elastic properties of the fluid are main determinants of the phase velocity. This finding is particularly significant in the context of /in vivo/ measurements, because it demonstrates that the subject-dependent properties of marrow may partly explain the inter-subject variability of speed of sound values.