Anisotropy of ultrasonic backscatter and attenuation from human calcaneus: implications for relative roles of absorption and scattering in determining attenuation

A Geometric Model of Ultrasound Backscatter to Describe Microstructural Anisotropy of Tissue

Methods to assess ultrasound backscatter anisotropy from clinical array transducers have recently been developed. However, they do not provide information about the anisotropy of microstructural features of the specimens. This work develops a simple geometric model, referred to as the secant model, of backscatter coefficient anisotropy. Specifically, we evaluate anisotropy of the frequency dependence of the backscatter coefficient parameterized in terms of effective scatterer size. We assess the model in phantoms with known scattering sources and in a skeletal muscle, a well-known anisotropic tissue. We demonstrate that the secant model can determine the orientation of the anisotropic scatterers, as well as accurately determining effective scatterer sizes, and it may classify isotropic versus anisotropic scatterers. The secant model may find utility in monitoring disease progression as well as characterizing normal tissue architectures.

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.

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.

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.

Anisotropy and Spatial Heterogeneity in Quantitative Ultrasound Parameters: Relevance to the Study of the Human Cervix

Imaging biomarkers based on quantitative ultrasound can offer valuable information about properties that inform tissue function and behavior such as microstructural organization (e.g., collagen alignment) and viscoelasticity (i.e., compliance). For example, the cervix feels softer as its microstructure remodels during pregnancy, an increase in compliance that can be objectively quantified with shear wave speed and therefore shear wave speed estimation is a potential biomarker of cervical remodeling. Other proposed biomarkers include parameters derived from the backscattered echo signal, such as attenuation and backscattered power loss, because such parameters can provide insight into tissue microstructural alignment and organization. Of these, attenuation values for the pregnant cervix have been reported, but large estimate variance reduces their clinical value. That said, parameter estimates based on the backscattered echo signal may be incorrect if assumptions they rely on, such as tissue isotropy and homogeneity, are violated. For that reason, we explored backscatter and attenuation parameters as potential biomarkers of cervical remodeling via careful investigation of the assumptions of isotropy and homogeneity in cervical tissue. Specifically, we estimated the angle- and spatial-dependence of parameters of backscattered power and acoustic attenuation in the ex vivo human cervix, using the reference phantom method and electronic steering of the ultrasound beam. We found that estimates are anisotropic and spatially heterogeneous, presumably because the tissue itself is anisotropic and heterogeneous. We conclude that appropriate interpretation of imaging biomarkers of cervical remodeling must account for tissue anisotropy and heterogeneity.

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

Quantifying Backscatter Anisotropy Using the Reference Phantom Method

Acoustic properties can be exploited to infer and evaluate tissue microstructure. However, common assumptions are that the medium of interest is homogeneous and isotropic, and that its underlying physical properties cause diffuse scattering. In this paper, we describe how we developed and tested novel parameters designed to address isotropy/anisotropy in backscattered echo signal power in complex biological tissues. Specifically, we explored isotropy/anisotropy in backscattered power in isotropic phantoms (spherical glass beads), an anisotropic phantom (dialysis phantom with rodlike fibers), and an in vivo human tissue with well-described anisotropy (bicep muscle). Our approach uses the reference phantom method to compensate for system transfer and diffraction losses when electronically beamsteering a linear array transducer. We define three parameters to quantify the presence and orientation of anisotropic scatterers, as well as address magnitude of anisotropy. We found that these parameters can detect and sense the degree of anisotropy in backscatter in both phantoms and bicep muscle. Bias of the summary anisotropy parameters, induced through a speed of sound mismatch of sample media and reference phantom, was less than 0.2 dB if the speed of sound was within ±20 m/s of the sample media. In summary, these new parameters may be useful for testing the assumption of isotropy as well as providing more detailed information about the underlying microstructural sources of backscatter in complex biological tissues.

Microstructural characterization of trabecular bone using ultrasonic backscattering and diffusion parameters

Finite differences time domain methods were utilized to simulate ultrasound propagation and scattering in anisotropic trabecular bone structures obtained from high resolution Computed Tomography (CT). The backscattered signals were collected and the incoherent contribution was extracted. The diffusion constant was calculated for propagations along and across the main direction of anisotropy, and was used to characterize the anisotropy of the trabecular microstructures. In anisotropic structures, the diffusion constant was significantly different in both directions, and the anisotropy of the diffusion constant was strongly correlated to the structural anisotropy measured on the CT images. These results indicate that metrics based on diffusion can be used to quantify the anisotropy of complex structures such as trabecular bone.

Solid volume fraction estimation of bone:marrow replica models using ultrasound transit time spectroscopy

The acceptance of broadband ultrasound attenuation (BUA) for the assessment of osteoporosis suffers from a limited understanding of both ultrasound wave propagation through cancellous bone and its exact dependence upon the material and structural properties. It has recently been proposed that ultrasound wave propagation in cancellous bone may be described by a concept of parallel sonic rays; the transit time of each ray defined by the proportion of bone and marrow propagated. A Transit Time Spectrum (TTS) describes the proportion of sonic rays having a particular transit time, effectively describing the lateral inhomogeneity of transit times over the surface aperture of the receive ultrasound transducer. The aim of this study was to test the hypothesis that the solid volume fraction (SVF) of simplified bone:marrow replica models may be reliably estimated from the corresponding ultrasound transit time spectrum. Transit time spectra were derived via digital deconvolution of the experimentally measured input and output ultrasonic signals, and compared to predicted TTS based on the parallel sonic ray concept, demonstrating agreement in both position and amplitude of spectral peaks. Solid volume fraction was calculated from the TTS; agreement between true (geometric calculation) with predicted (computer simulation) and experimentally-derived values were R(2)=99.9% and R(2)=97.3% respectively. It is therefore envisaged that ultrasound transit time spectroscopy (UTTS) offers the potential to reliably estimate bone mineral density and hence the established T-score parameter for clinical osteoporosis assessment.

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.

Enhanced correlation between quantitative ultrasound and structural and mechanical properties of bone using combined transmission-reflection measurement

Quantitative ultrasound (QUS) is capable of predicting the principal structural orientation of trabecular bone; this orientation is highly correlated with the mechanical strength of trabecular bone. Irregular shape of bone, however, would increase variation in such a prediction, especially under human in vivo measurement. This study was designed to combine transmission and reflection modes of QUS measurement to improve the prediction for the structural and mechanical properties of trabecular bone. QUS, mechanical testing, and micro computed tomography (μCT) scanning were performed on 24 trabecular bone cubes harvested from a bovine distal femur to obtain the mechanical and structural parameters. Transmission and reflection modes of QUS measurement in the transverse and frontal planes were performed in a confined 60° angle range with 5° increment. The QUS parameters, attenuation (ATT) and velocity (UV), obtained from transmission mode, were normalized to the specimen thickness acquired from reflection mode. Analysis of covariance showed that the combined transmission-reflection modes improved prediction for the structural and Young's modulus of bone in comparison to the traditional QUS measurement performed only in the medial-lateral orientation. In the transverse plane, significant improvement between QUS and μCT was found in ATT vs bone surface density (BS/BV) (p < 0.05), ATT vs trabecular thickness (Tb.Th) (p < 0.01), ATT vs degree of anisotropy (DA) (p < 0.05), UV vs trabecular bone number (Tb.N) (p < 0.05), and UV vs Tb.Th (p < 0.001). In the frontal plane, significant improvement was found in ATT vs structural model index (SMI) (p < 0.01), ATT vs bone volume fraction (BV/TV) (p < 0.01), ATT vs BS/BV (p < 0.001), ATT vs Tb.Th (p < 0.001), ATT vs DA (p < 0.001), and ATT vs modulus (p < 0.001), UV vs SMI (p < 0.01), UV vs BV/TV (p < 0.05), UV vs BS/BV (p < 0.05), UV vs Tb.Th (p < 0.01), UV vs trabecular spacing (p < 0.05), and UV vs modulus (p < 0.01). These data suggested that the combined transmission-reflection QUS method is capable of providing information more relevant to the structural and mechanical properties of trabecular bone.

Analysis of apparent integrated backscatter coefficient and backscattered spectral centroid shift in Calcaneus in vivo for the ultrasonic evaluation of osteoporosis

The purposes of our study were to evaluate the correlation among apparent integrated backscatter coefficient (AIB), spectral centroid shift (SCS) of ultrasonic backscatter signals and bone mineral density (BMD) and to examine the effectiveness of ultrasound variables as predictors of osteoporosis. A total of 1011 persons aged 21-80 y old were included. All study participants underwent BMD measurements of the lumbar spine (LSBMD) and the femoral neck (FNBMD). The participants also underwent calcaneal measurements to determine AIB and SCS with central frequencies of 3.5 (one transducer) and 5.0 MHz (the other transducer). AIB decreased with age and was positively correlated with BMD, while SCS increased with age and was negatively correlated with BMD. The correlation coefficient of SCS with LSBMD and FNBMD at 3.5 MHz was -0.72 and -0.70, respectively. The correlation coefficient at 5.0 MHz was -0.75 and -0.74, respectively. The correlation coefficient of AIB with LSBMD and FNBMD at 3.5 MHz was 0.65 and 0.63. The correlation coefficient at 5.0 MHz was 0.59 and 0.55, respectively. The correlation between SCS and BMD was significantly better than the correlation between AIB and BMD. Using receiver operating characteristic analysis, a significant difference was found between the areas under the curve for SCS and AIB at 3.5 MHz (0.781 vs. 0.715, respectively, p < 0.05), as well as at 5.0 MHz (0.782 vs. 0.709, respectively, p < 0.05). The optimum T-score threshold for SCS was -1.3 for both transducers. The sensitivity and specificity of SCS at 3.5 MHz and 5.0 MHz for the optimum threshold were 64%, 85%, 63% and 86%, respectively. In conclusion, the correlations among the ultrasound parameters and BMDs are strong. SCS performs better than AIB in differentiating patients with osteoporosis. Ultrasound variables may be taken into consideration as predictors of osteoporosis in the future considering its portability.

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.

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

Characterization of pseudarthrosis with ultrasound backscattered signals in rats

PURPOSE

To propose a novel model of pseudarthrosis in a small animal and to investigate the ability of backscatter parameters from ultrasound signals in differentiating normal bone from those ones with pseudarthrosis.

METHODS

Twelve Rattus norvegicus albinus free from pathogenic species (SPF) were randomly divided in two groups, with six animals each. In the Control group a surgical approach to the femur was made, followed by the synthesis of the muscle and skin layers. The Experimental group was submitted to an osteotomy of the femur and a vascularized flap of the fascia lata was interposed in the line of the fractured bone. Then the alignment and bone stabilization were accomplished, by using nylon stitch in U shape introduced in holes made in the proximal and distal fractured bone. Bone samples were scanned with ultrasound and signals were collected for each one to analyze the parameter Apparent Integrated Backscatter - AIB.

RESULTS

Radiological and anatomopathologic studies revealed the absence of bone consolidation with persistence of fiber-osteoid tissue. Values of the ultrasound parameter AIB from normal bones were statistically different from those with pseudarthrosis.

CONCLUSION

The experimental model was suitable for pseudarthrosis development in rats and the ultrasound backscatter parameters were able to identify such a bone disease in vitro.

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.

Ultrasound backscatter imaging provides frequency-dependent information on structure, composition and mechanical properties of human trabecular bone

The strength as well as the acoustic properties of trabecular bone are determined by its structure and composition. Consequently, tissue structure and compositional properties also affect the ultrasound propagation in bone. The diagnostic potential of ultrasound has not been fully exploited in clinical quantitative ultrasound devices. The aim of this study was to investigate the ability of quantitative ultrasound pulse-echo imaging, conducted over a broad range of frequencies (1 to 5 MHz), to predict the mechanics, composition and microstructure of trabecular bone. Ultrasound reflection and backscatter parameters correlated significantly with the ultimate strength of the trabecular bone and the bone volume fraction (r=0.76-0.90, n=20, p<0.01). Ultrasound backscatter associated significantly (independently of bone structure or mineral content) with the collagen content of the bone matrix (r=0.75, r(adjusted)=0.66, p<0.01). Interestingly, the applied ultrasound frequency seemed to relate the sensitivity of ultrasound backscatter to different properties of trabecular bone. At frequencies ranging from 1 to 3.5 MHz, the ultrasound backscatter associated significantly with the tissue mechanical and structural parameters. At 5MHz, the composition of the bone matrix was a more significant determinant of the measured backscatter. This study provides useful information for optimizing the use of pulse-echo measurements, and thereby further emphasizes the diagnostic potential of the ultrasound backscatter measurements of trabecular bone.

Effect of bone marrow on acoustic properties of trabecular bone--3D finite difference modeling study

The composition of bone marrow is influenced by many factors, such as age and diseases. The present numerical study investigates the contribution of marrow on the acoustic measurements of trabecular bone. Cylindrical bone samples (n = 11), extracted from three anatomical sites of human cadaver knees, were imaged with a high-resolution microtomography (microCT). Three-dimensional finite difference time domain (FDTD) models (Wave 3000 Pro 2.2, Cyberlogic Inc., NY, USA) were created using the segmented microCT images of each sample. First, we evaluated the effect of voxel size on the computer resource requirements, morphological parameters and acoustic simulations. Second, the effect of bone marrow on ultrasonic measurements was assessed. The simulations were repeated with two voxel sizes before and after substitution of bone marrow (i.e., fat) with water. The voxel size of the FDTD mesh controlled the fine structure of the modeled calcified matrix and significantly affected the simulation results. However, present simulations showed that the effect of bone marrow on ultrasound parameters can be reliably simulated with the applied voxel sizes of 72 and 90 microm. Ultrasound attenuation and speed were found (p < 0.01) to decrease and increase, respectively, when bone marrow was substituted with water. Moreover, reflection from the surface of the sample increased (p < 0.01) and backscatter from internal structures decreased (p < 0.01) after removal of marrow. The effect of bone marrow on the acoustic properties was stronger in samples with low bone volume fraction. The present results indicate that the amount and quality of bone marrow significantly influence the acoustic properties of trabecular bone. Possible interindividual differences in the composition of bone marrow may increase uncertainty in clinical ultrasound diagnostics of osteoporosis. Importantly, the effect is most significant in osteoporotic low-density bone.

Velocity dispersion in trabecular bone: influence of multiple scattering and of absorption

Speed of sound measurements are widely used clinically to assess bone strength. Trabecular bone is an attenuating composite material in which negative values of velocity dispersion have been measured, this behavior remaining poorly explained physically. The aim of this work is to describe the ultrasonic propagation in trabecular bone modeled by infinite cylinders immersed in a saturating matrix, and to derive the physical determinants of velocity dispersion. A homogenization model accounting for the coupling of multiple scattering and absorption phenomena allows the computation of phase velocity and of dispersion while varying bone properties. The present model is adapted from the generalized self-consistent method developed in the work of Yang and Mal [(1994). "Multiple-scattering of elastic waves in a fiber-reinforced composite," J. Mech. Phys. Solids 42, 1945-1968]. It predicts negative values of velocity dispersion, in agreement with experimental results obtained in phantoms mimicking trabecular bone. In trabecular bone, mostly negative and also positive values of velocity dispersion are predicted by the model, which span within the range of values measured experimentally. Scattering effects are responsible for the negative values of dispersion, whereas the frequency dependence of the absorption coefficient in bone marrow and/or in the trabeculae results in an increase in dispersion, which may then become positive.

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.

Frequency dependence of backscatter from thin, oblique, finite-length cylinders measured with a focused transducer-with applications in cancellous bone

A model is presented for the echo from a thin, oblique, finite-length cylinder. The echo is calculated from the line integral of the transducer directivity pattern along the cylinder axis. The model was validated with broadband pulse-echo measurements from (1) a perpendicular (to the ultrasound beam) nylon wire as a function of lateral displacement from the beam center, (2) a tilted nylon wire as a function of the angle of inclination relative to the ultrasound beam, and (3) a quasi-parallel-nylon-wire phantom, which mimicked the scattering properties of cancellous bone. The transducer directivity pattern (as a function of position and frequency) was measured with a membrane hydrophone. The model predicts an approximately cubic frequency dependence of backscatter coefficient from the phantom, as has been observed experimentally in cancellous bone. The model also predicts the relationship between the cylinder length and the exponent of a power law fit to backscatter coefficient versus frequency, which is 4 for very short (compared to a wavelength) cylinders and asymptotically approaches 3 for very long cylinders.

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.

Relationships of trabecular bone structure with quantitative ultrasound parameters: in vitro study on human proximal femur using transmission and backscatter measurements

The present study was designed to assess the relationships between QUS parameters and bone density or microarchitecture on samples of human femoral trabecular bone. The normalized slope of the frequency-dependent attenuation (nBUA), the speed of sound (SOS) and the broadband ultrasound backscatter coefficient (BUB) were measured on 37 specimens of pure trabecular bones removed from upper parts of fresh human femurs. Bone mineral density (BMD) was assessed using a clinical scanner. Finally, 8 mm diameter cylindrical cores were extracted from the specimens and their microarchitecture was reconstructed after synchrotron radiation microtomography experiments (isotropic resolution of 10 microm). A large number of microarchitectural parameters were computed, describing morphology, connectivity and geometry of the specimens. BMD correlated with all the microarchitectural parameters and the number of significant correlations was found among the architectural parameters themselves. All QUS parameters were significantly correlated to BMD (R=0.83 for nBUA, R=0.81 for SOS and R=0.69 for BUB) and to microarchitectural parameters (R=-0.79 between nBUA and Tb.Sp, R=-0.81 between SOS and Tb.Sp, R=-0.65 between BUB and BS/BV). Using multivariate model, it was found that microstructural parameters adds 10%, 19%, and 4% to the respective BMD alone contribution for the three variables BUA, SOS and BUB. Moreover, the RMSE was reduced by up to 50% for SOS, by up to 21% for nBUA and up to 11% when adding structural variables to BMD in explaining QUS results. Given the sample, which is not osteoporosis-enriched, the added contribution is quite substantial. The variability of SOS was indeed completely explained by a multivariate model including BMD and independent structural parameters (R(2)=0.94). The inverse problem on the data presented here has been addressed using simple and multiple linear regressions. It was shown that the predictions (in terms of R(2) or RMSE) of microarchitectural parameters was not enhanced when combining 2 or 3 QUS in multiple linear regressions compared to the prediction obtained with one QUS parameter alone. The best model was found for the prediction of Tb.Th() from BUA (R(2)=0.58, RMSE=17 microm). Given the high values of RMSE, these linear models appear of limited clinical value, suggesting that appropriate models have to be derived in order to solve the inverse problem. In this regard, a very interesting multivariate model was found for nBUA and BUB with Tb.Th and Tb.N, in agreement with single scattering theories by random medium. However, the source of residual variability of nBUA and BUB (15% and 45% respectively) remained unexplained.

The effect of phase cancellation on estimates of broadband ultrasound attenuation and backscatter coefficient in human calcaneus in vitro

Broadband ultrasound attenuation (BUA) is a clinically proven indicator of osteoporotic fracture risk. BUA measurements are typically performed in throughtransmission with single-element phase sensitive (PS) receivers and therefore can be compromised by phase cancellation artifact. Phase-insensitive (PI) receivers suppress phase cancellation artifact. To study the effect of phase cancellation on BUA measurements, through-transmission measurements were performed on 16 human calcaneus samples in vitro using a two-dimensional receiver array that enabled PS and PI BUA estimation. The means plus or minus standard deviations for BUA measurements were 22.1 +/- 15.8 dB/MHz (PS) and 17.6 +/- 7.2 dB/MHz (PI), suggesting that, on the average, approximately 20% of PS BUA values in vitro can be attributed to phase cancellation artifact. Therefore, although cortical plates are often regarded as the primary source of phase cancellation artifact, the heterogeneity of cancellous bone in the calcaneal interior may also be a significant source. Backscatter coefficient estimates in human calcaneus that are based on PS attenuation compensation overestimate 1) average magnitude of backscatter coefficient at 500 kHz by a factor of about 1.6 +/- 0.3 and 2) average exponent (n) of frequency dependence by about 0.34 +/- 0.12 (where backscatter coefficient is fit to a power law form proportional to frequency to the nth power).

Development of a method for ultrasound-guided placement of pedicle screws

Abstract-Many forms of spinal fusion involve the placement of long screws through the pedicles into the vertebral body. During the procedure, there is substantial risk of damage to vital neural and vascular structures due to the limited visibility of anatomic landmarks and high anatomic variability. As an alternative to current guidance systems, we have investigated the feasibility of performing ultrasound imaging through cancellous bone for the purpose of pedicle screw guidance. Quantitative ultrasonic characterization and A-mode imaging of seven defatted vertebral cancellous bone specimens was performed along the craniocaudal axis in water with unfocused, 1-MHz and 3.5- MHz broadband transducers. The center frequency attenuation increased considerably from 10.5 +/- 4.6 dB/cm at 1 MHz to 24.1 +/- 7.2 dB/cm at 3.5 MHz, while the speed of sound exhibited moderate positive dispersion, increasing from 1489 +/- 4.7 m/s at 1 MHz to 1494 +/- 4.2 m/s at 3.5 MHz. Despite the high attenuation and large specimen thickness (1.0-1.9 cm), A-mode imaging through cancellous bone to detect an aluminum reflector was possible in 83.2% and 70.1% of the cases at 1 MHz and 3.5 MHz, respectively. Specimen boundaries were identifiable with clinically sufficient average accuracy of 1.1 mm and 0.9 mm in the 1 MHz and 3.5 MHz A-mode images, respectively.

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.

Predictions of the modified Biot-Attenborough model for the dependence of phase velocity on porosity in cancellous bone

The modified Biot-Attenborough (MBA) model for acoustic wave propagation in porous media has been found useful to predict wave properties in cancellous bone. The present study is aimed at applying the MBA model to predict the dependence of phase velocity on porosity in cancellous bone. The MBA model predicts a phase velocity that decreases nonlinearly with porosity. The optimum values for input parameters of the MBA model, such as compressional speed c(m) of solid bone and phase velocity parameter s(2), were determined by comparing the predictions with previously published measurements in human calcaneus and bovine cancellous bone. The value of the phase velocity parameter s(2)=1.23 was obtained by curve fitting to the experimental data for 53 human calcaneus samples only, assuming a compressional speed c(m)=2500 m/s of solid bone. The root-mean-square error (RMSE) of the curve fit was 15.3m/s. The optimized value of s(2) for all 75 cancellous bone samples including 22 bovine samples was 1.42 with a value of 55 m/s for the RMSE of the curve fit. The latter fit was obtained by using of a value of c(m)=3200 m/s. Although the MBA model relies on the empirical parameters determined from experimental data, it is expected that the model can be usefully employed as a practical tool in the field of clinical ultrasonic bone assessment.

The effect of phase cancellation on estimates of calcaneal broadband ultrasound attenuation in vivo

Broadband ultrasonic attenuation (BUA) is a clinically-accepted measurement for prediction of osteoporotic fracture risk. Typical clinical BUA measurements are performed with phase-sensitive receivers and, therefore, can be affected by phase cancellation. In order to separate the effects of conventional attenuation (absorption plus scattering) from phase cancellation, BUA was measured on phantoms with acrylic wedge phase aberrators and on 73 women using both phase sensitive (PS) and phase insensitive (PI) reception. A clinical bone sonometer with a two-dimensional (2-D) receiver array was used. PI BUA measurements on phantoms with acrylic wedge phase aberrators were found to be far more resistant to phase cancellation than PS BUA measurements. In data from 73 women, means and standard deviations for BUA measurements were 81.4 +/- 21.4 dB/MHz (PS) and 67.2 +/- 9.7 dB/MHz (PI). The magnitude of the discrepancy between PS BUA and PI BUA tended to increase with bone mineral density (BMD).

Evaluation of ultrasonic scattering in human cancellous bone by using a binary mixture model

A weak scattering model based on small perturbations in a binary mixture is developed to estimate the ultrasonic scattering from human cancellous bone, which is modelled as a random isotropic continuum containing identical scatters. Ultrasonic scattering is determined by both velocity fluctuation and density fluctuation, when k(2)a(2) << 1 is satisfied. Two kinds of trabeculae thickness distributions, i.e. even distribution and Gauss distribution, are applied in the calculation of attenuation and backscattering. Frequency dependence of the backscatter coefficient is found to be Af(3.13) and Af(2.84) with the Gauss distribution and an even distribution, respectively. Both backscattering and attenuation change significantly against porosity for the case of high porosity. The predicted results are close to the measured ones from the literature. The errors of this theoretical model are also discussed in this paper.

Empirical angle-dependent Biot and MBA models for acoustic anisotropy in cancellous bone

The Biot and the modified Biot-Attenborough (MBA) models have been found useful to understand ultrasonic wave propagation in cancellous bone. However, neither of the models, as previously applied to cancellous bone, allows for the angular dependence of acoustic properties with direction. The present study aims to account for the acoustic anisotropy in cancellous bone, by introducing empirical angle-dependent input parameters, as defined for a highly oriented structure, into the Biot and the MBA models. The anisotropy of the angle-dependent Biot model is attributed to the variation in the elastic moduli of the skeletal frame with respect to the trabecular alignment. The angle-dependent MBA model employs a simple empirical way of using the parametric fit for the fast and the slow wave speeds. The angle-dependent models were used to predict both the fast and slow wave velocities as a function of propagation angle with respect to the trabecular alignment of cancellous bone. The predictions were compared with those of the Schoenberg model for anisotropy in cancellous bone and in vitro experimental measurements from the literature. The angle-dependent models successfully predicted the angular dependence of phase velocity of the fast wave with direction. The root-mean-square errors of the measured versus predicted fast wave velocities were 79.2 m s(-1) (angle-dependent Biot model) and 36.1 m s(-1) (angle-dependent MBA model). They also predicted the fact that the slow wave is nearly independent of propagation angle for angles about 50 degrees , but consistently underestimated the slow wave velocity with the root-mean-square errors of 187.2 m s(-1) (angle-dependent Biot model) and 240.8 m s(-1) (angle-dependent MBA model). The study indicates that the angle-dependent models reasonably replicate the acoustic anisotropy in cancellous bone.

Numerical simulation of the dependence of quantitative ultrasonic parameters on trabecular bone microarchitecture and elastic constants

Finite-difference numerical simulation of ultrasound propagation in complex media such as cancellous bone represents a fertile alternative to analytical approaches because it can manage the complex 3D bone structure by coupling the numerical computation with 3D numerical models of bone microarchitecture obtained from high-resolution imaging modalities. The objective of this work was to assess in silico the sensitivity of ultrasound parameters to controlled changes of microarchitecture and variation of elastic constants. The simulation software uses a finite-difference approach based on the Virieux numerical scheme. An incident plane wave was propagated through a volume of bone of approximately 5 x 5 x 8 mm(3). The volumes were reconstructed from high-resolution micro-computed tomography data. An iterative numerical scenario of "virtual osteoporosis" was implemented using a dedicated image processing algorithm in order to modify the initial 3D microstructures. Numerical computations of wave propagation were performed at each step of the process. The sensitivity to bone material properties was also tested by changing the elastic constants of bone tissue. Our results suggest that ultrasonic variables (slope of the frequency-dependent attenuation coefficient and speed of sound) are mostly influenced by bone volume fraction. However, material properties and structure also appear to play a role. The impact of modifications of the stiffness coefficients remained lower than the variability caused by structural variations. This study emphasizes the potential of numerical computations tools coupled to realistic 3D structures to elucidate the physical mechanisms of interaction between ultrasound and bone structure and to assess the sensitivity of ultrasound variables to different bone properties.

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.

Three-dimensional simulation of ultrasound propagation through trabecular bone structures measured by synchrotron microtomography

Three-dimensional numerical simulations of ultrasound transmission were performed through 31 trabecular bone samples measured by synchrotron microtomography. The synchrotron microtomography provided high resolution 3D mappings of bone structures, which were used as the input geometry in the simulation software developed in our laboratory. While absorption (i.e. the absorption of ultrasound through dissipative mechanisms) was not taken into account in the algorithm, the simulations reproduced major phenomena observed in real through-transmission experiments in trabecular bone. The simulated attenuation (i.e. the decrease of the transmitted ultrasonic energy) varies linearly with frequency in the MHz frequency range. Both the speed of sound (SOS) and the slope of the normalized frequency-dependent attenuation (nBUA) increase with the bone volume fraction. Twenty-five out of the thirty-one samples exhibited negative velocity dispersion. One sample was rotated to align the main orientation of the trabecular structure with the direction of ultrasonic propagation, leading to the observation of a fast and a slow wave. Coupling numerical simulation with real bone architecture therefore provides a powerful tool to investigate the physics of ultrasound propagation in trabecular structures. As an illustration, comparison between results obtained on bone modelled either as a fluid or a solid structure suggested the major role of mode conversion of the incident acoustic wave to shear waves in bone to explain the large contribution of scattering to the overall attenuation.

Recent developments in trabecular bone characterization using ultrasound

Currently available quantitative ultrasound technologies to assess cancellous bone are based on the measurements in transmission of speed of sound or slope of frequency-dependent attenuation (so called broadband ultrasonic attenuation). These two parameters are now considered as surrogate markers of site-matched bone mineral density. The ability of ultrasound techniques to provide non-bone mineral density-related bone properties (eg, microstructure) has not been clearly demonstrated yet. This is mainly because of two factors: a lack of understanding of ultrasound propagation with clear identification of the different underlying physical interactions; and the difficulty of performing experiments because of the limited sample size, the large number of statistical relationships to be tested with multiple variables, and the usual strong covariance observed between bone quantity and microarchitecture. The aim of this paper is to review the most recent development in the field of ultrasound characterization of trabecular bone. We present research work on ultrasound backscatter and how it could be used to estimate microarchitectural properties independently of bone quantity, and the first promising results obtained for the estimation of trabecular thickness. We then introduce numeric simulations of wave propagation through trabecular microarchitecture and show how it could contribute to elucidate and better characterize the physical underlying physics and result in more predictive models. These innovative acquisition schemes and the possibility of virtual experiments should altogether contribute to rapid advancement of ultrasonic bone characterization.

Speed of sound: relation to geometric characteristics of bone in children, adolescents, and adults

OBJECTIVES

To investigate the relation between volumetric bone mineral density (vBMD) and speed of sound (SOS).

STUDY DESIGN

Total and trabecular vBMD were measured by peripheral quantitative computed tomography at the forearm in a population of 216 individuals of a pediatric outpatient clinic. Moreover, SOS was measured by a quantitative ultrasound device (QUS) at the thumb, patella, and os calcis.

RESULTS

Linear regression analysis revealed that the prediction of SOS by vBMD is relatively weak (R2 < 0.1). Moreover, body height and measures of bone size have a stronger influence on SOS than vBMD. The influence of bone size on SOS also depends on the location of measurement (highest prediction of SOS by body height at patella with R2 = 0.56). Anthropometric characteristics have a stronger influence on SOS than measures of bone mineral density at the thumb and patella in comparison to os calcis (body height predicts SOS at os calcis, with R2 = 0.03). Conclusions QUS is not a suitable method to assess bone density. If QUS is applied for the assessment of bone development and of bone fracture risk, the measurement should be performed with consideration of anthropometric measurements.

Comparison of measurements of phase velocity in human calcaneus to Biot theory

Biot's theory for elastic propagation in porous media has previously been shown to be useful for modeling the dependence of phase velocity on porosity in bovine cancellous bone in vitro. In the present study, Biot's theory is applied to measurements of porosity-dependent phase velocity in 53 human calcanea in vitro. Porosity was measured using microcomputed tomography for some samples (n = 23) and estimated based on bone mineral densitometry for the remaining samples (n = 30). The phase velocity at 500 kHz was measured in a water tank using a through-transmission technique. Biot's theory performed well for the prediction of the dependence of sound speed on porosity. The trend was quasilinear, but both the theory and experiment show similar slight curvature. The root mean square error (RMSE) of predicted versus measured sound speed was 15.8 m/s.

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.

In vitro ultrasound measurement at the human femur

Measurements of broadband ultrasonic attenuation (BUA) represent an established means of skeletal status assessment in osteoporosis. Today, the skeletal site that is the most widely used clinically is the calcaneus (heel bone). However, we hypothesized that fracture risk predictions could be improved if ultrasound measurements were performed directly at the main site of fracture, e.g., the proximal femur. The goal of this paper is to demonstrate in vitro the feasibility of quantitative ultra sound (QUS) imaging at the upper part of fee femur, and to investigate the relationships of BUA to bone mineral density (BMD). Forty-four excise human femurs were measured in transmission with a pair of focused 0.5-MHz central frequency transducers. Two-dimensional scans were performed, and the radio frequency (rf) signals were recorded at each measurement point. A data acceptance criterion for region of interest (ROI) selection was established based on the linearity of the frequency-dependent attenuation. Five measurements with repositioning were performed on each sample to determine the reproducibility. Dual energy x-ray absorptiometry (DXA) scans have been performed on the samples for BMD measurements: Three ROIs were selected in the specimens: greater, trochanter intertrochanteric, and femoral neck regions. Coefficient of variations were in the range 1.6% to 2.5%. The determination coefficients (r(2)) between BUA and BMD in site-matched ROIs were 0.81, 0.78, and 0.73, respectively, for the greater trochanteric, intertrochanteric, and femoral neck regions. Our results are consistent with data previously shown at the calcaneus and demonstrate the feasibility of QUS measurements at the femur in vitro, with reasonable reproducibility.

Frequency and resolution dependence of the anisotropic impedance estimation in cortical bone using time-resolved scanning acoustic microscopy

The influences of frequency and spatial resolution on the anisotropic impedance estimation of cortical bone was investigated in the frequency range 25-100 MHz. A set of spherically focused transducers provided a spatial resolution in the range from 150 down to about 20 mum. Four embedded cortical bone samples (two male, two female, two donors aged <30 years, two donors aged >70 years) were cut with different orientations relative to the long axis of the femur (0-90 degrees ). From each section, impedance maps were acquired in the C-scan mode. Histogram evaluations showed a similar angular dependence with a characteristic off-axis maximum of the estimated impedance for all samples and frequencies. The impedance values obtained with the 25-MHz transducer were significantly lower than those obtained with the 50- and 100-MHz transducers. Morphological parameters of the macrostructure, for example, size and distribution of the haversian channels and the resulting porosity, were estimated from the high-resolution acoustic images. These structures appeared to have a significant influence on the measured properties of the bone matrix for the low-frequency and low-aperture measurements.

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.

Measurement of dependence of backscatter coefficient from cylinders on frequency and diameter using focused transducers--with applications in trabecular bone

A theory for the elastic scattering response from a cylinder insonified by a plane wave was previously derived by Faran. In the present paper, the empirical relationship between Faran's theory and measurements of backscatter coefficient from cylindrical targets using focused transducers is investigated. Experimental measurements of dependence of backscatter coefficient on frequency and diameter for nylon wires are reported. It is found that, under certain conditions (including weak, incoherent scattering), backscatter coefficient measurements from collections of cylindrical scatterers may be meaningfully compared with Faran's model predictions. At low frequencies, the theory and experimental measurements exhibit similar dependences on frequency and diameter, provided that the scatterers are not too densely packed. At higher frequencies, the fine structure of Faran's predictions becomes difficult to reproduce experimentally with a focused transducer. Implications regarding applications to characterization of trabecular bone are discussed.

Acoustic wave propagation in bovine cancellous bone: application of the Modified Biot-Attenborough model

Acoustic wave propagation in bovine cancellous bone is experimentally and theoretically investigated in the frequency range of 0.5-1 MHz. The phase velocity, attenuation coefficient, and broadband ultrasonic attenuation (BUA) of bovine cancellous bone are measured as functions of frequency and porosity. For theoretical estimation, the Modified Biot-Attenborough (MBA) model is employed with three new phenomenological parameters: the boundary condition, phase velocity, and impedance parameters. The MBA model is based on the idealization of cancellous bone as a nonrigid porous medium with circular cylindrical pores oriented normal to the surface. It is experimentally observed that the phase velocity is approximately nondispersive and the attenuation coefficient linearly increases with frequency. The MBA model predicts a slightly negative dispersion of phase velocity linearly with frequency and the nonlinear relationships of attenuation and BUA with porosity. The experimental results are in good agreement with the theoretical results estimated with the MBA model. It is expected that the MBA model can be usefully employed in the field of clinical bone assessment for the diagnosis of osteoporosis.

The dependence of ultrasonic backscatter on trabecular thickness in human calcaneus: theoretical and experimental results

Trabecular thickness within cancellous bone is an important determinant of osteoporotic fracture risk. Noninvasive assessment of trabecular thickness potentially could yield useful diagnostic information. Faran's theory of elastic scattering from a cylindrical object immersed in a fluid has been used to predict the dependence of ultrasonic backscatter on trabecular thickness. The theory predicts that, in the range of morphological and material properties expected for trabecular bone, the backscatter coefficient at 500 kHz should be approximately proportional to trabecular thickness to the power of 2.9. Experimental measurements of backscatter coefficient were performed on 43 human calcaneus samples in vitro. Mean trabecular thicknesses on the 43 samples were assessed using micro computed tomography (CT). A power law fit to the data showed that the backscatter coefficient empirically varied as trabecular thickness to the 2.8 power. The 95% confidence interval for this exponent was 1.7 to 3.9. The square of the correlation coefficient for the linear regression to the log transformed data was 0.40. This suggests that 40% of variations in backscatter may be attributed to variations in trabecular thickness. These results reinforce previous studies that offered validation for the Faran cylinder model for prediction of scattering properties of cancellous bone, and provide added evidence for the potential diagnostic utility of the backscatter measurement.

Prediction of frequency-dependent ultrasonic backscatter in cancellous bone using statistical weak scattering model

The goal of this study was to propose a model for the ultrasonic frequency-dependent backscatter coefficient in cancellous bone. This model allows us to address the inverse problem and to predict the mean trabecular thickness. A weak scattering model is used and the backscatter coefficient is expressed in terms of an autocorrelation function of the medium. Different autocorrelation functions (Gaussian, exponential and densely populated media) were used to compute the backscatter coefficient and comparison is made with experimental data for 19 specimens and for frequency ranging from 0.4 to 1.2 MHz. For each specimen, a nonlinear regression was performed and the mean trabecular thickness is estimated. Experimental data and theoretical predictions were averaged over the 19 specimens. A good agreement between experimental data and predictions was found for both the magnitude and the frequency-dependence of the backscatter coefficient. We also found a good agreement between the experimental mean trabecular thickness (Tb. Th = 130 +/- 6.5 micro m) derived from the analysis of bone 3-D microarchitecture using high-resolution microtomography and theoretical predictions (d(Gauss) = 140 +/- 10 micro m, d(exponential) = 153 +/- 12.5 micro m and d(dense) = 138 +/- 6.5 micro m). These results open interesting prospects for the estimation of the mean trabecular thickness from in vivo measurements.

Prediction of backscatter coefficient in trabecular bones using a numerical model of three-dimensional microstructure

A model of ultrasonic backscattering for cancellous bone saturated by water is proposed. This model assumes that scattering is caused by the solid trabeculae and describes the cancellous bone as a weak scattering medium. The backscatter coefficient is related to the spatial Fourier transform of bone microarchitecture and to the density and compressibility fluctuations between the solid trabeculae and the saturating fluid. The computations of the model make use of three-dimensional numerical images of bone microarchitecture, obtained by tomographic reconstructions with a 10 microm spatial resolution. With this model, the predictions of the frequency dependence and of the magnitude of the backscatter coefficient are reasonably accurate. The theoretical predictions are compared to experimental data obtained on 19 specimens. An accuracy error of approximately 1 dB was found (difference between the averaged experimental values and theoretical predictions). One limit of the model may come from inaccurate values of trabecular bone characteristics needed for the computations (density and longitudinal velocity), which are yet to be precisely determined for human trabecular bone. However, the model is only slightly sensitive to variations of bone material properties. It was found that an accuracy error of 2.2 dB at maximum resulted from inaccurate a priori values of bone material properties. A computation of the elastic mean free path in the medium suggests that multiple scattering plays a minor role in the working frequency bandwidth (0.4-1.2 MHz). It follows from these results that a weak scattering medium model may be appropriate to describe scattering from trabecular bone.

A Gaussian framework for modeling effects of frequency-dependent attenuation, frequency-dependent scattering, and gating

For a wide range of applications in medical ultrasound, power spectra of received signals are approximately Gaussian. It has been established previously that an ultrasound beam with a Gaussian spectrum propagating through a medium with linear attenuation remains Gaussian. In this paper, Gaussian transformations are derived to model the effects of scattering (according to a power law, as is commonly applicable in soft tissues, especially over limited frequency ranges) and gating (with a Hamming window, a commonly used gate function). These approximations are shown to be quite accurate even for relatively broad band systems with fractional bandwidths approaching 100%. The theory is validated by experiments in phantoms consisting of glass particles suspended in agar.