Cancellous bone fast and slow waves obtained with Bayesian probability theory correlate with porosity from computed tomography

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.

Bayesian inference of a human bone and biomaterials using ultrasonic transmitted signals

Ultrasonic techniques could be good candidates to aid the assessment of osteoporosis detection, due to their non-intrusiveness and low cost. While earlier studies made use of the measured ultrasonic phase velocity and attenuation inside the bone, very few have considered an inverse identification of both the intrinsic pore microstructure and the mechanical properties of the bone, based on Biot's model. The main purpose of this work is to present an in vitro methodology for bone identification, adopting a statistical Bayesian inference technique using ultrasonic transmitted signals, which allows the retrieval of the identified parameters and their uncertainty. In addition to the bone density, Young's modulus and Poisson's ratio, the bone pore microstructure parameters (porosity, tortuosity, and viscous length) are identified. These additional microstructural terms could improve the knowledge on the correlations between bone microstructure and bone diseases, since they provide more information on the trabecular structure. In general, the exact properties of the saturating fluid are unknown (bone marrow and blood in the case of bone study) so in this work, the fluid properties (water) are identified during the inference as a proof of concept.

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.

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

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

Fast decomposition of two ultrasound longitudinal waves in cancellous bone using a phase rotation parameter for bone quality assessment: Simulation study

Ultrasound signals that pass through cancellous bone may be considered to consist of two longitudinal waves, which are called fast and slow waves. Accurate decomposition of these fast and slow waves is considered to be highly beneficial in determination of the characteristics of cancellous bone. In the present study, a fast decomposition method using a wave transfer function with a phase rotation parameter was applied to received signals that have passed through bovine bone specimens with various bone volume to total volume (BV/TV) ratios in a simulation study, where the elastic finite-difference time-domain method is used and the ultrasound wave propagated parallel to the bone axes. The proposed method succeeded to decompose both fast and slow waves accurately; the normalized residual intensity was less than -19.5 dB when the specimen thickness ranged from 4 to 7 mm and the BV/TV value ranged from 0.144 to 0.226. There was a strong relationship between the phase rotation value and the BV/TV value. The ratio of the peak envelope amplitude of the decomposed fast wave to that of the slow wave increased monotonically with increasing BV/TV ratio, indicating the high performance of the proposed method in estimation of the BV/TV value in cancellous 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.

Experimental observation of ultrasound fast and slow waves through three-dimensional printed trabecular bone phantoms

In this paper, ultrasound measurements of 1:1 scale three-dimensional (3D) printed trabecular bone phantoms are reported. The micro-structure of a trabecular horse bone sample was obtained via synchrotron x-ray microtomography, converted to a 3D binary data set, and successfully 3D-printed at scale 1:1. Ultrasound through-transmission experiments were also performed through a highly anisotropic version of this structure, obtained by elongating the digitized structure prior to 3D printing. As in real anisotropic trabecular bone, both the fast and slow waves were observed. This illustrates the potential of stereolithography and the relevance of such bone phantoms for the study of ultrasound propagation in bone.

New perspectives in echographic diagnosis of osteoporosis on hip and spine

Currently, the accepted "gold standard" method for bone mineral density (BMD) measurement and osteoporosis diagnosis is dual-energy X-ray absorptiometry (DXA). However, actual DXA effectiveness is limited by several factors, including intrinsic accuracy uncertainties and possible errors in patient positioning and/or post-acquisition data analysis. DXA employment is also restricted by the typical issues related to ionizing radiation employment (high costs, need of dedicated structures and certified operators, unsuitability for population screenings). The only commercially-available alternative to DXA is represented by "quantitative ultrasound" (QUS) approaches, which are radiation-free, cheaper and portable, but they cannot be applied on the reference anatomical sites (lumbar spine and proximal femur). Therefore, their documented clinical usefulness is restricted to calcaneal applications on elderly patients (aged over 65 y), in combination with clinical risk factors and only for the identification of healthy subjects at low fracture risk. Literature-reported studies performed some QUS measurements on proximal femur, but their clinical translation is mostly hindered by intrinsic factors (e.g., device bulkiness). An innovative ultrasound methodology has been recently introduced, which performs a combined analysis of B-mode images and corresponding "raw" radiofrequency signals acquired during an echographic scan of the target reference anatomical site, providing two novel parameters: Osteoporosis Score and Fragility Score, indicative of BMD level and bone strength, respectively. This article will provide a brief review of the available systems for osteoporosis diagnosis in clinical routine contexts, followed by a synthesis of the most promising research results on the latest ultrasound developments for early osteoporosis diagnosis and fracture prevention.

Conventional, Bayesian, and Modified Prony's methods for characterizing fast and slow waves in equine cancellous bone

Conventional, Bayesian, and the modified least-squares Prony's plus curve-fitting (MLSP + CF) methods were applied to data acquired using 1 MHz center frequency, broadband transducers on a single equine cancellous bone specimen that was systematically shortened from 11.8 mm down to 0.5 mm for a total of 24 sample thicknesses. Due to overlapping fast and slow waves, conventional analysis methods were restricted to data from sample thicknesses ranging from 11.8 mm to 6.0 mm. In contrast, Bayesian and MLSP + CF methods successfully separated fast and slow waves and provided reliable estimates of the ultrasonic properties of fast and slow waves for sample thicknesses ranging from 11.8 mm down to 3.5 mm. Comparisons of the three methods were carried out for phase velocity at the center frequency and the slope of the attenuation coefficient for the fast and slow waves. Good agreement among the three methods was also observed for average signal loss at the center frequency. The Bayesian and MLSP + CF approaches were able to separate the fast and slow waves and provide good estimates of the fast and slow wave properties even when the two wave modes overlapped in both time and frequency domains making conventional analysis methods unreliable.

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.

Fast characterization of two ultrasound longitudinal waves in cancellous bone using an adaptive beamforming technique

The received signal in through-transmission ultrasound measurements of cancellous bone consists of two longitudinal waves, called the fast and slow waves. Analysis of these fast and slow waves may reveal characteristics of the cancellous bone that would be good indicators of osteoporosis. Because the two waves often overlap, decomposition of the received signal is an important problem in the characterization of bone quality. This study proposes a fast and accurate decomposition method based on the frequency domain interferometry imaging method with a modified wave transfer function that uses a phase rotation parameter. The proposed method accurately characterized the fast and slow waves in the experimental study, and the residual intensity, which was normalized with respect to the received signal intensity, was less than -20 dB over the bone specimen thickness range from 6 to 15 mm. In the simulation study, the residual intensity was less than -20 dB over the specimen thickness range from 3 to 8 mm. Decomposition of a single received signal takes only 5 s using a laptop personal computer with a single central processing unit. The proposed method has great potential to provide accurate and rapid measurements of indicators of osteoporosis in cancellous bone.

Fast and slow wave detection in bovine cancellous bone in vitro using bandlimited deconvolution and Prony's method

Fast and slow waves were detected in a bovine cancellous bone sample for thicknesses ranging from 7 to 12 mm using bandlimited deconvolution and the modified least-squares Prony's method with curve fitting (MLSP + CF). Bandlimited deconvolution consistently isolated two waves with linear-with-frequency attenuation coefficients as evidenced by high correlation coefficients between attenuation coefficient and frequency: 0.997 ± 0.002 (fast wave) and 0.986 ± 0.013 (slow wave) (mean ± standard deviation). Average root-mean-squared (RMS) differences between the two algorithms for phase velocities were 5 m/s (fast wave, 350 kHz) and 13 m/s (slow wave, 750 kHz). Average RMS differences for signal loss were 1.6 dB (fast wave, 350 kHz) and 0.4 dB (slow wave, 750 kHz). Phase velocities for thickness = 10 mm were 1726 m/s (fast wave, 350 kHz) and 1455 m/s (slow wave, 750 kHz). Results show support for the model of two waves with linear-with frequency attenuation, successful isolation of fast and slow waves, good agreement between bandlimited deconvolution and MLSP + CF as well as with a Bayesian algorithm, and potential variations of fast and/or slow wave properties with bone sample thickness.

Ultrasound to assess bone quality

Bone quality is determined by a variety of compositional, micro- and ultrastructural properties of the mineralized tissue matrix. In contrast to X-ray-based methods, the interaction of acoustic waves with bone tissue carries information about elastic and structural properties of the tissue. Quantitative ultrasound (QUS) methods represent powerful alternatives to ionizing x-ray based assessment of fracture risk. New in vivo applicable methods permit measurements of fracture-relevant properties, [eg, cortical thickness and stiffness at fragile anatomic regions (eg, the distal radius and the proximal femur)]. Experimentally, resonance ultrasound spectroscopy and acoustic microscopy can be used to assess the mesoscale stiffness tensor and elastic maps of the tissue matrix at microscale resolution, respectively. QUS methods, thus, currently represent the most promising approach for noninvasive assessment of components of fragility beyond bone mass and bone microstructure providing prospects for improved assessment of fracture risk.

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.

An experimental study on the ultrasonic wave propagation in cancellous bone: waveform changes during propagation

Wave propagation in a trabecular bone was experimentally investigated using an acoustic tube. For the purposes of this study, a cubic sample was gradually filed so the waveform change due to the sample thickness could be observed. The initial sample showed clear two-wave separation. As the sample became thinner, the fast and slow waves gradually overlapped. The apparent frequencies and amplitudes of the fast waves obtained from the time domain data decreased significantly for the smaller thicknesses. This indicates an increase in the apparent attenuation at the initial stage of the propagation. Next the authors investigated the distribution of the ultrasonic field after the transmission through the cancellous bone sample. In addition to a large aperture receiver, a needle-type ultrasonic transducer was used to observe the ultrasonic field. Within an area of the same size of the large transducer, the waveforms retrieved with the needle sensor exhibited high spatial variations; however, the averaged waveform in the plane was similar to the waveform obtained with the large aperture receiver. This indicates that the phase cancellation effect on the surface of the large aperture receiver can be one of the reasons for the strong apparent attenuation observed at the initial stages of the propagation.

Estimation of multipath transmission parameters for quantitative ultrasound measurements of bone

When applying quantitative ultrasound (QUS) measurements to bone for predicting osteoporotic fracture risk, the multipath transmission of sound waves frequently occurs. In the last 10 years, the interest in separating multipath QUS signals for their analysis awoke, and led to the introduction of several approaches. Here, we compare the performances of the two fastest algorithms proposed for QUS measurements of bone: the modified least-squares Prony method (MLSP), and the space alternating generalized expectation maximization algorithm (SAGE) applied in the frequency domain. In both approaches, the parameters of the transfer functions of the sound propagation paths are estimated. To provide an objective measure, we also analytically derive the Cramér-Rao lower bound of variances for any estimator and arbitrary transmit signals. In comparison with results of Monte Carlo simulations, this measure is used to evaluate both approaches regarding their accuracy and precision. Additionally, with simulations using typical QUS measurement settings, we illustrate the limitations of separating two superimposed waves for varying parameters with focus on their temporal separation. It is shown that for good SNRs around 100 dB, MLSP yields better results when two waves are very close. Additionally, the parameters of the smaller wave are more reliably estimated. If the SNR decreases, the parameter estimation with MLSP becomes biased and inefficient. Then, the robustness to noise of the SAGE clearly prevails. Because a clear influence of the interrelation between the wavelength of the ultrasound signals and their temporal separation is observable on the results, these findings can be transferred to QUS measurements at other sites. The choice of the suitable algorithm thus depends on the measurement conditions.

Estimation of fast and slow wave properties in cancellous bone using Prony's method and curve fitting

The presence of two longitudinal waves in poroelastic media is predicted by Biot's theory and has been confirmed experimentally in through-transmission measurements in cancellous bone. Estimation of attenuation coefficients and velocities of the two waves is challenging when the two waves overlap in time. The modified least squares Prony's (MLSP) method in conjuction with curve-fitting (MLSP + CF) is tested using simulations based on published values for fast and slow wave attenuation coefficients and velocities in cancellous bone from several studies in bovine femur, human femur, and human calcaneus. The search algorithm is accelerated by exploiting correlations among search parameters. The performance of the algorithm is evaluated as a function of signal-to-noise ratio (SNR). For a typical experimental SNR (40 dB), the root-mean-square errors (RMSEs) for one example (human femur) with fast and slow waves separated by approximately half of a pulse duration were 1 m/s (slow wave velocity), 4 m/s (fast wave velocity), 0.4 dB/cm MHz (slow wave attenuation slope), and 1.7 dB/cm MHz (fast wave attenuation slope). The MLSP + CF method is fast (requiring less than 2 s at SNR = 40 dB on a consumer-grade notebook computer) and is flexible with respect to the functional form of the parametric model for the transmission coefficient. The MLSP + CF method provides sufficient accuracy and precision for many applications such that experimental error is a greater limiting factor than estimation error.

Simulations of ultrasound propagation in random arrangements of elliptic scatterers: occurrence of two longitudinal waves

Ultrasound propagation in clusters of elliptic (two-dimensional) or ellipsoidal (three-dimensional) scatterers randomly distributed in a fluid is investigated numerically. The essential motivation for the present work is to gain a better understanding of ultrasound propagation in trabecular bone. Bone microstructure exhibits structural anisotropy and multiple wave scattering. Some phenomena remain partially unexplained, such as the propagation of two longitudinal waves. The objective of this study was to shed more light on the occurrence of these two waves, using finite-difference simulations on a model medium simpler than bone. Slabs of anisotropic, scattering media were randomly generated. The coherent wave was obtained through spatial and ensemble-averaging of the transmitted wavefields. When varying relevant medium parameters, four of them appeared to play a significant role for the observation of two waves: (i) the solid fraction, (ii) the direction of propagation relatively to the scatterers orientation, (iii) the ability of scatterers to support shear waves, and (iv) a continuity of the solid matrix along the propagation. These observations are consistent with the hypothesis that fast waves are guided by the locally plate/bar-like solid matrix. If confirmed, this interpretation could significantly help developing approaches for a better understanding of trabecular bone micro-architecture using ultrasound.