Propagation of two longitudinal waves in a cancellous bone with the closed pore boundary

Two-dimensional shape discrimination by sighted people using simulated virtual echoes

In this study, a new research method using psychoacoustic experiments and acoustic simulations is proposed for human echolocation research. A shape discrimination experiment was conducted for sighted people using pitch-converted virtual echoes from targets of dissimilar two-dimensional (2D) shapes. These echoes were simulated using a three-dimensional acoustic simulation based on a finite-difference time-domain method from Bossy, Talmat, and Laugier [(2004). J. Acoust. Soc. Am. 115, 2314-2324]. The experimental and simulation results suggest that the echo timbre and pitch determined based on the sound interference may be effective acoustic cues for 2D shape discrimination. The newly developed research method may lead to more efficient future studies of human echolocation.

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.

Can low-frequency guided waves at the tibia paired with machine learning differentiate between healthy and osteopenic/osteoporotic subjects? A pilot study

PURPOSE

Axial transmission quantitative acoustics (ax-QA) has shown to be a promising tool for assessing bone health and properties in a safe, inexpensive, and portable manner. This study investigated the efficacy of low-frequency ax-QA measured at the tibia, paired with a support vector machine (SVM) approach for combining multiple acoustic indicators, to diagnose osteoporosis as defined by bone mineral density.

METHODS

This pilot study measured 41 female subjects using ax-QA (flexural mode, 3 kHz) at the tibia and using dual X-ray absorptiometry (DXA) at the lumbar spine, femoral neck, and distal radius. For each location, a threshold classifier and SVM were trained to differentiate between healthy and non-healthy subjects based on the phase velocity at different frequencies. Receiver Operating Characteristics and area under curve values (AUC) were used to assess the classifiers' performances for various thresholds and class-weights.

RESULTS

The SVM outperformed the threshold classifier for all three bone locations at low false positive rates. While differentiation between healthy and non-healthy bone states was poor for the spine (AUC: 0.56 ± 0.04), good to moderate performances were observed for the radius (AUC: 0.83 ± 0.03) and hip (AUC: 0.71 ± 0.04).

CONCLUSIONS

Low-frequency ax-QA has demonstrated potential for complementing DXA in screening for osteoporosis at the radius and hip. Through further addition of acoustic indicators ax-QA could provide a diagnostic alternative in third-world countries, and bring bone health screening and monitoring into the hands of clinicians and general health practitioners everywhere.

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

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

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.

Two-wave behavior under various conditions of transition area from cancellous bone to cortical bone

The two-wave phenomenon, the wave separation of a single ultrasonic pulse in cancellous bone, is expected to be a useful tool for the diagnosis of osteoporosis. However, because actual bone has a complicated structure, precise studies on the effect of transition conditions between cortical and cancellous parts are required. This study investigated how the transition condition influenced the two-wave generation using three-dimensional X-ray CT images of an equine radius and a three-dimensional simulation technique. As a result, any changes in the boundary between cortical part and trabecular part, which gives the actual complex structure of bone, did not eliminate the generation of either the primary wave or the secondary wave at least in the condition of clear trabecular alignment. The results led us to the possibility of using the two-wave phenomenon in a diagnostic system for osteoporosis in cases of a complex boundary.

Measurements of ultrasound velocity and attenuation in numerical anisotropic porous media compared to Biot's and multiple scattering models

This article quantitatively investigates ultrasound propagation in numerical anisotropic porous media with finite-difference simulations in 3D. The propagation media consist of clusters of ellipsoidal scatterers randomly distributed in water, mimicking the anisotropic structure of cancellous bone. Velocities and attenuation coefficients of the ensemble-averaged transmitted wave (also known as the coherent wave) are measured in various configurations. As in real cancellous bone, one or two longitudinal modes emerge, depending on the micro-structure. The results are confronted with two standard theoretical approaches: Biot's theory, usually invoked in porous media, and the Independent Scattering Approximation (ISA), a classical first-order approach of multiple scattering theory. On the one hand, when only one longitudinal wave is observed, it is found that at porosities higher than 90% the ISA successfully predicts the attenuation coefficient (unlike Biot's theory), as well as the existence of negative dispersion. On the other hand, the ISA is not well suited to study two-wave propagation, unlike Biot's model, at least as far as wave speeds are concerned. No free fitting parameters were used for the application of Biot's theory. Finally we investigate the phase-shift between waves in the fluid and the solid structure, and compare them to Biot's predictions of in-phase and out-of-phase motions.

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.

Multichannel instantaneous frequency analysis of ultrasound propagating in cancellous bone

An ultrasonic pulse propagating in cancellous bone can be separated into two waves depending on the condition of the specimen. These two waves, which are called the fast wave and the slow wave, provide important information for the diagnosis of osteoporosis. The present study proposes to utilize a signal processing method that extracts the instantaneous frequency (IF) of waveforms from multiple spectral channels. The instantaneous frequency was expected to be able to show detailed time-frequency properties of ultrasonic waves being transmitted through cancellous bone. The employed method, termed the multichannel instantaneous frequency (MCIF) method, showed robustness against background noise as compared to the IF that was directly derived from the original waveform. The extracted IF revealed that the frequency of the fast wave was affected by both the propagation distance within the specimen and the bone density, independently. On the other hand, the alternation of the center frequency of the originally transmitted wave did not produce proportional changes in the extracted IF values of the fast waves, suggesting that the fast wave IF mainly reflected the thickness of the specimens. These findings may provide the possibility of obtaining a more precise diagnosis of osteoporosis.

Numerical investigation of reflection properties of fast and slow longitudinal waves in cancellous bone

The basic reflection properties of fast and slow waves propagating in cancellous bone were numerically investigated using finite-difference time-domain (FDTD) simulations with 18 microcomputed tomographic (µCT) models of bovine cancellous bone. The simulated results showed that both reflection coefficients of the fast and slow waves linearly increased with porosity.

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.

Slow and fast ultrasonic wave detection improvement in human trabecular bones using Golay code modulation

The identification of fast and slow waves propagating through trabecular bone is a challenging task due to temporal wave overlap combined with the high attenuation of the fast wave in the presence of noise. However, it can provide valuable information about bone integrity and become a means for monitoring osteoporosis. The objective of this work is to apply different coded excitation methods for this purpose. The results for single-sine cycle pulse, Golay code, and chirp excitations are compared. It is shown that Golay code is superior to the other techniques due to its signal enhancement while exhibiting excellent resolution without the ambiguity of sidelobes.

Two-wave propagation imaging to evaluate the structure of cancellous bone

The two-wave phenomenon reflects not only bone mass but also the complex bone structure of cancellous bone. We propose a new simple imaging technique based on the two-wave phenomenon for investigating the anisotropic structure of cancellous bone. A cylindrical specimen of cancellous bone was obtained from a bovine femur. The structure (alignment of trabeculae) of the specimen was obtained from 3-D X-ray micro computed tomography imaging. Using a conventional ultrasonic pulse technique, we rotated the receiver around the specimen to investigate the ultrasonic fields after propagation within the specimen. The ultrasonic propagation image clearly showed the effect of the bone structure. We found that the fast wave showed apparent refraction, whereas the slow wave did not. Fast-wave propagation imaging is thus a simple and convenient technique for easy interpretation of the anisotropic structure. This imaging technique has the potential to become a powerful tool to investigate the structure of trabeculae during in vivo measurements.