Interlaboratory comparison of ultrasonic backscatter, attenuation, and speed measurements

Variation of ultrasound image lateral spectrum with assumed speed of sound and true scatterer density

One method of estimating sound speed in diagnostic ultrasound imaging consists of choosing the speed of sound that generates the sharpest image, as evaluated by the lateral frequency spectrum of the squared B-mode image. In the current work, simulated and experimental data on a typical (47 mm aperture, 3.3-10.0 MHz response) linear array transducer are used to investigate the accuracy of this method. A range of candidate speeds of sound (1240-1740 m/s) was used, with a true speed of sound of 1490 m/s in simulations and 1488 m/s in experiments. Simulations of single point scatterers and two interfering point scatterers at various locations with respect to each other gave estimate errors of 0.0-2.0%. Simulations and experiments of scatterer distributions with a mean scatterer spacing of at least 0.5 mm gave estimate errors of 0.1-4.0%. In the case of lower scatterer spacing, the speed of sound estimates become unreliable due to a decrease in contrast of the sharpness measure between different candidate speeds of sound. This suggests that in estimating speed of sound in tissue, the region of interest should be dominated by a few, sparsely spaced scatterers. Conversely, the decreasing sensitivity of the sharpness measure to speed of sound errors for higher scatterer concentrations suggests a potential method for estimating mean scatterer spacing.

Reference characterisation of sound speed and attenuation of the IEC agar-based tissue-mimicking material up to a frequency of 60 MHz

To support the development of clinical applications of high-frequency ultrasound, appropriate tissue-mimicking materials (TMMs) are required whose acoustic properties have been measured using validated techniques. This paper describes the characterisation of the sound speed (phase velocity) and attenuation coefficient of the International Electrotechnical Commission (IEC) agar-based TMM over the frequency range 1 to 60 MHz. Measurements implemented a broadband through-transmission substitution immersion technique over two overlapping frequency ranges, with co-axially aligned 50 MHz centre-frequency transducers employed for characterisation above 15 MHz. In keeping with usual practice employed within the technical literature, thin acoustic windows (membranes) made of 12-μm-thick Mylar protected the TMM from water damage. Various important sources of uncertainty that could compromise measurement accuracy have been identified and evaluated through a combination of experimental studies and modelling. These include TMM sample thickness, measured both manually and acoustically, and the influence of interfacial losses that, even for thin protective membranes, are significant at the frequencies of interest. In agreement with previous reports, the attenuation coefficient of the IEC TMM exhibited non-linear frequency dependence, particularly above 20 MHz, yielding a value of 0.93 ± 0.04 dB cm(-1) MHz(-1) at 60 MHz, derived at 21 ± 0.5°C. For the first time, phase velocity, measured with an estimated uncertainty of ±3.1 m s(-1), has been found to be dispersive over this extended frequency range, increasing from 1541 m s(-1) at 1 MHz to 1547 m s(-1) at 60 MHz. This work will help standardise acoustic property measurements, and establishes a reference measurement capability for TMMs underpinning clinical applications at elevated frequencies.

An improved tissue-mimicking polyacrylamide hydrogel phantom for visualizing thermal lesions with high-intensity focused ultrasound

A recipe was created to improve the tissue-mimicking (TM) bovine serum albumin (BSA) polyacrylamide hydrogel (PAG) reported in our previous study (Choi MJ, Guntur SR, Lee KI, Paeng DG, Coleman AJ. Ultrasound Med Biol 2013; 29:439-448). In that work, the concentration of acrylamide in TM BSA PAG was increased to make its attenuation coefficient the same as that of a tissue. However, this increase made the PAG stiffer and less homogeneous. In addition, the increase in acrylamide caused a significant increase in temperature over the denaturation threshold of BSA during polymerization, which required forced cooling so that the PAG did not become opaque at room temperature after polymerization. To eliminate those shortcomings, we substituted the increased acrylamide with a viscous polysaccharide liquid (corn syrup). The concentration of corn syrup was optimized to 20% (w/v, tested in the volume of 50 mL), so that the acoustic properties of the PAG would be close to those of human liver. The improved TM (iTM) BSA PAG constructed in this study had a speed of sound of 1588 ± 9 m/s, an attenuation coefficient of 0.51 ± 0.06 dB cm(-1) at 1 MHz and a backscattering coefficient of 0.22 ± 0.09 × 10(-3) sr(-1) cm(-1) MHz(-1). The density and acoustic impedance were 1057 kg/m(3) and 1.68 MRayl, respectively, and the non-linear parameter (B/A) was 5.9 ± 0.3. The thermal, optical and mechanical properties were almost the same as those of the BSA PAG (Lafon et al.2005). Experimental verification indicated that the thermal lesions visualized in the proposed iTM BSA PAG by high-intensity focused ultrasound were highly reproducible. In conclusion, iTM BSA PAG was proven to eliminate TM BSA PAG shortcomings effectively and is expected to be a promising test phantom for clinical high-intensity focused ultrasound device.

Experimental evaluation of spectral-based quantitative ultrasound imaging using plane wave compounding

Quantitative ultrasound (QUS) based on backscatter coefficient (BSC) estimation has shown potential for tissue characterization. Beamforming using plane wave compounding has advantages for echographic, Doppler, and elastographic imaging; however, to date, plane wave compounding has not been experimentally evaluated for the purpose of BSC estimation. In this study, two BSC-derived parameters (i.e., the BSC midband fit and intercept) were estimated from experimental data obtained using compound plane wave beamforming. For comparison, QUS parameters were also estimated from data obtained using both fixed focus and dynamic receive beamforming. An ultrasound imaging system equipped with a 9-MHz center frequency, 64-element array was used to collect data up to a depth of 45 mm. Two gelatin phantoms with randomly distributed 20-μm inclusions with a homogeneous scatterer concentration and a two-region scatterer concentration were used for assessing the precision and lateral resolution of QUS imaging, respectively. The use of plane wave compounding resulted in accurate QUS estimation (i.e., bias in the BSC parameters of less than 2 dB) and relatively constant lateral resolution (i.e., BSC midband fit 10% to 90% rise distance ranging between 1.0 and 1.5 mm) throughout a 45 mm field of view. Although both fixed focus and dynamic receive beamforming provided the same performance around the focal depth, the reduction in SNR away from the focus resulted in a reduced field of view in the homogeneous phantom (i.e., only 28 mm). The lateral resolution also degraded away from the focus, with up to a 2-fold and 10-fold increase in the rise distance at 20 mm beyond the focal depth for dynamic receive and fixed focus beamforming, respectively. These results suggest that plane wave compounding has the potential to improve the performance of spectral-based quantitative ultrasound over other conventional beamforming strategies.

Quantitative ultrasound imaging for monitoring in situ high-intensity focused ultrasound exposure

Quantitative ultrasound (QUS) imaging is hypothesized to map temperature elevations induced in tissue with high spatial and temporal resolution. To test this hypothesis, QUS techniques were examined to monitor high-intensity focused ultrasound (HIFU) exposure of tissue. In situ experiments were conducted on mammary adenocarcinoma tumors grown in rats and lesions were formed using a HIFU system. A thermocouple was inserted into the tumor to provide estimates of temperature at one location. Backscattered time-domain waveforms from the tissue during exposure were recorded using a clinical ultrasonic imaging system. Backscatter coefficients were estimated using a reference phantom technique. Two parameters were estimated from the backscatter coefficient (effective scatterer diameter (ESD) and effective acoustic concentration (EAC). The changes in the average parameters in the regions corresponding to the HIFU focus over time were correlated to the temperature readings from the thermocouple. The changes in the EAC parameter were consistently correlated to temperature during both heating and cooling of the tumors. The changes in the ESD did not have a consistent trend with temperature. The mean ESD and EAC before exposure were 120 ± 16 μm and 32 ± 3 dB/cm3, respectively, and changed to 144 ± 9 μm and 51 ± 7 dB/cm3, respectively, just before the last HIFU pulse was delivered to the tissue. After the tissue cooled down to 37 °C, the mean ESD and EAC were 126 ± 8 μm and 35 ± 4 dB/cm3, respectively. Peak temperature in the range of 50-60 °C was recorded by a thermocouple placed just behind the tumor. These results suggest that QUS techniques have the potential to be used for non-invasive monitoring of HIFU exposure.

Tissue-mimicking gel phantoms for thermal therapy studies

Tissue-mimicking phantoms that are currently available for routine biomedical applications may not be suitable for high-temperature experiments or calibration of thermal modalities. Therefore, design and fabrication of customized thermal phantoms with tailored properties are necessary for thermal therapy studies. A multitude of thermal phantoms have been developed in liquid, solid, and gel forms to simulate biological tissues in thermal therapy experiments. This article is an attempt to outline the various materials and techniques used to prepare thermal phantoms in the gel state. The relevant thermal, electrical, acoustic, and optical properties of these phantoms are presented in detail and the benefits and shortcomings of each type are discussed. This review could assist the researchers in the selection of appropriate phantom recipes for their in vitro study of thermal modalities and highlight the limitations of current phantom recipes that remain to be addressed in further studies.

Development of novel imaging probe for optical/acoustic radiation imaging (OARI)

PURPOSE

Optical/acoustic radiation imaging (OARI) is a novel imaging modality being developed to interrogate the optical and mechanical properties of soft tissues. OARI uses acoustic radiation force to generate displacement in soft tissue. Optical images before and after the application of the force are used to generate displacement maps that provide information about the mechanical properties of the tissue under interrogation. Since the images are optical images, they also represent the optical properties of the tissue as well. In this paper, the authors present the first imaging probe that uses acoustic radiation force in conjunction with optical coherence tomography (OCT) to provide information about the optical and mechanical properties of tissues to assist in the diagnosis and staging of epithelial cancers, and in particular bladder cancer.

METHODS

The OARI prototype probe consisted of an OCT probe encased in a plastic sheath, a miniaturized transducer glued to a plastic holder, both of which were encased in a 10 cm stainless steel tube with an inner diameter of 10 mm. The transducer delivered an acoustic intensity of 18 W/cm(2) and the OCT probe had a spatial resolution of approximately 10-20 μm. The tube was filled with deionized water for acoustic coupling and covered by a low density polyethylene cap. The OARI probe was characterized and tested on bladder wall phantoms. The phantoms possessed Young's moduli ranging from 10.2 to 12 kPa, mass density of 1.05 g/cm(3), acoustic attenuation coefficient of 0.66 dB/cm MHz, speed of sound of 1591 m/s, and optical scattering coefficient of 1.80 mm(-1). Finite element model (FEM) theoretical simulations were performed to assess the performance of the OARI probe.

RESULTS

The authors obtained displacements of 9.4, 8.7, and 3.4 μm for the 3%, 4%, and 5% bladder wall phantoms, respectively. This shows that the probe is capable of generating optical images, and also has the ability to generate and track displacements in tissue. This will provide information about the optical and mechanical properties of the tissue to assist in epithelial cancer detection. The corresponding theoretical FEM displacement was 5.8, 5.4, and 5.0 μm for the 3%, 4%, and 5% phantoms, respectively. Deviation between OARI displacement and FEM displacement is due to the resolution of the crosscorrelation algorithm used to track the displacement. To the authors' knowledge, this is the first probe that successfully combines OCT with a source of acoustic radiation force.

CONCLUSIONS

The OARI probe has the ability to provide information about the mechanical and optical properties of phantoms and soft tissue. This could prove useful in early epithelial cancer detection. Because the probe is 10 mm in diameter, it is currently only useful for skin and oral applications. The probe would have to be reduced in size to make it applicable for cancer detection in other internal sites. Future work will focus on utilizing phase-sensitive optical coherence elastography to obtain the resulting OARI displacements, improving the resolution of the probe, and enable physicians to better evaluate the mechanical properties of soft tissues.

Low-echo sphere phantoms and methods for assessing imaging performance of medical ultrasound scanners

Tissue-mimicking phantoms and software for quantifying the ability of human observers to detect small low-echo spheres as a function of depth have been developed. Detectability is related to the imager's ability to delineate the boundary of a 3-D object such as a spiculated tumor. The phantoms accommodate a broad range of transducer shapes and sizes. Three phantoms are described: one with 2-mm-diameter spheres (for higher frequencies), one with 3.2-mm-diameter spheres (for lower frequencies) and one with 4-mm-diameter spheres (for lower frequencies). The spheres are randomly distributed in each phantom. The attenuation coefficients of spheres and surroundings are nearly identical; thus, compromising shadowing or enhancement distal to spheres does not occur. Reproducibility results are given for pairs of independent data sets involving eight different combinations of scanner, transducer and console settings. The following comparison results are also reported: (i) only the selected frequency differs; (ii) transducers and scan parameters are nearly the same but manufacturers differ; (iii) ordinary B-scanning, spatial compounding and tissue harmonic imaging are addressed. The phantoms and software promise to be valuable tools for scanning system and setup comparisons and for acceptance testing.

Signal to noise ratio comparisons for ultrasound attenuation slope estimation algorithms

PURPOSE

Attenuation imaging has a promising role in the detection of tissue abnormalities. The authors have previously compared three different frequency domain ultrasound attenuation estimation methods, for accuracy and bias. The mean estimated attenuation value in a region of interest has been the determining factor of how well a method performs; however, the noise level has not been quantified for attenuation estimated using different methods.

METHODS

The authors compare three different frequency domain ultrasound attenuation estimation methods [the reference phantom method (RPM), the centroid downshift method (CEN), and the hybrid method (HYB)] using the signal to noise ratio (SNR) metric. Both simulated and experimental tissue-mimicking phantoms are used in the performance comparison study, evaluating the impact of the variation in acoustical properties.

RESULTS

For attenuation estimation in a tissue-mimicking phantom with a known attenuation coefficient of 0.5 dB/cm/MHz, all the three methods estimated the attenuation coefficient to be ≈ 0.49 dB/cm/MHz for a transmit center frequency of 6 MHz, however, the signal to noise ratio obtained was found to be 8.5, 5.7, and 2.2 for the HYB, RPM, and CEN methods, respectively. These results demonstrate the need for the SNR metric in the comparison of different algorithms and to evaluate the impact of varying different ultrasound system and tissue parameters.

CONCLUSIONS

In this paper, the authors demonstrate that although the estimated mean attenuation value with a region of interest may be closely estimated using different methods, the signal to noise ratio obtained of the estimates can vary significantly. The centroid downshift method presented with the lowest signal-to-noise ratio of the methods compared. The hybrid method was the least susceptible to changes in the acoustical properties and provided unbiased attenuation coefficient estimates with the highest signal-to-noise ratios.

Evaluation of the impact of backscatter intensity variations on ultrasound attenuation estimation

PURPOSE

Quantitative ultrasound based approaches such as attenuation slope estimation can be used to determine underlying tissue properties and eventually used as a supplemental diagnostic technique to B-mode imaging. The authors investigate the impact of backscatter intensity and frequency dependence variations on the attenuation slope estimation accuracy.

METHODS

The authors compare three frequency domain based attenuation slope estimation algorithms, namely, a spectral difference method, the reference phantom method, and two spectral shift methods: a hybrid method and centroid downshift method. Both the reference phantom and hybrid method use a tissue-mimicking phantom with well-defined acoustic properties to reduce system dependencies and diffraction effects. The normalized power spectral ratio obtained is then filtered by a Gaussian filter centered at the transmit center frequency in the hybrid method. A spectral shift method is then used to estimate the attenuation coefficient from the normalized and filtered spectrum. The centroid downshift method utilizes the shift in power spectrum toward lower frequencies with depth. Numerical phantoms that incorporate variations in the backscatter intensity from -3 to 3 dB, by varying the scatterer number density and variations in the scatterer diameters ranging from 10 to 100 μm are simulated. Experimental tissue mimicking phantoms with three different scatterer diameter ranges (5-40, 75-90, and 125-150 μm) are also used to evaluate the accuracy of the estimation methods.

RESULTS

The reference phantom method provided accurate results when the acoustical properties of the reference and the sample are well matched. Underestimation occurs when the reference phantom possessed a higher sound speed than the sample, and overestimation occurs when the reference phantom had a lower sound speed than the sample. The centroid downshift method depends significantly on the bandwidth of the power spectrum, which in turn depends on the frequency dependence of the backscattering. The hybrid method was the least susceptible to changes in the sample's acoustic properties and provided the lowest standard deviation in the numerical simulations and experimental evaluations.

CONCLUSIONS

No significant variations in the estimation accuracy of the attenuation coefficient were observed with an increase in the scatterer number density in the simulated numerical phantoms for the three methods. Changes in the scatterer diameters, which result in different frequency dependence of backscatter, do not significantly affect attenuation slope estimation with the reference phantom and hybrid approaches. The centroid method is sensitive to variations in the scatterer diameter due to the frequency shift introduced in the power spectrum.

A multimodality vascular imaging phantom of an abdominal aortic aneurysm with a visible thrombus

PURPOSE

With the continuous development of new stent grafts and implantation techniques, it has now become technically feasible to treat abdominal aortic aneurysms (AAA) with challenging anatomy using endovascular repair with standard, fenestrated, or branched stent-grafts. In vitro experimentations are very useful to improve stent-graft design and conformability or imaging guidance for stent-graft delivery or follow-up. Vascular replicas also help to better understand the limitation of endovascular approaches in challenging anatomy and possibly improve surgical planning or training by practicing high risk clinical procedures in the laboratory to improve outcomes in the operating room. Most AAA phantoms available have a very basic anatomy, which is not representative of the clinical reality. This paper presents a method of fabrication of a realistic AAA phantom with a visible thrombus, as well as some mechanical properties characterizing such phantom.

METHODS

A realistic AAA geometry replica of a real patient anatomy taken from a multidetector computed tomography (CT) scan was manufactured. To demonstrate the multimodality imaging capability of this new phantom with a thrombus visible in magnetic resonance (MR) angiography, CT angiography (CTA), digital subtraction angiography (DSA), and ultrasound, image acquisitions with all these modalities were performed by using standard clinical protocols. Potential use of this phantom for stent deployment was also tested. A rheometer allowed defining hyperelastic and viscoelastic properties of phantom materials.

RESULTS

MR imaging measurements of SNR and CNR values on T1 and T2-weighted sequences and MR angiography indicated reasonable agreement with published values of AAA thrombus and abdominal components in vivo. X-ray absorption also lay within normal ranges of AAA patients and was representative of findings observed on CTA, fluoroscopy, and DSA. Ultrasound propagation speeds for developed materials were also in concordance with the literature for vascular and abdominal tissues.

CONCLUSIONS

The mimicked abdominal tissues, AAA wall, and surrounding thrombus were developed to match imaging features of in vivo MR, CT, and ultrasound examinations. This phantom should be of value for image calibration, segmentation, and testing of endovascular devices for AAA endovascular repair.

Comparison between experimental and 2-D numerical studies of multiple scattering in Inconel600 by means of array probes

Ultrasonic non-destructive testing of polycrystalline structures can be disturbed by scattering at grain boundaries. Understanding and modeling this so-called "structural noise" is crucial for characterization as well as detection purposes. Structural noise can be considered as a fingerprint of the material under investigation, since it contains information about its microstructure. The interpretation of experimental data necessitates an accurate comprehension of complex phenomena that occur in multiple scattering media and thus robust scattering models. In particular, numerical models can offer the opportunity to realize parametrical studies on controlled microstructures. However, the ability of the model to simulate wave propagation in complex media must be validated. In that perspective, the main objective of the present work is to evaluate the ability of the finite-element code ATHENA 2D to reproduce typical features of multiple wave scattering in the context of ultrasonic non-destructive evaluation, with an array of sources and receivers. Experiments were carried out with a 64-element array, around 2 MHz. The sample was a mock-up of Inconel600 exhibiting a coarse grain structure with a known grain size distribution. The numerical model of this microstructure is based on Voronoi diagrams. Two physical parameters were used to compare numerical and experimental data: the coherent backscattering peak, and the singular value distribution of the array response matrix. Though the simulations are 2-D, a good agreement was found between simulated and experimental data.

A tissue mimicking polyacrylamide hydrogel phantom for visualizing thermal lesions generated by high intensity focused ultrasound

An optically transparent tissue-mimicking (TM) phantom whose acoustic properties are close to those of tissue was constructed for visualizing therapeutic effects by high intensity focused ultrasound (HIFU). The TM phantom was designed to improve a widely used standard bovine serum albumin (BSA) polyacrylamide hydrogel (PAG), which attenuated ultrasound far less than tissue and, unlike tissue, did not scatter ultrasound. A modified recipe has been proposed in the study by adding scattering glass beads with diameters of 40-80 μm (0.002% w/v) and by raising the concentration of acrylamide (30% v/v). The TM BSA-PAG constructed has an acoustic impedance of 1.67 MRayls, a speed of sound of 1576 m/s, an attenuation coefficient of 0.52 dB/cm at 1 MHz, a backscattering coefficient of 0.242 × 10(-3) 1/sr/cm at 1 MHz and a nonlinear parameter (B/A) of 5.7. These parameters are close to those of liver. The thermal and optical properties are almost the same as the standard BSA-PAG. The characteristic features of the thermal lesions by HIFU were observed to be more accurately visualized in the TM BSA-PAG than in the standard BSA-PAG. In conclusion, the proposed TM BSA-PAG acoustically mimics tissue better than the standard BSA-PAG and is expected to be preferentially used for assuring if a clinical HIFU device produces the thermal lesion as planned.

Comparison of ultrasound attenuation and backscatter estimates in layered tissue-mimicking phantoms among three clinical scanners

Backscatter and attenuation coefficient estimates are needed in many quantitative ultrasound strategies. In clinical applications, these parameters may not be easily obtained because of variations in scattering by tissues overlying a region of interest (ROI). The goal of this study is to assess the accuracy of backscatter and attenuation estimates for regions distal to nonuniform layers of tissue-mimicking materials. In addition, this work compares results of these estimates for "layered" phantoms scanned using different clinical ultrasound machines. Two tissue-mimicking phantoms were constructed, each exhibiting depth-dependent variations in attenuation or backscatter. The phantoms were scanned with three ultrasound imaging systems, acquiring radio frequency echo data for offline analysis. The attenuation coefficient and the backscatter coefficient (BSC) for sections of the phantoms were estimated using the reference phantom method. Properties of each layer were also measured with laboratory techniques on test samples manufactured during the construction of the phantom. Estimates of the attenuation coefficient versus frequency slope, α(0), using backscatter data from the different systems agreed to within 0.24 dB/cm-MHz. Bias in the α(0) estimates varied with the location of the ROI. BSC estimates for phantom sections whose locations ranged from 0 to 7 cm from the transducer agreed among the different systems and with theoretical predictions, with a mean bias error of 1.01 dB over the used bandwidths. This study demonstrates that attenuation and BSCs can be accurately estimated in layered inhomogeneous media using pulse-echo data from clinical imaging systems.

Time domain attenuation estimation method from ultrasonic backscattered signals

Ultrasonic attenuation is important not only as a parameter for characterizing tissue but also for compensating other parameters that are used to classify tissues. Several techniques have been explored for estimating ultrasonic attenuation from backscattered signals. In the present study, a technique is developed to estimate the local ultrasonic attenuation coefficient by analyzing the time domain backscattered signal. The proposed method incorporates an objective function that combines the diffraction pattern of the source/receiver with the attenuation slope in an integral equation. The technique was assessed through simulations and validated through experiments with a tissue mimicking phantom and fresh rabbit liver samples. The attenuation values estimated using the proposed technique were compared with the attenuation estimated using insertion loss measurements. For a data block size of 15 pulse lengths axially and 15 beamwidths laterally, the mean attenuation estimates from the tissue mimicking phantoms were within 10% of the estimates using insertion loss measurements. With a data block size of 20 pulse lengths axially and 20 beamwidths laterally, the error in the attenuation values estimated from the liver samples were within 10% of the attenuation values estimated from the insertion loss measurements.

Changes in ultrasonic properties of liver tissue in vitro during heating-cooling cycle concomitant with thermal coagulation

The present work considers the ultrasonic properties of porcine liver tissue in vitro measured during heating concomitant with thermal coagulation followed by natural cooling, so as to provide information about changes in the ultrasonic properties of the tissue after thermal coagulation. The excised liver samples were heated in a degassed water bath up to 75°C and naturally cooled down to 30°C. The tissue was observed to begin thermally coagulating at temperatures lower than 75°C. The ultrasonic parameters considered include the speed of sound, the attenuation coefficient, the backscatter coefficient and the nonlinear parameter of B/A. They were more sensitive to temperature when heating than during natural cooling. All of the parameters were shown to rise significantly on completion of the heating-cooling cycle. At 35°C after thermal coagulation, the B/A value was increased by 96%, the attenuation and backscatter coefficients were increased by 50%∼68% and 33%∼37%, respectively, in the typical frequency ranges of 3 MHz∼5 MHz used for ultrasonic imaging and the speed of sound was increased by 1.4%. The results of this study added to the evidence that tissue characterization, in particular, based on the B/A could be valuable for ultrasonically imaging the thermal lesions following high-intensity focused ultrasound (HIFU) surgery.

Temperature dependent ultrasonic characterization of biological media

Quantitative ultrasound (QUS) is an imaging technique that can be used to quantify tissue microstructure giving rise to scattered ultrasound. Other ultrasonic properties, e.g., sound speed and attenuation, of tissues have been estimated versus temperature elevation and found to have a dependence with temperature. Therefore, it is hypothesized that QUS parameters may be sensitive to changes in tissue microstructure due to temperature elevation. Ultrasonic backscatter experiments were performed on tissue-mimicking phantoms and freshly excised rabbit and beef liver samples. The phantoms were made of agar and contained either mouse mammary carcinoma cells (4T1) or chinese hamster ovary cells (CHO) as scatterers. All scatterers were uniformly distributed spatially at random throughout the phantoms. All the samples were scanned using a 20-MHz single-element f/3 transducer. Quantitative ultrasound parameters were estimated from the samples versus increases in temperature from 37 °C to 50 °C in 1 °C increments. Two QUS parameters were estimated from the backscatter coefficient [effective scatterer diameter (ESD) and effective acoustic concentration (EAC)] using a spherical Gaussian scattering model. Significant increases in ESD and decreases in EAC of 20%-40% were observed in the samples over the range of temperatures examined. The results of this study indicate that QUS parameters are sensitive to changes in temperature.

Theoretical and phantom based investigation of the impact of sound speed and backscatter variations on attenuation slope estimation

Quantitative ultrasound features such as the attenuation slope, sound speed and scatterer size, have been utilized to evaluate pathological variations in soft tissues such as the liver and breast. However, the impact of variations in the sound speed and backscatter due to underlying fat content or fibrotic changes, on the attenuation slope has not been addressed. Both numerical and acoustically uniform tissue-mimicking experimental phantoms are used to demonstrate the impact of sound speed variations on attenuation slope using clinical real-time ultrasound scanners equipped with linear array transducers. Radiofrequency data at center frequencies of 4 and 5 MHz are acquired for the experimental and numerical phantoms respectively. Numerical phantom sound speeds between 1480 and 1600 m/s in increments of 20 m/s for attenuation coefficients of 0.3, 0.4, 0.5, 0.6, and 0.7 dB/cm/MHz are simulated. Variations in the attenuation slope when the backscatter intensity of the sample is equal, 3 dB higher, and 3 dB lower than the reference is also evaluated. The sound speed for the experimental tissue-mimicking phantoms were 1500, 1540, 1560 and 1580 m/s respectively, with an attenuation coefficient of 0.5 dB/cm/MHz. Radiofrequency data is processed using three different attenuation estimation algorithms, i.e. the reference phantom, centroid downshift, and a hybrid method. In both numerical and experimental phantoms our results indicate a bias in attenuation slope estimates when the reference phantom sound speed is higher (overestimation) or lower (underestimation) than that of the sample. This bias is introduced via a small spectral shift in the normalized power spectra of the reference and sample with different sound speeds. The hybrid method provides the best estimation performance, especially for sample attenuation coefficient values lower than that of the reference phantom. The performance of all the methods deteriorates when the attenuation coefficient of the reference phantom is lower than that of the sample. In addition, the hybrid method is the least sensitive to sample backscatter intensity variations.

Absolute backscatter coefficient estimates of tissue-mimicking phantoms in the 5-50 MHz frequency range

Absolute backscatter coefficients in tissue-mimicking phantoms were experimentally determined in the 5-50 MHz frequency range using a broadband technique. A focused broadband transducer from a commercial research system, the VisualSonics Vevo 770, was used with two tissue-mimicking phantoms. The phantoms differed regarding the thin layers covering their surfaces to prevent desiccation and regarding glass bead concentrations and diameter distributions. Ultrasound scanning of these phantoms was performed through the thin layer. To avoid signal saturation, the power spectra obtained from the backscattered radio frequency signals were calibrated by using the signal from a liquid planar reflector, a water-brominated hydrocarbon interface with acoustic impedance close to that of water. Experimental values of absolute backscatter coefficients were compared with those predicted by the Faran scattering model over the frequency range 5-50 MHz. The mean percent difference and standard deviation was 54% ± 45% for the phantom with a mean glass bead diameter of 5.40 μm and was 47% ± 28% for the phantom with 5.16 μm mean diameter beads.

Development and characterization of a tissue-mimicking material for high-intensity focused ultrasound

A tissue-mimicking material (TMM) for the acoustic and thermal characterization of high-intensity focused ultrasound (HIFU) devices has been developed. The material is a high-temperature hydrogel matrix (gellan gum) combined with different sizes of aluminum oxide particles and other chemicals. The ultrasonic properties (attenuation coefficient, speed of sound, acoustical impedance, and the thermal conductivity and diffusivity) were characterized as a function of temperature from 20 to 70°C. The backscatter coefficient and nonlinearity parameter B/A were measured at room temperature. Importantly, the attenuation coefficient has essentially linear frequency dependence, as is the case for most mammalian tissues at 37°C. The mean value is 0.64f(0.95) dB·cm(-1) at 20°C, based on measurements from 2 to 8 MHz. Most of the other relevant physical parameters are also close to the reported values, although backscatter signals are low compared with typical human soft tissues. Repeatable and consistent temperature elevations of 40°C were produced under 20-s HIFU exposures in the TMM. This TMM is appropriate for developing standardized dosimetry techniques, validating numerical models, and determining the safety and efficacy of HIFU devices.

Characterization of the spatial resolution of different high-frequency imaging systems using a novel anechoic-sphere phantom

The spatial resolution of high-frequency ultrasound (HFU, >20 MHz) imaging systems is usually determined using wires perpendicular to the beam. Recently, two tissue-mimicking phantoms (TMPs) were developed to estimate three-dimensional (3-D) resolution. Each TMP consists of nine 1-cm-wide slabs of tissue-mimicking material containing randomly distributed anechoic spheres. All anechoic spheres in one slab have the same dimensions, and their diameter is increased from 0.1 mm in the first slab to 1.09 mm in the last. The scattering background for one set of slabs was fabricated using 3.5-μm glass beads; the second set used 6.4-μm glass beads. The ability of a HFU system to detect these spheres against a speckle background provides a realistic estimation of its 3-D spatial resolution. In the present study, these TMPs were used with HFU systems using single-element transducers, linear arrays, and annular arrays. The TMPs were immersed in water and each slab was scanned using two commercial imaging systems and a custom HFU system based on a 5-element annular array. The annular array had a nominal center frequency of 40 MHz, a focal length of 12 mm, and a total aperture of 6 mm. A synthetic-focusing algorithm was used to form images with an increased depth-of-field. The penetration depth was increased by using a linear-chirp signal spanning 15 to 65 MHz over 4 μs. Results obtained with the custom system were compared with those of the commercial systems (40-MHz probes) in terms of sphere detection, i.e., 3-D spatial resolution, and contrast-to-noise ratio (CNR). Resulting B-mode images indicated that only the linear-array transducer failed to clearly resolve the 0.2-mm spheres, which showed that the 3-D spatial resolution of the single-element and annular-array transducers was superior to that of the linear array. The single-element transducer could only detect these spheres over a narrow 1.5 mm depth-of-field, whereas the annular array was able to detect them to depths of at least 7 mm. For any size of the anechoic spheres, the annular array excited by a chirp-coded signal provided images of the highest contrast, with a maximum CNR of 1.8 at the focus, compared with 1.3 when using impulse excitation and 1.6 with the single-element transducer and linear array. This imaging configuration also provided CNRs above 1.2 over a wide depth range of 8 mm, whereas CNRs would quickly drop below 1 outside the focal zone of the other configurations.

On the estimation of backscatter coefficients using single-element focused transducers

The ultimate goal of quantitative ultrasound (QUS) imaging methods based on backscatter coefficient (BSC) estimates is to obtain system-independent structural information about samples. In the current study, three BSC estimation methods were compared and evaluated using the same backscattered pressure datasets in order to assess their consistency. BSC estimates were obtained from two phantoms with embedded glass spheres and compared to theoretical BSCs calculated using size distributions estimated using optical microscopy. Effective scatterer diameter and concentration estimates of the glass spheres were also obtained from the estimated BSCs. One estimation method needed to be compensated by more than an order of magnitude in amplitude in order to produce BSCs comparable to the other two methods. All calibration methods introduced different frequency-dependent effects, which could have noticeable effects on the bias of QUS estimates derived from experimental BSCs. Although in most cases the experimental QUS estimates obtained with all three methods were observed to differ by less than 10%, larger differences are expected depending on both the pressure focusing gain of the transducer (proportional to the ratio of the square of the aperture radius to the product of the wavelength and focal length) and ka range used in the estimation.

Exvivo ultrasound attenuation coefficient for human cervical and uterine tissue from 5 to 10 MHz

Attenuation estimation and imaging in the cervix has been utilized to evaluate the onset of cervical ripening during pregnancy. This feature has also been utilized for the acoustic characterization of leiomyomas and myometrial tissue. In this paper, we present direct narrowband substitution measurement values of the variation in the ultrasonic attenuation coefficient in ex vivo human uterine and cervical tissue, in the 5-10 MHz frequency range. At 5 MHz, the attenuation coefficient values are similar for the different orientations of uterine tissue with values of 4.1-4.2 dB/cm, 5.1 dB/cm for the leiomyoma, and 6.3 dB/cm for the cervix. As the frequency increases, the attenuation coefficient values increase and are also spread out, with a value of approximately 12.6 dB/cm for the uterus (both parallel and perpendicular), 16.0 for the leiomyoma, and 26.8 dB/cm for the cervix at 10 MHz. The attenuation coefficient measured increases monotonically over the frequency range measured following a power law.

Anechoic sphere phantoms for estimating 3-D resolution of very-high-frequency ultrasound scanners

Two phantoms have been constructed for assessing performance of high-frequency ultrasound imagers. They also allow for periodic quality assurance tests and training technicians in the use of higher-frequency scanners. The phantoms contain eight blocks of tissue-mimicking material; each block contains a spatially random distribution of suitably small anechoic spheres having a small distribution of diameters. The eight mean sphere diameters are distributed from 0.10 to 1.09 mm. The two phantoms differ primarily in terms of the frequency dependence of the backscatter coefficient of the background material. Because spheres have no preferred orientation, all three (spatial) dimensions of resolution contribute to sphere detection on an equal basis; thus, the resolution is termed 3-D. Two high-frequency scanners are compared. One employs single-element (fixed focus) transducers (25 and 55 MHz), and the other employs variable focus linear arrays (20, 30, and 40 MHz). The depth range for detection of spheres of each size is determined corresponding to determination of 3-D resolution as a function of depth. As expected, the single-element transducers are severely limited in useful imaging depth ranges compared with the linear arrays. In this preliminary report, only one human observer analyzed images.

Egg white as a blood coagulation surrogate

Egg white, a protein-containing solution, is characterized as a blood coagulation surrogate for the acoustical and thermal evaluation of therapeutic ultrasound, especially high intensity focused ultrasound (HIFU) devices. Physical properties, including coagulation temperature, frequency dependent attenuation, sound speed, viscosity, and thermal properties, were measured as a function of temperature (20-95 degrees C). Thermal coagulation and attenuation (5-12 and 1 MHz) of cow blood, pig blood, and human blood also were assessed and compared with egg white. For a 30 s thermal exposure, both egg white and blood samples (3 mm thickness) started to denature at 65 degrees C and coagulate into an elastic gel at 85 degrees C. The attenuation of egg white was found to be similar to that of the blood samples, having values of 0.23f(1.09), 1.58f(0.61), and 2.7f(0.5) dB/cm at 20, 75, and 95 degrees C, respectively. This significant attenuation increase with temperature was determined to be caused mainly by bubble cavity formation. The other temperature-dependent parameters are also similar to the reported values for blood. These properties make egg white a potentially useful bench testing tool for the safety and efficacy evaluation of therapeutic ultrasound devices.

Tissue-mimicking bladder wall phantoms for evaluating acoustic radiation force-optical coherence elastography systems

PURPOSE

Acoustic radiation force-optical coherence elastography (ARF-OCE) systems are novel imaging systems that have the potential to simultaneously quantify and characterize the optical and mechanical properties of in vivo tissues. This article presents the construction of bladder wall phantoms for use in ARF-OCE systems. Mechanical, acoustic, and optical properties are reported and compared to published values for the urinary bladder.

METHODS

The phantom consisted of 0.2000 +/- 0.0089 and 6.0000 +/- 0.2830 microm polystyrene microspheres (Polysciences Inc., Warrington, PA, Catalog Nos. 07304 and 07312), 7.5 +/- 1.5 microm copolymer microspheres composed of acrylonitrile and vinylidene chloride, (Expancel, Duluth, GA, Catalog No. 461 DU 20), and bovine serum albumin within a gelatin matrix. Young's modulus was measured by successive compression of the phantom and obtaining the slope of the resulting force-displacement data. Acoustic measurements were performed using the transmission method. The phantoms were submerged in a water bath and placed between transmitting and receiving 13 mm diameter unfocused transducers operating at a frequency of 3.5 MHz. A MATLAB algorithm to extract the optical scattering coefficient from optical coherence tomography (OCT) images of the phantom was used.

RESULTS

The phantoms possess a Young's modulus of 17.12 +/- 2.72 kPa, a mass density of 1.05 +/- 0.02 g/cm3, an acoustic attenuation coefficient of 0.66 +/- 0.08 dB/cm/MHz, a speed of sound of 1591 +/- 8.76 m/s, and an optical scattering coefficient of 1.80 +/- 0.23 mm(-1). Ultrasound and OCT images of the bladder wall phantom are presented.

CONCLUSIONS

A material that mimics the mechanical, optical, and acoustic properties of healthy bladder wall has been developed. This tissue-mimicking bladder wall phantom was developed as a control tool to investigate the feasibility of using ARF-OCE to detect the mechanical and optical changes that may be indicative of the onset or development of cancer in the urinary bladder. By following the methods used in this article, phantoms matching the optical, acoustic, and mechanical properties of other biological tissues can also be constructed.

Interlaboratory comparison of backscatter coefficient estimates for tissue-mimicking phantoms

Ultrasonic backscatter is useful for characterizing tissues and several groups have reported methods for estimating backscattering properties. Previous interlaboratory comparisons have been made to test the ability to accurately estimate the backscatter coefficient (BSC) by different laboratories around the world. Results of these comparisons showed variability in BSC estimates but were acquired only for a relatively narrow frequency range, and, most importantly, lacked reference to any independent predictions from scattering theory. The goal of this study was to compare Faran-scattering-theory predictions with cooperatively-measured backscatter coefficients for low-attenuating and tissue-like attenuating phantoms containing glass sphere scatterers of different sizes for which BSCs can independently be predicted. Ultrasonic backscatter measurementswere made for frequencies from 1 to 12 MHz. Backscatter coefficients were estimated using two different planar-reflector techniques at two laboratories for two groups of phantoms. Excellent agreement was observed between BSC estimates from both laboratories. In addition, good agreement with the predictions of Faran's theory was obtained, with average fractional (bias) errors ranging from 8-14%. This interlaboratory comparison demonstrates the ability to accurately estimate parameters derived from the BSC, including an effective scatterer size and the acoustic concentration, both of which may prove useful for diagnostic applications of ultrasound tissue characterization.

A novel coded excitation scheme to improve spatial and contrast resolution of quantitative ultrasound imaging

Quantitative ultrasound (QUS) imaging techniques based on ultrasonic backscatter have been used successfully to diagnose and monitor disease. A method for improving the contrast and axial resolution of QUS parametric images by using the resolution enhancement compression (REC) technique is proposed. Resolution enhancement compression is a coded excitation and pulse compression technique that enhances the -6-dB bandwidth of an ultrasonic imaging system. The objective of this study was to combine REC with QUS (REC-QUS) and evaluate and compare improvements in scatterer diameter estimates obtained using the REC technique to conventional pulsing methods. Simulations and experimental measurements were conducted with a single-element transducer (f/4) having a center frequency of 10 MHz and a -6-dB bandwidth of 80%. Using REC, the -6-dB bandwidth was enhanced to 155%. Images for both simulation and experimental measurements contained a signal-to-noise ratio of 28 dB. In simulations, to monitor the improvements in contrast a software phantom with a cylindrical lesion was evaluated. In experimental measurements, tissue-mimicking phantoms that contained glass spheres with different scatterer diameters were evaluated. Estimates of average scatterer diameter in the simulations and experiments were obtained by comparing the normalized backscattered power spectra to theory over the -6-dB bandwidth for both conventional pulsing and REC. Improvements in REC-QUS over conventional QUS were quantified through estimate bias and standard deviation, contrast-to-noise ratio, and histogram analysis of QUS parametric images. Overall, a 51% increase in contrast and a 60% decrease in the standard deviation of average scatterer diameter estimates were obtained during simulations, while a reduction of 34% to 71% was obtained in the standard deviation of average scatterer diameter for the experimental results.

Noninvasive determination of in situ heating rate using kHz acoustic emissions and focused ultrasound

For high-intensity focused ultrasound (HIFU) to be widely applicable in the clinic, robust methods of treatment planning, guidance and delivery need to be developed. These technologies would greatly benefit if patient specific tissue parameters could be provided as inputs so that the treatment planning and monitoring schemes are customized and tailored on a case by case basis. A noninvasive method of estimating the local in situ acoustic heating rate using the heat transfer equation (HTE) and applying novel signal processing techniques is presented in this article. The heating rate is obtained by experimentally measuring the time required to raise the temperature of the therapeutic focus from a baseline temperature to boiling (here assumed to be 100 degrees C for aqueous media) and then solving the heat transfer equation iteratively to find the heating rate that results in the onset of boiling. The onset of boiling is noninvasively detected by measuring the time instant of onset of acoustic emissions in the audible frequency range due to violent collapse of bubbles. In vitro experiments performed in a tissue mimicking alginate phantom and excised turkey breast muscle tissue demonstrate that the noninvasive estimates of heating rate are in good agreement with those obtained independently using established methods. The results show potential for the applicability of these techniques in therapy planning and monitoring for therapeutic dose optimization using real-time acoustic feedback.

Copolymer-in-oil phantom materials for elastography

Phantoms that mimic mechanical and acoustic properties of soft biological tissues are essential to elasticity imaging investigation and to elastography device characterization. Several materials including agar/gelatin, polyvinyl alcohol and polyacrylamide gels have been used successfully in the past to produce tissue phantoms, as reported in the literature. However, it is difficult to find a phantom material with a wide range of stiffness, good stability over time and high resistance to rupture. We aim at developing and testing a new copolymer-in-oil phantom material for elastography. The phantom is composed of a mixture of copolymer, mineral oil and additives for acoustic scattering. The mechanical properties of phantoms were evaluated with a mechanical test instrument and an ultrasound-based elastography technique. The acoustic properties were investigated using a through-transmission water-substituting method. We showed that copolymer-in-oil phantoms are stable over time. Their mechanical and acoustic properties mimic those of most soft tissues: the Young's modulus ranges from 2.2-150 kPa, the attenuation coefficient from 0.4-4.0 dB.cm(-1) and the ultrasound speed from 1420-1464 m/s. Their density is equal to 0.90 +/- 0.04 g/cm3. The results suggest that copolymer-in-oil phantoms are attractive materials for elastography.

Monitoring of cell death in epithelial cells using high frequency ultrasound spectroscopy

Spectral and wavelet analyses were performed on ultrasound radiofrequency (RF) data collected from centrifuged cell samples containing HEp-2 cells after induction of apoptosis by exposure to camptothecin. Samples were imaged at several time points after drug exposure using high-frequency ultrasound in the range from 10-60 MHz. A 20-MHz transducer with a f-number of 2.35 and a 40-MHz transducer with a f-number of 3 were used for collecting the RF data. Normalized power spectra were computed from the backscattered ultrasound signals within a region-of-interest (ROI) for further analysis. Spectral slopes, integrated backscatter coefficients (IBCs) and wavelet parameters were estimated as a function of treatment time to monitor acoustic property changes during apoptosis. Changes in spectral parameters were detected starting six hours after treatment and coincided with changes in corresponding histology. Throughout the course of chemotherapy, variation in estimates of the spectral slope of up to 35% were observed. During the treatment, IBCs increased by 400% compared with estimates obtained from the control samples. Changes in spectral parameters are hypothesized to be linked to structural cell changes during apoptosis. In addition, the sensitivity of a wavelet-based analysis to the ultrasonic assessment of cellular changes was investigated. Results of the wavelet analysis showed variations similar to the spectral parameters. Where values of the spectral slope decreased, estimates of the scaling factors increased. Because wavelet analysis preserves the signal-time localization, its application will be potentially beneficial for assessing treatment responses in vivo. The current study contributes toward the development of a non-invasive method for monitoring apoptosis as a measure of the success of chemotherapeutic treatment of cancer.

Anthropomorphic phantoms for assessment of strain imaging methods involving saline-infused sonohysterography

Two anthropomorphic uterine phantoms were developed that allow assessment and comparison of strain imaging systems adapted for use with saline-infused sonohysterography (SIS). Tissue-mimicking (TM) materials consist of dispersions of safflower oil in gelatin. TM fibroids are stiffer than the TM myometrium/cervix, and TM polyps are softer. The first uterine phantom has 3-mm-diameter TM fibroids distributed randomly in TM myometrium. The second uterine phantom has a 5-mm and 8-mm spherical TM fibroid, in addition to a 5-mm spherical and a 12.5-mm-long (medicine capsule-shaped) TM endometrial polyp protruding into the endometrial cavity; also, a 10-mm spherical TM fibroid projects from the serosal surface. Strain images using the first phantom show the stiffer 3-mm TM fibroids in the myometrium. Results from the second uterine phantom show that, as expected, parts of inclusions projecting into the uterine cavity will appear very stiff, whether they are stiff or soft. Results from both phantoms show that although there is a five-fold difference in the Young's moduli values, there is not a significant difference in the strain in the transition from the TM myometrium to the TM fat. These phantoms allow for realistic comparison and evolution of SIS strain imaging techniques and can aid clinical personnel to develop skills for SIS strain imaging.

Development and characterization of a blood mimicking fluid for high intensity focused ultrasound

A blood mimicking fluid (BMF) has been developed for the acoustic and thermal characterizations of high intensity focused ultrasound (HIFU) ablation devices. The BMF is based on a degassed and de-ionized water solution dispersed with low density polyethylene microspheres, nylon particles, gellan gum, and glycerol. A broad range of physical parameters, including attenuation coefficient, speed of sound, viscosity, thermal conductivity, and diffusivity, were characterized as a function of temperature (20-70 degrees C). The nonlinear parameter B/A and backscatter coefficient were also measured at room temperature. Importantly, the attenuation coefficient is linearly proportional to the frequency (2-8 MHz) with a slope of about 0.2 dB cm(-1) MHz(-1) in the 20-70 degrees C range as in the case of human blood. Furthermore, sound speed and bloodlike backscattering indicate the usefulness of the BMF for ultrasound flow imaging and ultrasound-guided HIFU applications. Most of the other temperature-dependent physical parameters are also close to the reported values in human blood. These properties make it a unique HIFU research tool for developing standardized exposimetry techniques, validating numerical models, and determining the safety and efficacy of HIFU ablation devices.

Noninvasive measurement of local thermal diffusivity using backscattered ultrasound and focused ultrasound heating

Previously, noninvasive methods of estimating local tissue thermal and acoustic properties using backscattered ultrasound have been proposed in the literature. In this article, a noninvasive method of estimating local thermal diffusivity in situ during focused ultrasound heating using beamformed acoustic backscatter data and applying novel signal processing techniques is developed. A high intensity focused ultrasound (HIFU) transducer operating at subablative intensities is employed to create a brief local temperature rise of no more than 10 degrees C. Beamformed radio-frequency (RF) data are collected during heating and cooling using a clinical ultrasound scanner. Measurements of the time-varying "acoustic strain", that is, spatiotemporal variations in the RF echo shifts induced by the temperature related sound speed changes, are related to a solution of the heat transfer equation to estimate the thermal diffusivity in the heated zone. Numerical simulations and experiments performed in vitro in tissue mimicking phantoms and excised turkey breast muscle tissue demonstrate agreement between the ultrasound derived thermal diffusivity estimates and independent estimates made by a traditional hot-wire technique. The new noninvasive ultrasonic method has potential applications in thermal therapy planning and monitoring, physiological monitoring and as a means of noninvasive tissue characterization.

High frequency ultrasound tissue characterization and acoustic microscopy of intracellular changes

The objective of this work is to investigate changes in the acoustic properties of cells when exposed to chemotherapy for monitoring treatment response. High frequency ultrasound spectroscopy (10-60 MHz) and scanning acoustic microscopy (0.9 GHz) were performed on HeLa cells (Ackermann et al. 1954, Masters 2002) that were exposed to the chemotherapeutic agent cisplatin. Ultrasonic radio-frequency data were acquired from pellets containing HeLa cells after exposure to cisplatin to induce apoptosis. Scanning acoustic and laser fluorescence microscopy images were recorded from single HeLa cells exposed to the same drug. Data acquisition in both cases was performed at several time points throughout the chemotherapeutic treatment for up to 27 h. In the high frequency ultrasound investigation, normalized power spectra were calculated within a region-of-interest. A 20 MHz transducer (f-number 2.35) and a 40 MHz transducer (f-number 3) were used for the data collection in the high frequency ultrasound experiments. The backscatter coefficients, integrated backscatter coefficients, mid-band fit and spectral slope were computed as a function of treatment time to monitor acoustical property changes during apoptosis. Acoustic attenuation was measured using the spectral substitution technique at all time points. Spectral parameter changes were detected after 12 h of exposure and coincided with the initiation of cell damage as assessed by optical microscopy. Integrated backscatter coefficients increased by over 100% between 0 h and 24 h of treatment, with small changes in the associated attenuation ( approximately 0.1 dB/[MHz cm]). Acoustic microscopy was performed at 0.9 GHz frequency. The cell structure was imaged using staining in laser fluorescence microscopy. All cells showed excellent correspondence between the locations of apoptotic nuclear condensation observed in optical imaging and changes in attenuation contrast in acoustic microscopy images. The time after drug exposure at which such changes occurred in the optical images were coincident with the time of changes detected in the acoustic microscopy images and the high frequency ultrasound experiments.

The acoustic properties, centered on 20 MHZ, of an IEC agar-based tissue-mimicking material and its temperature, frequency and age dependence

The purpose of this study was to characterize the ultrasonic properties of agar-based tissue-mimicking materials (TMMs) at ultrasound frequencies centered around 20 MHz. The TMM acoustic properties measured are the amplitude attenuation coefficient alpha (dB cm(-1)MHz(-1)), the speed of sound (ms(-1)) and the backscattered power spectral density (distribution of power per unit frequency normalized to the total received power) characteristics of spectral slope (dB MHz(-1)), y-axis intercept (dB) and reflected power (dB). The acoustic properties are measured over a temperature range of 22 to 37 degrees C. An intercomparison of results between two independent ultrasound measurement laboratories is also presented. A longitudinal study of the acoustic properties over a period of two years is also detailed, and the effect of water immersion on the acoustic properties of TMM is measured. In addition, the physical parameters of mass density rho (kg m(-3)) and specific heat capacity C (J kg(-1) K(-1)) are included. The measurement techniques used were based on the substitution technique using both broadband and narrowband pulses centered on 20 MHz. Both the attenuation coefficient and speed of sound (both group and phase) showed good agreement with the expected values of 0.5 dB cm(-1) MHz(-1) and 1540 ms(-1), respectively, with average values over the three-year period of 0.49 dBcm(-1)MHz1 (SD +/- 0.05) and 1540.9 ms(-1) (SD +/- 8.7). These results also showed agreement between the two independent measurement laboratories. Speed of sound and attenuation coefficient were shown to change with temperature with rates of + 2.1 m s(-1) degrees C(-1) and -0.005 dB cm(-1) MHz(-1) degrees C(-1), respectively. Attenuation changed linearly with frequency at the high frequency range of 17 to 23 MHz, and speed of sound was found to be independent of frequency in this range. The spectral slope of relative backscattered power for the material increased with frequency at typically 1.5 dB MHz(-1). This compared favorably with theoretical spectral slope values, calculated for a variety of scatterer sizes, albeit at a lower frequency range. It is also noticed that, on extrapolation back to lower frequencies, the backscatter is comparable with that measured at 7 MHz. Overall, this non-commercial agar-based TMM is shown to perform as expected at the higher frequency range of 17 to 23 MHz and is seen to retain its acoustic properties of attenuation and speed of sound over a three-year period.

Research interface on a programmable ultrasound scanner

MOTIVATION

Commercial ultrasound machines in the past did not provide the ultrasound researchers access to raw ultrasound data. Lack of this ability has impeded evaluation and clinical testing of novel ultrasound algorithms and applications.

OBJECTIVES

Recently, we developed a flexible ultrasound back-end where all the processing for the conventional ultrasound modes, such as B, M, color flow and spectral Doppler, was performed in software. The back-end has been incorporated into a commercial ultrasound machine, the Hitachi HiVision 5500. The goal of this work is to develop an ultrasound research interface on the back-end for acquiring raw ultrasound data from the machine.

METHODS

The research interface has been designed as a software module on the ultrasound back-end. To increase the amount of raw ultrasound data that can be spooled in the limited memory available on the back-end, we have developed a method that can losslessly compress the ultrasound data in real time.

RESULTS AND DISCUSSION

The raw ultrasound data could be obtained in any conventional ultrasound mode, including duplex and triplex modes. Furthermore, use of the research interface does not decrease the frame rate or otherwise affect the clinical usability of the machine. The lossless compression of the ultrasound data in real time can increase the amount of data spooled by approximately 2.3 times, thus allowing more than 6s of raw ultrasound data to be acquired in all the modes. The interface has been used not only for early testing of new ideas with in vitro data from phantoms, but also for acquiring in vivo data for fine-tuning ultrasound applications and conducting clinical studies. We present several examples of how newer ultrasound applications, such as elastography, vibration imaging and 3D imaging, have benefited from this research interface. Since the research interface is entirely implemented in software, it can be deployed on existing HiVision 5500 ultrasound machines and may be easily upgraded in the future.

CONCLUSIONS

The developed research interface can aid researchers in the rapid testing and clinical evaluation of new ultrasound algorithms and applications. Additionally, we believe that our approach would be applicable to designing research interfaces on other ultrasound machines.

Improved scatterer size estimation using backscatter coefficient measurements with coded excitation and pulse compression

Scatterer size estimates from ultrasonic backscatter coefficient measurements have been used to differentiate diseased tissue from normal. A low echo signal-to-noise ratio (eSNR) leads to increased bias and variance in scatterer size estimates. One way to improve the eSNR is to use coded excitation (CE). The normalized backscatter coefficient was measured from three tissue-mimicking phantoms by using CE and conventional pulsing (CP) techniques. The three phantoms contained randomly spaced glass beads with median diameters of 30, 45, and 82 mum, respectively. Measurements were made with two weakly focused, single-element transducers (f(0)=5 MHz and f(0)=10 MHz). For CE, a linear frequency modulated chirp with a time bandwidth product of 40 was used and pulse compression was accomplished by the use of a Wiener filter. Preliminary results indicated that improved estimation bias versus penetration depth was obtained by using CE compared to CP. The depth of penetration, where the accuracy of scatterer diameter estimates (absolute divergence <25%) were obtained with the 10 MHz transducer, was increased up to 50% by using CE versus CP techniques. In addition, for a majority of the phantoms, the increase in eSNR from CE resulted in a modest reduction in estimate variance versus depth of penetration.

In vivo ultrasonic attenuation slope estimates for detecting cervical ripening in rats: Preliminary results

To effectively postpone preterm birth, cervical ripening needs to be detected and delayed. As the cervix ripens, the spacing between the collagen fibers increases and fills with water, hyaluronan, decorin, and enzymes suggesting that the ultrasonic attenuation of the cervix should decrease. The decrease in ultrasonic attenuation may be detectable, leading to an effective means of detecting cervical ripening. Herein, the traditional attenuation slope-estimation algorithm based on measuring the downshift in center frequency of the ultrasonic backscattered signal with propagation depth was modified and applied to the cervix of rats. The modified algorithm was verified using computer simulations and an ex vivo tissue sample before being evaluated in in vivo animal studies. Spherically-focused f/3 transducers with 33-MHz center frequencies and with 9-mm focal lengths were used in both the simulations and experiments. The accuracy was better than 15% in the simulations, and the attenuation slope of the cervix in the ex vivo experiment was 2.6+/-0.6 dB/cm-MHz, which is comparable to 2.5+/-0.4 dB/cm-MHz measured using a through-transmission insertion loss technique. For the in vivo experiments, a statistically significant effect of ultrasonic attenuation with gestational age was not observed. The large variances in the in vivo results were most likely due to the natural variation in attenuation for biological tissue between animals.

Non-invasive ultrasound-based temperature imaging for monitoring radiofrequency heating—phantom results

Minimally invasive therapies (such as radiofrequency ablation) are becoming more commonly used in the United States for the treatment of hepatocellular carcinomas and liver metastases. Unfortunately, these procedures suffer from high recurrence rates of hepatocellular carcinoma ( approximately 34-55%) or metastases following ablation therapy. The ability to perform real-time temperature imaging while a patient is undergoing radiofrequency ablation could provide a significant reduction in these recurrence rates. In this paper, we demonstrate the feasibility of ultrasound-based temperature imaging on a tissue-mimicking phantom undergoing radiofrequency heating. Ultrasound echo signals undergo time shifts with increasing temperature, which are tracked using 2D correlation-based speckle tracking methods. Time shifts or displacements in the echo signal are accumulated, and the gradient of these time shifts are related to changes in the temperature of the tissue-mimicking phantom material using a calibration curve generated from experimental data. A tissue-mimicking phantom was developed that can undergo repeated radiofrequency heating procedures. Both sound speed and thermal expansion changes of the tissue-mimicking material were measured experimentally and utilized to generate the calibration curve relating temperature to the displacement gradient. Temperature maps were obtained, and specific regions-of-interest on the temperature maps were compared to invasive temperatures obtained using fiber-optic temperature probes at the same location. Temperature elevation during a radiofrequency ablation procedure on the phantom was successfully tracked to within +/-0.5 degrees C.

Anisotropy of the backscatter coefficient of formalin-fixed ovine myocardium

The objective of this study was to measure the backscatter coefficient of formalin-fixed myocardial tissue as a function of angle of insonification relative to the myocardial fiber direction. Backscatter measurements were performed on eight cylindrical formalin-fixed lamb myocardial specimens and compensated for attenuation and diffraction effects to determine the backscatter coefficient. The backscatter coefficient at 5 MHz was found to be maximum for insonification perpendicular to the predominant myofiber orientation and minimum for parallel insonification, with values of (17+/-14) and (1.2+/-0.7) x 10(-4) cm(-1) sr(-1) (mean+/-standard deviation), respectively.

The Ultrasonix 500RP: a commercial ultrasound research interface

Unlike researchers in magnetic resonance imaging who have considerable access to high level tools and to data at a very basic level on their scanners, those involved with ultrasound have found little in the way of meaningful and widespread access to even the most basic echo signals in their clinical systems. Interest has emerged, however, in ultrasound research interfaces on commercial scanners to provide access to raw ultrasound data and control of basic research functions. This paper describes initial experience gained on one such ultrasound system. The Ultrasonix 500RP system provides research access to the data at multiple points in the signal processing chain and allows control over most imaging parameters. The Ultrasonix system allows for three methods of research control. One is implemented along with the standard clinical imaging software using "mouseover" screens on the periphery of the application window. These screens are configured by the user to display various signal processing variables, which can be modified in real time. Second, the system can be controlled via a user-written remote control client application interacting through the clinical exam software. Lastly, the user can write a complete application which initializes the basic ultrasound module but need not use the Ultrasonix clinical exam software. All of the modes can be done locally on the scanner itself or via a network, and are based on software developed in C++ with libraries supplied with the scanner. Two examples are presented in this paper from the evaluation of the system in "real world" applications. Measurements of absolute backscatter coefficients and attenuation coefficients versus frequency are shown and elastograms utilizing spatial compounding are described.

An Automated Measurement System of Ultrasonic Properties of Tofu and Acorn Curd (Dotori Muk)

Tofu was suggested as Tissue Mimicking Materials (TMMs) due to its low cost, availability, gelatinous nature, and shape flexibility. However, the similarity in acoustic properties of the different packets of Tofu and acorn curd (Muk) is critical for TMMs but not fully investigated yet. In order to improve the accuracy for time-consuming and multi-stepped measurements of acoustic properties, an automated measurement system was developed. Using the automated system, the acoustic properties Tofu and Muk were measured to know whether they are similar among the different samples. In conclusion, the acoustic properties were measured to be similar among the different packets of Tofu and Muk, suggesting them as potential TMMs.

Anthropomorphic Breast Phantoms for Qualification of Investigators for ACRIN Protocol 6666

The purpose of this study was to evaluate various ultrasonic properties of breast phantoms developed for use in qualifying investigators for participation in the American College of Radiology Imaging Network (ACRIN) protocol 6666, "Screening Breast Ultrasound in High-Risk Women." Specifically, a tool was sought to consistently measure the performance of radiology personnel in detecting and characterizing lesions similar to those expected with screening breast ultrasonography (US). The phantoms are equivalent to one another except for the randomization of positions of 14 of the 17 simulated lesions. The lesions differ in depth and ultrasonic properties. Representative values of propagation speed, attenuation, relative echogenicity, and mass density are reported for all tissue-mimicking components. Beam refraction occurs at the interface between the subcutaneous fat layer and the glandular parenchyma and can result in beam distortion artifacts similar to those encountered in clinical breast US.

Characterizing acoustic attenuation of homogeneous media using focused impulsive acoustic radiation force

A new method to characterize a material's attenuation using acoustic radiation force is proposed. Comparison of displacement magnitudes generated in a homogeneous material by acoustic radiation force excitations can be used to estimate the material's attenuation when the excitations are applied over a range of focal depths while maintaining a constant lateral focal configuration. Acoustic attenuations are related to the inverse of the excitation focal depth that yields the greatest focal zone displacement for this protocol. Experimental studies in calibrated tissue-mimicking phantoms are presented to demonstrate the feasibility of this method. Attenuations ranging from 0.3-1.5 dB/cm/MHz were characterized over excitation focal depths ranging from 5-30 mm, with an accuracy of 0.1 +/- 0.15 dB/cm/MHz. As currently implemented, this method is limited to characterizing materials that have homogeneous material properties and acoustic attenuations. This method for characterizing acoustic attenuation can be performed using conventional diagnostic scanners without any additional hardware and could also be performed concurrently with acoustic radiation force-based imaging modalities to generate images of mechanical properties and attenuation that are spatially co-registered with B-mode images.

Stability of heterogeneous elastography phantoms made from oil dispersions in aqueous gels

A set of five tissue-mimicking phantoms with cylindrical inclusions were produced for assessing long-term stability of geometry and elastic properties and assessing accuracy of determination of elastic properties. The base aqueous materials were either gelatin or a mixture of agar and gelatin. Stiffness was controlled by selection of the volume percent consisting of microscopic safflower oil droplets. Cylinder diameters remained unchanged within 1% or 2% over many months. Strain ratios from elastograms of the phantoms were stable over many months, implying that elastic contrasts were also stable. Test samples, called production samples, for measurement of Young's moduli were made at the time of manufacture of each phantom and were stored separately from one another. Each production sample was homogeneous and consisted of either inclusion material or background material. For all five phantoms, it was found that the elastic contrast computed using Young's modulus values determined using the production samples accurately represented the true elastic contrasts in the corresponding phantom. This finding was established by the fact that the (true) elastic contrasts determined using samples excised from the phantoms themselves agreed with the elastic contrasts obtained using the homogeneous production samples.