Using ultrasound CBE imaging without echo shift compensation for temperature estimation

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.

Detection of microwave ablation coagulation areas using ultrasound Nakagami imaging based on Gaussian pyramid decomposition: A feasibility study

In this paper, we explored the feasibility of using ultrasound Nakagami-m parametric imaging based on Gaussian pyramid decomposition (GPD) to detect microwave ablation coagulation areas. Monte Carlo simulation and phantom simulation results demonstrated that a 2-layer GPD model was sufficient to achieve the same m parameter estimation accuracy, smoothness and resolution as 3-layer and 4-layer. The performances of GPD, moment-based estimator (MBE) and window-modulated compounding (WMC) algorithms were compared in terms of parameter estimation, smoothness, resolution and contrast-to-noise (CNR). Results showed that the m parameter estimation obtained by GPD algorithm was better than that of MBE and WMC algorithms except the small window size (27 × 5). When using a window size of >3 pulse lengths, GPD algorithm could achieve better smoothness and CNR than MBE and WMC algorithms, but there was a certain loss of axial resolution. The computation time of GPD algorithm was less than that of WMC algorithm, while about 2.24 times that of MBE algorithm. Experimental results of porcine liver microwave ablation ex vivo (n = 20) illustrated that the average areas under the operating characteristic curve (AUCs) of Nakagami m_GPD, m_MBE and m_WMC parametric imaging and homodyned-K (HK) α and k parametric imaging to detect coagulation areas were significantly improved by polynomial approximation (PAX). Kruskal-Wallis test showed that the accuracy of coagulation area detection obtained by PAX imaging of m_GPD parameter had no significant difference with that of m_MBE, m_WMC, HK_α and HK_k parameters. This preliminary study suggested that Nakagami imaging based on GPD algorithm may have the potential to detect microwave ablation coagulation areas.

Ultrasound single-phase CBE imaging for monitoring radiofrequency ablation of the liver tumor: A preliminary clinical validation

Radiofrequency ablation (RFA) is an alternative treatment for early-stage hepatocellular carcinoma (HCC). The production of gas bubbles by RFA indicates threshold temperature of tissue necrosis and results in changes in backscattered energy (CBE) when ultrasound monitors RFA. In this study, ultrasound single-phase CBE imaging was used as a means of monitoring RFA of the liver tumor by analyzing the backscattering of ultrasound from gas bubbles in the liver. A total of 19 HCC patients were enrolled in the study. An ultrasound system was used during RFA to monitor the ablation process and acquire raw image data consisting of backscattered signals for single-phase CBE imaging. On the basis of single-phase CBE imaging, the area corresponding to the range of gas bubbles was compared with the tumor sizes and ablation zones estimated from computed tomography. During RFA, ultrasound single-phase CBE imaging enabled improved visualization of gas bubbles. Measured gas bubble areas by CBE were related to tumor size (the Spearman correlation coefficient r_s = 0.86; p < 0.05); less dependent on the ablation zone. Approximately 95% of the data fell within the limits of agreement in Bland-Altman plots, and 58% of the data fell within the 95% CI. This study suggests that single-phase CBE imaging provides information about liver tumor size because of the abundant vessels in liver tumors that promote the generation of gas bubbles, which serve as natural contrast agents in RFAs to enhance ultrasound backscattering. Ultrasound single-phase CBE imaging may allow clinicians to determine if the required minimum RFA efficacy level is reached by assessing gas bubbles in the liver tumors.

High contrast ultrasonic method with multi-spatiotemporal compounding for monitoring catheter-based ultrasound thermal therapy: Development and Ex Vivo Evaluations

OBJECTIVE

Changes in ultrasound backscatter energy (CBE) imaging can monitor thermal therapy. Catheter-based ultrasound (CBUS) can treat deep tumors with precise spatial control of energy deposition and ablation zones, of which CBE estimation can be limited by low contrast and robustness due to small or inconsistent changes in ultrasound data. This study develops a multi-spatiotemporal compounding CBE (MST-CBE) imaging approach for monitoring specific to CBUS thermal therapy.

METHODS

Ex vivo thermal ablations were performed with stereotactic positioning of a 180 directional CBUS applicator, temperature monitoring probes, endorectal US probe, and subsequent lesion sectioning and measurement. Five frames of raw radiofrequency data were acquired throughout in 15s intervals. Using window-by-window estimation methods, absolute and positive components of MST-CBE images at each point were obtained by the compounding ratio of squared envelope data within an increasing spatial size in each short-time window.

RESULTS

Compared with conventional US, Nakagami, and CBE imaging, the detection contrast and robustness quantified by tissue-modification-ratio improved by 37.24.7 (p<0.001), 37.55.2 (p<0.001), and 6.44.0 dB (p<0.05) in the MST-CBE imaging, respectively. Correlation coefficient and bias between cross-sectional dimensions of the ablation zones measured in tissue sections and estimated from MST-CBE were up to 0.91 (p<0.001) and -0.02 mm2, respectively.

CONCLUSION

The MST-CBE approach can monitor the detailed changes within target tissues and effectively characterize the dimensions of the ablation zone during CBUS energy deposition.

SIGNIFICANCE

The MST-CBE approach could be practical for improved accuracy and contrast of monitoring and evaluation for CBUS thermal therapy.

Frequency-domain CBE imaging for ultrasound localization of the HIFU focal spot: a feasibility study

High-intensity focused ultrasound (HIFU) is a well-accepted tool for noninvasive thermal therapy. To control the quality of HIFU treatment, the focal spot generated in tissues must be localized. Ultrasound imaging can monitor heated regions; in particular, the change in backscattered energy (CBE) allows parametric imaging to visualize thermal information in the tissue. Conventional CBE imaging constructed in the spatial domain may be easily affected by noises when the HIFU focal spot is visualized. This study proposes frequency-domain CBE imaging to improve noise tolerance and image contrast in HIFU focal spot monitoring. Phantom experiments were performed in a temperature-controlled environment. HIFU of 2.12 MHz was applied to the phantoms, during which a clinical scanner equipped with a 3-MHz convex array transducer was used to collect raw image data consisting of backscattered signals for B-mode, spatial-, and frequency-domain CBE imaging. Concurrently, temperature changes were measured at the focal spot using a thermocouple for comparison with CBE values by calculating the correlation coefficient r. To further analyze CBE image contrast levels, a contrast factor was introduced, and an independent t-test was performed to calculate the probability value p. Experimental results showed that frequency-domain CBE imaging performed well in thermal distribution visualization, enabling quantitative detection of temperature changes. The CBE value calculated in the frequency domain also correlated strongly with that obtained using the conventional spatial-domain approach (r = 0.97). In particular, compared with the image obtained through the conventional method, the contrast of the CBE image obtained using the method based on frequency-domain analysis increased by 2.5-fold (4 dB; p < 0.05). Frequency-domain computations may constitute a new strategy when ultrasound CBE imaging is used to localize the focal spot in HIFU treatment planning.

Real-Time Radiofrequency Ablation Lesion Depth Estimation Using Multi-frequency Impedance With a Deep Neural Network and Tree-Based Ensembles

OBJECTIVE

Design and optimization of statistical models for use in methods for estimating radiofrequency ablation (RFA) lesion depths in soft real-time performance.

METHODS

Using tissue multi-frequency complex electrical impedance data collected from a low-cost embedded system, a deep neural network (NN) and tree-based ensembles (TEs) were trained for estimating the RFA lesion depth via regression.

RESULTS

Addition of frequency sweep data, previous depth data, and previous RF power state data boosted accuracy of the statistical models. The root mean square errors were 2 mm for NN and 0.5 mm for TEs for previous statistical models and the root mean square errors were 0.4 mm for NN and 0.04 mm for TEs for the statistical models presented in this paper. Simulation ablation performance showed a mean difference against physical measurements of 0.5 ±0.2 mm for the NN-based depth estimation method and 0.7 ±0.4 mm for the TE-based depth estimation method.

CONCLUSION

The results show that multi-frequency data significantly improves the depth estimation performance of the statistical models.

SIGNIFICANCE

The RFA lesion depth estimation methods presented in this work achieve millimeter-resolution accuracy with soft real-time performance on an ARMv7-based embedded system for potential translation to clinical RFA technologies.

On the Feasibility of Ultrasound Imaging Enrichment by Medium-Temperature Changes

We describe the concept of a new imaging modality based on the tracking and dynamic modeling of local intensity changes (ICs) observed in conventional ultrasound images collected during a medium-temperature change. We computed the pixel-by-pixel IC from averaged B-mode images that exhibited different behaviors with varying temperature resulting from changes in the speed of sound, which consequently induce changes in the backscattered energy. Moreover, for each pixel, a first-order polynomial model was adjusted to the different temperature-dependent ICs. The representation of the polynomial angular parameter in 2D pixel space was used to obtain a parametric image. The results obtained by simulations and with real B-mode images indicated that this new ultrasound imaging modality was able to enhance the contrast and highlight structures that were poorly visible or even undetected in conventional images. A temperature change of 3°C was found to be sufficient to generate appropriate images with the proposed method. In addition, if a temperature change of 6°C was considered, the thermal dose, measured as the cumulative number of equivalent minutes at 43°C (CEM43°C), was 2.4 CEM43°C, which is a value that is considered safe according to the literature. We provide a proof-of-concept of a new imaging modality that opens new opportunities for the enhancement of ultrasound images and consequently contributes to improvements in ultrasound-based diagnoses. Our approach is based on images returned by commercial ultrasound scanners. Therefore, it can be implemented in any ultrasound system and is independent of specific ultrasound hardware and software data acquisition characteristics.

Ultrasound single-phase CBE imaging for monitoring radiofrequency ablation

Radiofrequency (RF) ablation (RFA) is the most commonly used minimally invasive procedure for thermal ablation of liver tumors. Ultrasound not only provides real-time feedback of the electrode location for RFA guidance but also enables visualization of the tissue temperature. Changes in backscattered energy (CBE) have been widely applied to ultrasound temperature imaging for assessing thermal ablation. Pilot studies have revealed that significant shadowing features appear in CBE imaging and are caused by the electrode and RFA-induced gas bubbles. To resolve this problem, the current study proposed ultrasound single-phase CBE imaging based on positive CBE values. An in vitro model with tissue samples derived from the porcine tenderloin was used to validate the proposed method. During RFA with various electrode lengths, ultrasound scans of tissue samples were obtained using a clinical ultrasound scanner equipped with a convex array transducer of 3 MHz. Raw image data comprising 256 scan lines of backscattered RF signals were acquired for B-mode, conventional CBE, and single-phase CBE imaging by using the proposed algorithmic scheme. The ablation sizes estimated using CBE imaging and gross examinations were compared to calculate the correlation coefficient. The experimental results indicated that single-phase CBE imaging largely suppressed artificial CBE information in the shadowed region. Moreover, compared with conventional CBE imaging, single-phase CBE imaging provided a more accurate estimation of ablation sizes (the correlation coefficient was higher than 0.8).

Adaptive ultrasound temperature imaging for monitoring radiofrequency ablation

Radiofrequency ablation (RFA) has been widely used as an alternative treatment modality for liver tumors. Monitoring the temperature distribution in the tissue during RFA is required to assess the thermal dosage. Ultrasound temperature imaging based on the detection of echo time shifts has received the most attention in the past decade. The coefficient k, connecting the temperature change and the echo time shift, is a medium-dependent parameter used to describe the confounding effects of changes in the speed of sound and thermal expansion as temperature increases. The current algorithm of temperature estimate based on echo time shift detection typically uses a constant k, resulting in estimation errors when ablation temperatures are higher than 50°C. This study proposes an adaptive-k algorithm that enables the automatic adjustment of the coefficient k during ultrasound temperature monitoring of RFA. To verify the proposed algorithm, RFA experiments on in vitro porcine liver samples (total n = 15) were performed using ablation powers of 10, 15, and 20 W. During RFA, a clinical ultrasound system equipped with a 7.5-MHz linear transducer was used to collect backscattered signals for ultrasound temperature imaging using the constant- and adaptive-k algorithms. Concurrently, an infrared imaging system and thermocouples were used to measure surface temperature distribution of the sample and internal ablation temperatures for comparisons with ultrasound estimates. Experimental results demonstrated that the proposed adaptive-k method improved the performance in visualizing the temperature distribution. In particular, the estimation errors were also reduced even when the temperature of the tissue is higher than 50°C. The proposed adaptive-k ultrasound temperature imaging strategy has potential to serve as a thermal dosage evaluation tool for monitoring high-temperature RFA.

Adaptive ultrasound temperature imaging for monitoring radiofrequency ablation

Radiofrequency ablation (RFA) has been widely used as an alternative treatment modality for liver tumors. Monitoring the temperature distribution in the tissue during RFA is required to assess the thermal dosage. Ultrasound temperature imaging based on the detection of echo time shifts has received the most attention in the past decade. The coefficient k, connecting the temperature change and the echo time shift, is a medium-dependent parameter used to describe the confounding effects of changes in the speed of sound and thermal expansion as temperature increases. The current algorithm of temperature estimate based on echo time shift detection typically uses a constant k, resulting in estimation errors when ablation temperatures are higher than 50°C. This study proposes an adaptive-k algorithm that enables the automatic adjustment of the coefficient k during ultrasound temperature monitoring of RFA. To verify the proposed algorithm, RFA experiments on in vitro porcine liver samples (total n = 15) were performed using ablation powers of 10, 15, and 20 W. During RFA, a clinical ultrasound system equipped with a 7.5-MHz linear transducer was used to collect backscattered signals for ultrasound temperature imaging using the constant- and adaptive-k algorithms. Concurrently, an infrared imaging system and thermocouples were used to measure surface temperature distribution of the sample and internal ablation temperatures for comparisons with ultrasound estimates. Experimental results demonstrated that the proposed adaptive-k method improved the performance in visualizing the temperature distribution. In particular, the estimation errors were also reduced even when the temperature of the tissue is higher than 50°C. The proposed adaptive-k ultrasound temperature imaging strategy has potential to serve as a thermal dosage evaluation tool for monitoring high-temperature RFA.

Temperature monitoring by channel data delays: Feasibility based on estimated delays magnitude for cardiac ablation

Ultrasound thermometry is based on measuring tissue temperature by its impact on ultrasound wave propagation. This study focuses on the use of transducer array channel data (not beamformed) and examines how a layer of increased velocity (heat induced) affects the travel-times of the ultrasound backscatter signal. Based on geometric considerations, a new equation was derived for the change in time delay as a function of temperature change. The resulting expression provides insight into the key factors that link change in temperature to change in travel time. It shows that velocity enters in combination with heating geometry: complementary information is needed to compute velocity from the changes in travel time. Using the bio-heat equation as a second source of information in the derived expressions, the feasibility of monitoring the temperature increase during cardiac ablation therapy using channel data was investigated. For an intra-cardiac (ICE) probe, using this "time delay error approach" would not be feasible, while for a trans-esophageal array transducer (TEE) transducer it might be feasible.

Feasibility of Using Ultrasonic Nakagami Imaging for Monitoring Microwave-Induced Thermal Lesion in Ex Vivo Porcine Liver

The feasibility of using ultrasonic Nakagami imaging to evaluate thermal lesions induced by microwave ablation (MWA) in ex vivo porcine liver was explored. Dynamic changes in echo amplitudes and Nakagami parameters in the region of the MWA-induced thermal lesion, as well as the contrast-to-noise ratio (CNR) between the MWA-induced thermal lesion and the surrounding normal tissue, were calculated simultaneously during the MWA procedure. After MWA exposure, a bright hyper-echoic region appeared in ultrasonic B-mode and Nakagami parameter images as an indicator of the thermal lesion. Mean values of the Nakagami parameter in the thermal lesion region increased to 0.58, 0.71 and 0.91 after 1, 3 and 5 min of MVA. There were no significant differences in envelope amplitudes in the thermal lesion region among ultrasonic B-mode images obtained after different durations of MWA. Unlike ultrasonic B-mode images, Nakagami images were less affected by the shadow effect in monitoring of MWA exposure, and a fairly complete hyper-echoic region was observed in the Nakagami image. The mean value of the Nakagami parameter increased from approximately 0.47 to 0.82 during MWA exposure. At the end of the postablation stage, the mean value of the Nakagami parameter decreased to 0.55 and was higher than that before MWA exposure. CNR values calculated for Nakagami parameter images increased from 0.13 to approximately 0.61 during MWA and then decreased to 0.26 at the end of the post-ablation stage. The corresponding CNR values calculated for ultrasonic B-mode images were 0.24, 0.42 and 0.17. This preliminary study on ex vivo porcine liver suggested that Nakagami imaging have potential use in evaluating the formation of MWA-induced thermal lesions. Further in vivo studies are needed to evaluate the potential application.

Calibration and Evaluation of Ultrasound Thermography Using Infrared Imaging

Real-time monitoring of the spatiotemporal evolution of tissue temperature is important to ensure safe and effective treatment in thermal therapies including hyperthermia and thermal ablation. Ultrasound thermography has been proposed as a non-invasive technique for temperature measurement, and accurate calibration of the temperature-dependent ultrasound signal changes against temperature is required. Here we report a method that uses infrared thermography for calibration and validation of ultrasound thermography. Using phantoms and cardiac tissue specimens subjected to high-intensity focused ultrasound heating, we simultaneously acquired ultrasound and infrared imaging data from the same surface plane of a sample. The commonly used echo time shift-based method was chosen to compute ultrasound thermometry. We first correlated the ultrasound echo time shifts with infrared-measured temperatures for material-dependent calibration and found that the calibration coefficient was positive for fat-mimicking phantom (1.49 ± 0.27) but negative for tissue-mimicking phantom (-0.59 ± 0.08) and cardiac tissue (-0.69 ± 0.18°C-mm/ns). We then obtained the estimation error of the ultrasound thermometry by comparing against the infrared-measured temperature and revealed that the error increased with decreased size of the heated region. Consistent with previous findings, the echo time shifts were no longer linearly dependent on temperature beyond 45°C-50°C in cardiac tissues. Unlike previous studies in which thermocouples or water bath techniques were used to evaluate the performance of ultrasound thermography, our results indicate that high-resolution infrared thermography is a useful tool that can be applied to evaluate and understand the limitations of ultrasound thermography methods.

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.

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.

Feasibility of non-invasive temperature estimation by the assessment of the average gray-level content of B-mode images

This paper assesses the potential of the average gray-level (AVGL) from ultrasonographic (B-mode) images to estimate temperature changes in time and space in a non-invasive way. Experiments were conducted involving a homogeneous bovine muscle sample, and temperature variations were induced by an automatic temperature regulated water bath, and by therapeutic ultrasound. B-mode images and temperatures were recorded simultaneously. After data collection, regions of interest (ROIs) were defined, and the average gray-level variation computed. For the selected ROIs, the AVGL-Temperature relation were determined and studied. Based on uniformly distributed image partitions, two-dimensional temperature maps were developed for homogeneous regions. The color-coded temperature estimates were first obtained from an AVGL-Temperature relation extracted from a specific partition (where temperature was independently measured by a thermocouple), and then extended to the other partitions. This procedure aimed to analyze the AVGL sensitivity to changes not only in time but also in space. Linear and quadratic relations were obtained depending on the heating modality. We found that the AVGL-Temperature relation is reproducible over successive heating and cooling cycles. One important result was that the AVGL-Temperature relations extracted from one region might be used to estimate temperature in other regions (errors inferior to 0.5 °C) when therapeutic ultrasound was applied as a heating source. Based on this result, two-dimensional temperature maps were developed when the samples were heated in the water bath and also by therapeutic ultrasound. The maps were obtained based on a linear relation for the water bath heating, and based on a quadratic model for the therapeutic ultrasound heating. The maps for the water bath experiment reproduce an acceptable heating/cooling pattern, and for the therapeutic ultrasound heating experiment, the maps seem to reproduce temperature profiles consistent with the pressure field of the transducer, and in agreement with temperature maps developed by COMSOL®MultiPhysics simulations.

An approach for the visualization of temperature distribution in tissues according to changes in ultrasonic backscattered energy

Previous studies developed ultrasound temperature-imaging methods based on changes in backscattered energy (CBE) to monitor variations in temperature during hyperthermia. In conventional CBE imaging, tracking and compensation of the echo shift due to temperature increase need to be done. Moreover, the CBE image does not enable visualization of the temperature distribution in tissues during nonuniform heating, which limits its clinical application in guidance of tissue ablation treatment. In this study, we investigated a CBE imaging method based on the sliding window technique and the polynomial approximation of the integrated CBE (ICBE_pa image) to overcome the difficulties of conventional CBE imaging. We conducted experiments with tissue samples of pork tenderloin ablated by microwave irradiation to validate the feasibility of the proposed method. During ablation, the raw backscattered signals were acquired using an ultrasound scanner for B-mode and ICBE_pa imaging. The experimental results showed that the proposed ICBE_pa image can visualize the temperature distribution in a tissue with a very good contrast. Moreover, tracking and compensation of the echo shift were not necessary when using the ICBE_pa image to visualize the temperature profile. The experimental findings suggested that the ICBE_pa image, a new CBE imaging method, has a great potential in CBE-based imaging of hyperthermia and other thermal therapies.

Ultrasonic evaluation of microwave-induced thermal lesions based on wavelet analysis of mean scatterer spacing

The microwave ablation has become an important manner for tumor treatment. In this paper, we proposed a new method for evaluation of microwave-induced thermal lesions using the wavelet analysis of the mean scatterer spacing (MSS). First, the ultrasonic radiofrequency (RF) data of normal and coagulated porcine liver tissues was collected through the temperature-controlled water bath heating experiments. The convex array ultrasound probe with a center frequency of 3.5 MHz was used. Second, the wavelet analysis was used to compute the MSS of normal and coagulated porcine liver tissues, respectively. Finally, the microwave-induced thermal lesions were detected based on the differences in the MSS between normal tissues and thermal lesions. Eighteen cases of microwave ablation experiments and 20 cases of water bath heating experiments were conducted on fresh porcine liver samples. The MSS of normal porcine liver tissues was 1.15±0.12 mm, and the MSS of coagulated porcine liver tissues was 0.93±0.07 mm. Six cases of thermal lesions were compared between the MSS-detected area and the caliper-measured area, and the MSS-detected area had an error of (13.55±5.29) %. The experimental results indicated that the proposed method could be used in preliminary detection and evaluation of microwave-induced thermal lesions.