An investigation of the use of transmission ultrasound to measure acoustic attenuation changes in thermal therapy

Ultrasonic Attenuation of an Agar, Silicon Dioxide, and Evaporated Milk Gel Phantom

BACKGROUND

It has been demonstrated that agar-based gel phantoms can emulate the acoustic parameters of real tissues and are the most commonly used tissue-mimicking materials for high-intensity focused ultrasound applications. The following study presents ultrasonic attenuation measurements of agar-based phantoms with different concentrations of additives (percent of agar, silicon dioxide and evaporated milk) in an effort of matching the material's acoustic property as close as possible to human tissues.

METHODS

Nine different agar-based phantoms with various amounts of agar, silicon dioxide, and evaporated milk were prepared. Attenuation measurements of the samples were conducted using the through-transmission immersion techniques.

RESULTS

The ultrasonic attenuation coefficient of the agar-based phantoms varied in the range of 0.30-1.49 dB/cm-MHz. The attenuation was found to increase in proportion to the concentration of agar and evaporated milk. Silicon dioxide was found to significantly contribute to the attenuation coefficient up to 4% weight to volume (w/v) concentration.

CONCLUSION

The acoustic attenuation coefficient of agar-based phantoms can be adjusted according to the tissue of interest in the range of animal and human tissues by the proper selection of agar, silicon dioxide, and evaporated milk.

Estimating the ultrasound attenuation coefficient using complementary Golay codes

Accurate evaluation of ultrasonic wave attenuation is important in many medical applications of ultrasound. The aim of this work is to present a thorough analysis of the effectiveness of using complementary Golay coded sequences (CGCS) during the evaluation of ultrasound attenuation in tissue-like materials, especially at greater depths or at high attenuation. In order to compare the results of the attenuation measurement with the use of CGCS transmission and a short two sine cycles pulse, ultrasound backscattered from medium with predefined attenuation of 0.3, 0.7 and 2 dB/[MHz × cm] were simulated. Also for the same transmission signals, measurements of ultrasound echoes scattered in the tissue phantom with an attenuation of 0.5 dB/[MHz × cm] were performed. In the case of numerically simulated data, for the CGCS excitation, the maximum depth for which the attenuation was correctly determined increased from 55 mm to 80 mm for the 0.7 dB/[MHz × cm] phantom and from 20 mm to 50 mm for the 2 dB/[MHz × cm] phantom compared to excitation of the transducer with a short two sine cycles pulse. When the measurement data obtained using the tissue phantom was used to estimate the attenuation coefficient, the relative error was determined to be 6% and 16% for the depths of 10 mm and 40 mm for the short two sine cycles pulse excitation, respectively. Corresponding values for CGCS excitation and considered depths were 2% and 4%. The use of CGSC sequence during attenuation measurements increases measurement accuracy and can improve medical diagnostic techniques.

Radiation-force-based estimation of acoustic attenuation using harmonic motion imaging (HMI) in phantoms and in vitro livers before and after HIFU ablation

Acoustic attenuation represents the energy loss of the propagating wave through biological tissues and plays a significant role in both therapeutic and diagnostic ultrasound applications. Estimation of acoustic attenuation remains challenging but critical for tissue characterization. In this study, an attenuation estimation approach was developed using the radiation-force-based method of harmonic motion imaging (HMI). 2D tissue displacement maps were acquired by moving the transducer in a raster-scan format. A linear regression model was applied on the logarithm of the HMI displacements at different depths in order to estimate the acoustic attenuation. Commercially available phantoms with known attenuations (n = 5) and in vitro canine livers (n = 3) were tested, as well as HIFU lesions in in vitro canine livers (n = 5). Results demonstrated that attenuations obtained from the phantoms showed a good correlation (R² = 0.976) with the independently obtained values reported by the manufacturer with an estimation error (compared to the values independently measured) varying within the range of 15-35%. The estimated attenuation in the in vitro canine livers was equal to 0.32   ±   0.03 dB cm(-1) MHz(-1), which is in good agreement with the existing literature. The attenuation in HIFU lesions was found to be higher (0.58   ±   0.06 dB cm(-1) MHz(-1)) than that in normal tissues, also in agreement with the results from previous publications. Future potential applications of the proposed method include estimation of attenuation in pathological tissues before and after thermal ablation.

Scatterer number density considerations in reference phantom-based attenuation estimation

Attenuation estimation and imaging have the potential to be a valuable tool for tissue characterization, particularly for indicating the extent of thermal ablation therapy in the liver. Often the performance of attenuation estimation algorithms is characterized with numerical simulations or tissue-mimicking phantoms containing a high scatterer number density (SND). This ensures an ultrasound signal with a Rayleigh distributed envelope and a signal-to-noise ratio (SNR) approaching 1.91. However, biological tissue often fails to exhibit Rayleigh scattering statistics. For example, across 1647 regions of interest in five ex vivo bovine livers, we obtained an envelope SNR of 1.10 ± 0.12 when the tissue was imaged with the VFX 9L4 linear array transducer at a center frequency of 6.0 MHz on a Siemens S2000 scanner. In this article, we examine attenuation estimation in numerical phantoms, tissue-mimicking phantoms with variable SNDs and ex vivo bovine liver before and after thermal coagulation. We find that reference phantom-based attenuation estimation is robust to small deviations from Rayleigh statistics. However, in tissue with low SNDs, large deviations in envelope SNR from 1.91 lead to subsequently large increases in attenuation estimation variance. At the same time, low SND is not found to be a significant source of bias in the attenuation estimate. For example, we find that the standard deviation of attenuation slope estimates increases from 0.07 to 0.25 dB/cm-MHz as the envelope SNR decreases from 1.78 to 1.01 when estimating attenuation slope in tissue-mimicking phantoms with a large estimation kernel size (16 mm axially × 15 mm laterally). Meanwhile, the bias in the attenuation slope estimates is found to be negligible (<0.01 dB/cm-MHz). We also compare results obtained with reference phantom-based attenuation estimates in ex vivo bovine liver and thermally coagulated bovine liver.

Mean scatterer spacing estimation in normal and thermally coagulated ex vivo bovine liver

The liver has been hypothesized to have a unique arrangement of microvasculature that presents as an arrangement of quasiperiodic scatterers to an interrogating ultrasound pulse. The mean scatterer spacing (MSS) of these quasiperiodic scatterers has been proposed as a useful quantitative ultrasound biomarker for characterizing liver tissue. Thermal ablation is an increasingly popular method for treating hepatic tumors, and ultrasonic imaging approaches for delineating the extent of thermal ablation are in high demand. In this work, we examine the distribution of estimated MSS in thermally coagulated bovine liver and normal untreated bovine liver ex vivo. We estimate MSS by detecting local maxima in the spectral coherence function of radio frequency echoes from a clinical transducer, the Siemens VFX 9L4 transducer operating on an S2000 scanner. We find that normal untreated bovine liver was characterized by an MSS of approximately 1.3 mm. We examined regions of interest 12 mm wide laterally, and ranging from 12 mm to 18 mm axially, in 2 mm increments. Over these parameters, the mode of the MSS estimates was between 1.25 and 1.37 mm. On the other hand, estimation of MSS in thermally coagulated liver tissue yields a distribution of MSS estimates whose mode varied between 0.45 and 1.0 mm when examining regions of interest over the same sizes. We demonstrate that the estimated MSS in thermally coagulated liver favors small spacings because the randomly positioned scatterers in this tissue are better modeled as aperiodic scatterers. The submillimeter spacings result from the fact that this was the most probable spacing to be estimated if the discretely sampled spectral coherence function was a uniformly random two-dimensional function.

Optimization of the algorithms for estimating the ultrasonic attenuation along the propagation path

In this study, we perform statistical analysis on two methods used to estimate the total ultrasound attenuation along the propagation path from the surface of the transducer to a region of interest at a particular depth; namely, the spectral-fit method and the multiple-filter method. We derive mathematical equations for the bias and variance in the attenuation estimates as a function of region of interest (ROI) size, imaging system bandwidth, and number of independent Gaussian filters (for the multiple filter method). We use numerical simulations to validate the mathematical equations and compare the two algorithms. The results show that the variance in the total attenuation coefficient estimates obtained with the two methods are comparable, and that the estimates are unbiased. For the multiple filter method, the optimal number of Gaussian filters is two.

Investigation of power and frequency for 3D conformal MRI-controlled transurethral ultrasound therapy with a dual frequency multi-element transducer

Transurethral ultrasound therapy uses real-time magnetic resonance (MR) temperature feedback to enable the 3D control of thermal therapy accurately in a region within the prostate. Previous canine studies showed the feasibility of this method in vivo. The aim of this study was to reduce the procedure time, while maintaining targeting accuracy, by investigating new combinations of treatment parameters. Simulations and validation experiments in gel phantoms were used, with a collection of nine 3D realistic target prostate boundaries obtained from previous preclinical studies, where multi-slice MR images were acquired with the transurethral device in place. Acoustic power and rotation rate were varied based on temperature feedback at the prostate boundary. Maximum acoustic power and rotation rate were optimised interdependently, as a function of prostate radius and transducer operating frequency. The concept of dual frequency transducers was studied, using the fundamental frequency or the third harmonic component depending on the prostate radius. Numerical modelling enabled assessment of the effects of several acoustic parameters on treatment outcomes. The range of treatable prostate radii extended with increasing power, and tended to narrow with decreasing frequency. Reducing the frequency from 8 MHz to 4 MHz or increasing the surface acoustic power from 10 to 20 W/cm(2) led to treatment times shorter by up to 50% under appropriate conditions. A dual frequency configuration of 4/12 MHz with 20 W/cm(2) ultrasound intensity exposure can treat entire prostates up to 40 cm(3) in volume within 30 min. The interdependence between power and frequency may, however, require integrating multi-parametric functions in the controller for future optimisations.

Estimating the total ultrasound attenuation along the propagation path by using a reference phantom

In this study, an algorithm previously developed for estimating the total ultrasonic attenuation along the propagation path from the surface of the transducer to a region of interest (ROI) in tissue, was modified to make it more practical for use in clinical settings. Specifically, the algorithm was re-derived for when a tissue mimicking phantom rather than a planar reflector is used to obtain the reference power spectrum. The reference power spectrum is needed to compensate for the transfer function of the transmitted pulse, the transfer function of transducer, and the diffraction effects that result from focusing/beam forming. The modified algorithm was tested on simulated radio frequency (RF) echo lines obtained from two samples that have different scatterer sizes and different attenuation coefficient slopes, one of which was used as a reference. The mean and standard deviation of the percent errors in the attenuation coefficient estimates (ACEs) were less than 5% and 10%, respectively, for ROIs that contain more than 10 pulse lengths and more than 25 independent echo lines. The proposed algorithm was also tested on two tissue mimicking phantoms that have attenuation coefficient slopes of 0.7 dB/cm-MHz and 0.5 dB/cm-MHz respectively, the latter being the reference phantom. When a single element spherically focused source was used, the mean and standard deviation of the percent errors in the ACEs were less than 5% and 10% respectively for windows that contain more than 10 pulse lengths and more than 17 independent echo lines. When a clinical array transducer was used, the mean and standard deviation of the percent errors in the ACEs were less than 5% and 25%, respectively, for windows that contain more than 12 pulse lengths and more than 45 independent echo lines.

Real-time full-field photoacoustic imaging using an ultrasonic camera

A photoacoustic imaging system that incorporates a commercial ultrasonic camera for real-time imaging of two-dimensional (2-D) projection planes in tissue at video rate (30 Hz) is presented. The system uses a Q-switched frequency-doubled Nd:YAG pulsed laser for photoacoustic generation. The ultrasonic camera consists of a 2-D 12 x 12 mm CCD chip with 120 x 120 piezoelectric sensing elements used for detecting the photoacoustic pressure distribution radiated from the target. An ultrasonic lens system is placed in front of the chip to collect the incoming photoacoustic waves, providing the ability for focusing and imaging at different depths. Compared with other existing photoacoustic imaging techniques, the camera-based system is attractive because it is relatively inexpensive and compact, and it can be tailored for real-time clinical imaging applications. Experimental results detailing the real-time photoacoustic imaging of rubber strings and buried absorbing targets in chicken breast tissue are presented, and the spatial resolution of the system is quantified.