Relationships of Ultrasonic Backscatter With Bone Densities and Microstructure in Bovine Cancellous Bone

Full-Waveform Inversion Imaging of Cortical Bone Using Phased Array Tomography

Classic ultrasound bone imaging modalities usually demand either a prior knowledge or an advanced estimation on speed of sound (SoS), which not only renders to a burdensome imaging process but also supplies a limited resolution. To overcome these drawbacks, this article proposed a frequency-domain full-waveform inversion (FDFWI) modality using phased array tomography for high-accuracy cortical bone imaging. A transmission scenario of ultrasound wave in 2-D space was presented in the frequency domain to simulate the forward wavefield propagation. Iterations in the inversion process were performed by matching the simulation wavefield to the experimental one from low to high discrete frequency points. Moreover, the association between the maximum initial frequency and the initial SoS model was explored to prevent the occurrence of cycle-skipping phenomenon, which could lead to the outcomes being trapped in local minima. The feasibility and effectiveness of the proposed imaging scheme were testified by simulation, phantom, and ex-vivo studies, with mean relative errors of cortical part being 3.18%, 8.71%, and 9.36%, respectively. It is verified that the proposed FDFWI method is an effective way for parametric imaging of cortical bone without any prior knowledge of sound speed.

Backscatter measurement of cancellous bone using the ultrasound transit time spectroscopy

Recently, ultrasound transit time spectroscopy (UTTS) was proposed as a promising method for bone quantitative ultrasound measurement. Studies have showed that UTTS could estimate the bone volume fraction and other trabecular bone structure in ultrasonic through-transmission measurements. The goal of this study was to explore the feasibility of UTTS to be adapted in ultrasonic backscatter measurement and further evaluate the performance of backscattered ultrasound transit time spectrum (BS-UTTS) in the measurement of cancellous bone density and structure. First, taking ultrasonic attenuation into account, the concept of BS-UTTS was verified on ultrasonic backscatter signals simulated from a set of scatterers with different positions and intensities. Then, in vitro backscatter measurements were performed on 26 bovine cancellous bone specimens. After a logarithmic compression of the BS-UTTS, a linear fitting of the log-compressed BS-UTTS versus ultrasonic propagated distance was performed and the slope and intercept of the fitted line for BS-UTTS were determined. The associations between BS-UTTS parameters and cancellous bone features were analyzed using simple linear regression. The results showed that the BS-UTTS could make an accurate deconvolution of the backscatter signal and predict the position and intensity of the simulated scatterers eliminating phase interference, even the simulated backscatter signal was with a relatively low signal-to-noise ratio. With varied positions and intensities of the scatterers, the slope of the fitted line for the log-compressed BS-UTTS versus ultrasonic propagated distance (i.e., slope of BS-UTTS for short) yield a high agreement (r2 = 99.84%-99.96%) with ultrasonic attenuation in simulated backscatter signal. Compared with the high-density cancellous bone, the low-density specimen showed more abundant backscatter impulse response in the BS-UTTS. The slope of BS-UTTS yield a significant correlation with bone mineral density (r = 0.87; p < 0.001), BV/TV (r = 0.87; p < 0.001), and cancellous bone microstructures (r up to 0.87; p < 0.05). The intercept of BS-UTTS was also significantly correlated with bone densities (r = -0.87; p < 0.001) and trabecular structures (|r|=0.43-0.80; p < 0.05). However, the slope of the BS-UTTS underestimated attenuation when measurements were performed experimentally. In addition, a significant non-linear relationship was observed between the measured attenuation and the attenuation estimated by the slope of the BS-UTTS. This study demonstrated that the UTTS method could be adapted to ultrasonic backscatter measurement of cancellous bone. The derived slope and intercept of BS-UTTS could be used in the measurement of bone density and microstructure. The backscattered ultrasound transit time spectroscopy might have potential in the diagnosis of osteoporosis in the clinic.

The Protocol of Ultrasonic Backscatter Measurements of Musculoskeletal Properties

This study aims to introduce the protocol for ultrasonic backscatter measurements of musculoskeletal properties based on a novel ultrasonic backscatter bone diagnostic (UBBD) instrument. Dual-energy X-ray absorptiometry (DXA) can be adopted to measure bone mineral density (BMD) in the hip, spine, legs and the whole body. The muscle and fat mass in the legs and the whole body can be also calculated by DXA body composition analysis. Based on the proposed protocol for backscatter measurements by UBBD, ultrasonic backscatter signals can be measured in vivo, deriving three backscatter parameters [apparent integral backscatter (AIB), backscatter signal peak amplitude (BSPA) and the corresponding arrival time (BSPT)]. AIB may provide important diagnostic information about bone properties. BSPA and BSPT may be important indicators of muscle and fat properties. The standardized backscatter measurement protocol of the UBBD instrument may have the potential to evaluate musculoskeletal characteristics, providing help for promoting the application of the backscatter technique in the clinical diagnosis of musculoskeletal disorders (MSDs), such as osteoporosis and muscular atrophy.

Conical shell X-ray beam tomosynthesis and micro-computed tomography for microarchitectural characterisation

Bone quality is commonly used to diagnose bone diseases such as osteoporosis, with many studies focusing on microarchitecture for fracture prediction. In this study a bovine distal femur was imaged using both micro-computed tomography (µCT) and tomosynthesis using focal construct geometry (FCG) for comparison of microarchitectural parameters. Six regions of interest (ROIs) were compared between the two imaging modalities, with both global and adaptive methods used to binarize the images. FCG images were downsampled to the same pixel size as the µCT images. Bone morphometrics were determined using BoneJ, for each imaging modality, binarization technique and ROI. Bone area/total area was found to have few significant differences between FCG and µCT (p < 0.05 for two of six ROIs). Fractal Dimension had only one significant difference (p < 0.05 for one of six ROIs) between µCT and downsampled FCG (where pixel size was equalized). Trabecular thickness and trabecular spacing were observed to follow trends as observed for the corresponding µCT images, although many absolute values were significantly different (p < 0.05 for between one and six ROIs depending on image types used). This study demonstrates the utility of tomosynthesis for measurement of microarchitectural morphometrics.

Ultrasonic Through-Transmission Measurements of Human Musculoskeletal and Fat Properties

The study described here was aimed at investigating the feasibility of using the ultrasonic through-transmission technique to estimate human musculoskeletal and fat properties. Five hundred eighty-two volunteers were assessed by dual-energy X-ray absorptiometry (DXA) and ultrasonic transmission techniques. Bone mineral density (BMD), muscle and fat mass were measured for both legs and the whole body. Hip BMD and spine BMD were also measured. Ultrasonic transmission measurements were performed on the heel, and the measured parameters were broadband ultrasound attenuation (BUA), speed of sound (SOS), ultrasonic stiffness index (SI), T-score and Z-score, which were significantly correlated with all measured BMDs. The optimal correlation was observed between SI and left-leg BMD (p < 0.001) before and after adjustment for age, sex and body mass index (BMI). The linear and partial correlation analyses revealed that BUA and SOS were closely associated with muscle and fat mass, respectively. Multiple regressions revealed that muscle and fat mass significantly contributed to the prediction of transmission parameters, explaining up to 17.83% (p < 0.001) variance independently of BMD. The results suggest that the ultrasonic through-transmission technique could help in the clinical diagnosis of skeletal and muscular system diseases.

In Vivo Comparison of Backscatter Techniques for Ultrasonic Bone Assessment at the Femoral Neck

Ultrasonic techniques are being developed to detect changes in cancellous bone caused by osteoporosis. The goal of this study was to test the relative in vivo performance of eight backscatter parameters developed over the last several years for ultrasonic bone assessment: apparent integrated backscatter (AIB), frequency slope of apparent backscatter (FSAB), frequency intercept of apparent backscatter (FIAB), normalized mean of the backscatter difference (nMBD), normalized slope of the backscatter difference (nSBD), normalized intercept of the backscatter difference (nIBD), normalized backscatter amplitude ratio (nBAR) and backscatter amplitude decay constant (BADC). Backscatter measurements were performed on the left and right femoral necks of 80 adult volunteers (age = 25 ± 11 y) using an imaging system equipped with a convex array transducer. For comparison, additional ultrasonic measurements were performed at the left and right heel using a commercially available heel-bone ultrasonometer that measured the stiffness index. Six of the eight backscatter parameters (all but nSBD and nIBD) exhibited similar and highly significant (p < 0.000001) left-right correlations (0.51 ≤ R ≤ 0.68), indicating sensitivity to naturally occurring variations in bone tissue. Left-right correlations for the stiffness index measured at the heel (R = 0.75) were not significantly better than those produced by AIB, FSAB and FIAB. The short-term precisions of AIB, nMBD, nBAR and BADC (7.8%-11.7%) were comparable to that of the stiffness index measured with the heel-bone ultrasonometer (7.5%).

Index-Rotated Fast Ultrasound Imaging of Cortical Bone Based on Predicted Velocity Model

Due to the significant acoustic impedance contrast at cortical boundaries, highly inside attenuation, and the unknown sound velocity distribution, accurate ultrasound cortical bone imaging remains a challenge, especially for the traditional pulse-echo modalities using unique sound velocity. Moreover, the large amounts of data recorded by multielement probe results in a relatively time-consuming reconstruction process. To overcome these limitations, this article proposed an index-rotated fast ultrasound imaging method based on predicted velocity model (IR-FUI-VP) for cortical cross section ultrasound tomography (UST) imaging, utilizing ray-tracing synthetic aperture (RTSA). In virtue of ring probe, the sound velocity model was predicted in advance using bent-ray inversion (BRI). With the predicted velocity model, index-rotated fast ultrasound imaging (IR-FUI) was further applied to image the cortical cross sections in the sectors corresponding to the dynamic apertures (DAs) and ring center. The final result was merged by all sector images. One cortical bone phantom and two ex vivo bovine femurs were utilized to demonstrate the performance of the proposed method. Compared to the conventional synthetic aperture (SA) imaging, the method can not only accurately image the outer cortical boundary but also precisely reconstruct the inner cortical surface. The mean relative errors of the predicted sound velocity in the region of interest (ROI) were all smaller than 7%, and the mean errors of cortical thickness are all less than 0.31 mm. The reconstructed images of bovine femurs were in good agreement with the reference images scanned by micro-computed tomography ( μ CT) with respect to the morphology and thickness. The speed of IR-FUI is about 3.73 times faster than the traditional SA. It is proved that the proposed IR-FUI-VP-based UST is an effective way for fast and accurate cortical bone imaging.

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.

The feasibility study of the transmission mode photoacoustic measurement of human calcaneus bone in vivo

The photoacoustic (PA) technique is uniquely positioned for biomedical applications primarily due to its ability to visualize optical absorption contrast in deep tissue at ultrasound resolution. In this work, via both three-dimensional (3D) numerical simulations and in vivo experiments on human subjects, we investigated the possibility of PA measurement of human calcaneus bones in vivo in a non-invasive manner, as well as its feasibility to differentiate osteoporosis patients from normal subjects. The results from the simulations and the experiments both demonstrated that, when one side of the heel is illuminated by laser with light fluence under the ANSI safety limit, the PA signal generated in the human calcaneus bone can be detected by an ultrasonic transducer at the other side of the heel (i.e. transmission mode). Quantitative power spectral analyses of the calcaneus bone PA signals were also conducted, demonstrating that the microarchitectural changes in calcaneus bone due to osteoporosis can be detected, as reflected by enhanced high frequency components in detected PA bone signal. Further statistical analysis of the experimental results from 10 osteoporosis patients and 10 healthy volunteers showed that the weighted frequency as a quantified PA spectral parameter can differentiate the two subject groups with statistical significance.

Fourier-Domain Ultrasonic Imaging of Cortical Bone Based on Velocity Distribution Inversion

There is a significant acoustic impedance contrast between the cortical bone and the surrounding soft tissue, resulting in difficulty for ultrasound penetration into bone tissue with high frequency. It is challenging for the conventional pulse-echo modalities to give accurate cortical bone images using uniform sound velocity model. To overcome these limitations, an ultrasound imaging method called full-matrix Fourier-domain synthetic aperture based on velocity inversion (FM-FDSA-VI) was developed to provide accurate cortical bone images. The dual linear arrays were located on the upper and lower sides of the imaging region. After full-matrix acquisition with two identical linear array probes facing each other, travel-time inversion was used to estimate the velocity distribution in advance. Then, full-matrix Fourier-domain synthetic aperture (FM-FDSA) imaging based on the estimated velocity model was applied twice to image the cortical bone, utilizing the data acquired from top and bottom linear array, respectively. Finally, to further improve the image quality, the two images were merged to give the ultimate result. The performance of the method was verified by two simulated models and two bone phantoms (i.e., regular and irregular hollow bone phantom). The mean relative errors of estimated sound velocity in the region-of-interest (ROI) are all below 12%, and the mean errors of cortical section thickness are all less than 0.3 mm. Compared to the conventional synthetic aperture (SA) imaging, the FM-FDSA-VI method is able to accurately image cortical bone with respect to the structure. Moreover, the result of irregular bone phantom was close to the image scanned by microcomputed tomography ( μ CT) in terms of macro geometry and thickness. It is demonstrated that the proposed FM-FDSA-VI method is an efficient way for cortical bone ultrasonic imaging.

Ultrasonic Backscatter Measurements of Human Cortical and Trabecular Bone Densities in a Head-Down Bed-Rest Study

This study aims to investigate the feasibility of quantitative ultrasonic backscatter in evaluating human cortical and trabecular bone densities in vivo based on a head-down-tilt bed rest study, with 36 participants tested through 90 d of bed rest and 180 d of recovery. Backscatter measurements were performed using an ultrasonic backscatter bone diagnostic instrument. Backscatter parameters were calculated with a dynamic signal-of-interest method, which was proposed to ensure the same ultrasonic interrogated volume in cortical and trabecular bones. The backscatter parameters exhibited significant correlations with site-matched bone densities provided by high-resolution peripheral quantitative computed tomography (0.33 < |R| < 0.72, p < 0.05). Some bone densities and backscatter parameters exhibited significant changes after the 90-d bed rest. The proposed method can be used to characterize bone densities, and the portable ultrasonic backscatter bone diagnostic device might be used to non-invasively reveal mean bone loss (across a group of people) after long-term bed rest and microgravity conditions of spaceflight missions.

Three-dimensional-printed replica models of bone for experimentally decoupling trabecular bone properties contribution to ultrasound propagation parameters

A detailed investigation of the relationship between ultrasonic (US) properties and trabecular bone microstructure is difficult because of the great variability in the bone loss process. The aim of this work was twofold. First, to verify by compressive tests that the three-dimensional (3D)-printer is able to produce precisely and repeatedly "bone replica models" of different size and density. Following, replicas of the original specimens with two different polymers and thinned trabeculae models were used to investigate US properties (speed of sound, SOS, and backscatter coefficient), aiming to deconvolute the influence of material properties on ultrasound characteristics. The results revealed that matrix material properties influence only the magnitude of the backscatter coefficient, whereas the characteristic undulated patterns are related to the trabecular structure. Simulation of perforation and thinning of cancellous bone, associated with bone loss, showed that SOS and mechanical properties were reduced perfectly linearly with apparent density when structure deteriorated. The 3D-printed bone replicas have the potential to enable systematic investigations of the influence of structure on both acoustical and mechanical properties and evaluate changes caused by bone loss. The development of replicas from materials with properties close to those of bone will permit quantitative conclusions for trabecular bone.

Ultrasonic Bone Assessment Using the Backscatter Amplitude Decay Constant

Ultrasonic backscatter techniques are being developed to detect changes in bone caused by osteoporosis. The present study introduces a new technique that measures the exponential decay in the amplitude of the backscatter signal quantified by a parameter called the backscatter amplitude decay constant (BADC). Measurements were performed on 54 specimens of cancellous bone from 14 human femurs using a 3.5-MHz transducer. Six methods were tested to determine BADC. The recommended method measures the time slope of the natural log of the rectified signal. Measured values of BADC ranged from approximately 0.1 μs^-1 to 0.6 μs^-1. Moderate to strong correlations (Spearman's ρ >0.7) were found between BADC and the density and microstructural characteristics of the specimens determined using X-ray microcomputed tomography. The results of this study suggest that BADC may be able to detect changes in the density and microstructure of cancellous bone caused by osteoporosis and other diseases.

A Combined Ultrasonic Backscatter Parameter for Bone Status Evaluation in Neonates

Metabolic bone disease (MBD) is one of the major complications of prematurity. Ultrasonic backscatter technique has the potential to be a portable and noninvasive method for early diagnosis of MBD. This study firstly applied CAS to neonates, which was defined as a linear combination of the apparent integrated backscatter coefficient (AIB) and spectral centroid shift (SCS). The objective was to evaluate the feasibility of ultrasonic backscatter technique for assessing neonatal bone health using AIB, SCS, and CAS. Ultrasonic backscatter measurements at 3.5 MHz, 5.0 MHz, and 7.5 MHz were performed on a total of 505 newborns within 48 hours after birth. The values of backscatter parameters were calculated and compared among gestational age groups. Correlations between backscatter parameters, gestational age, anthropometric indices, and biochemical markers were analyzed. The optimal predicting models for CAS were determined. The results showed term infants had lower SCS and higher AIB and CAS than preterm infants. Gestational age and anthropometric indices were negatively correlated with SCS (|r| = 0.45 – 0.57, P < 0.001), and positively correlated with AIB (|r| = 0.36 – 0.60, P < 0.001) and CAS (|r| = 0.56 – 0.69, P < 0.001). Biochemical markers yielded weak or nonsignificant correlations with backscatter parameters. CAS had relatively stronger correlations with the neonatal variables than AIB and SCS. At 3.5 MHz and 5.0 MHz, only gestational age (P < 0.001) independently contributed to the measurements of CAS, and could explain up to 40.5% – 44.3% of CAS variation. At 7.5 MHz, the combination of gestational age (P < 0.001), head circumference (P = 0.002), and serum calcium (P = 0.037) explained up to 40.3% of CAS variation. This study suggested ultrasonic backscatter technique was feasible to evaluate neonatal bone status. CAS was a promising parameter to provide more information about bone health than AIB or SCS alone.

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.

Ultrasonic Backscatter Difference Measurement of Bone Health in Preterm and Term Newborns

Metabolic bone disease of prematurity remains a significant problem for preterm infants. Quantitative ultrasound (QUS) has potential as a non-invasive tool for assessing bone health of newborns. The aim of this study was to assess bone health in preterm and term newborns using ultrasonic backscatter difference measurement. This study analyzed a total of 493 neonates, including 239 full-term infants (gestational age [GA] >37 wk), 201 preterm I infants (GA: 32-37 wk) and 53 extreme preterm II infants (GA <32 wk). Ultrasonic backscatter measurements were performed on the calcaneus of infants at birth, and the normalized mean of the backscatter difference spectrum (nMBD) was calculated as an ultrasonic index of neonatal bone status. Simple and multiple linear regressions were performed to determine the association of ultrasonic nMBD with GA, anthropometric characteristics and biochemical markers. Statistically significant differences in GA, anthropometric characteristics (birth weight, birth length [BL], birth head circumference and body mass index [BMI]) and biochemical markers (alkaline phosphatase, serum calcium and serum phosphate) were observed among preterm and term infants. The nMBD for term infants (median = 3.72 dB/μs, interquartile range [IR] = 1.95 dB/μs) was significantly higher than that for preterm I infants (median = 1.95 dB/μs, IR = 3.12 dB/μs), which was, in turn, significantly higher than that for preterm II infants (median = 0.19 dB/μs, IR = 3.50 dB/μs). The nMBD yielded moderate correlations (ρ = 0.57-0.62, p < 0.001) with GA and anthropometric characteristics and weak correlations (|ρ| = 0.08-0.21, p < 0.001 or not significant) with biochemical markers. Multivariate regressions revealed that only BL (p = 0.002) and BMI (p = 0.032) yielded significantly independent contributions to the nMBD measurement, and combinations of BL and BMI could explain up to 42% of the variation of nMBD in newborn infants. This study found that ultrasonic backscatter difference measurement might be helpful in bone health evaluation in preterm and term newborns. The utility of ultrasonic backscatter measurement in diagnosis of metabolic bone disease in infants should be investigated further.