Interlaboratory comparison of ultrasonic backscatter, attenuation, and speed measurements

In a study involving 10 different sites, independent results of measurements of ultrasonic properties on equivalent tissue-mimicking samples are reported and compared. The properties measured were propagation speed, attenuation coefficients, and backscatter coefficients. Reasonably good agreement exists for attenuation coefficients, but less satisfactory results were found for propagation speeds. As anticipated, agreement was not impressive in the case of backscatter coefficients. Results for four sites agreed rather well in both absolute values and frequency dependence, and results from other sites were lower by as much as an order of magnitude. The study is valuable for laboratories doing quantitative studies.

Ultrasonic tissue characterization with integrated backscatter. Acute myocardial ischemia, reperfusion, and stunned myocardium in patients

We have previously shown in studies of experimental animals that myocardium exhibits a cardiac cycle-dependent variation of integrated backscatter that reflects regional myocardial contractile performance and that is blunted promptly after arterial occlusion and recovers after reperfusion. To define the clinical utility of ultrasonic tissue characterization with integrated backscatter for detection of acute myocardial infarction and reperfusion, 21 patients (14 men and seven women) were studied in the cardiac care unit within the first 24 hours (mean time, 11.3 hours; range, 3.5-23.8 hours) after the onset of symptoms indicative of acute myocardial infarction with conventional two-dimensional and M-mode echocardiography and with analysis of integrated backscatter. The magnitude of cyclic variation of integrated backscatter was measured from several sites within acute infarct regions and normal regions remote from the infarct zone for each patient. The average magnitude of cyclic variation among all patients (n = 21) was 4.8 +/- 0.5 dB in normal regions compared with 0.8 +/- 0.3 dB in infarct regions (p less than 0.05) within the first 24 hours after the onset of symptoms. Among the patients who had two studies, 15 (mean, 7.1 days; range, 2-31 days for second study) underwent coronary arteriography to define vessel patency. In patients with vessels with documented patency (n = 10), the magnitude of cyclic variation in infarct regions increased over time from 1.3 +/- 0.6 to 2.5 +/- 0.5 dB from the initial to final study (p less than 0.05). Patients with occluded infarct-related arteries (n = 5) exhibited no significant recovery of cyclic variation (0.3 +/- 0.3-0.6 +/- 0.3 dB). A blinded analysis of standard two-dimensional echocardiographic images revealed no significant recovery of wall thickening in either group over the same time intervals. Ultrasonic tissue characterization promptly detects acute myocardial infarction and may delineate potential beneficial effects of coronary artery reperfusion manifest by restoration of cyclic variation of integrated backscatter in the presence of severe wall motion abnormalities.

Ultrasound integrated backscatter tissue characterization of remote myocardial infarction in human subjects

To determine whether quantitative ultrasound tissue characterization differentiates normal myocardial regions from segments of remote infarction, 32 consecutive patients with a diagnosis of previous myocardial infarction were evaluated. Images were obtained in real time with a modified two-dimensional ultrasound system capable of providing continuous signals in proportion to the logarithm of integrated backscatter along each A line. In 15 patients, adequate parasternal long-axis images that delineated both normal and infarct segments were obtained with standard time-gain compensation. Image data were analyzed to yield both magnitude and delay (electrocardiographic R wave to nadir normalized for the QT interval) of the cyclic variation of backscatter. Cyclic variation was present in 55 of 56 normal myocardial sites, averaging (mean +/- SEM) 3.2 +/- 0.2 dB in magnitude and exhibiting a mean normalized delay of 0.87 +/- 0.03. The magnitude of cyclic variation in infarct segments was significantly reduced to 1.1 +/- 0.2 dB (42 sites), and the delay was markedly increased to 1.47 +/- 0.12 (21 sites) (p less than 0.0001 for both). In 20 of 42 infarct sites, no cyclic variation was detectable. Thus, ultrasound tissue characterization quantitatively differentiated infarct segments from normal myocardium in patients with remote myocardial infarction.

Interlaboratory comparison of ultrasonic backscatter coefficient measurements from 2 to 9 MHz

OBJECTIVE

As are the attenuation coefficient and sound speed, the backscatter coefficient is a fundamental ultrasonic property that has been used to characterize many tissues. Unfortunately, there is currently far less standardization for the ultrasonic backscatter measurement than for the other two, as evidenced by a previous American Institute of Ultrasound in Medicine (AIUM)-sponsored interlaboratory comparison of ultrasonic backscatter, attenuation, and speed measurements (J Ultrasound Med 1999; 18:615-631). To explore reasons for these disparities, the AIUM Endowment for Education and Research recently supported this second interlaboratory comparison, which extends the upper limit of the frequency range from 7 to 9 MHz.

METHODS

Eleven laboratories were provided with standard test objects designed and manufactured at the University of Wisconsin (Madison, WI). Each laboratory was asked to perform ultrasonic measurements of sound speed, attenuation coefficients, and backscatter coefficients. Each laboratory was blinded to the values of the ultrasonic properties of the test objects at the time the measurements were performed.

RESULTS

Eight of the 11 laboratories submitted results. The range of variation of absolute magnitude of backscatter coefficient measurements was about 2 orders of magnitude. If the results of 1 outlier laboratory are excluded, then the range is reduced to about 1 order of magnitude. Agreement regarding frequency dependence of backscatter was better than reported in the previous interlaboratory comparison. For example, when scatterers were small compared with the ultrasonic wavelength, experimental frequency-dependent backscatter coefficient data obtained by the participating laboratories were usually consistent with the expected Rayleigh scattering behavior (proportional to frequency to the fourth power).

CONCLUSIONS

Greater standardization of backscatter measurement methods is needed. Measurements of frequency dependence of backscatter are more consistent than measurements of absolute magnitude.

Early identification with ultrasonic integrated backscatter of viable but stunned myocardium in dogs

It has been shown that canine and human hearts exhibit a cardiac cycle-dependent variation of integrated backscatter (cyclic variation) that reflects intrinsic regional contractile performance. To determine whether ultrasound tissue characterization can identify viable though stunned myocardium before recovery of regional wall thickening, transient ischemic injury was produced in eight open chest dogs for 15 min followed by reperfusion for 2 h. Cyclic variation and wall thickening were measured before ischemia, at 15 min after the onset of ischemia and at selected intervals after the onset of reperfusion from multiple sites within the ischemic zone with a novel combined two-dimensional and M-mode acquisition system. Cyclic variation and wall thickening were computed from digitized M-mode integrated backscatter images with an algorithm developed and validated for this purpose. Magnitude and "delay" of cyclic variation and wall thickening were compared. Delay represents the degree of synchrony of regional cyclic variation or wall thickening with global ventricular mechanical systole. Baseline cyclic variation and wall thickening magnitudes were 3.8 +/- 0.2 dB and 37 +/- 1.4%, respectively. With ischemia, cyclic variation and wall thickening decreased to 1.7 +/- 0.2 dB and 17 +/- 2%, respectively (p less than 0.05, compared with baseline). Cyclic variation recovered to baseline levels within 20 min after reperfusion (3.3 +/- 0.4 dB, p = NS). Wall thickening remained depressed for 2 h after the onset of reperfusion (23 +/- 2%, p less than 0.05 compared with baseline). Delay of cyclic variation in a unitless ratio expressed as delay (in milliseconds) divided by the QT interval (in milliseconds) increased from 0.87 +/- 0.03 at baseline to 1.10 +/- 0.12 with ischemia, a change consistent with mild asynchrony, and returned to baseline (0.95 +/- 0.07, p = NS compared with baseline) within 20 min after reperfusion. Delay of wall thickening was 0.88 +/- 0.02 at baseline, increased to 1.23 +/- 0.09 with ischemia and remained significantly increased 2 h after reperfusion (1.07 +/- 0.05, p less than 0.05 compared with baseline). Recovery time constants for cyclic variation and wall thickening with reperfusion reflected earlier restoration of cyclic variation (8.1 min) than of wall thickening (420.5 min). Thus, cyclic variation recovers before wall thickening with reperfusion. Its analysis appears to provide a useful index of the presence of viable and potentially salvageable tissue in regions of stunned myocardium that is independent of wall thickening.

Measurements of phase velocity and group velocity in human calcaneus

Ultrasonic velocity in calcaneus correlates highly with bone mineral density, which is a good predictor of osteoporotic fracture risk. Several commercial bone sonometers perform a velocity measurement based on the transit time of a broadband pulse to assess skeletal status. This approach is somewhat problematic, however, because several authors have reported ambiguities in measurements in calcaneus. Phase velocity is an alternative that may be less dependent on device spectral characteristics. In addition, dispersion (the frequency-dependence of phase velocity) is a fundamental property worth investigating to increase understanding of interaction between ultrasound and bone. To compare two group-velocity measurement methods and one phase-velocity measurement method, a polycarbonate sample (for method validation) and 24 human calcanei were investigated in vitro. Phase velocity in calcaneus at 500 kHz was 1511 m/s +/- 30 m/s (mean +/- standard deviation). Average phase velocity decreased approximately linearly with frequency (-18 m/s MHz). The two group velocity measurements were comparable and tended to be slightly lower than phase velocity. The magnitude of dispersion showed little correlation with bone mineral density.

The effects of frequency-dependent attenuation and dispersion on sound speed measurements: applications in human trabecular bone

Sound speed may be measured by comparing the transit time of a broadband ultrasonic pulse transmitted through an object with that transmitted through a reference water path. If the speed of sound in water and the thickness of the sample are known, the speed of sound in the object may be computed. To measure the transit time differential, a marker such as a zero-crossing, may be used. A sound speed difference between the object and water shifts all markers backward or forward. Frequency-dependent attenuation and dispersion may alter the spectral characteristics of the waveform, thereby distorting the locations of markers and introducing variations in sound-speed estimates. Theory is derived to correct for this distortion for Gaussian pulses propagating through linearly attenuating, weakly dispersive media. The theory is validated using numerical analysis, measurements on a tissue mimicking phantom, and on 24 human calcaneus samples in vitro. Variations in soft tissue-like media are generally not exceptionally large for most applications but can be substantial, particularly for high bandwidth pulses propagating through media with high attenuation coefficients. At 500 kHz, variations in velocity estimates in bone can be very substantial, on the order of 40 to 50 m/s because of the high attenuation coefficient of bone. In trabecular bone, the effects of frequency-dependent attenuation are considerable, and the effects of dispersion are negligible.

Assessment of bone density using ultrasonic backscatter

The goal of this project was to investigate the utility of ultrasonic backscatter for the assessment of bone status. Ultrasound offers a low-cost, portable, nonionizing alternative or complement to common X-ray- or radioisotope (gamma ray)-based methods of bone densitometry. Ultrasonic backscatter may provide useful information not revealed by ultrasonic attenuation and sound-speed densitometers. Backscatter is sensitive to microstructural variations in acoustic impedance and should therefore provide information regarding architecture (which is related to fracture risk), as well as density. Ultrasonic backscatter at 2.25 MHz and CT bone densitometric data have been acquired from 10 healthy human volunteers. The degree of correlation between CT and ultrasonic backscatter is high (r = 0.87, p < 0.001). The envelope signal-to-noise ratio was 1.81 +/- 0.08 (mean +/- standard deviation). This suggests that the number of scatterers per resolution cell is large, the radiofrequency signal approximately obeys circular Gaussian statistics, and the envelope obeys Rayleigh statistics. These results indicate promise for ultrasonic backscatter as a substitute for or an adjunct to other ultrasonic measurements (attenuation and sound speed) and X-ray measurements for the assessment of bone status.

Ultrasonic scattering from cancellous bone: a review

This paper reviews theory, measurements, and computer simulations of scattering from cancellous bone reported by many laboratories. Three theoretical models (binary mixture, Faran cylinder, and weak scattering) for scattering from cancellous bone have demonstrated some consistency with measurements of backscatter. Backscatter is moderately correlated with bone mineral density in human calcaneus in vitro (r(2) = 0.66 - 0.68). Backscatter varies approximately as frequency cubed and trabecular thickness cubed in human calcaneus and femur in vitro. Backscatter from human calcaneus and bovine tibia exhibits substantial anisotropy. So far, backscatter has demonstrated only modest clinical utility. Computer simulation models have helped to elucidate mechanisms underlying scattering from cancellous bones.

Application of autoregressive spectral analysis to cepstral estimation of mean scatterer spacing

The problem of estimation of mean scatterer spacing in an object containing regularly spaced structures is addressed. An autoregressive (AR) spectral estimation method is compared with a conventional fast Fourier transform (FFT)-based approach for this task. Regularly spaced structures produce a periodicity in the power spectrum of ultrasonic backscatter. This periodicity is manifested as a peak in the cepstrum. A phantom was constructed for comparison of the two methods. It contained regularly spaced nylon filaments. It also contained randomly positioned glass spheres that produced incoherent backscatter. In an experiment in which this target was interrogated using broadband ultrasound, the AR spectral estimate offered considerable improvement over the FFT when the analysis gate length was on the order of the structural dimension. Advantages included improved resolution, reduction in bias and variance of scatterer spacing estimates, and greater resistance to ringing artifacts. Data were also acquired from human liver in vivo. AR spectral estimates on human data exhibited a decreased dependence on gate length. These results offer promise for enhanced spatial resolution and accuracy in ultrasonic tissue characterization and nondestructive evaluation of materials.

Biologically relevant photoacoustic imaging phantoms with tunable optical and acoustic properties

Established medical imaging technologies such as magnetic resonance imaging and computed tomography rely on well-validated tissue-simulating phantoms for standardized testing of device image quality. The availability of high-quality phantoms for optical-acoustic diagnostics such as photoacoustic tomography (PAT) will facilitate standardization and clinical translation of these emerging approaches. Materials used in prior PAT phantoms do not provide a suitable combination of long-term stability and realistic acoustic and optical properties. Therefore, we have investigated the use of custom polyvinyl chloride plastisol (PVCP) formulations for imaging phantoms and identified a dual-plasticizer approach that provides biologically relevant ranges of relevant properties. Speed of sound and acoustic attenuation were determined over a frequency range of 4 to 9 MHz and optical absorption and scattering over a wavelength range of 400 to 1100 nm. We present characterization of several PVCP formulations, including one designed to mimic breast tissue. This material is used to construct a phantom comprised of an array of cylindrical, hemoglobin-filled inclusions for evaluation of penetration depth. Measurements with a custom near-infrared PAT imager provide quantitative and qualitative comparisons of phantom and tissue images. Results indicate that our PVCP material is uniquely suitable for PAT system image quality evaluation and may provide a practical tool for device validation and intercomparison.

Ultrasonic attenuation in human calcaneus from 0.2 to 1.7 MHz

Ultrasonic attenuation has been demonstrated to be a useful measurement in the diagnosis of osteoporosis. Most studies have employed ultrasound in a range of frequencies from about 200 kHz-300 kHz to 600 kHz-1 MHz, and many have assumed a linear dependence of attenuation on frequency. In order to investigate the attenuation properties of human calcaneus at higher frequencies, 16 defatted human calcanea were interrogated in vitro using two matched pairs of transducers with center frequencies of 500 kHz and 2.25 MHz. The linear dependence of attenuation on frequency seems to extend up to at least 1.7 MHz. The correlation between attenuation coefficient and frequency from 400 kHz to 1.7 MHz was r = 0.999 (95% confidence interval, CI, = 0.998-1.00). The measurements suggest that some deviations from linear frequency dependence of attenuation may occur at lower frequencies (below 400 kHz), however.

Conditionally Increased Acoustic Pressures in Nonfetal Diagnostic Ultrasound Examinations Without Contrast Agents: A Preliminary Assessment

The mechanical index (MI) has been used by the US Food and Drug Administration (FDA) since 1992 for regulatory decisions regarding the acoustic output of diagnostic ultrasound equipment. Its formula is based on predictions of acoustic cavitation under specific conditions. Since its implementation over 2 decades ago, new imaging modes have been developed that employ unique beam sequences exploiting higher-order acoustic phenomena, and, concurrently, studies of the bioeffects of ultrasound under a range of imaging scenarios have been conducted. In 2012, the American Institute of Ultrasound in Medicine Technical Standards Committee convened a working group of its Output Standards Subcommittee to examine and report on the potential risks and benefits of the use of conditionally increased acoustic pressures (CIP) under specific diagnostic imaging scenarios. The term "conditionally" is included to indicate that CIP would be considered on a per-patient basis for the duration required to obtain the necessary diagnostic information. This document is a result of that effort. In summary, a fundamental assumption in the MI calculation is the presence of a preexisting gas body. For tissues not known to contain preexisting gas bodies, based on theoretical predications and experimentally reported cavitation thresholds, we find this assumption to be invalid. We thus conclude that exceeding the recommended maximum MI level given in the FDA guidance could be warranted without concern for increased risk of cavitation in these tissues. However, there is limited literature assessing the potential clinical benefit of exceeding the MI guidelines in these tissues. The report proposes a 3-tiered approach for CIP that follows the model for employing elevated output in magnetic resonance imaging and concludes with summary recommendations to facilitate Institutional Review Board (IRB)-monitored clinical studies investigating CIP in specific tissues.

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.

Relationships of ultrasonic backscatter with ultrasonic attenuation, sound speed and bone mineral density in human calcaneus

Ultrasonic attenuation and sound speed have been investigated in trabecular bone by numerous authors. Ultrasonic backscatter has received much less attention. To investigate relationships among these three ultrasonic parameters and bone mineral density (BMD), 30 defatted human calcanei were investigated in vitro. Normalized broadband ultrasonic attenuation (nBUA), sound speed (SOS), and logarithm of ultrasonic backscatter coefficient (LBC) were measured. Bone mineral density was assessed using single-beam dual energy x-ray absorptiometry (DEXA). The correlation coefficients of least squares linear regressions of the three individual ultrasound (US) parameters with BMD were 0.84 (nBUA), 0.84 (SOS) and 0.79 (LBC). The 95% confidence intervals for the correlation coefficients were 0. 69-0.92 (nBUA), 0.68-0.92 (SOS) and 0.60-0.90 (LBC). The correlations among pairs of US variables ranged from 0.63-0.79. Variations in nBUA accounted for r(2) = 62% of the variations in LBC. Variations in SOS accounted for r(2) = 40% of the variations in LBC. These results suggest that ultrasonic backscattering properties may contain substantial information not already contained in nBUA and SOS. A multiple regression model including all three US variables was somewhat more predictive of BMD than a model including only nBUA and SOS.

The dependence of ultrasonic backscatter on trabecular thickness in human calcaneus: theoretical and experimental results

Trabecular thickness within cancellous bone is an important determinant of osteoporotic fracture risk. Noninvasive assessment of trabecular thickness potentially could yield useful diagnostic information. Faran's theory of elastic scattering from a cylindrical object immersed in a fluid has been used to predict the dependence of ultrasonic backscatter on trabecular thickness. The theory predicts that, in the range of morphological and material properties expected for trabecular bone, the backscatter coefficient at 500 kHz should be approximately proportional to trabecular thickness to the power of 2.9. Experimental measurements of backscatter coefficient were performed on 43 human calcaneus samples in vitro. Mean trabecular thicknesses on the 43 samples were assessed using micro computed tomography (CT). A power law fit to the data showed that the backscatter coefficient empirically varied as trabecular thickness to the 2.8 power. The 95% confidence interval for this exponent was 1.7 to 3.9. The square of the correlation coefficient for the linear regression to the log transformed data was 0.40. This suggests that 40% of variations in backscatter may be attributed to variations in trabecular thickness. These results reinforce previous studies that offered validation for the Faran cylinder model for prediction of scattering properties of cancellous bone, and provide added evidence for the potential diagnostic utility of the backscatter measurement.

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.

Interference between wave modes may contribute to the apparent negative dispersion observed in cancellous bone

Previous work has shown that ultrasonic waves propagating through cancellous bone often exhibit a linear-with-frequency attenuation coefficient, but a decrease in phase velocity with frequency (negative dispersion) that is inconsistent with the causality-imposed Kramers-Kronig relations. In the current study, interfering wave modes similar to those observed in bone are shown to potentially contribute to the observed negative dispersion. Biot theory, the modified Biot-Attenborogh model, and experimental results are used to aid in simulating multiple-mode wave propagation through cancellous bone. Simulations entail constructing individual wave modes exhibiting a positive dispersion using plausible velocities and amplitudes, and then summing the individual modes to create mixed-mode output wave forms. Results of the simulations indicate that mixed-mode wave forms can exhibit negative dispersion when analyzed conventionally under the assumption that only one wave is present, even when the individual interfering waves exhibit positive dispersions in accordance with the Kramers-Kronig relations. Furthermore, negative dispersion is observed when little or no visual evidence of interference exists in the time-domain data. Understanding the mechanisms responsible for the observed negative dispersion could aid in determining the true material properties of cancellous bone, as opposed to the apparent properties measured using conventional data analysis techniques.

Measurements of ultrasonic backscatter coefficients in human liver and kidney in vivo

Ultrasonic backscatter coefficients, in the range of 2.0-4.0 MHz, were measured in normal human livers and kidneys in vivo. In liver, data were acquired and analyzed from 15 normal volunteers and 19 patients with hepatitis. No significant difference between normal and chronic hepatitis was found. The power-law fit to the backscatter coefficient in normal liver as a function of frequency was eta(f) = 4.5 x 10(-5) f1.6 cm-1 Str-1. This is comparable to that measured by other investigators in in vitro preparations of human and animal liver and to that measured by two other teams of investigators in in vivo human liver. In kidney, data were acquired from 11 normal volunteers. The power-law fit to the backscatter coefficient in normal kidney was eta (f) = 2.3 x 10(-5) f2.1 cm-1 Str-1. This is in the range of that measured by other investigators in in vitro preparations of human and animal kidney. In order to assess the system dependence of in vivo abdominal organ backscatter coefficients, measurements were performed using two different ultrasonic data-acquisition systems. The two systems exhibited close agreement.

A stratified model to predict dispersion in trabecular bone

Frequency-dependent phase velocity (dispersion) has previously been measured in trabecular bone by several groups. In contrast to most biologic tissues, phase velocity in trabecular bone tends to decrease with frequency. A stratified model, consisting of alternating layers of bone and marrow (in vivo) or water (in vitro), has been employed in an attempt to explain this phenomenon. Frequency-dependent phase velocity was measured from 300 to 700 kHz in 1) phantoms consisting of regularly spaced thin parallel layers of polystyrene sheets in water and 2) 30 calcaneus samples in vitro. For the polystyrene phantoms, the agreement between theory and experiment was good. For the calcaneus samples, the model has some limited usefulness (uncertainty of about 5%) in predicting average phase velocity. More importantly, the model seems to perform consistently well for predicting the frequency dependence of phase velocity in calcaneus.

Characterization of trabecular bone using the backscattered spectral centroid shift

Ultrasonic attenuation in bone in vivo is generally measured using a through-transmission method at the calcaneus. Although attenuation in calcaneus has been demonstrated to be a useful predictor for osteoporotic fracture risk, measurements at other clinically important sites, such as hip and spine, could potentially contain additional useful diagnostic information. Through-transmission measurements may not be feasible at these sites due to complex bone shapes and the increased amount of intervening soft tissue. Centroid shift from the backscattered signal is an index of attenuation slope and has been used previously to characterize soft tissues. In this paper, centroid shift from signals backscattered from 30 trabecular bone samples in vitro were measured. Attenuation slope also was measured using a through-transmission method. The correlation coefficient between centroid shift and attenuation slope was -0.71. The 95% confidence interval was (-0.86, -0.47). These results suggest that the backscattered spectral centroid shift may contain useful diagnostic information potentially applicable to hip and spine.

Comparison of measurements of phase velocity in human calcaneus to Biot theory

Biot's theory for elastic propagation in porous media has previously been shown to be useful for modeling the dependence of phase velocity on porosity in bovine cancellous bone in vitro. In the present study, Biot's theory is applied to measurements of porosity-dependent phase velocity in 53 human calcanea in vitro. Porosity was measured using microcomputed tomography for some samples (n = 23) and estimated based on bone mineral densitometry for the remaining samples (n = 30). The phase velocity at 500 kHz was measured in a water tank using a through-transmission technique. Biot's theory performed well for the prediction of the dependence of sound speed on porosity. The trend was quasilinear, but both the theory and experiment show similar slight curvature. The root mean square error (RMSE) of predicted versus measured sound speed was 15.8 m/s.

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.

Radiological Society of North America/Quantitative Imaging Biomarker Alliance Shear Wave Speed Bias Quantification in Elastic and Viscoelastic Phantoms

OBJECTIVES

To quantify the bias of shear wave speed (SWS) measurements between different commercial ultrasonic shear elasticity systems and a magnetic resonance elastography (MRE) system in elastic and viscoelastic phantoms.

METHODS

Two elastic phantoms, representing healthy through fibrotic liver, were measured with 5 different ultrasound platforms, and 3 viscoelastic phantoms, representing healthy through fibrotic liver tissue, were measured with 12 different ultrasound platforms. Measurements were performed with different systems at different sites, at 3 focal depths, and with different appraisers. The SWS bias across the systems was quantified as a function of the system, site, focal depth, and appraiser. A single MRE research system was also used to characterize these phantoms using discrete frequencies from 60 to 500 Hz.

RESULTS

The SWS from different systems had mean difference 95% confidence intervals of ±0.145 m/s (±9.6%) across both elastic phantoms and ± 0.340 m/s (±15.3%) across the viscoelastic phantoms. The focal depth and appraiser were less significant sources of SWS variability than the system and site. Magnetic resonance elastography best matched the ultrasonic SWS in the viscoelastic phantoms using a 140 Hz source but had a - 0.27 ± 0.027-m/s (-12.2% ± 1.2%) bias when using the clinically implemented 60-Hz vibration source.

CONCLUSIONS

Shear wave speed reconstruction across different manufacturer systems is more consistent in elastic than viscoelastic phantoms, with a mean difference bias of < ±10% in all cases. Magnetic resonance elastographic measurements in the elastic and viscoelastic phantoms best match the ultrasound systems with a 140-Hz excitation but have a significant negative bias operating at 60 Hz. This study establishes a foundation for meaningful comparison of SWS measurements made with different platforms.