Influence of microarchitecture alterations on ultrasonic backscattering in an experimental simulation of bovine cancellous bone aging

An experimental model which can simulate physical changes that occur during aging was developed in order to evaluate the effects of change of mineral content and microstructure on ultrasonic properties of bovine cancellous bone. Timed immersion in hydrochloric acid was used to selectively alter the mineral content. Scanning electron microscopy and histological staining of the acid-treated trabeculae demonstrated a heterogeneous structure consisting of a mineralized core and a demineralized layer. The presence of organic matrix contributed very little to normalized broadband ultrasound attenuation (nBUA) and speed of sound. All three ultrasonic parameters, speed of sound, nBUA and backscatter coefficient, were sensitive to changes in apparent density of bovine cancellous bone. A two-component model utilizing a combination of two autocorrelation functions (a densely populated model and a spherical distribution) was used to approximate the backscatter coefficient. The predicted attenuation due to scattering constituted a significant part of the measured total attenuation (due to both scattering and absorption mechanisms) for bovine cancellous bone. Linear regression, performed between trabecular thickness values and estimated from the model correlation lengths, showed significant linear correlation, with R(2)=0.81 before and R(2)=0.80 after demineralization. The accuracy of estimation was found to increase with trabecular thickness.     

Characterization of dense bovine cancellous bone tissue microstructure by ultrasonic backscattering using weak scattering models

A weak scattering model was proposed for the ultrasonic frequency-dependent backscatter in dense bovine cancellous bone, using two autocorrelation functions to describe the medium: one with discrete homogeneities (spherical distribution of equal spheres) and another, which considers tissue as an inhomogeneous continuum (densely populated medium). The inverse problem to estimate trabecular thickness of bone tissue has been addressed. A combination of the two autocorrelation functions was required to closely approximate the backscatter from bovine bone with various microarchitecture, given that the shape of trabeculae ranges from a rodlike to a platelike shape. Because of the large variation in trabecular thickness, both at an intraspecimen and an interspecimen level, thickness distributions for individual trabeculae for each bone specimen were obtained, and dominant trabecular sizes were determined. Comparison of backscatter measurements to theoretical predictions indicated that there were more than one dominant trabecular sizes that scatter sound for most specimens. Linear regression, performed between dominant trabecular thickness and estimated correlation length, showed significant linear correlation (R(2)=0.81). Attenuation due to scattering by a continuous distribution of scatterers was predicted to be linear over a frequency range from 0.3 to 0.9 MHz, suggesting a possibility that scattering may be a significant source of attenuation.     

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

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

Theoretical and experimental study of the ultrasonic attenuation in bovine cancellous bone

In this study a theoretical approach for the estimation of ultrasonic attenuation is proposed. The approach combines two models which take into account both absorption and scattering. Attenuation due to absorption is studied by using the Biot’s analytical model whereas that due to scattering is described by means of a generalized weak scattering model which is formulated for binary mixtures. The scattering model takes account of the density fluctuation of the porous medium in addition to the propagation velocity fluctuation. For the calculation of the attenuation coefficient due to absorption, experimental values have been used to link size of pores to porosity. The theoretical results have been compared with experimental data obtained on bovine cancellous bone samples filled with water. Using an immersion acoustic transmission method, the ultrasonic attenuation has been measured at a frequency range between 0.1 and 1.0 MHz for 12 bovine cancellous bone samples with a porosity range between 40% and 70%. The prediction of attenuation with this model appears to correspond more closely to its experimentally observed behavior. This study indicates that scattering is the predominant mechanism which is responsible for attenuation in trabecular bone. Furthermore, it shows that the density fluctuations contribute significantly to the phenomenon of attenuation and cannot thus be neglected.     

骨小梁材料特性对超声背散射信号的影响

基于时域有限差分法(FDTD)建立了松质骨的超声背散射仿真系统,研究了骨小梁材料特性对超声背散射信号的影响。首次得到松质骨中的超声背散射系数(BSC)和积分背散射系数(IBC)随骨小梁材料参数(密度、拉梅常数、黏度系数及声阻抗系数)的变化关系。研究结果表明,IBC随骨小梁密度的增加而增加;BSC和IBC随拉梅常数的增加而增加、随第一黏度系数的增加而近似线性地减小,第二黏度的变化对背散射信号的影响很小;背散射参数随阻抗系数的增加而减小。说明松质骨中的超声背散射特性不仅受骨矿密度(BMD)和骨微结构的影响,还与骨小梁的材料参数密切相关。研究结果有利于理解松质骨中超声的背散射特性,对松质骨骨质状况的评价有一定帮助。      

Effect of trabecular bone material properties on ultrasonic backscattering signals

In this study, ultrasonic backscattering signals in cancellous bones were obtained by finite difference time domain （FDTD） simulations, and the effect of trabecular material properties on these signals was analyzed. The backscatter coefficient （BSC） and integrated backscatter coefficient （IBC） were numerically investigated for varying trabecular bone material properties, including density, Lame coefficients, viscosities, and resistance coefficients. The results show that the BSC is a complex function of trabecular bone density, and the IBC increases as density increases. The BSC and IBC increase with the first and second Lame coefficients. While not very sensitive to the second viscosity of the trabeculae, the BSC and IBC decrease as the first viscosity and resistance coefficients increase. The results demonstrate that, in addition to bone mineral density （BMD） and microarchitecture, trabecular material properties significantly influence ultrasonic backseattering signals in cancellous bones. This research furthers the understanding of ultrasonic backscattering in cancellous bones and the characterization of cancellous bone status.     

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

Although bone sonometry has been demonstrated to be useful in the diagnosis of osteoporosis, much remains to be learned about the processes governing the interactions between ultrasound and bone. In order to investigate these processes, ultrasonic attenuation and backscatter in two orientations were measured in 43 human calcaneal specimens in vitro at 500 kHz. In the mediolateral (ML) orientation, the ultrasound propagation direction is approximately perpendicular to the trabecular axes. In the anteroposterior (AP) orientation, a wide range of angles between the ultrasound propagation direction and trabecular axes is encountered. Average attenuation slope was 18% greater while average backscatter coefficient was 50% lower in the AP orientation compared with the ML orientation. Backscatter coefficient in both orientations approximately conformed to a cubic dependence on frequency, consistent with a previously reported model. These results support the idea that absorption is a greater component of attenuation than scattering in human calcaneal trabecular bone.     

Reproducibility of the quantitative assessment of cartilage morphology and trabecular bone structure with magnetic resonance imaging at 7 T

To assess the reproducibility of quantitative measurements of cartilage morphology and trabecular bone structure of the knee at 7 T, high-resolution sagittal spoiled gradient-echo images and high-resolution axial fully refocused steady-state free-precession (SSFP) images from six healthy volunteers were acquired with a 7-T scanner. The subjects were repositioned between repeated scans to test the reproducibility of the measurements. The reproducibility of each measurement was evaluated using the coefficient(s) of variation (CV). The computed CV were 1.13% and 1.55% for cartilage thickness and cartilage volume, respectively, and were 2.86%, 1.07%, 2.27% and 3.30% for apparent bone volume over total volume fraction (app.BV/TV), apparent trabecular number (app.Tb.N), apparent trabecular separation (app.Tb.Sp) and apparent trabecular thickness (app.Tb.Th), respectively. The results demonstrate that quantitative assessment of cartilage morphology and trabecular bone structure is reproducible at 7 T and motivates future musculoskeletal applications seeking the high-field strength's superior signal-to-noise ratio.     

Frequency dependence of ultrasonic backscatter from human trabecular bone: theory and experiment

A model describing the frequency dependence of backscatter coefficient from trabecular bone is presented. Scattering is assumed to originate from the surfaces of trabeculae, which are modeled as long thin cylinders with radii small compared with the ultrasonic wavelength. Experimental ultrasonic measurements at 500 kHz, 1 MHz, and 2.25 MHz from a wire target and from trabecular bone samples from human calcaneus in vitro are reported. In both cases, measurements are in good agreement with theory. For mediolateral insonification of calcaneus at low frequencies, including the typical diagnostic range (near 500 kHz), backscatter coefficient is proportional to frequency cubed. At higher frequencies, the frequency response flattens out. The data also suggest that at diagnostic frequencies, multiple scattering effects on the average are relatively small for the samples investigated. Finally, at diagnostic frequencies, the data suggest that absorption is likely to be a larger component of attenuation than scattering.     

Statistics of ultrasonic scatterer size estimation with a reference phantom

A theoretical expression for the variance of scatterer size estimates is derived for a modified least squares size estimator used in conjunction with a reference phantom method for backscatter coefficient measurement. A Gaussian spatial autocorrelation function is assumed. Simulations and phantom experiments were performed to verify the results for backscatter and size variances. The dependence of size estimate errors upon free experimental parameters is explored. Implications of the findings for the optimization of scatterer size estimation are discussed. The utility of scatterer size parametric imaging is examined through the signal to noise ratio comparison with standard ultrasonic B-mode imaging.     

X射线显微CT用于大鼠骨小梁结构分析的研究

应用X射线显微CT(X-μCT)对正常及骨质疏松大鼠的骨小梁结构进行了分析,并与骨组织形态计量法的测量值进行了比较,探讨了X射线光谱技术在骨结构分析中的应用。实验对大鼠骨样品进行X-μCT扫描,扫描条件为 80 kVp,80 μA,360°旋转,帧平均4帧,角度增益 0.4°,分辨率14 μm。三维重建并分析了骨小梁结构,结构参数包括骨体积分数(BV/TV)、骨小梁厚度(Tb.Th)、骨小梁数量(Tb.N)以及骨小梁间隔(Tb.Sp)。结果表明,采用X-μCT分析不同组大鼠的骨小梁结构参数值之间存在显著差异(P<0.05),测定值与传统骨组织形态计量法的测定值显著相关,其中胫骨骨小梁BV/TV,Tb.Th,Tb.N,Tb.Sp的相关系数r分别为0.984,0.960,0.995,0.988,腰椎骨小梁各结构参数的相关系数分别为0.938,0.968,0.877,0.951。因此,X-μCT可以较好地呈现并区分正常骨组织、骨质疏松骨组织以及经雌激素治疗后骨组织的微观结构,可以实现对骨小梁结构参数的分析测定,与骨组织形态计量法相比是一种更精确、立体、快速且无损测量骨微结构和评价骨质量的方法。     

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

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

Relationships of quantitative ultrasound parameters with cancellous bone microstructure in human calcaneus in vitro

Ultrasound parameters (attenuation, phase velocity, and backscatter), bone mineral density (BMD), and microarchitectural features were measured on 29 human cancellous calcaneus samples in vitro. Regression analysis was performed to predict ultrasound parameters from BMD and microarchitectural features. The best univariate predictors of the ultrasound parameters were the indexes of bone quantity: BMD and bone volume fraction (BV/TV). The most predictive univariate models for attenuation, phase velocity, and backscatter coefficient yielded adjusted squared correlation coefficients of 0.69-0.73. Multiple regression models yielded adjusted correlation coefficients of 0.74-0.83. Therefore attenuation, phase velocity, and backscatter are primarily determined by bone quantity, but multiple regression models based on bone quantity plus microarchitectural features achieve slightly better predictive performance than models based on bone quantity alone.     

Short-TE projection reconstruction NMR microscopy of trabecular bone.

The aim of this study was to assess the potential of projection reconstruction (PR) NMR microscopy in the quantitative evaluation of trabecular bone architecture. Short-TE PR spin-echo microimages were acquired at 7.05 T on normal bone explants. The main structural parameters such as bone volume fraction (BVF), trabecular thickness (Tb.Th.) and trabecular separation (Tb.Sp.) were obtained from the 3D microimages using the method of directed secants. Quantitative structural data were then compared with those derived from conventional spin-echo microimages. Our study indicates that projection reconstruction NMR microscopy promises to be more accurate than the conventional FTI method in the analysis of trabecular bone.     

Study on the mean trabecular bone spacing based on fundamental frequency estimation method

Ultrasonic backscatter signals from the cancellous bone can be used to diagnose osteoporosis effectively due to its ability to provide the information of bone microstructure. Mean trabecular bone spacing(MTBS)is one of the important parameters for characterization of bone microstructure.This paper proposed a MTBS estimating method based on the fundamental frequency estimation,which was applied to backscatter signals from simulations, and in vitro bovine trabeculae.The estimated MTBS were compared with those of simplified inverse filter tracking(SIFT)algorithm and autoregressive(AR)cepstrum method.The results demonstrated that the proposed method is very robust for the MTBS estimation with more precise estimates and smaller estimated variance in the presence of a small signal-to-noise ratio (SNR),and a large scattering strength ratio of diffuse scatterers to regular ones.     

Fundamental precision limitations for measurements of frequency dependence of backscatter: applications in tissue-mimicking phantoms and trabecular bone.

Various models for ultrasonic scattering from trabecular bone have been proposed. They may be evaluated to a certain extent by comparison with experimental measurements. In order to appreciate limitations of these comparisons, it is important to understand measurement precision. In this article, an approach proposed by Lizzi and co-workers is adapted to model precision of estimates of frequency-dependent backscatter for scattering targets (such as trabecular bone) that contain many scatterers per resolution cell. This approach predicts uncertainties in backscatter due to the random nature of the interference of echoes from individual scatterers as they are summed at the receiver. The model is validated in experiments on a soft-tissue-mimicking phantom and on 24 human calcaneus samples interrogated in vitro. It is found that while random interference effects only partially explain measured variations in the magnitude of backscatter, they are virtually entirely responsible for observed variations in the frequency dependence (exponent of a power law fit) of backscatter.     

超声背散射法评价松质骨状况的研究

超声背散射法评价松质骨状况及诊断骨质疏松症是近年来医学超声领域内的研究热点之一,现已取得了显著的进展。本文将介绍近年来超声背散射法及其参量评价松质骨状况的研究进展,并分析超声背散射相关参量频谱质心偏移量（SCS）和平均骨小梁间距（TbSp）与骨矿密度（BMD）的相关性。研究结果表明,超声背散射参量与BMD有较高的相关性。最后提出了将来研究中需要努力的方向。     

Frequency-dependent attenuation and backscatter coefficients in bovine trabecular bone from 0.2 to 1.2 MHz

The frequency-dependent attenuation and backscatter coefficients were measured in 25 bovine femoral trabecular bone samples from 0.2 to 1.2 MHz. When the average attenuation coefficient was fitted to a nonlinear power law α(f)=α(0)+α(1)f(n), the exponent n was found to be 1.65. In contrast, the average backscatter coefficient was fitted to a power law η(f)=η(1)f(n) and the exponent n was measured as 3.25. The apparent bone density was significantly correlated with the parameter α(1) (0.2-0.7 MHz: r = 0.852, 0.6-1.2 MHz: r = 0.832) as well as the backscatter coefficient (0.5 MHz: r = 0.751, 1.0 MHz: r = 0.808).     

Statistics of the envelope of ultrasonic backscatter from human trabecular bone

The paper describes the investigations intended to compare the results of experimental measurements of backscattering properties of the trabecular bone with the results of computer simulations. Ultrasonic RF echoes were collected using two bone scanners operating at 0.58 and 1.3 MHz. The simulations of the backscattered RF echoes were performed using the scattering model of the trabecular bone that consisted of cylindrical and spherical elements uniformly distributed in water-like medium. For each measured or simulated RF backscatter the statistical properties of the signal envelope were determined. Experimental results suggest deviations of the backscattering properties from the Rayleigh distribution. The results of simulation suggest that deviation from Rayleigh distribution depends on the variation of trabeculae diameters and the number of thin trabeculae. Experimentally determined deviations corresponded well to the deviations calculated from simulated echoes assuming trabeculae thickness variation equaled to the earlier published histomorphometric study results.     

Effects of risedronate in a rat model of osteopenia due to orchidectomy and disuse: densitometric, histomorphometric and microtomographic studies

Dual energy X-ray absorptiometry (DXA), histomorphometry and X-ray microtomography (microCT) were used to assess effects of risedronate and testosterone in a combined rat model of orchidectomy (ORX) and local paralysis induced by botulinum neurotoxin (BTX). Four groups of mature rats were studied for 1 month: SHAM operated; ORX and right hindlimb immobilization (BTX); ORX+BTX+risedronate or testosterone. Changes in bone and body composition were measured by DXA (BMC, lean and fat mass), histomorphometry (BV/TV(2D), Tb.Th and microarchitectural parameters) and microCT (BV/TV(3D), SMI and cortical parameters). ORX and BTX had additive effects on bone loss since differences were maximized on the immobilized bone. The decrease in BMC on the tibial metaphysis reached -33.6% vs. -11.3% in the non-immobilized limb. BV/TV and Tb.N decreased and Tb.Sp increased in both hindlimbs whereas Tb.Th was significantly lower only in the immobilized limb. Decrease of tibial cortical area and thickness was greater in the immobilized limb. Risedronate prevented BMC, BV/TV and architecture loss but not reduction in Tb.Th. Cortical bone was preserved only in the non-immobilized limb. Testosterone was unable to prevent trabecular and cortical bone loss, but it prevents loss of whole body lean mass. In conclusion, ORX and BTX resulted in additive effects on bone loss. Risedronate had protective effects on trabecular bone loss but was less effective on cortical bone.