Inverse problems in cancellous bone: estimation of the ultrasonic properties of fast and slow waves using Bayesian probability theory

What kind of waves are measured in trabecular bone?

The paper discusses the fundamental mechanisms underlying the interaction between ultrasound and trabecular bone, which is considered a two-phase material. When fluid-saturated cancellous bone is interrogated by ultrasound, in some cases, one or two wave modes are observed. Many authors claim that these waves correspond to the fast and slow waves predicted by Biot's theory of elastic wave propagation in fluid-saturated porous media. Within our analysis of the physical conditions, predictions of the existing two-phase models of the propagation of ultrasonic waves in the material as well as numerical simulations for fluid-saturated trabecular bone were performed. On the basis of the theoretical results (from numerical studies) and arguments presented in this paper, we aimed to answer the question of whether two waves observed in ultrasonic wave transmission studies can be interpreted as the fast and slow waves predicted by Biot's theory.

Mechanisms of Interaction of Ultrasound With Cancellous Bone: A Review

Ultrasound is now a clinically accepted modality in the management of osteoporosis. The most common commercial clinical devices assess fracture risk from measurements of attenuation and sound speed in cancellous bone. This review discusses fundamental mechanisms underlying the interaction between ultrasound and cancellous bone. Because of its two-phase structure (mineralized trabecular network embedded in soft tissue-marrow), its anisotropy, and its inhomogeneity, cancellous bone is more difficult to characterize than most soft tissues. Experimental data for the dependencies of attenuation, sound speed, dispersion, and scattering on ultrasound frequency, bone mineral density, composition, microstructure, and mechanical properties are presented. The relative roles of absorption, scattering, and phase cancellation in determining attenuation measurements in vitro and in vivo are delineated. Common speed of sound metrics, which entail measurements of transit times of pulse leading edges (to avoid multipath interference), are greatly influenced by attenuation, dispersion, and system properties, including center frequency and bandwidth. However, a theoretical model has been shown to be effective for correction for these confounding factors in vitro and in vivo. Theoretical and phantom models are presented to elucidate why cancellous bone exhibits negative dispersion, unlike soft tissue, which exhibits positive dispersion. Signal processing methods are presented for separating "fast" and "slow" waves (predicted by poroelasticity theory and supported in cancellous bone) even when the two waves overlap in time and frequency domains. Models to explain dependencies of scattering on frequency and mean trabecular thickness are presented and compared with measurements. Anisotropy, the effect of the fluid filler medium (marrow in vivo or water in vitro), phantoms, computational modeling of ultrasound propagation, acoustic microscopy, and nonlinear properties in cancellous bone are also discussed.

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

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

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.

Predicting Structural Properties of Cortical Bone by Combining Ultrasonic Attenuation and an Artificial Neural Network (ANN): 2-D FDTD Study

The goal of this paper is to predict the micro-architectural parameters of cortical bone such as pore diameter ( ϕ ) and porosity ( v ) from ultrasound attenuation measurements using an artificial neural network (ANN). Slices from a 3-D CT scan of human femur are obtained. The micro-architectural parameters of porosity such as average pore size and porosity are calculated using image processing. When ultrasound waves propagate in porous structures, attenuation is observed due to scattering. Two-dimensional finite-difference time-domain simulations are carried out to obtain frequency dependent attenuation in those 2D structures. An artificial neural network (ANN) is then trained with the input feature vector as the frequency dependent attenuation and output as pore diameter ( ϕ ) and porosity ( v ) . The ANN is composed of one input layer, 3 hidden layers and one output layer, all of which are fully connected. 340 attenuation data sets were acquired and trained over 2000 epochs with a batch size of 32. Data was split into train, validation and test. It was observed that the ANN predicted the micro-architectural parameters of the cortical bone with high accuracies and low losses with a minimum R2 (goodness of fit) value of 0.95. ANN approaches could potentially help inform the solution of inverse-problems to retrieve bone porosity from ultrasound measurements. Ultimately, those inverse-problems could be used for the non-invasive diagnosis and monitoring of osteoporosis.

Low Variance Estimation of Backscatter Quantitative Ultrasound Parameters Using Dynamic Programming

One of the main limitations of ultrasound imaging is that image quality and interpretation depend on the skill of the user and the experience of the clinician. Quantitative ultrasound (QUS) methods provide objective, system-independent estimates of tissue properties, such as acoustic attenuation and backscattering properties of tissue, which are valuable as objective tools for both diagnosis and intervention. Accurate and precise estimation of these properties requires correct compensation for intervening tissue attenuation. Prior attempts to estimate intervening-tissue attenuation based on minimizing cost functions that compared backscattered echo data to models have resulted in limited precision and accuracy. To overcome these limitations, in this paper, we incorporate the prior information of piecewise continuity of QUS parameters as a regularization term into our cost function. We further propose to calculate this cost function using dynamic programming (DP), a computationally efficient optimization algorithm that finds the global optimum. Our results on tissue-mimicking phantoms show that DP substantially outperforms a published least squares method in terms of both estimation bias and variance.

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.

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.

Attenuated Fractional Wave Equations With Anisotropy

This paper develops new fractional calculus models for wave propagation. These models permit a different attenuation index in each coordinate to fully capture the anisotropic nature of wave propagation in complex media. Analytical expressions that describe power law attenuation and anomalous dispersion in each direction are derived for these fractional calculus models.

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.

A deconvolution method for deriving the transit time spectrum for ultrasound propagation through cancellous bone replica models

The acceptance of broadband ultrasound attenuation for the assessment of osteoporosis suffers from a limited understanding of ultrasound wave propagation through cancellous bone. It has recently been proposed that the ultrasound wave propagation can be described by a concept of parallel sonic rays. This concept approximates the detected transmission signal to be the superposition of all sonic rays that travel directly from transmitting to receiving transducer. The transit time of each ray is defined by the proportion of bone and marrow propagated. An ultrasound transit time spectrum describes the proportion of sonic rays having a particular transit time, effectively describing lateral inhomogeneity of transit times over the surface of the receiving ultrasound transducer. The aim of this study was to provide a proof of concept that a transit time spectrum may be derived from digital deconvolution of input and output ultrasound signals. We have applied the active-set method deconvolution algorithm to determine the ultrasound transit time spectra in the three orthogonal directions of four cancellous bone replica samples and have compared experimental data with the prediction from the computer simulation. The agreement between experimental and predicted ultrasound transit time spectrum analyses derived from Bland–Altman analysis ranged from 92% to 99%, thereby supporting the concept of parallel sonic rays for ultrasound propagation in cancellous bone. In addition to further validation of the parallel sonic ray concept, this technique offers the opportunity to consider quantitative characterisation of the material and structural properties of cancellous bone, not previously available utilising ultrasound.

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.

Nested sampling applied in Bayesian room-acoustics decay analysis

Room-acoustic energy decays often exhibit single-rate or multiple-rate characteristics in a wide variety of rooms/halls. Both the energy decay order and decay parameter estimation are of practical significance in architectural acoustics applications, representing two different levels of Bayesian probabilistic inference. This paper discusses a model-based sound energy decay analysis within a Bayesian framework utilizing the nested sampling algorithm. The nested sampling algorithm is specifically developed to evaluate the Bayesian evidence required for determining the energy decay order with decay parameter estimates as a secondary result. Taking the energy decay analysis in architectural acoustics as an example, this paper demonstrates that two different levels of inference, decay model-selection and decay parameter estimation, can be cohesively accomplished by the nested sampling algorithm.

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

A Bayesian probability theory approach for separating overlapping ultrasonic fast and slow waves in cancellous bone has been previously introduced. The goals of this study were to investigate whether the fast and slow waves obtained from Bayesian separation of an apparently single mode signal individually correlate with porosity and to isolate the fast and slow waves from medial-lateral insonification of the calcaneus. The Bayesian technique was applied to trabecular bone data from eight human calcanei insonified in the medial-lateral direction. The phase velocity, slope of attenuation (nBUA), and amplitude were determined for both the fast and slow waves. The porosity was assessed by micro-computed tomography (microCT) and ranged from 78.7% to 94.1%. The method successfully separated the fast and slow waves from medial-lateral insonification of the calcaneus. The phase velocity for both the fast and slow wave modes showed an inverse correlation with porosity (R(2) = 0.73 and R(2) = 0.86, respectively). The slope of attenuation for both wave modes also had a negative correlation with porosity (fast wave: R(2) = 0.73, slow wave: R(2) = 0.53). The fast wave amplitude decreased with increasing porosity (R(2) = 0.66). Conversely, the slow wave amplitude modestly increased with increasing porosity (R(2) = 0.39).

In vitro estimation of fast and slow wave parameters of thin trabecular bone using space-alternating generalized expectation-maximization algorithm

In testing cancellous bone using ultrasound, two types of longitudinal Biot's waves are observed in the received signal. These are known as fast and slow waves and their appearance depend on the alignment of bone trabeculae in the propagation path and the thickness of the specimen under test (SUT). They can be used as an effective tool for the diagnosis of osteoporosis because wave propagation behavior depends on the bone structure. However, the identification of these waves in the received signal can be difficult to achieve. In this study, ultrasonic wave propagation in a 4mm thick bovine cancellous bone in the direction parallel to the trabecular alignment is considered. The observed Biot's fast and slow longitudinal waves are superimposed; which makes it difficult to extract any information from the received signal. These two waves can be separated using the space alternating generalized expectation maximization (SAGE) algorithm. The latter has been used mainly in speech processing. In this new approach, parameters such as, arrival time, center frequency, bandwidth, amplitude, phase and velocity of each wave are estimated. The B-Scan images and its associated A-scans obtained through simulations using Biot's finite-difference time-domain (FDTD) method are validated experimentally using a thin bone sample obtained from the femoral-head of a 30 months old bovine.

Determining attenuation properties of interfering fast and slow ultrasonic waves in cancellous bone

Previous studies have shown that interference between fast waves and slow waves can lead to observed negative dispersion in cancellous bone. In this study, the effects of overlapping fast and slow waves on measurements of the apparent attenuation as a function of propagation distance are investigated along with methods of analysis used to determine the attenuation properties. Two methods are applied to simulated data that were generated based on experimentally acquired signals taken from a bovine specimen. The first method uses a time-domain approach that was dictated by constraints imposed by the partial overlap of fast and slow waves. The second method uses a frequency-domain log-spectral subtraction technique on the separated fast and slow waves. Applying the time-domain analysis to the broadband data yields apparent attenuation behavior that is larger in the early stages of propagation and decreases as the wave travels deeper. In contrast, performing frequency-domain analysis on the separated fast waves and slow waves results in attenuation coefficients that are independent of propagation distance. Results suggest that features arising from the analysis of overlapping two-mode data may represent an alternate explanation for the previously reported apparent dependence on propagation distance of the attenuation coefficient of cancellous bone.

Numerical investigation of ultrasound refraction caused by oblique orientation of trabecular network in cancellous bone

Ultrasound propagation through cancellous bone can be greatly affected by the trabecular structure. In the present study, the ultrasound propagation for the oblique orientation of the trabecular network was numerically investigated using 3-D finite-difference time-domain (FDTD) simulations. The models of cancellous bone were reconstructed from X-ray microcomputed tomographic (μCT) images of a bovine bone. Cancellous bone models with various orientations of the trabecular network were realized by cutting the μCT images rotated from 0 to 90°. Ultrasound waveforms propagating through these cancellous bone models were simulated while changing the receiving position. The refraction of the ultrasound wave for the oblique angle of the main orientation was investigated on the basis of the variation in the arrival time and peak amplitude. As the propagation direction approached the direction parallel to the main orientation, the arrival time of the first peak became less and the peak amplitude became smaller. This means that the wave of the first peak, which corresponded to a fast wave, propagated in the direction perpendicular to the main orientation. In addition, a strong correlation between the first-peak amplitude and the arrival time was observed in the porosity range of 0.68 to 0.85, in which the slope of the amplitude with respect to time increased linearly with porosity.