Ultrasound monitoring of in vitro radio frequency ablation by echo decorrelation imaging

Decorrelation Time Mapping as an Analysis Tool for Nanobubble-Based Contrast Enhanced Ultrasound Imaging

Nanobubbles (NBs; ~100-500 nm diameter) are preclinical ultrasound (US) contrast agents that expand applications of contrast enhanced US (CEUS). Due to their sub-micron size, high particle density, and deformable shell, NBs in pathological states of heightened vascular permeability (e.g. in tumors) extravasate, enabling applications not possible with microbubbles (~1000-10,000 nm diameter). A method that can separate intravascular versus extravascular NB signal is needed as an imaging biomarker for improved tumor detection. We present a demonstration of decorrelation time (DT) mapping for enhanced tumor NB-CEUS imaging. In vitro models validated the sensitivity of DT to agent motion. Prostate cancer mouse models validated in vivo imaging potential and sensitivity to cancerous tissue. Our findings show that DT is inversely related to NB motion, offering enhanced detail of NB dynamics in tumors, and highlighting the heterogeneity of the tumor environment. Average DT was high in tumor regions (~9 s) compared to surrounding normal tissue (~1 s) with higher sensitivity to tumor tissue compared to other mapping techniques. Molecular NB targeting to tumors further extended DT (11 s) over non-targeted NBs (6 s), demonstrating sensitivity to NB adherence. From DT mapping of in vivo NB dynamics we demonstrate the heterogeneity of tumor tissue while quantifying extravascular NB kinetics and delineating intra-tumoral vasculature. This new NB-CEUS-based biomarker can be powerful in molecular US imaging, with improved sensitivity and specificity to diseased tissue and potential for use as an estimator of vascular permeability and the enhanced permeability and retention (EPR) effect in tumors.

Ultrasound k-nearest neighbor entropy imaging: Theory, algorithm, and applications

Ultrasound information entropy is a flexible approach for analyzing ultrasound backscattering. Shannon entropy imaging based on probability distribution histograms (PDHs) has been implemented as a promising method for tissue characterization and diagnosis. However, the bin number affects the stability of entropy estimation. In this study, we introduced the k-nearest neighbor (KNN) algorithm to estimate entropy values and proposed ultrasound KNN entropy imaging. The proposed KNN estimator leveraged the Euclidean distance between data samples, rather than the histogram bins by conventional PDH estimators. We also proposed cumulative relative entropy (CRE) imaging to analyze time-series radiofrequency signals and applied it to monitor thermal lesions induced by microwave ablation (MWA). Computer simulation phantom experiments were conducted to validate and compare the performance of the proposed KNN entropy imaging, the conventional PDH entropy imaging, and Nakagami-m parametric imaging in detecting the variations of scatterer densities and visualizing inclusions. Clinical data of breast lesions were analyzed, and porcine liver MWA experiments ex vivo were conducted to validate the performance of KNN entropy imaging in classifying benign and malignant breast tumors and monitoring thermal lesions, respectively. Compared with PDH, the entropy estimation based on KNN was less affected by the tuning parameters. KNN entropy imaging was more sensitive to changes in scatterer densities and performed better visualizable capability than typical Shannon entropy (TSE) and Nakagami-m parametric imaging. Among different imaging methods, KNN-based Shannon entropy (KSE) imaging achieved the higher accuracy in classification of benign and malignant breast tumors and KNN-based CRE imaging had larger lesion-to-normal contrast when monitoring the ablated areas during MWA at different powers and treatment durations. Ultrasound KNN entropy imaging is a potential quantitative ultrasound approach for tissue characterization.

Monitoring of thermal lesions in ultrasound using fully convolutional neural networks: A preclinical study

Accurate monitoring of thermal ablation regions is an important guarantee for successful ablation treatment, which mainly depends on the subjective judgment of radiologists in current clinical practice. This work innovatively applied fully convolutional neural networks (FCNs) for detection and monitoring of thermal ablation regions in ultrasound (US) and comprehensively compared the performance of VGG16-FCN, U-Net, UNet++, Attention U-Net, MultiResUNet, and ResUNet, which have shown outstanding performance in medical image segmentation. The input of the models was US echo envelope data backscattered from the ablated regions. Excised porcine liver ablation dataset and clinical liver tumors ablation dataset were respectively used to evaluate the prediction ability of the models. With 1000 excised porcine liver ablation samples for training and 200 samples for testing, the UNet++ achieves both the highest Dice score (DSC) of 0.7824 ± 0.1098 and the best Hausdorff distance (HD) of 2.70 ± 1.38 mm. Additionally, considering potential clinical usage, we also tested the model generalizability by training on the excised dataset and testing on the clinical data, in which we obtained the performance with the highest DSC obtained by the ResUNet and the best HD by the UNet++. Our comparative study suggests that both UNet++ and ResUNet have relatively outstanding segmentation performance among all compared models, which are potential candidates for automatic segmentation of thermal ablation regions in US during clinical ablation treatment.

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.

Three-dimensional echo decorrelation monitoring of radiofrequency ablation in ex vivo bovine liver

Three-dimensional (3D) echo decorrelation imaging was investigated for monitoring radiofrequency ablation (RFA) in ex vivo bovine liver. RFA experiments (N = 14) were imaged by 3D ultrasound using a matrix array, with in-phase and quadrature complex echo volumes acquired about every 11 s. Tissue specimens were then frozen at -80 °C, sectioned, and semi-automatically segmented. Receiver operating characteristic (ROC) curves were constructed for assessing ablation prediction performance of 3D echo decorrelation with three potential normalization approaches, as well as 3D integrated backscatter (IBS). ROC analysis indicated that 3D echo decorrelation imaging is potentially a good predictor of local RFA, with the best prediction performance observed for globally normalized decorrelation. Tissue temperatures, recorded by four thermocouples integrated into the RFA probe, showed good correspondence with spatially averaged decorrelation and statistically significant but weak correlation with measured echo decorrelation at the same spatial locations. In tests predicting ablation zones using a weighted K-means clustering approach, echo decorrelation performed better than IBS, with smaller root mean square volume errors and higher Dice coefficients relative to measured ablation zones. These results suggest that 3D echo decorrelation and IBS imaging are capable of real-time monitoring of thermal ablation, with potential application to clinical treatment of liver tumors.

A Review of Imaging Methods to Assess Ultrasound-Mediated Ablation

Ultrasound ablation techniques are minimally invasive alternatives to surgical resection and have rapidly increased in use. The response of tissue to HIFU ablation differs based on the relative contributions of thermal and mechanical effects, which can be varied to achieve optimal ablation parameters for a given tissue type and location. In tumor ablation, similar to surgical resection, it is desirable to include a safety margin of ablated tissue around the entirety of the tumor. A factor in optimizing ablative techniques is minimizing the recurrence rate, which can be due to incomplete ablation of the target tissue. Further, combining focal ablation with immunotherapy is likely to be key for effective treatment of metastatic cancer, and therefore characterizing the impact of ablation on the tumor microenvironment will be important. Thus, visualization and quantification of the extent of ablation is an integral component of ablative procedures. The aim of this review article is to describe the radiological findings after ultrasound ablation across multiple imaging modalities. This review presents readers with a general overview of the current and emerging imaging methods to assess the efficacy of ultrasound ablative treatments.

Safety of Irreversible Electroporation Ablation of the Pancreas

OBJECTIVES

To evaluate the safety of irreversible electroporation (IRE) on swine pancreatic tissue including its effects on peripancreatic vessels, bile ducts, and bowel.

METHODS

Eighteen Yorkshire pigs underwent IRE ablation of the pancreas successfully and without clinical complications. Contrast-enhanced computed tomography angiography and laboratory studies before the IRE ablation with follow-up computed tomography angiography, laboratory testing, and pathological examination up to 4 weeks postablation were performed.

RESULTS

In a subset of cases, anatomic peripancreatic vessel narrowing was seen by 1 week postablation, persisting at 4 weeks postablation, without apparent functional impairment of blood flow. Laboratory studies revealed elevated amylase and lipase at 24 hours post-IRE, suggestive of acute pancreatitis, which normalized by 4 weeks post-IRE. There was extensive pancreatic tissue damage 24 hours after IRE with infiltration of immune cells, which was gradually replaced by fibrotic tissue. Ductal regeneration without loss of pancreatic acinar tissue and glandular function was observed at 1 and 4 weeks postablation.

CONCLUSIONS

In our study, we demonstrated and confirmed the safety and minimal complications of IRE ablation in the pancreas and its surrounding vital structures. These results show the potential of IRE as an alternative treatment modality in patients with pancreatic cancer, especially those with locally advanced disease.

Effect of Overpressure on Acoustic Emissions and Treated Tissue Histology in ex Vivo Bulk Ultrasound Ablation

Bulk ultrasound ablation is a thermal therapy approach in which tissue is heated by unfocused or weakly focused sonication (average intensities on the order of 100 W/cm²) to achieve coagulative necrosis within a few minutes exposure time. Assessing the role of bubble activity, including acoustic cavitation and tissue vaporization, in bulk ultrasound ablation may help in making bulk ultrasound ablation safer and more effective for clinical applications. Here, two series of ex vivo ablation trials were conducted to investigate the role of bubble activity and tissue vaporization in bulk ultrasound ablation. Fresh bovine liver tissue was ablated with unfocused, continuous-wave ultrasound using ultrasound image-ablate arrays sonicating at 31 W/cm² (0.9 MPa amplitude) for either 20 min at a frequency of 3.1 MHz or 10 min at 4.8 MHz. Tissue specimens were maintained at a static overpressure of either 0.52 or 1.2 MPa to suppress bubble activity and tissue vaporization or at atmospheric pressure for control groups. A passive cavitation detector was used to record subharmonic (1.55 or 2.4 MHz), broadband (1.2-1.5 MHz) and low-frequency (5-20 kHz) acoustic emissions. Treated tissue was stained with 2% triphenyl tetrazolium chloride to evaluate thermal lesion dimensions. Subharmonic emissions were significantly reduced in overpressure groups compared with control groups. Correlations observed between acoustic emissions and lesion dimensions were significant and positive for the 3.1-MHz series, but significant and negative for the 4.8-MHz series. The results indicate that for bulk ultrasound ablation, where both acoustic cavitation and tissue vaporization are possible, bubble activity can enhance ablation in the absence of tissue vaporization, but can reduce thermal lesion dimensions in the presence of vaporization.

In vivo ultrasound thermal ablation control using echo decorrelation imaging in rabbit liver and VX2 tumor

The utility of echo decorrelation imaging feedback for real-time control of in vivo ultrasound thermal ablation was assessed in rabbit liver with VX2 tumor. High-intensity focused ultrasound (HIFU) and unfocused (bulk) ablation were performed using 5 MHz linear image-ablate arrays. Treatments comprised up to nine lower-power sonications, followed by up to nine higher-power sonications, ceasing when the average cumulative echo decorrelation within a control region of interest exceeded a predefined threshold (- 2.3, log10-scaled echo decorrelation per millisecond, corresponding to 90% specificity for tumor ablation prediction in previous in vivo experiments). This threshold was exceeded in all cases for both HIFU (N = 12) and bulk (N = 8) ablation. Controlled HIFU trials achieved a significantly higher average ablation rate compared to comparable ablation trials without image-based control, reported previously. Both controlled HIFU and bulk ablation trials required significantly less treatment time than these previous uncontrolled trials. Prediction of local liver and VX2 tumor ablation using echo decorrelation was tested using receiver operator characteristic curve analysis, showing prediction capability statistically equivalent to uncontrolled trials. Compared to uncontrolled trials, controlled trials resulted in smaller thermal ablation regions and higher contrast between echo decorrelation in treated vs. untreated regions. These results indicate that control using echo decorrelation imaging may reduce treatment duration and increase treatment reliability for in vivo thermal ablation.

Application of Ultrasound Thermal Imaging for Monitoring Laser Ablation in Ex Vivo Cardiac Tissue

BACKGROUND AND OBJECTIVE

Laser ablation can be used to treat atrial fibrillation by thermally isolating pulmonary veins. In this study, we evaluated the feasibility of high-resolution (<1 mm) ultrasound thermal imaging to monitor spatial temperature distribution during laser ablation on ex vivo cardiac tissue.

STUDY DESIGN/MATERIALS AND METHODS

Laser ablation (808 nm) was performed on five porcine cardiac tissue samples. A thermocouple was used to measure the interstitial tissue temperature during the laser ablation process. Tissue-strain-based ultrasound thermal imaging was conducted to monitor the spatial distribution of the temperature in the cardiac tissue. The tissue temperature was estimated from the time shifts of ultrasound signals owing to the changes in the speed of sound and was compared with the measured temperature. The temperature estimation coefficient k of porcine cardiac tissue was calculated from the estimated thermal strain and the measured temperature. The degree of tissue coagulation (temperatures > 50°C) was derived from the estimated temperature and was compared with that of the tested cardiac tissue.

RESULTS

The estimated tissue temperature using strain-based ultrasound thermal imaging at a depth of 1 mm agreed with thermocouple measurements. During the 30-second period of the laser ablation process, the estimated tissue temperature increased from 25 to 70°C at a depth of 0.1 mm, while the estimated temperature at a depth of 1 mm increased up to 46°C. Owing to the uncertainty of the coefficient k, the k value of the porcine cardiac tissue varied from 160 to 220°C with temperature changes of up to 20°C. The estimated coagulation region in the ultrasound thermal imaging was 20% wider (+0.6 mm) but 9% shallower (-0.1 mm) than the measured region of the ablated porcine cardiac tissue.

CONCLUSIONS

The current study demonstrated the feasibility of temperature monitoring with the use of ultrasound thermal imaging during the laser ablation on ex vivo porcine cardiac tissue. The high-resolution ultrasound thermal imaging could map the spatial distribution of the tissue temperature. The proposed method can be used to monitor the temperature and thermal coagulation to achieve effective laser ablation for atrial fibrillation. Lasers Surg. Med. © 2019 Wiley Periodicals, Inc.

Monitoring Microwave Ablation Using Ultrasound Echo Decorrelation Imaging: An ex vivo Study

In this study, a microwave-induced ablation zone (thermal lesion) monitoring method based on ultrasound echo decorrelation imaging was proposed. A total of 15 cases of ex vivo porcine liver microwave ablation (MWA) experiments were carried out. Ultrasound radiofrequency (RF) signals at different times during MWA were acquired using a commercial clinical ultrasound scanner with a 7.5-MHz linear-array transducer. Instantaneous and cumulative echo decorrelation images of two adjacent frames of RF data were calculated. Polynomial approximation images were obtained on the basis of the thresholded cumulative echo decorrelation images. Experimental results showed that the instantaneous echo decorrelation images outperformed conventional B-mode images in monitoring microwave-induced thermal lesions. Using gross pathology measurements as the reference standard, the estimation of thermal lesions using the polynomial approximation images yielded an average accuracy of 88.60%. We concluded that instantaneous ultrasound echo decorrelation imaging is capable of monitoring the change of thermal lesions during MWA, and cumulative ultrasound echo decorrelation imaging and polynomial approximation imaging are feasible for quantitatively depicting thermal lesions.

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).

Optimized Echo Decorrelation Imaging Feedback for Bulk Ultrasound Ablation Control

Feasibility of controlling bulk ultrasound (US) thermal ablation using echo decorrelation imaging was investigated in ex vivo bovine liver. The first of two ablation and control procedures used a sequence of constant-intensity sonication cycles, ceased when the minimum echo decorrelation within a control region of interest (ROI) exceeded a predetermined threshold. The second procedure used a variable-intensity sonication sequence, with spatially averaged decorrelation as the stopping criterion. US exposures and echo decorrelation imaging were performed by a linear image-ablate array. Based on preliminary experiments, control ROIs and thresholds for the minimum-decorrelation and average-decorrelation criteria were specified. Controlled trials for the minimum-decorrelation and average-decorrelation criteria were compared with uncontrolled trials employing 9 or 18 cycles of matching sonication sequences. Lesion dimensions, treatment times, ablation rates, and areas under receiver operating characteristic curves were statistically compared. Successfully controlled trials using both criteria required significantly shorter treatment times than corresponding 18-cycle treatments, with better ablation prediction performance than uncontrolled 9-cycle and 18-cycle treatments. Either control approach resulted in greater ablation rate than corresponding 9-cycle or 18-cycle uncontrolled approaches. A post hoc analysis studied the effect of exchanging control criteria between the two series of controlled experiments. For either group, the average time needed to exceed the alternative decorrelation threshold approximately matched the average duration of successfully controlled experimental trials. These results indicate that either approach, using minimum-decorrelation or average-decorrelation criteria, is feasible for control of bulk US ablation. In addition, use of a variable-intensity sonication sequence for bulk US thermal ablation can result in larger ablated regions compared to constant-intensity sonication sequences.

Real-Time Spatiotemporal Control of High-Intensity Focused Ultrasound Thermal Ablation Using Echo Decorrelation Imaging in ex Vivo Bovine Liver

The ability to control high-intensity focused ultrasound (HIFU) thermal ablation using echo decorrelation imaging feedback was evaluated in ex vivo bovine liver. Sonications were automatically ceased when the minimum cumulative echo decorrelation within the region of interest exceeded an ablation control threshold, determined from preliminary experiments as -2.7 (log-scaled decorrelation per millisecond), corresponding to 90% specificity for local ablation prediction. Controlled HIFU thermal ablation experiments were compared with uncontrolled experiments employing two, five or nine sonication cycles. Means and standard errors of the lesion width, area and depth, as well as receiver operating characteristic curves testing ablation prediction performance, were computed for each group. Controlled trials exhibited significantly smaller average lesion area, width and treatment time than five-cycle or nine-cycle uncontrolled trials and also had significantly greater prediction capability than two-cycle uncontrolled trials. These results suggest echo decorrelation imaging is an effective approach to real-time HIFU ablation control.

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.

Echo Decorrelation Imaging of Rabbit Liver and VX2 Tumor during In Vivo Ultrasound Ablation

In open surgical procedures, image-ablate ultrasound arrays performed thermal ablation and imaging on rabbit liver lobes with implanted VX2 tumor. Treatments included unfocused (bulk ultrasound ablation, N = 10) and focused (high-intensity focused ultrasound ablation, N = 13) exposure conditions. Echo decorrelation and integrated backscatter images were formed from pulse-echo data recorded during rest periods after each therapy pulse. Echo decorrelation images were corrected for artifacts using decorrelation measured prior to ablation. Ablation prediction performance was assessed using receiver operating characteristic curves. Results revealed significantly increased echo decorrelation and integrated backscatter in both ablated liver and ablated tumor relative to unablated tissue, with larger differences observed in liver than in tumor. For receiver operating characteristic curves computed from all ablation exposures, both echo decorrelation and integrated backscatter predicted liver and tumor ablation with statistically significant success, and echo decorrelation was significantly better as a predictor of liver ablation. These results indicate echo decorrelation imaging is a successful predictor of local thermal ablation in both normal liver and tumor tissue, with potential for real-time therapy monitoring.

Detection of tissue coagulation by decorrelation of ultrasonic echo signals in cavitation-enhanced high-intensity focused ultrasound treatment

BACKGROUND

A noninvasive technique to monitor thermal lesion formation is necessary to ensure the accuracy and safety of high-intensity focused ultrasound (HIFU) treatment. The purpose of this study is to ultrasonically detect the tissue change due to thermal coagulation in the HIFU treatment enhanced by cavitation microbubbles.

METHODS

An ultrasound imaging probe transmitted plane waves at a center frequency of 4.5 MHz. Ultrasonic radio-frequency (RF) echo signals during HIFU exposure at a frequency of 1.2 MHz were acquired. Cross-correlation coefficients were calculated between in-phase and quadrature (IQ) data of two B-mode images with an interval time of 50 and 500 ms for the estimation of the region of cavitation and coagulation, respectively. Pathological examination of the coagulated tissue was also performed to compare with the corresponding ultrasonically detected coagulation region.

RESULTS

The distribution of minimum hold cross-correlation coefficient between two sets of IQ data with 50-ms intervals was compared with a pulse inversion (PI) image. The regions with low cross-correlation coefficients approximately corresponded to those with high brightness in the PI image. The regions with low cross-correlation coefficients in 500-ms intervals showed a good agreement with those with significant change in histology.

CONCLUSIONS

The results show that the regions of coagulation and cavitation could be ultrasonically detected as those with low cross-correlation coefficients between RF frames with certain intervals. This method will contribute to improve the safety and accuracy of the HIFU treatment enhanced by cavitation microbubbles.

Dependence of ultrasound echo decorrelation on local tissue temperature during ex vivo radiofrequency ablation

This study investigates echo decorrelation imaging, an ultrasound method for thermal ablation monitoring. The effect of tissue temperature on the mapped echo decorrelation parameter was assessed in radiofrequency ablation experiments performed on ex vivo bovine liver tissue. Echo decorrelation maps were compared with corresponding tissue temperatures simulated using the finite element method. For both echo decorrelation imaging and integrated backscatter imaging, the mapped tissue parameters correlated significantly but weakly with local tissue temperature. Receiver operating characteristic (ROC) curves were used to assess the ability of echo decorrelation and integrated backscatter to predict tissue temperature greater than 40, 60, and 80 °C. Significantly higher area under the ROC curve (AUROC) values were obtained for prediction of tissue temperatures greater than 40, 60, and 80 °C using echo decorrelation imaging (AUROC = 0.871, 0.948 and 0.966) compared to integrated backscatter imaging (AUROC = 0.865, 0.877 and 0.832).

Optimization of tissue physical parameters for accurate temperature estimation from finite-element simulation of radiofrequency ablation

Computational finite element models are commonly used for the simulation of radiofrequency ablation (RFA) treatments. However, the accuracy of these simulations is limited by the lack of precise knowledge of tissue parameters. In this technical note, an inverse solver based on the unscented Kalman filter (UKF) is proposed to optimize values for specific heat, thermal conductivity, and electrical conductivity resulting in accurately simulated temperature elevations. A total of 15 RFA treatments were performed on ex vivo bovine liver tissue. For each RFA treatment, 15 finite-element simulations were performed using a set of deterministically chosen tissue parameters to estimate the mean and variance of the resulting tissue ablation. The UKF was implemented as an inverse solver to recover the specific heat, thermal conductivity, and electrical conductivity corresponding to the measured area of the ablated tissue region, as determined from gross tissue histology. These tissue parameters were then employed in the finite element model to simulate the position- and time-dependent tissue temperature. Results show good agreement between simulated and measured temperature.

Catheter-based ultrasound technology for image-guided thermal therapy: current technology and applications

Catheter-based ultrasound (CBUS) is applied to deliver minimally invasive thermal therapy to solid cancer tumours, benign tissue growth, vascular disease, and tissue remodelling. Compared to other energy modalities used in catheter-based surgical interventions, unique features of ultrasound result in conformable and precise energy delivery with high selectivity, fast treatment times, and larger treatment volumes. We present a concise review of CBUS technology being currently utilized in animal and clinical studies or being developed for future applications. CBUS devices have been categorised into interstitial, endoluminal and endovascular/cardiac applications. Basic applicator designs, site-specific evaluations and possible treatment applications have been discussed in brief. Particular emphasis has been given to ablation studies that incorporate image guidance for applicator placement, therapy monitoring, feedback control, and post-procedure assessment. Examples of devices included here span the entire spectrum of the development cycle from preliminary simulation-based design studies to implementation in clinical investigations. The use of CBUS under image guidance has the potential for significantly improving precision and applicability of thermal therapy delivery.

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.

Analysis of tissue changes, measurement system effects, and motion artifacts in echo decorrelation imaging

Echo decorrelation imaging, a method for mapping ablation-induced ultrasound echo changes, is analyzed. Local echo decorrelation is shown to approximate the decoherence spectrum of tissue reflectivity. Effects of the ultrasound measurement system, echo signal windowing, electronic noise, and tissue motion on echo decorrelation images are determined theoretically, leading to a method for reduction of motion and noise artifacts. Theoretical analysis is validated by simulations and experiments. Simulated decoherence of the scattering medium was recovered with root-mean-square error less than 10% with accuracy dependent on the correlation window size. Motion-induced decorrelation measured in an ex vivo pubovisceral muscle model showed similar trends to theoretical motion-induced decorrelation for a 2.1 MHz curvilinear array with decorrelation approaching unity for 3-4 mm elevational displacement or 1-1.6 mm range displacement. For in vivo imaging of porcine liver by a 7 MHz linear array, theoretical decorrelation computed using image-based motion estimates correlated significantly with measured decorrelation (r = 0.931, N = 10). Echo decorrelation artifacts incurred during in vivo radiofrequency ablation in the same porcine liver were effectively compensated based on the theoretical echo decorrelation model and measured pre-treatment decorrelation. These results demonstrate the potential of echo decorrelation imaging for quantification of heat-induced changes to the scattering tissue medium during thermal ablation.

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.

Prediction and suppression of HIFU-induced vessel rupture using passive cavitation detection in an ex vivo model

Background

Occlusion of blood vessels using high-intensity focused ultrasound (HIFU) is a potential treatment for arteriovenous malformations and other neurovascular disorders. However, attempting HIFU-induced vessel occlusion can also cause vessel rupture, resulting in hemorrhage. Possible rupture mechanisms include mechanical effects of acoustic cavitation and heating of the vessel wall.

Methods

HIFU exposures were performed on 18 ex vivo porcine femoral arteries with simultaneous passive cavitation detection. Vessels were insonified by a 3.3-MHz focused source with spatial-peak, temporal-peak focal intensity of 15,690–24,430 W/cm² (peak negative-pressure range 10.92–12.52 MPa) and a 50% duty cycle for durations up to 5 min. Time-dependent acoustic emissions were recorded by an unfocused passive cavitation detector and quantified within low-frequency (10–30 kHz), broadband (0.3–1.1 MHz), and subharmonic (1.65 MHz) bands. Vessel rupture was detected by inline metering of saline flow, recorded throughout each treatment. Recorded emissions were grouped into ‘pre-rupture’ (0–10 s prior to measured point of vessel rupture) and ‘intact-vessel’ (>10 s prior to measured point of vessel rupture) emissions. Receiver operating characteristic curve analysis was used to assess the ability of emissions within each frequency band to predict vessel rupture.

Based on these measurements associating acoustic emissions with vessel rupture, a real-time feedback control module was implemented to monitor acoustic emissions during HIFU treatment and adjust the ultrasound intensity, with the goal of maximizing acoustic power delivered to the vessel while avoiding rupture. This feedback control approach was tested on 10 paired HIFU exposures of porcine femoral and subclavian arteries, in which the focal intensity was stepwise increased from 9,117 W/cm² spatial-peak temporal-peak (SPTP) to a maximum of 21,980 W/cm², with power modulated based on the measured subharmonic emission amplitude. Time to rupture was compared between these feedback-controlled trials and paired controller-inactive trials using a paired Wilcoxon signed-rank test.

Results

Subharmonic emissions were found to be the most predictive of vessel rupture (areas under the receiver operating characteristic curve (AUROC) = 0.757, p < 10^-16) compared to low-frequency (AUROC = 0.657, p < 10^-11) and broadband (AUROC = 0.729, p < 10^-16) emissions. An independent-sample t test comparing pre-rupture to intact-vessel emissions revealed a statistically significant difference between the two groups for broadband and subharmonic emissions (p < 10^-3), but not for low-frequency emissions (p = 0.058).

In a one-sided paired Wilcoxon signed-rank test, activation of the control module was shown to increase the time to vessel rupture (T_- = 8, p = 0.0244, N = 10). In one-sided paired t tests, activation of the control module was shown to cause no significant difference in time-averaged focal intensity (t = 0.362, p = 0.363, N = 10), but was shown to cause delivery of significantly greater total acoustic energy (t = 2.037, p = 0.0361, N = 10).

Conclusions

These results suggest that acoustic cavitation plays an important role in HIFU-induced vessel rupture. In HIFU treatments for vessel occlusion, passive monitoring of acoustic emissions may be useful in avoiding hemorrhage due to vessel rupture, as shown in the rupture suppression experiments.

Considering angle selection when using ultrasound electrode displacement elastography to evaluate radiofrequency ablation of tissues

Percutaneous radiofrequency ablation (RFA) is a minimally invasive treatment to thermally destroy tumors. Ultrasound-based electrode-displacement elastography is an emerging technique for evaluating the region of RFA-induced lesions. The angle between the imaging probe and the RFA electrode can influence electrode-displacement elastography when visualizing the ablation zone. We explored the angle effect on electrode-displacement elastography to measure the ablation zone. Phantoms embedded with meatballs were fabricated and then ablated using an RFA system to simulate RFA-induced lesions. For each phantom, a commercial ultrasound scanner with a 7.5 MHz linear probe was used to acquire raw image data at different angles, ranging from 30° to 90° at increments of 10°, to construct electrode-displacement images and facilitate comparisons with tissue section images. The results revealed that the ablation regions detected using electrode-displacement elastography were highly correlated with those from tissue section images when the angle was between 30° and 60°. However, the boundaries of lesions were difficult to distinguish, when the angle was larger than 60°. The experimental findings suggest that angle selection should be considered to achieve reliable electrode-displacement elastography to describe ablation zones.

In vivo thermal ablation monitoring using ultrasound echo decorrelation imaging

Previous work indicated that ultrasound echo decorrelation imaging can track and quantify changes in echo signals to predict thermal damage during in vitro radiofrequency ablation (RFA). In the in vivo studies reported here, the feasibility of using echo decorrelation imaging as a treatment monitoring tool was assessed. RFA was performed on normal swine liver (N = 5), and ultrasound ablation using image-ablate arrays was performed on rabbit liver implanted with VX2 tumors (N = 2). Echo decorrelation and integrated backscatter were computed from Hilbert transformed pulse-echo data acquired during RFA and ultrasound ablation treatments. Receiver operating characteristic (ROC) curves were employed to assess the ability of echo decorrelation imaging and integrated backscatter to predict ablation. Area under the ROC curves (AUROC) was determined for RFA and ultrasound ablation using echo decorrelation imaging. Ablation was predicted more accurately using echo decorrelation imaging (AUROC = 0.832 and 0.776 for RFA and ultrasound ablation, respectively) than using integrated backscatter (AUROC = 0.734 and 0.494).

Real-time monitoring of HIFU treatment using pulse inversion

Ultrasound (US) imaging is widely used for the real-time guidance of high-intensity focused ultrasound (HIFU) treatment at a relatively low cost. However, ultrasound image guided HIFU (USgHIFU) is limited in the real-time monitoring of HIFU treatment due to the large amplitude HIFU signals received by the US imaging transducer. The amplitude of the HIFU scattered signal is generally much higher than the amplitude of the pulse-echo signal received by the imaging transducer. This creates an interference pattern obscuring the image of the tissue. As such, it is difficult to monitor lesion location. This paper proposes a real-time monitoring method to be performed concurrently with the HIFU insonation, but without HIFU interference, which allows for the improvement of treatment accuracy and safety in USgHIFU. The proposed method utilizes the physical properties of pulse inversion which is capable of removing the fundamental and odd harmonic components of the HIFU interference. Therefore, it is possible to secure the desired spectral bandwidth used to construct US images for HIFU treatment monitoring. The performance of the proposed method was evaluated through experiments with both a bovine serum albumin phantom and a chicken breast. The results demonstrated that the proposed method is capable of providing interference-free US images, thus successfully allowing for US imaging during HIFU treatment.

2pBAb5. Validation of three-dimensional strain tracking by volumetric ultrasound image correlation in a pubovisceral muscle model

Little is understood about the biomechanical changes leading to pelvic floor disorders such as stress urinary incontinence. In order to measure regional biomechanical properties of the pelvic floor muscles in vivo, a three dimensional (3D) strain tracking technique employing correlation of volumetric ultrasound images has been implemented. In this technique, local 3D displacements are determined as a function of applied stress and then converted to strain maps. To validate this approach, an in vitro model of the pubovisceral muscle, with a hemispherical indenter emulating the downward stress caused by intra-abdominal pressure, was constructed. Volumetric B-scan images were recorded as a function of indenter displacement while muscle strain was measured independently by a sonomicrometry system (Sonometrics). Local strains were computed by ultrasound image correlation and compared with sonomicrometry-measured strains to assess strain tracking accuracy. Image correlation by maximizing an exponential likelihood function was found more reliable than the Pearson correlation coefficient. Strain accuracy was dependent on sizes of the subvolumes used for image correlation, relative to characteristic speckle length scales of the images. Decorrelation of echo signals was mapped as a function of indenter displacement and local tissue orientation. Strain measurement accuracy was weakly related to local echo decorrelation.

High-frequency rapid B-mode ultrasound imaging for real-time monitoring of lesion formation and gas body activity during high-intensity focused ultrasound ablation

The goal of this study was to examine the ability of high-frame-rate, high-resolution imaging to monitor tissue necrosis and gas-body activities formed during high-intensity focused ultrasound (HIFU) application. Ex vivo porcine cardiac tissue specimens (n = 24) were treated with HIFU exposure (4.33 MHz, 77 to 130 Hz pulse repetition frequency (PRF), 25 to 50% duty cycle, 0.2 to 1 s, 2600 W/cm(2)). RF data from B-mode ultrasound imaging were obtained before, during, and after HIFU exposure at a frame rate ranging from 77 to 130 Hz using an ultrasound imaging system with a center frequency of 55 MHz. The time history of changes in the integrated backscatter (IBS), calibrated spectral parameters, and echo-decorrelation parameters of the RF data were assessed for lesion identification by comparison against gross sections. Temporal maximum IBS with +12 dB threshold achieved the best identification with a receiver-operating characteristic (ROC) curve area of 0.96. Frame-to-frame echo decorrelation identified and tracked transient gas-body activities. Macroscopic (millimeter-sized) cavities formed when the estimated initial expansion rate of gas bodies (rate of expansion in lateral-to-beam direction) crossed 0.8 mm/s. Together, these assessments provide a method for monitoring spatiotemporal evolution of lesion and gas-body activity and for predicting macroscopic cavity formation.

High-frequency ultrasound m-mode imaging for identifying lesion and bubble activity during high-intensity focused ultrasound ablation

Effective real-time monitoring of high-intensity focused ultrasound (HIFU) ablation is important for application of HIFU technology in interventional electrophysiology. This study investigated rapid, high-frequency M-mode ultrasound imaging for monitoring spatiotemporal changes during HIFU application. HIFU (4.33 MHz, 1 kHz PRF, 50% duty cycle, 1 s, 2600‒6100 W/cm²) was applied to ex vivo porcine cardiac tissue specimens with a confocally and perpendicularly aligned high-frequency imaging system (Visualsonics Vevo 770, 55 MHz center frequency). Radio-frequency (RF) data from M-mode imaging (1 kHz PRF, 2 s × 7 mm) was acquired before, during and after HIFU treatment (n = 12). Among several strategies, the temporal maximum integrated backscatter with a threshold of +12 dB change showed the best results for identifying final lesion width (receiver-operating characteristic curve area 0.91 ± 0.04, accuracy 85 ± 8%, compared with macroscopic images of lesions). A criterion based on a line-to-line decorrelation coefficient is proposed for identification of transient gas bodies.

Influence of temperature variations on the average grayscale of B-mode images of HIFU-induced lesions

The aim of this work is to understand how the variations in grayscale values of B-mode ultrasound (US) images can be used as an approach for non-invasive temperature estimation. In order to obtain real-time monitoring of HIFU treatment, an US imaging system and HIFU were synchronized. Images were acquired using an electronic convex imaging probe. An 8% BSA tissue-mimicking polyacrylamide gel was used for the experiments. First, the HIFU power was set to 10 W. The application of HIFU resulted in the appearance of hyperechoic regions that were used to place a thermocouple tip at the focal spot by means of US imaging. Afterwards, the power was gradually increased up to 40 W for 4 min. The temperature sampling rate was set to 5 Hz. For each temperature sample the synchronization system captured one frame. The region of interest (ROI) was manually selected and a set of morphological operations were implemented in order to obtain the shape and size of the thermal lesion. From here, the average grayscale (AGS) and area of the thermal lesion were calculated to assess temperature quantification. The AGS parameter showed a maximum correlation coefficient of 0.6626 as a function of temperature whereas the thermal lesion appeared. In contrast, area values showed a greater correlation coefficient of 0.9122. In conclusion, temperature shows a non-linear behavior with respect to the parameters estimated due to the nature of the thermal lesion formation by HIFU exposure.

Advanced hepatic ablation technique for creating complete cell death: irreversible electroporation

PURPOSE

To evaluate the effectiveness of irreversible electroporation (IRE) in hepatic tissue ablation and the radiologic-pathologic correlation of IRE-induced cell death.

MATERIALS AND METHODS

On approval of the animal research committee, 16 Yorkshire pigs underwent ultrasonography (US)-guided IRE of normal liver. A total of 55 ablation zones were created, which were imaged with US, magnetic resonance (MR) imaging, and computed tomography (CT) and evaluated with immunohistochemical analysis, including hematoxylin-eosin (H-E), Von Kossa, and von Willibrand factor (vWF) staining and terminal deoxynucleotidyl transferase dUTP nick end labeling (TUNEL) assay.

RESULTS

At gross section examination, the mean diameter of the ablation zones was 33.5 mm + or - 3.0 (standard deviation) and was achieved in 6.9 minutes (mean total procedure time per ablation), with a mean difference of 2.5 mm + or - 3.6 between US and gross section measurements (r = 0.804). No complications were seen in any of the 16 animals. IRE ablation zones were well characterized with US, CT, and MR imaging, and real-time monitoring was feasible with US. H-E, Von Kossa, and vWF staining showed complete cell death, with a sharply demarcated treatment area. Bile ducts and vessels were completely preserved. Areas of complete cell death were stained positive for apoptotic markers (TUNEL, BCL-2 oncoprotein), suggesting involvement of the apoptotic process in the pathophysiology of cell death caused by IRE.

CONCLUSION

In an animal model, IRE proved to be a fast, safe, and potent ablative method, causing complete tissue death by means of apoptosis. Cell death is seen with full preservation of periablative zone structures, including blood vessels, bile ducts, and neighboring nonablated tissues.

MR-guided liver tumor ablation employing open high-field 1.0T MRI for image-guided brachytherapy

OBJECTIVE

To determine the feasibility and safety of image-guided brachytherapy employing a modified open high-field MR system.

METHODS

This is a follow-up study of a development project enabling technologies for interventional use of 1.0T open MRI. Modifications included coils and in-bore visualization, fluoroscopic sequences and user interfaces. We recruited 104 patients with 224 liver malignancies to receive MR-guided brachytherapy. Interventions were performed >20 min after Gd-EOB-DTPA. We recorded interventional parameters including the intervention time (from acquisition of the first scout until the final sequence for brachytherapy treatment planning). Two reviewers assessed MR-fluoroscopic images in comparison to plain CT as used in CT intervention, applying a rating scale of 1-10. Statistical analysis included Friedman and Kendall's W tests.

RESULTS

We employed freehand puncture with interactive dynamic imaging for navigation. Technical success rate was 218 complete ablations in 224 tumours (97%). The median intervention time was 61 min. We recorded no adverse events related to the use of MRI. No major complications occurred. The rate of minor complications was 4%. Local control at 3 months was 96%. Superiority of MR-fluoroscopic, Gd-EOB-DTPA-enhanced images over plain CT was highly significant (P < 0.001).

CONCLUSION

MR-guided brachytherapy employing open high-field MRI is feasible and safe.