Ultrasound temperature estimation based on probability variation of backscatter data

Scatterer spacing based on Gabor atoms matched from harmonic ultrasound echoes to improve assessment of microwave-induced thermal lesions

BACKGROUND

Microwave ablation (MWA) is a minimally invasive alternative for the treatment of unresectable liver tumors. To verify the effectiveness and safety of MWA, it is critical to measure the temperature variation and assess the regions of the microwave-induced thermal lesions.

PURPOSE

Recent studies have indicated that the locations of optimally matched Gabor atoms (LOMGA) from ultrasound radiofrequency (RF) echo signals allow accurate and stable scatterer spacing estimation. Herein, a harmonic-based LOMGA method is proposed to estimate the scatterer spacing for improving the assessment of microwave-induced thermal lesions.

METHODS

The mean scatterer spacing (MSS) is estimated via the LOMGA method incorporating the selection of concise atoms from separated second-harmonic RF echo signals with the pulse-inversion algorithm for thermal lesion evaluation. In vitro experiments, 10 fresh porcine liver samples were ablated at different time nodes during the ablation period, and 200 sets of second-harmonic and fundamental RF echo signals were randomly selected from the regions of interest in the coagulated liver samples for MSS estimation. The means and standard deviations of the MSSs, as well as the linear regression for the mean MSSs, were calculated from fundamental and second-harmonic signals for comparison and evaluation, the receiver operating characteristic (ROC) curves for the 200 sets of fundamental-based and harmonic-based MSS estimates from the 10 liver samples at five pairs of adjacent time nodes were calculated, and one-way analysis of variance (ANOVA) tests were performed for the five pairs of adjacent time nodes. The fundamental and harmonic-based p-values in the ANOVA tests and the areas under the ROC curves (AUCs) were calculated to statistically analyze the differences in the MSSs between adjacent time nodes.

RESULTS

The harmonic-based increments in the intensity variation and coherent components were larger than the fundamental-based increments with the increasing ablation time. The harmonic-based MSSs from the 10 liver samples at five pairs of adjacent time nodes were found to be highly statistically significant (p < 0.01). Thus, the harmonic-based MSSs had greater variations. Compared with the fundamental-based results, for the five preset ST values, the average increment in the harmonic-based mean slopes was 69.22% and the average decrement in the mean standard deviations was 11.67% for the linear-fitting MSS results, and the results were statistically significant (p < 0.05).

CONCLUSION

Harmonic-based MSSs are more sensitive and robust to variations in coagulated tissues, which is advantageous for the assessment of microwave-induced thermal lesions.

Appropriate indicator of modeling error for threshold-based model selection in statistics-based ultrasound tissue characterization

Analysis of the envelope statistics of ultrasound echo signals contributes to quantitative tissue characterization in medical ultrasound. Many probability distribution model functions have been studied, and the model function that should be used for tissue characterization depends on the type of disease, even in the same organ. Thus, an appropriate model selection is important for an accurate diagnosis. In this study, we aimed to select a model using threshold processing for modeling errors instead of a simple selection by minimizing the modeling error. For this purpose, we compared several indicators of modeling errors using random number simulations, ultrasonic simulation, and phantom experiment. The results validated that the Mahalanobis distance of moments is an appropriate indicator because it enables the use of a constant threshold value, regardless of the type of model function and data length.

Temperature Monitoring During Microwave Hyperthermia Based on BP-Nakagami Distribution

OBJECTIVE

The purpose of this study is to accurately monitor temperature during microwave hyperthermia. We propose a temperature estimation model BP-Nakagami based on neural network for Nakagami distribution.

METHODS

In this work, we designed the microwave hyperthermia experiment of fresh ex vivo pork tissue and phantom, collected ultrasonic backscatter data at different temperatures, modeled these data using Nakagami distribution, and calculated Nakagami distribution parameter m. A neural network model was built to train the relationship between Nakagami distribution parameter m and temperature, and a BP-Nakagami temperature model with good fitting was obtained. The temperature model is used to draw the two-dimensional temperature distribution map of biological tissues in microwave hyperthermia. Finally, the temperature estimated by the model is compared with the temperature measured by thermocouples.

RESULTS

The error between the temperature estimated by the temperature model and the temperature measured by the thermocouple is within 1°C in the range of 25°C-50°C for ex vivo pork tissue, and the error between the temperature estimated by the temperature model and the temperature measured by the thermocouple is within 0.5°C in the range of 25°C-50°C for phantom.

CONCLUSIONS

The results show that the temperature estimation model proposed by us is an effective model for monitoring the internal temperature change of biological tissues.

Quantification of limitations in statistical analysis of ultrasound echo envelope amplitudes

Ultrasound echo envelope statistics have been widely studied for quantitative tissue characterization. In ultrasound measurements, the size of the region-of-interest (ROI) is limited by several factors, such as the locality of the tissue characteristics and the depth dependence of the acoustic field of the ultrasound beam. In this case, the evaluated echo envelope statistics vary even when the envelope amplitudes follow the same population without any noise. In this study, the statistical variance of the moments caused by this finite number of samples was quantified based on the central limit theorem and the law of error propagation. The proposed principles were validated by random number simulation and used to quantify the statistical variance of Nakagami parameter estimation. Finally, the effective number of independent samples in an ultrasonic measurement was quantified based on the relationship between the ROI size and the ultrasound spatial resolution.

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.

Correlations between B-mode ultrasound image texture features and tissue temperatures in hyperthermia

PURPOSE

The noninvasive monitoring of mild hyperthermia or thermal ablation is important to guarantee therapeutic safety and efficacy. The potential of ultrasound B-mode image texture features in monitoring temperature or coagulation zones studied in this article.

MATERIALS AND METHODS

The experiments carried out on eighteen in vitro porcine liver samples heated from 20°C to 60°C in the water bath. The ultrasound radiofrequency signal at different temperature collected to reconstruct B-mode ultrasound images. The texture features based on gray level histogram (GLH), gray level co-occurrence matrix (GLCM), and gray level-gradient co-occurrence matrix (GGCM) extracted, respectively. Accordingly, we analyze the correlations between these texture features and temperature based on the experiment results.

RESULTS

The results showed that five texture feature parameters closely related to temperature, including mean gray scale of GLH, homogeneity of GLCM, hybrid entropy, inverse difference moment, and correlation of GGCM. Some of these feature parameters have correlation coefficients larger than 0.9 within the temperature range of 20°C to 60°C.

CONCLUSIONS

The above-mentioned five feature parameters expected to apply for noninvasive monitoring of MH or TA.

Unsupervised deep learning based approach to temperature monitoring in focused ultrasound treatment

Temperature monitoring in ultrasound (US) imaging is important for various medical treatments, such as high-intensity focused US (HIFU) therapy or hyperthermia. In this work, we present a deep learning based approach to temperature monitoring based on radio-frequency (RF) US data. We used Siamese neural networks in an unsupervised way to spatially compare RF data collected at different time points of the heating process. The Siamese model consisted of two identical networks initially trained on a large set of simulated RF data to assess tissue backscattering properties. To illustrate our approach, we experimented with a tissue-mimicking phantom and an ex-vivo tissue sample, which were both heated with a HIFU transducer. During the experiments, we collected RF data with a regular US scanner. To determine spatiotemporal variations in temperature distribution within the samples, we extracted small 2D patches of RF data and compared them with the Siamese network. Our method achieved good performance in determining the spatiotemporal distribution of temperature during heating. Compared with the temperature monitoring based on the change in radio-frequency signal backscattered energy parameter, our method provided more smooth spatial parametric maps and did not generate ripple artifacts. The proposed approach, when fully developed, might be used for US based temperature monitoring of tissues.

Ultrasonic Thermal Monitoring of the Brain Using Golay-Coded Excitations-Feasibility Study

Thermal monitoring during focused ultrasound (FUS) transcranial procedures is mandatory and commonly performed by MRI. Transcranial ultrasonic thermal monitoring is an attractive alternative. Furthermore, using the therapeutic FUS transducer itself for this task is highly desirable. Nonetheless, such application is challenged by massive skull-induced signal attenuation and aberrations. This study examined the feasibility of implementing the Golay-coded excitations (CoE) for temperature monitoring in bovine brain samples in the range of 35 °C-43 °C (hyperthermia). Feasibility was assessed using computer simulations, water-based phantoms, and ex vivo bovine brain white-matter samples. The samples were gradually heated to about 45 °C and sonicated during cool down with a 1-MHz therapeutic FUS implementing Golay CoE. Initially, a calibration curve correlating the normalized time-of-flight (TOF) changes and the temperature was generated. Next, a bovine bone was positioned between the FUS and the brain samples, and the scanning process was repeated for different fresh samples. The calibration curve was then used as a mean for estimating the temperature, which was compared to thermocouple measurements. The simulations demonstrated a substantial improvement in signal-to-noise ratio (SNR) and suggested that the implementation of 4-bit sequences is advantageous. The experimental measurements with bone demonstrated good temperature estimation with an average absolute error for the water phantoms and brains of 1.46 °C ± 1.22 °C and 1.23 °C ± 0.99 °C, respectively. In conclusion, a novel noninvasive method utilizing the Golay CoE for ultrasonic thermal monitoring using a therapeutic FUS transducer is introduced. This method can lead to the development of an acoustic tool for brain thermal monitoring.

Statistical Analysis of Ultrasonic Scattered Echoes Enables the Non-invasive Measurement of Temperature Elevations inside Tumor Tissue during Oncological Hyperthermia

Non-invasive monitoring of temperature elevations inside tumor tissue is imperative for the oncological thermotherapy known as hyperthermia. In the present study, two cancer patients, one with a developing right renal cell carcinoma and the other with pseudomyxoma peritonei, underwent hyperthermia. The two patients were irradiated with radiofrequency current for 40 min during hyperthermia. We report the results of our clinical trial study in which the temperature increases inside the tumor tissues of patients with right renal cell carcinoma and pseudomyxoma peritonei induced by radiofrequency current irradiation for 40 min could be detected by statistical analysis of ultrasonic scattered echoes. The Nakagami shape parameter m varies depending on the temperature of the medium. We calculated the Nakagami shape parameter m by statistical analysis of the ultrasonic echoes scattered from the tumor tissues. The temperature elevations inside the tumor tissues were expressed as increases in brightness on 2-D hot-scale maps of the specific parameter α_mod, indicating the absolute values of the percentage changes in m values. In the α_mod map for each tumor tissue, the brightness clearly increased with treatment time. In quantitative analysis, the mean values of α_mod were calculated. The mean value of α_mod for the right renal cell carcinoma increased to 1.35 dB with increasing treatment time, and the mean value of α_mod for pseudomyxoma peritonei increased to 1.74 with treatment time. The increase in both α_mod brightness and the mean value of α_mod implied temperature elevations inside the tumor tissues induced by the radiofrequency current; thus, the acoustic method is promising for monitoring temperature elevations inside tumor tissues during hyperthermia.

A study on understanding the physical mechanism of change in ultrasonic envelope statistical property during temperature elevation

PURPOSE

Our previous studies demonstrate that the variation in ultrasonic envelope statistics is correlated with the temperature change inside scattering media. This variation is identified as the change in the scatterer structure during thermal expansion or contraction. However, no specific evidence has been verified to date. This study numerically reproduces the change in the scatterer distribution during thermal expansion or contraction using finite element simulations and also investigates how the situation is altered by different material properties.

METHODS

The material properties of a linear elastic solid depend on the thermal expansion coefficient, thermal conductivity, specific heat, and initial scatterer number density. Three-dimensional displacements, calculated in the simulation, were sequentially used to update the positions of the randomly distributed scatterers. Ultrasound signals from the scatterer distribution were generated by simulating a 7.5-MHz linear array transducer whose specifications were the same as those in the experimental measurements of several phantoms and excised porcine livers. To represent the change in the envelope statistical feature, the absolute value of the ratio change in the logarithmic Nakagami (NA) parameter, Δ m , at each time was calculated as a value normalized with the initial NA parameter.

RESULTS

The change in the scatterer number density relates to the volume change during temperature elevation. The magnitude of the Δ m shift against the temperature change increases depending on the higher thermal expansion coefficient. In contrast, the relationship between Δ m and the scatterer number density is similar with any material property. Additionally, the changes in Δ m obtained by several experimental phantoms with low to high scatterer number densities are comparable with the numerical simulation results.

CONCLUSIONS

The change in Δ m is indirectly related to the change in the scatterer number density owing to the volume change during thermal expansion or contraction.

Deep impact of superficial skin inking: acoustic analysis of underlying tissue

Background: Skin tattoos are a common decoration, but profound scientific study whether the presence of a skin tattoo alters the acoustic response from superficial tissue, and therefore from underlying tissue, was previously lacking. Any image aberrations caused by tattoo presence may have been thought negligible, yet empirically found artifacts in brightness-mode images of tattooed skin suggest otherwise. This study investigated the nature of these artifacts theoretically and experimentally in extremely simplified cases of perfectly flat and homogenous layered media and in tattooed pork.

Methods: Theory was derived for computing the acoustic response from horizontally and vertically layered media containing a thin inked layer. Experiments were performed in vitro. Artificial and pork skin were tattooed, attached to phantom material, and sonicated with a 13–6-MHz probe. The speed of sound of these materials was determined, and the perceived refraction angles was measured.

Results: The measured speeds of sound of tattooed materials were higher than those of their uninked counterparts. The presence of tattoo ink was found to have increased the linear acoustic attenuation by 1 dB/cm. This value is negligible for typical tattoos of only few millimeters. The perceived critical refraction angles of adjacent materials could be detected, and their corresponding speeds of sound were quantified. These coincided with values derived from theory.

Conclusion: The ratio of speeds of sound of adjacent materials was shown to create distinct highlights in brightness-mode images. The artifacts observed in in vitro and in vivo brightness-mode scans were explained from near-vertical transitions between areas of different sound speed. This is the first study correlating so-called critical refraction highlighting with speed-of-sound information. In addition, it was found that phantom material is a room-temperature acoustic alternative for experiments on live human skin. In summary, the presence of superficial tattoos has a small but quantifiable effect on the acoustic response from deeper tissues.

Method for estimating average grey-level's measurement uncertainty from ultrasound images for non-invasive estimation of temperature in different tissue types

The objective of this work is to assess, on metrological basis, the average grey-levels (AVGL) calculated from B-Mode images for estimating temperature variations non-invasively in different kinds of tissues. Thermal medicine includes several thermal therapies, being hyperthermia the most noted and well known. Recently, efforts have been made to understand the benefits of ultrasound hyperthermia at mild temperature levels, i.e., between 39 °C and 41 °C. Moreover, the best practices on ultrasound bio-effects research have been encouraged by recommending that temperature rise in the region of interest should be measured even when a thermal mechanism is not being tested. In this work, the average grey-levels (AVGL) calculated from B-Mode images were assessed for non-invasive temperature estimation in a porcine tissue sample containing two different tissue types, fat and muscle, with temperature varying from 35 °C to 41 °C. The sample was continuously imaged with an ultrasound scanner, and simultaneously the temperature was measured. The achieved results were assessed under the light of the measurement uncertainty in order to allow comparability among different ultrasound thermometry methods. The highest expanded uncertainty of estimating temperature variation using AVGL was determined as 0.68 °C.

Temperature elevation in tissue detected in vivo based on statistical analysis of ultrasonic scattered echoes

It is demanded to monitor temperature in tissue during oncological hyperthermia therapy. In the present study, we non-invasively measured the temperature elevation inside the abdominal cavity and tumour tissue of a living rat induced by capacitive-coupled radiofrequency heating. In the analysis of ultrasound scattered echoes, the Nakagami shape parameter m in each region of interest was estimated at each temperature. The Nakagami shape parameter m has temperature dependence; hence, the temperature increase inside tissue specimens can be detected with the m values. By carrying out in vivo experiments, we visualized the temperature increase inside the abdominal cavity and tumour tissue of living rats using two-dimensional hot-scale images indicating the absolute values of the ratio changes of the m values. In both the abdominal cavity and tumour tissue, the brightness in the hot-scale images clearly increased with increasing temperature. The increases in brightness in the hot-scale images imply the temperature elevations inside the abdominal cavity and tumour tissue of the living rats. The study results prove that the acoustic method we proposed is a promising method for monitoring changes in the internal temperature of the human body under hyperthermia treatment.

Speed of Sound and Attenuation Temperature Dependence of Bovine Brain: Ex Vivo Study

OBJECTIVES

Brain treatments using focused ultrasound (FUS) offer a new range of noninvasive transcranial therapies. The acoustic energy deposition during these procedures may induce a temperature elevation in the tissue; therefore, noninvasive thermal monitoring is essential. Magnetic resonance imaging is the current adopted monitoring modality, but its high operational costs and limited availability may hinder the accessibility to FUS treatments. Aiming at the development of a thermometric ultrasound (US) method for the brain, the specific objective of this investigation was to study the acoustic thermal response of the speed of sound (SOS) and attenuation coefficient (AC) of different brain tissues: namely white matter (WM) and cortical matter.

METHODS

Sixteen ex vivo bovine brain samples were investigated. These included 7 WM and 9 cortical matter samples. The samples were gradually heated to about 45°C and then repeatedly scanned while cooling using a computerized US system in the through-transmission mode. The temperature was simultaneously registered with thermocouples. From the scans, the normalized SOS and AC for both tissues were calculated.

RESULTS

The results demonstrated a characteristic cooldown temporal behavior for the normalized AC and SOS curves, which were related to the temperature. The SOS curves enabled clear differentiation between the tissue types but depicted more scattered trajectories for the WM tissue. As for the AC curves, the WM depicted a linear behavior in relation to the temperature. However, both tissue types had rather similar temperature patterns.

CONCLUSIONS

These findings may contribute to the development of a US temperature-monitoring method during FUS procedures.

Non-invasive Estimation of Temperature during Physiotherapeutic Ultrasound Application Using the Average Gray-Level Content of B-Mode Images: A Metrological Approach

Healing therapies that make use of ultrasound are based on raising the temperature in biological tissue. However, it is not possible to heal impaired tissue by applying a high dose of ultrasound. The temperature of the tissue is ultimately the physical quantity that has to be assessed to minimize the risk of undesired injury. Invasive temperature measurement techniques are easy to use, despite the fact that they are detrimental to human well being. Another approach to assessing a rise in tissue temperature is to derive the material's general response to temperature variations from ultrasonic parameters. In this article, a method for evaluating temperature variations is described. The method is based on the analytical study of an ultrasonic image, in which gray-level variations are correlated to the temperature variations in a tissue-mimicking material. The physical assumption is that temperature variations induce wave propagation changes modifying the backscattered ultrasound signal, which are expressed in the ultrasonographic images. For a temperature variation of about 15°C, the expanded uncertainty for a coverage probability of 0.95 was found to be 2.5°C in the heating regime and 1.9°C in the cooling regime. It is possible to use the model proposed in this article in a straightforward manner to monitor temperature variation during a physiotherapeutic ultrasound application, provided the tissue-mimicking material approach is transferred to actual biological tissue. The novelty of such approach resides in the metrology-based investigation outlined here, as well as in its ease of reproducibility.

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.

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.

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.

Evaluation of thrombolysis by using ultrasonic imaging: an in vitro study

The hematocrit of a thrombus is a key factor associated with the susceptibility to thrombolysis. Ultrasonic imaging is currently the first-line screening tool for thrombus examinations. Different hematocrits result in different acoustical structures of thrombi, which alter the behavior of ultrasonic backscattering. This study explored the relationships among thrombolytic efficiencies, hematocrits, and ultrasonic parameters (the echo intensity and backscattered statistics). Porcine thrombi with different hematocrits, ranging from 0% to 50%, were induced in vitro. An ultrasonic scanner was used to scan thrombi and acquire raw image data for B-mode (echo intensity measurements) and Nakagami imaging (backscattered statistics analysis). Experiments on thrombolysis were performed using urokinase to explore the effect of the hematocrit on thrombolytic efficiency. Results showed that the weight loss ratio of thrombi exponentially decreased as the hematocrit increased from 0% to 50%. Compared with the echo intensity obtained from the conventional B-scan, the Nakagami parameter predicts the weight loss ratio, increasing from 0.6 to 1.2 as the weight loss ratio decreased from 0.67 to 0.26. The current findings suggest that using Nakagami imaging characterizing thrombi provides information of backscattered statistics, which may be associated with the thrombolytic efficiency.

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.

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.

Sensitivity enhanced nanothermal sensors for photoacoustic temperature mapping

Photoacoustic imaging can be used to guide and validate the therapeutic outcome of nanoparticle-mediated photothermal therapy through its ability to visualize the delivery of nanoparticle contrast agents, image the temperature distribution inside living tissue, and confirm tissue coagulation. In this image-guided process, temperature mapping plays a critical role for thermal dosage control. Therefore, developing a sensitive and accurate photoacoustic technique to quantitatively measure the temperature distribution during thermal therapy is essential. In this study, we investigated and demonstrated that silica-coated gold nanorods, can provide a multi-fold improvement in sensitivity of the photoacoustic temperature mapping compared to gold nanorods without silica coating, and serve as a nanothermal sensor to accurately and quantitatively visualize temperature distributions during photothermal therapy.

Relationship between ultrasound backscattered statistics and the concentration of fatty droplets in livers: an animal study

Ultrasound grayscale B-mode imaging is the most frequently used modality for examining fatty liver. Different concentrations and arrangements of fatty droplets in the liver may produce different statistical distributions of ultrasound backscatter signals, which may be treated as a useful clue for assessing the stage of fatty liver. To verify this point, we investigate the relationship between changes in backscattered statistics and the concentration of fatty droplets in the liver. Fatty liver was induced in rats fed a methionine-choline-deficient diet. Livers were excised from rats for in vitro ultrasound scanning using a single-element transducer. The envelopes of the acquired raw ultrasound signals were used for the analysis of the backscattered statistics by ultrasound Nakagami parametric imaging, which has been shown as a reliable tool to model the statistical distribution of ultrasound backscatter signals. Histological analyses and the measurements of triglyceride and cholesterol in the rat liver were conducted for comparison with the Nakagami parameter. Results show that the ultrasound Nakagami parameter has an excellent correlation with the concentration of fatty droplets, demonstrating that ultrasound backscatter statistics depend on the degree of fatty liver in rats.

Postmastoidectomy effusion measurement using a delay-line ultrasound transducer: cadaver experiments

Mastoidectomy is a surgical procedure for removing mastoid air cells and has been used widely to establish drainage and cleaning infections from the mastoid bone. Computed tomography (CT) and magnetic resonance imaging (MRI) are the primary instruments used to clinically detect mastoid effusion following surgery. This study examines the feasibility of ultrasonography to detect postsurgical mastoid effusion. In vitro ultrasound measurements were conducted on 10 cadavers. For each sample, mastoidectomy was performed, and saline and whole blood samples were injected into the surgical mastoid cavity to simulate different effusion properties. A 2.25-MHz delay-line ultrasound transducer was used to obtain ultrasound backscattered data from the mastoid. The data were used to analyze echo intensity and the backscattered statistics by estimating the scaling and Nakagami parameters of the Nakagami distribution model. The results show that the scaling and Nakagami parameters can successfully detect mastoid effusion. Specifically, the Nakagami parameter is capable of characterizing the mastoid effusion properties. This indicates that ultrasound measurement based on a combination of delay-line transducer and Nakagami parameter estimation is a potential real-time diagnostic tool for evaluating mastoid effusion following mastoidectomy without radiation exposure.