Nakagami imaging for detecting thermal lesions induced by high-intensity focused ultrasound in tissue

Therapeutic ultrasound transducer technology and monitoring techniques: a review with clinical examples

The exponential growth of therapeutic ultrasound applications demonstrates the power of the technology to leverage the combinations of transducer technology and treatment monitoring techniques to effectively control the preferred bioeffect to elicit the desired clinical effect.Objective: This review provides an overview of the most commonly used bioeffects in therapeutic ultrasound and describes existing transducer technologies and monitoring techniques to ensure treatment safety and efficacy.Methods and materials: Literature reviews were conducted to identify key choices that essential in terms of transducer design, treatment parameters and procedure monitoring for therapeutic ultrasound applications. Effective combinations of these options are illustrated through descriptions of several clinical indications, including uterine fibroids, prostate disease, liver cancer, and brain cancer, that have been successful in leveraging therapeutic ultrasound to provide effective patient treatments.Results: Despite technological constraints, there are multiple ways to achieve a desired bioeffect with therapeutic ultrasound in a target tissue. Visualizations of the interplay of monitoring modality, bioeffect, and applied acoustic parameters are presented that demonstrate the interconnectedness of the field of therapeutic ultrasound. While the clinical indications explored in this review are at different points in the clinical evaluation path, based on the ever expanding research being conducted in preclinical realms, it is clear that additional clinical applications of therapeutic ultrasound that utilize a myriad of bioeffects will continue to grow and improve in the coming years.Conclusions: Therapeutic ultrasound will continue to improve in the next decades as the combination of transducer technology and treatment monitoring techniques will continue to evolve and be translated in clinical settings, leading to more personalized and efficient therapeutic ultrasound mediated therapies.

Decorrelated compounding of synthetic aperture ultrasound imaging to detect low contrast thermal lesions induced by focused ultrasound

PURPOSE

Decorrelated Compounding (DC) for synthetic aperture ultrasound can reduce speckle variation in images, suggesting enhanced detectability of low-contrast targets in tissue including thermal lesions produced by focused ultrasound (FUS). The DC imaging method has primarily been investigated in simulation and in phantom studies. This work investigates the feasibility of the DC method in monitoring thermal therapy via image guidance and non-invasive thermometry based on the change in backscattered energy (CBE).

METHODS

Ex vivo porcine tissue was exposed to FUS exposures at acoustic powers of 5 W and 1 W, with peak pressure amplitudes of 0.64 MPa and 0.27 MPa respectively. During FUS exposure, RF echo data frames was acquired using a 7.8 MHz linear array probe and a Verasonics Vantage^TM ultrasound scanner (Verasonics Inc., Redmond, WA). RF echo data was taken to produce B-mode images, as reference images. Synthetic aperture RF echo data was also acquired and processed using delay-and-sum (DAS), a combination of spatial and frequency compounding referred to as Traditional Compounding (TC), and the proposed DC imaging methods. Image quality was assessed using the contrast-to-noise ratio (CNR) at the FUS beam focus, and the speckle SNR (sSNR) of the background region as preliminary metrics. A calibrated thermocouple was placed near the FUS beam focus for temperature measurements and calibrations using the CBE method.

RESULTS

The DC imaging method significantly improved image quality to detect low contrast thermal lesions in treated ex vivo porcine tissue in comparison to other imaging methods. In comparison to B-mode imaging, the lesion CNR measured using the DC imaging was shown to improve up to a factor of approximately 5.5. The corresponding sSNR improved by a factor of approximately 4.2 in comparison to B-mode imaging. CBE calculation using the DC imaging method yielded more precise measurements of the backscattered energy compared to other imaging methods studied.

CONCLUSIONS

The despeckling performance of the DC imaging method significantly improves the lesion CNR in comparison to B-mode imaging. This suggests that the proposed method can detect low-contrast thermal lesions induced by FUS therapy that are not detectable using standard B-mode imaging. Furthermore, the signal change at the focal point were more precisely measured by DC imaging, and the signal change in response to FUS exposure follows the temperature profile more closely than changes measured using B-mode, as well as synthetic aperture DAS and TC images. These suggest that DC imaging can potentially be used with the CBE method to improve non-invasive thermometry.

A Review of Quantitative Ultrasound-Based Approaches to Thermometry and Ablation Zone Identification Over the Past Decade

Percutaneous thermal therapy is an important clinical treatment method for some solid tumors. It is critical to use effective image visualization techniques to monitor the therapy process in real time because precise control of the therapeutic zone directly affects the prognosis of tumor treatment. Ultrasound is used in thermal therapy monitoring because of its real-time, non-invasive, non-ionizing radiation, and low-cost characteristics. This paper presents a review of nine quantitative ultrasound-based methods for thermal therapy monitoring and their advances over the last decade since 2011. These methods were analyzed and compared with respect to two applications: ultrasonic thermometry and ablation zone identification. The advantages and limitations of these methods were compared and discussed, and future developments were suggested.

Nakagami imaging with related distributions for advanced thermogram pseudocolorization

Pseudocoloring algorithms embedded in the software of thermal cameras gradually colorize original intensity thermograms generated by detecting temperatures and contrast. Maximum and minimum based algorithms, however, executed by thresholding, applied to intensity thermograms for revealing and coloring the outliers instead. Although the common pseudocoloring protocols employed for general purposes may provide crucial information on the superficial contrast between radiation emitted by various sources; their common kernel is not sufficient for detecting and differentiating high radiated regions from surrounding areas, which is mandatory for recognition of abnormalities. Therefore, we propose novel imaging methodology based on Nakagami and related distributions, including gamma, Rayleigh, Weibull, chi-square and exponential, for enhancing thermal images and also for creating adequate discrimination. We initially define the boundaries of tumor and surrounding area in a synthetically generated breast thermogram already diagnosed as retroareolar tumor. Using Nakagami and transformations supported by mathematical foundations, we conducted several experiments to find the discrimination factor of the pseudocoloring techniques by calculating difference of average contrast between the tumor and the surrounding area. The performance is greatly encouraging that we achieved considerably better discrimination factor, designated for this study, up to 106.80 compared to the results of existing built-in pseudocolorization methods computed as 11.56.

Ultrasonic Nakagami Imaging of High-intensity Focused Ultrasound-induced Thermal Lesions in Porcine Livers: Ex Vivo Study

High-intensity focused ultrasound (HIFU) has demonstrated the capacity to be used for local thermal ablation in clinical surgery; however, relying solely on conventional ultrasound B-mode imaging to monitor HIFU thermal ablation and determine ablation levels remains a challenge. Here, we experimentally demonstrate the ability to use Nakagami imaging to monitor HIFU-induced thermal lesions in porcine livers ex vivo. Ultrasonic Nakagami imaging has been proven to be able to characterize tissues with different scatterer concentrations and distributions. The pathological sections from HIFU thermally ablated porcine liver tissues reveal that normal and denatured tissues significantly differ in scatterer concentration and distribution. Therefore, we believe that Nakagami imaging can be used to monitor thermal ablation by tracing Nakagami parameter changes in liver tissues. The ex vivo porcine liver experiments were performed using a homemade HIFU device synchronized with a commercial diagnostic ultrasound scanner to obtain the ultrasound envelope data before and after thermal ablation. These data were used to evaluate the performance of thermal lesion characterization using Nakagami imaging and were compared with those derived from conventional B-mode imaging. Experimental results showed that Nakagami imaging can be used to identify thermal lesions, which are difficult to visualize using conventional B-mode imaging because there is no apparent bubble formation. In cases with apparent bubble formation, Nakagami imaging could provide a more accurate estimation of lesion size and position. In addition, the Nakagami imaging algorithm is characterized by low computational complexity, which means it can be easily integrated as postprocessing for existing array imaging systems.

High-Intensity Focused Ultrasound Thermal Lesion Detection Using Entropy Imaging of Ultrasound Radio Frequency Signal Time Series

Background:

During the past few decades, high-intensity focused ultrasound (HIFU) modality has been gaining surging interest in various therapeutic applications such as non- or minimally-invasive cancer treatment. Among other attributes, robust and real-time HIFU treatment monitoring and lesion detection have become essential issues for successful clinical acceptance of the modality. More recently, ultrasound radio frequency (RF) time series imaging has been studied by a number of researchers.

Materials and Methods:

The objective of this study is to investigate the applicability of entropy parameter of RF time series of ultrasound backscattered signals, a. k. a. Entropy imaging, toward HIFU thermal lesion detection. To this end, five fresh ex vivo porcine muscle tissue samples were exposed to HIFU exposures with total acoustic powers ranging from 30 to 110 Watts. The contrast-to-speckle ratio (CSR) values of the entropy images and their corresponding B-mode images of pre-, during- and post-HIFU exposure for each acoustic power were calculated.

Results:

The novelty of this study is the use of Entropy parameter on ultrasound RF time series for the first time. Statistically significant differences were obtained between the CSR values for the B mode and entropy images at various acoustic powers. In case of 110 Watt, a CSR value 3.4 times higher than B-mode images was accomplished using the proposed method. Furthermore, the proposed method is compared with the scaling parameter of Nakagami imaging and same data which are used in this study.

Conclusion:

Entropy has the potential for using as an imaging parameter for differentiating lesions in HIFU surgery.

Three-dimensional Visualization of Ultrasound Backscatter Statistics by Window-modulated Compounding Nakagami Imaging

In this study, the window-modulated compounding (WMC) technique was integrated into three-dimensional (3D) ultrasound Nakagami imaging for improving the spatial visualization of backscatter statistics. A 3D WMC Nakagami image was produced by summing and averaging a number of 3D Nakagami images (number of frames denoted as N) formed using sliding cubes with varying side lengths ranging from 1 to N times the transducer pulse. To evaluate the performance of the proposed 3D WMC Nakagami imaging method, agar phantoms with scatterer concentrations ranging from 2 to 64 scatterers/mm³ were made, and six stages of fatty liver (zero, one, two, four, six, and eight weeks) were induced in rats by methionine-choline-deficient diets (three rats for each stage, total n = 18). A mechanical scanning system with a 5-MHz focused single-element transducer was used for ultrasound radiofrequency data acquisition. The experimental results showed that 3D WMC Nakagami imaging was able to characterize different scatterer concentrations. Backscatter statistics were visualized with various numbers of frames; N = 5 reduced the estimation error of 3D WMC Nakagami imaging in visualizing the backscatter statistics. Compared with conventional 3D Nakagami imaging, 3D WMC Nakagami imaging improved the image smoothness without significant image resolution degradation, and it can thus be used for describing different stages of fatty liver in rats.

Determining temperature distribution in tissue in the focal plane of the high (>100 W/cm(2)) intensity focused ultrasound beam using phase shift of ultrasound echoes

In therapeutic applications of High Intensity Focused Ultrasound (HIFU) the guidance of the HIFU beam and especially its focal plane is of crucial importance. This guidance is needed to appropriately target the focal plane and hence the whole focal volume inside the tumor tissue prior to thermo-ablative treatment and beginning of tissue necrosis. This is currently done using Magnetic Resonance Imaging that is relatively expensive. In this study an ultrasound method, which calculates the variations of speed of sound in the locally heated tissue volume by analyzing the phase shifts of echo-signals received by an ultrasound scanner from this very volume is presented. To improve spatial resolution of B-mode imaging and minimize the uncertainty of temperature estimation the acoustic signals were transmitted and received by 8 MHz linear phased array employing Synthetic Transmit Aperture (STA) technique. Initially, the validity of the algorithm developed was verified experimentally in a tissue-mimicking phantom heated from 20.6 to 48.6 °C. Subsequently, the method was tested using a pork loin sample heated locally by a 2 MHz pulsed HIFU beam with focal intensity ISATA of 129 W/cm(2). The temperature calibration of 2D maps of changes in the sound velocity induced by heating was performed by comparison of the algorithm-determined changes in the sound velocity with the temperatures measured by thermocouples located in the heated tissue volume. The method developed enabled ultrasound temperature imaging of the heated tissue volume from the very inception of heating with the contrast-to-noise ratio of 3.5-12 dB in the temperature range 21-56 °C. Concurrently performed, conventional B-mode imaging revealed CNR close to zero dB until the temperature reached 50 °C causing necrosis. The data presented suggest that the proposed method could offer an alternative to MRI-guided temperature imaging for prediction of the location and extent of the thermal lesion prior to applying the final HIFU treatment.

Nakagami-based total variation method for speckle reduction in thyroid ultrasound images

A good statistical model is necessary for the reduction in speckle noise. The Nakagami model is more general than the Rayleigh distribution for statistical modeling of speckle in ultrasound images. In this article, the Nakagami-based noise removal method is presented to enhance thyroid ultrasound images and to improve clinical diagnosis. The statistics of log-compressed image are derived from the Nakagami distribution following a maximum a posteriori estimation framework. The minimization problem is solved by optimizing an augmented Lagrange and Chambolle's projection method. The proposed method is evaluated on both artificial speckle-simulated and real ultrasound images. The experimental findings reveal the superiority of the proposed method both quantitatively and qualitatively in comparison with other speckle reduction methods reported in the literature. The proposed method yields an average signal-to-noise ratio gain of more than 2.16 dB over the non-convex regularizer-based speckle noise removal method, 3.83 dB over the Aubert-Aujol model, 1.71 dB over the Shi-Osher model and 3.21 dB over the Rudin-Lions-Osher model on speckle-simulated synthetic images. Furthermore, visual evaluation of the despeckled images shows that the proposed method suppresses speckle noise well while preserving the textures and fine details.

Effects of Estimators on Ultrasound Nakagami Imaging in Visualizing the Change in the Backscattered Statistics from a Rayleigh Distribution to a Pre-Rayleigh Distribution

Ultrasound Nakagami imaging has recently attracted interest as an imaging technique for analyzing envelope statistics. Because the presence of structures has a strong effect on estimation of the Nakagami parameter, previous studies have indicated that Nakagami imaging should be used specifically for characterization of soft tissues with fewer structures, such as liver tissues. Typically, changes in the properties of the liver parenchyma cause the backscattered statistics to transform from a Rayleigh distribution to a pre-Rayleigh distribution, and this transformation can be visualized using a Nakagami imaging technique. However, different estimators result in different estimated values; thus, the performance of a Nakagami image may depend on the type of estimator used. This study explored the effects of various estimators on ultrasound Nakagami imaging to describe the backscattered statistics as they change from a Rayleigh distribution to a pre-Rayleigh distribution. Simulations and clinical measurements involving patients with liver fibrosis (n = 85) yielded image data that were used to construct B-mode and conventional Nakagami images based on the moment estimator (denoted as mINV images) and maximum-likelihood estimator (denoted as mML images). In addition, novel window-modulated compounding Nakagami images based on the moment estimator (denoted as mWMC images) were also obtained. The means and standard deviations of the Nakagami parameters were examined as a function of the backscattered statistics. The experimental results indicate that the mINV, mML and mWMC images enabled quantitative visualization of the change in backscattered statistics from a Rayleigh distribution to a pre-Rayleigh distribution. Importantly, the mWMC image is superior to both mINV and mML images because it simultaneously realizes sensitive detection of the backscattered statistics and a reduction of estimation variance for image smoothness improvement. We therefore recommend using mWMC image as a novel strategy in Nakagami imaging technique for liver tissue characterization.

Thermometry and ablation monitoring with ultrasound

In this review we present the current status of ultrasound thermometry and ablation monitoring, with emphasis on the diverse approaches published in the literature and with an eye on which methods are closest to clinical reality. It is hoped that this review will serve as a guide to the expansion of sonographic methods for treatment monitoring and thermometry since the last brief review in 2007.

Monitoring radiofrequency ablation using real-time ultrasound Nakagami imaging combined with frequency and temporal compounding techniques

Gas bubbles induced during the radiofrequency ablation (RFA) of tissues can affect the detection of ablation zones (necrosis zone or thermal lesion) during ultrasound elastography. To resolve this problem, our previous study proposed ultrasound Nakagami imaging for detecting thermal-induced bubble formation to evaluate ablation zones. To prepare for future applications, this study (i) created a novel algorithmic scheme based on the frequency and temporal compounding of Nakagami imaging for enhanced ablation zone visualization, (ii) integrated the proposed algorithm into a clinical scanner to develop a real-time Nakagami imaging system for monitoring RFA, and (iii) investigated the applicability of Nakagami imaging to various types of tissues. The performance of the real-time Nakagami imaging system in visualizing RFA-induced ablation zones was validated by measuring porcine liver (n = 18) and muscle tissues (n = 6). The experimental results showed that the proposed algorithm can operate on a standard clinical ultrasound scanner to monitor RFA in real time. The Nakagami imaging system effectively monitors RFA-induced ablation zones in liver tissues. However, because tissue properties differ, the system cannot visualize ablation zones in muscle fibers. In the future, real-time Nakagami imaging should be focused on the RFA of the liver and is suggested as an alternative monitoring tool when advanced elastography is unavailable or substantial bubbles exist in the ablation zone.

Real-time monitoring of high-intensity focused ultrasound thermal therapy using the manifold learning method

High-intensity focused ultrasound (HIFU) induces thermal lesions by increasing the tissue temperature in a tight focal region. The main ultrasound imaging techniques currently used to monitor HIFU treatment are standard pulse-echo B-mode ultrasound imaging, ultrasound temperature estimation and elastography-based methods. The present study was carried out on ex vivo animal tissue samples, in which backscattered radiofrequency (RF) signals were acquired in real time at time instances before, during and after HIFU treatment. The manifold learning algorithm, a non-linear dimensionality reduction method, was applied to RF signals whichconstruct B-mode images to detect the HIFU-induced changes among the image frames obtained during HIFU treatment. In this approach, the embedded non-linear information in the region of interest of sequential images is represented in a 2-D manifold with the Isomap algorithm, and each image is depicted as a point on the reconstructed manifold. Four distinct regions are chosen in the manifold corresponding to the four phases of HIFU treatment (before HIFU treatment, during HIFU treatment, immediately after HIFU treatment and 10-min after HIFU treatment). It was found that disorganization of the points is achieved by increasing the acoustic power, and if the thermal lesion has been formed, the regions of points related to pre- and post-HIFU significantly differ. Moreover, the manifold embedding was repeated on 2-D moving windows in RF data envelopes related to pre- and post-HIFU exposure data frames. It was concluded that if mean values of the points related to pre- and post-exposure frames in the reconstructed manifold are estimated, and if the Euclidean distance between these two mean values is calculated and the sliding window is moved and this procedure is repeated for the whole image, a new image based on the Euclidean distance can be formed in which the HIFU thermal lesion is detectable.

Picture the future: emerging imaging modalities

The past century has witnessed accelerated development in imaging modalities. Better anatomical visualisation and improved data analysis have improved survival rates. Through emerging functional, molecular and structural imaging modalities, better anatomical visualisation has been extended to cellular and molecular detail, improving diagnosis and management of diseases. This article reviews the advances made in emerging imaging modalities as well as their potential applications in targeted therapy.

A survey of ultrasound elastography approaches to percutaneous ablation monitoring

Percutaneous thermal ablation has been widely used as a minimally invasive treatment for tumors. Treatment monitoring is essential for preventing complications while ensuring treatment efficacy. Mechanical testing measurements on tissue reveal that tissue stiffness increases with temperature and ablation duration. Different types of imaging methods can be used to monitor ablation procedures, including temperature or thermal strain imaging, strain imaging, modulus imaging, and shear modulus imaging. Ultrasound elastography demonstrates the potential to become the primary imaging modality for monitoring percutaneous ablation. This review briefly presented the state-of-the-art ultrasound elastography approaches for monitoring radiofrequency ablation and microwave ablation. These techniques were divided into four groups: quasi-static elastography, acoustic radiation force elastography, sonoelastography, and applicator motion elastography. Their advantages and limitations were compared and discussed. Future developments were proposed with respect to heat-induced bubbles, tissue inhomogeneities, respiratory motion, three-dimensional monitoring, multi-parametric monitoring, real-time monitoring, experimental data center for percutaneous ablation, and microwave ablation monitoring.