Dual energy computed tomography thermometry during hepatic microwave ablation in an ex-vivo porcine model

In Vivo Temperature Mapping during CT-Guided Thermal Ablations: Proof of Feasibility Using Spectral CT

This study aimed to perform temperature mapping during clinical ablations by leveraging the known relationship between computed tomography (CT) signal and temperature. First, the spectral CT signal was characterized as a function of temperature for soft tissue-mimicking and fat-mimicking materials across a wide range of temperatures in both frozen and heated states. Spectral CT images were acquired while fiber optic thermal sensors continuously recorded local temperatures, which were correlated with electron density changes. Substantially different material-specific thermal expansion coefficients and temperature dependencies were measured and well matched reported values. In patient images, effective atomic number images were used to segment soft tissue, fat, and bone. Electron density-temperature relationships for fat and soft tissue were used to generate temperature maps during microwave ablations and cryoablations. In vivo temperature maps correlated well with interventional radiologist expectations and demonstrated smooth transitions at fat/soft tissue interfaces.

Reproducible spectral CT thermometry with liver-mimicking phantoms for image-guided thermal ablation

Objectives. Evaluate the reproducibility, temperature tolerance, and radiation dose requirements of spectral CT thermometry in tissue-mimicking phantoms to establish its utility for non-invasive temperature monitoring of thermal ablations.Methods. Three liver mimicking phantoms embedded with temperature sensors were individually scanned with a dual-layer spectral CT at different radiation dose levels during heating (35 °C-80 °C). Physical density maps were reconstructed from spectral results using varying reconstruction parameters. Thermal volumetric expansion was then measured at each temperature sensor every 5 °C in order to establish a correlation between physical density and temperature. Linear regressions were applied based on thermal volumetric expansion for each phantom, and coefficient of variation for fit parameters was calculated to characterize reproducibility of spectral CT thermometry. Additionally, temperature tolerance was determined to evaluate effects of acquisition and reconstruction parameters. The resulting minimum radiation dose to meet the clinical temperature accuracy requirement was determined for each slice thickness with and without additional denoising.Results. Thermal volumetric expansion was robustly replicated in all three phantoms, with a correlation coefficient variation of only 0.43%. Similarly, the coefficient of variation for the slope and intercept were 9.6% and 0.08%, respectively, indicating reproducibility of the spectral CT thermometry. Temperature tolerance ranged from 2 °C to 23 °C, decreasing with increased radiation dose, slice thickness, and iterative reconstruction level. To meet the clinical requirement for temperature tolerance, the minimum required radiation dose ranged from 20, 30, and 57 mGy for slice thickness of 2, 3, and 5 mm, respectively, but was reduced to 2 mGy with additional denoising.Conclusions. Spectral CT thermometry demonstrated reproducibility across three liver-mimicking phantoms and illustrated the clinical requirement for temperature tolerance can be met for different slice thicknesses. The reproducibility and temperature accuracy of spectral CT thermometry enable its clinical application for non-invasive temperature monitoring of thermal ablation.

Reproducible spectral CT thermometry with liver-mimicking phantoms for image-guided thermal ablation

Objectives

Evaluate the reproducibility, temperature sensitivity, and radiation dose requirements of spectral CT thermometry in tissue-mimicking phantoms to establish its utility for non-invasive temperature monitoring of thermal ablations.

Materials and Methods

Three liver mimicking phantoms embedded with temperature sensors were individually scanned with a dual-layer spectral CT at different radiation dose levels during heating and cooling (35 to 80 °C). Physical density maps were reconstructed from spectral results using a range of reconstruction parameters. Thermal volumetric expansion was then measured at each temperature sensor every 5°C in order to establish a correlation between physical density and temperature. Linear regressions were applied based on thermal volumetric expansion for each phantom, and coefficient of variation for fit parameters was calculated to characterize reproducibility of spectral CT thermometry. Additionally, temperature sensitivity was determined to evaluate the effect of acquisition parameters, reconstruction parameters, and image denoising. The resulting minimum radiation dose to meet the clinical temperature sensitivity requirement was determined for each slice thickness, both with and without additional denoising.

Results

Thermal volumetric expansion was robustly replicated in all three phantoms, with a correlation coefficient variation of only 0.43%. Similarly, the coefficient of variation for the slope and intercept were 9.6% and 0.08%, respectively, indicating reproducibility of the spectral CT thermometry. Temperature sensitivity ranged from 2 to 23 °C, decreasing with increased radiation dose, slice thickness, and iterative reconstruction level. To meet the clinical requirement for temperature sensitivity, the minimum required radiation dose ranged from 20, 30, and 57 mGy for slice thickness of 2, 3, and 5 mm, respectively, but was reduced to 2 mGy with additional denoising.

Conclusions

Spectral CT thermometry demonstrated reproducibility across three liver-mimicking phantoms and illustrated the clinical requirement for temperature sensitivity can be met for different slice thicknesses. Moreover, additional denoising enables the use of more clinically relevant radiation doses, facilitating the clinical translation of spectral CT thermometry. The reproducibility and temperature accuracy of spectral CT thermometry enable its clinical application for non-invasive temperature monitoring of thermal ablation.

Non-invasive mass and temperature quantifications with spectral CT

Spectral CT has been increasingly implemented clinically for its better characterization and quantification of materials through its multi-energy results. It also facilitates calculation of physical density, allowing for non-invasive mass measurements and temperature evaluations by manipulating the definition of physical density and thermal volumetric expansion, respectively. To develop spectral physical density quantifications, original and parametrized Alvarez-Macovski model and electron density-physical density model were validated with a phantom. The best physical density model was then implemented on clinical spectral CT scans of ex vivo bovine muscle to determine the accuracy and effect of acquisition parameters on mass measurements. In addition, the relationship between physical density and changes in temperature was evaluated by scanning and subjecting the tissue to a range of temperatures. The parametrized Alvarez-Macovski model performed best in both model development and validation with errors within ± 0.02 g/mL. No effect from acquisition parameters was observed in mass measurements, which demonstrated accuracy with a maximum percent error of 0.34%. Furthermore, physical density was strongly correlated (R of 0.9781) to temperature changes through thermal volumetric expansion. Accurate and precise spectral physical density quantifications enable non-invasive mass measurements for pathological detection and temperature evaluation for thermal therapy monitoring in interventional oncology.

Predicting the coagulation volume induced by microwave ablation of hepatocellular carcinoma: the role of deposited energy, ex-vivo bovine liver charts and central hyperdense area on post-treatment CT

PURPOSE

To study the correlation between the overall coagulation zone (A) attained in percutaneous microwave ablation (MWA) of hepatocellular carcinomas (HCC) and: (1) the hyperdense zone (C) visible in the central part of zone A on post-treatment unenhanced CT scans; (2) the deposited energy; (3) the coagulation zones observed on ex-vivo bovine liver.

MATERIALS AND METHODS

The post-procedural computed tomography (CT) scans of HCCs treated with a single energy deployment through the same 2450 MHz MWA system were retrospectively analyzed, retrieving the dimensions of A and C zones and the deposited energy (E). Ex-vivo bovine liver MWA with the same system were performed and analyzed to determine the same quantities by gross-pathologic examination and CT imaging.

RESULTS

A total of 101 HCC treatments were analyzed. The average coagulation volumes increased linearly with deposited energy (1.11 cc/kJ, R² = 0.90, 4.2 kJ ≤ E ≤ 48 kJ), similarly to ex-vivo findings (1.38 cc/kJ, R² =0.97, 7.2 kJ ≤ E ≤ 144 kJ). The long axis (L) and short axis (D) of zones A and C held a fairly constant ratio both in-vivo (L_C/L_A=0.43 ± 0.13; D_C/D_A=0.42 ± 0.10) and ex-vivo (L_C/L_A = 0.49 ± 0.07; D_C/D_A = 0.28 ± 0.06).

CONCLUSIONS

The average dimensions of the ablation zone induced by the considered system on HCC increase linearly with the deposited energy and are fairly well predicted by the corresponding ex-vivo dimensions. The ratio between each linear dimension of A and C zones was found to be roughly constant over a large deposited energy span, both ex-vivo and in-vivo.

Temperature Monitoring in Hyperthermia Treatments of Bone Tumors: State-of-the-Art and Future Challenges

Bone metastases and osteoid osteoma (OO) have a high incidence in patients facing primary lesions in many organs. Radiotherapy has long been the standard choice for these patients, performed as stand-alone or in conjunction with surgery. However, the needs of these patients have never been fully met, especially in the ones with low life expectancy, where treatments devoted to pain reduction are pivotal. New techniques as hyperthermia treatments (HTs) are emerging to reduce the associated pain of bone metastases and OO. Temperature monitoring during HTs may significantly improve the clinical outcomes since the amount of thermal injury depends on the tissue temperature and the exposure time. This is particularly relevant in bone tumors due to the adjacent vulnerable structures (e.g., spinal cord and nerve roots). In this Review, we focus on the potential of temperature monitoring on HT of bone cancer. Preclinical and clinical studies have been proposed and are underway to investigate the use of different thermometric techniques in this scenario. We review these studies, the principle of work of the thermometric techniques used in HTs, their strengths, weaknesses, and pitfalls, as well as the strategies and the potential of improving the HTs outcomes.

CT-based thermometry with virtual monoenergetic images by dual-energy of fat, muscle and bone using FBP, iterative and deep learning-based reconstruction

Objectives

The aim of this study was to evaluate the sensitivity of CT-based thermometry for clinical applications regarding a three-component tissue phantom of fat, muscle and bone. Virtual monoenergetic images (VMI) by dual-energy measurements and conventional polychromatic 120-kVp images with modern reconstruction algorithms adaptive statistical iterative reconstruction-Volume (ASIR-V) and deep learning image reconstruction (DLIR) were compared.

Methods

A temperature-regulating water circuit system was developed for the systematic evaluation of the correlation between temperature and Hounsfield units (HU). The measurements were performed on a Revolution CT with gemstone spectral imaging technology (GSI). Complementary measurements were performed without GSI (voltage 120 kVp, current 130–545 mA). The measured object was a tissue equivalent phantom in a temperature range of 18 to 50°C. The evaluation was carried out for VMI at 40 to 140 keV and polychromatic 120-kVp images.

Results

The regression analysis showed a significant inverse linear dependency between temperature and average HU regardless of ASIR-V and DLIR. VMI show a higher temperature sensitivity compared to polychromatic images. The temperature sensitivities were 1.25 HU/°C (120 kVp) and 1.35 HU/°C (VMI at 140 keV) for fat, 0.38 HU/°C (120 kVp) and 0.47 HU/°C (VMI at 40 keV) for muscle and 1.15 HU/°C (120 kVp) and 3.58 HU/°C (VMI at 50 keV) for bone.

Conclusions

Dual-energy with VMI enables a higher temperature sensitivity for fat, muscle and bone. The reconstruction with ASIR-V and DLIR has no significant influence on CT-based thermometry, which opens up the potential of drastic dose reductions.

Key Points

• Virtual monoenergetic images (VMI) enable a higher temperature sensitivity for fat (8%), muscle (24%) and bone (211%) compared to conventional polychromatic 120-kVp images.

• With VMI, there are parameters, e.g. monoenergy and reconstruction kernel, to modulate the temperature sensitivity. In contrast, there are no parameters to influence the temperature sensitivity for conventional polychromatic 120-kVp images.

• The application of adaptive statistical iterative reconstruction-Volume (ASIR-V) and deep learning–based image reconstruction (DLIR) has no effect on CT-based thermometry, opening up the potential of drastic dose reductions in clinical applications.

Hyperspectral image-based analysis of thermal damage for ex-vivo bovine liver utilizing radiofrequency ablation

BACKGROUND & OBJECTIVE

Thermal ablation is the predominant methodology to treat liver tumors for segregating patients who are not permitted to have surgical intervention. However, noticing or predicting the size of the thermal strategies is a challenging endeavor. We aim to analyze the effects of ablation district volume following radiofrequency ablation (RFA) of ex-vivo liver exploiting a custom Hyperspectral Imaging (HSI) system.

MATERIALS AND METHODS

RFA was conducted on the ex-vivo bovine liver at focal and peripheral blood vessel sites and observed by Custom HSI system, which has been designed to assess the exactness and proficiency using visible and near-infrared wavelengths region for tissue thermal effect. The experiment comprised up to ten trials with RFA. The experiment was carried out in two stages to assess the percentage of the thermal effect on the investigated sample superficially and for the side penetration effect. Measuring the diffuse reflectance (Ŗ_d) of the sample to identify the spectral reflectance shift which could differentiate between normal and ablated tissue exploiting the designed cross-correlation algorithm for monitoring of thermal ablation.

RESULTS

Determination of the diffuse reflection (Ŗ_d) spectral signature responses from normal, thermal effected, and thermal ablation regions of the investigated liver sample. Where the ideal wavelength range at (600-640 nm) could discriminate between these different regions. Then, exploited the converted RGB image of the HS liver tissue after RFA for more validations which shows that the optimum wavelength for differentiation at (530-560 nm and 600-640 nm). Finally, applying statistical analysis to validate our results presenting that wavelength 600 nm had the highest standard deviation (δ) to differentiate between various thermally affected regions regarding the normal tissue and wavelength 640 nm shows the highest (δ) to differentiate between the ablated and normal regions.

CONCLUSION

The designed and implemented medical imaging system incorporated the hyperspectral camera capabilities with the associate cross-correlation algorithm that could successfully distinguish between the ablated and thermally affected regions to assist the surgery during the tumor therapy.

Computed Tomography Thermography for Ablation Zone Prediction in Microwave Ablation and Cryoablation: Advantages and Challenges in an Ex Vivo Porcine Liver Model

PURPOSE

The aim of this study was to investigate the diagnostic accuracy of computed tomography (CT) for the prediction of ablation zones from microwave ablation (MWA) and cryoablation (CA) in an ex vivo porcine liver model.

METHODS

Sequential (30 seconds) CT scans were acquired during and after MWA and CA in an ex vivo porcine liver model. We generated 120-kVp equivalent reconstructions of generic dual-energy CT data sets, and comprehensive region-of-interest measurements were statistically correlated with invasive temperature monitoring using Pearson correlation coefficient. Binary logistic regression was performed for prediction of successful ablation.

RESULTS

With the use of pooled data from 6 lesions in 2 separate experiments, correlation analysis of attenuation in Hounsfield units (HU) and temperature yielded r = -0.79 [confidence interval (CI), -0.85 to -0.71] for MWA and r = 0.62 (CI, 0.55 to 0.67) for CA.For MWA, there was a linear association between attenuation and temperature up to 75°C; thus, linear regression yielded a slope of -2.00 HU/°C (95% CI, -1.58 to -2.41). For CA, a linear association between attenuation and temperature was observed in the cooling phase with a slope of 2.11 HU/°C (95% CI, 1.79 to 2.58). In MWA treatment, binary logistic regression separated less than 70°C and greater than 70°C with 89.2% accuracy. Within the ice ball, temperatures above and below -20°C were distinguished with 65.3% accuracy.

CONCLUSIONS

Our experiments reveal several difficulties in predicting ablation zone temperature from CT attenuation. Microwave ablation leads to gas production in the tissue, which degrades the accuracy of noninvasive temperature measurement, especially at higher temperatures. In CA, CT thermometry is limited by ice ball formation, which leads to homogeneous attenuation, nearly independent of temperature. Further research is needed to define the role of CT thermography in ablation zone monitoring in liver malignancies.

Evaluation of the thermal sensitivity of porcine liver in CT-guided cryoablation: an initial study

PURPOSE

To evaluate computed tomography (CT)-based thermometry in cryoablation, the thermal sensitivity of an ex-vivo porcine liver was determined in an initial study design.

METHODS

The CT-guided cryoablation was performed in three porcine liver samples over a period of 10 min. Fiber optic temperature probes were positioned parallel to the shaft of the cryoprobe in an axial slice orientation. During ablation, temperature measurements were performed simultaneously with CT imaging at 5 s intervals. On the CT images, the average CT number was calculated for a region of interest of 3 × 3 pixels just below the tip of each temperature probe. A linear regression analysis was performed using eleven data sets to determine the dependence of the CT number on the temperature.

RESULTS

With decreasing temperature, an increasing hypodense area around the tip of the cryoprobe was observed on the CT images and decreasing values of the CT number were determined. Starting at a temperature of - 40°C a linear relation between the CT number and the temperature was determined and a thermal sensitivity of 0.95 HU/°C (R²  = 0.73) was obtained. The thermal sensitivity was used to calculate color-coded temperature maps. The calculated temperature distribution corresponds quantitatively to the increasing hypodense area.

CONCLUSIONS

A noninvasive CT-based temperature determination during cryoablation in a normal ex vivo porcine liver is feasible. A thermal sensitivity of 0.95 HU/°C was determined by linear regression analysis. A color-coded map of the temperature distribution was presented.

In vivo noninvasive three-dimensional (3D) assessment of microwave thermal ablation zone using non-contrast-enhanced x-ray CT

PURPOSE

To develop an image processing methodology for noninvasive three-dimensional (3D) quantification of microwave thermal ablation zones in vivo using x-ray computed tomography (CT) imaging without injection of a contrast enhancing material.

METHODS

Six microwave (MW) thermal ablation procedures were performed in three pigs. The ablations were performed with a constant heating duration of 8 min and power level of 30 W. During the procedure images from sixty 1 mm thick slices were acquired every 30 s. At the end of all ablation procedures for each pig, a contrast-enhanced scan was acquired for reference. Special algorithms for addressing challenges stemming from the 3D in vivo setup and processing the acquired images were prepared. The algorithms first rearranged the data to account for the oblique needle orientation and for breathing motion. Then, the gray level variance changes were analyzed, and optical flow analysis was applied to the treated volume in order to obtain the ablation contours and reconstruct the ablation zone in 3D. The analysis also included a special correction algorithm for eliminating artifacts caused by proximal major blood vessels and blood flow. Finally, 3D reference reconstructions from the contrast-enhanced scan were obtained for quantitative comparison.

RESULTS

For four ablations located >3 mm from a large blood vessel, the mean dice similarity coefficient (DSC) and the mean absolute radial discrepancy between the contours obtained from the reference contrast-enhanced images and the contours produced by the algorithm were 0.82 ± 0.03 and 1.92 ± 1.47 mm, respectively. In two cases of ablation adjacent to large blood vessels, the average DSC and discrepancy were: 0.67 ± 0.6 and 2.96 ± 2.15 mm, respectively. The addition of the special correction algorithm utilizing blood vessels mapping improved the mean DSC and the mean absolute discrepancy to 0.85 ± 0.02 and 1.19 ± 1.00 mm, respectively.

CONCLUSIONS

The developed algorithms provide highly accurate detailed contours in vivo (average error < 2.5 mm) and cope well with the challenges listed above. Clinical implementation of the developed methodology could potentially provide real time noninvasive 3D accurate monitoring of MW thermal ablation in-vivo, provided that the radiation dose can be reduced.

Liver microwave ablation: a systematic review of various FDA-approved systems

OBJECTIVES

The aim of the present study is to analyze preclinical and clinical data on the performance of the currently US Food and Drug Administration (FDA)-approved microwave ablation (MWA) systems.

METHODS

A review of the literature, published between January 1, 2005, and December 31, 2016, on seven FDA-approved MWA systems, was conducted. Ratio of ablation zone volume to applied energy R(AZ:E) and sphericity indices were calculated for ex vivo and in vivo experiments.

RESULTS

Thirty-four studies with ex vivo, in vivo, and clinical data were summarized. In total, 14 studies reporting data on ablation zone volume and applied energy were included for comparison R(AZ:E). A significant correlation between volume and energy was found for the ex vivo experiments (r = 0.85, p < 0.001) in contrast to the in vivo experiments (r = 0.54, p = 0.27).

CONCLUSION

Manufacturers' algorithms on microwave ablation zone sizes are based on preclinical animal experiments with normal liver parenchyma. Clinical data reporting on ablation zone volume in relation to applied energy and sphericity index during MWA are scarce and require more adequate reporting of MWA data.

KEY POINTS

• Clinical data reporting on the ablation zone volume in relation to applied energy during microwave ablation are scarce. • Manufacturers' algorithms on microwave ablation zone sizes are based on preclinical animal experiments with normal liver parenchyma. • Preclinical data do not predict actual clinical ablation zone volumes in patients with liver tumors.

Microwave thermal ablation using CT-scanner for predicting the variation of ablated region over time: advantages and limitations

This study aims at investigating in real-time the structural and dynamical changes occurring in an ex vivo tissue during a microwave thermal ablation (MTA) procedure. The experimental set-up was based on ex vivo liver tissue inserted in a dedicated box, in which 3 fibre-optic (FO) temperature probes were introduced to measure the temperature increase over time. Computed tomography (CT) imaging technique was exploited to experimentally study in real-time the Hounsfield Units (HU) modification occurring during MTA. The collected image data were processed with a dedicated MATLAB tool, developed to analyse the FO positions and HU modifications from the CT images acquired over time before and during the ablation procedures. The radial position of a FO temperature probe (rFO) and the value of HU in the region of interest (ROI) containing the probe (HUo), along with the corresponding value of HU in the contralateral ROI with respect to the MTA antenna applicator (HUc), were determined and registered over time during and after the MTA procedure. Six experiments were conducted to confirm results. The correlation between temperature and the above listed predictors was investigated using univariate and multivariate analysis. At the multivariate analysis, the time, rFO and HUc resulted significant predictive factors of the logarithm of measured temperature. The correlation between predicted and measured temperatures was 0.934 (p   <  0.001). The developed tool allows identifying and registering the image-based parameters useful for predicting the temperature variation over time in each investigated voxel by taking into consideration the HU variation.