Sensitivity of CT-based Thermometry for Clinical Applications

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

Authors: Andreas Heinrich¹ & Sebastian Schenkl² & David Buckreus¹ & Felix V. Güttler¹ & Ulf K-M. Teichgräber¹

¹ Department of Radiology, Jena University Hospital – Friedrich Schiller University, Am Klinikum 1, 07747 Jena, Germany

² Institute of Forensic Medicine, Jena University Hospital – Friedrich Schiller University, Am Klinikum 1, 07747 Jena, Germany

Abstract

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 ofdrastic dosereductionsin clinical applications.

Keywords: Computed tomography, Thermometry, Dual-energy scanned projection, Deep learning

Abbreviations

ASIR-V	Adaptive statistical iterative reconstructionVolume
CI	Confidence interval
DLIR	Deep learning image reconstruction
FBP	Filtered back projection
GSI	Gemstone spectral imaging
VMI	Virtual monoenergetic images

Introduction

Radiodensity depends on the temperature of the tissue, since the density decreases if the temperature of the tissue increases. The changes in tissue temperature can therefore be measured indirectly via the computed tomography (CT) numbers or Hounsfield units (HU), also called CTbased thermometry. However, the sensitivity (change of HU per °C) depends on the volumetric thermal expansion coefficient of the specific material [1]. The temperature sensitivity was reported to range between −0.40 and −0.29 HU/°C for water [1], between −0.52 and −0.60 HU/°C for porcine liver [2, 3], and between −0.45 and −0.43 HU/°C for muscle [4, 5]. Possible applications of CT-based thermometry would be the monitoring of tissue temperatures during high-frequency treatments or microwave ablation of tumors in the liver and kidneys [2, 3, 6–9]. Despite these clear relationships, the method is rarely used in clinical practice so far [10], because the repeatability of quantitative CT numerical measurements has not been guaranteed and a higher required ionizing radiation, e.g. for thermometry due to more frequent scans [11]. Further research is required to systematically analyze the influencing factors on CT-based thermometry [12, 13], so as demonstrate if better levels of thermal-spatial resolution could be achieved, along with lower patient radiation doses [1]. With recent technical advances in CT (Fig. 1), such as dual-energy and iterative reconstruction, e.g. adaptive statistical iterative reconstruction-Volume (ASIR-V) [14, 15], or deep learning image reconstruction (DLIR) [16], the question arises whether these advances have improved the quality of CT-based thermometry. DLIR represents the latest stage of development in CT image reconstruction and promises significant improvements in image quality and radiation dose reduction compared to ASIR-V and filtered back projection (FBP) [16–18]. A dual-energy CT with a gemstone spectral imaging (GSI) detector measures the raw data of an object during a continuous energy change between low (80 kVp) and high (140 kVp) voltage [19]. This generates two data records with different attenuation values based on the corresponding energy levels, which allow the reconstruction of material decomposition images and virtual monoenergetic images (VMI). VMI demonstrated significantly better signal-to-noise ratio (SNR), contrast-to-noise ratio (CNR) and subjective score compared to conventional 120-kVp polychromatic images [20].

The objective of this study was to evaluate CT-based thermometry for fat, muscle and bone using state-of-the-art imaging technology. This includes the dose reduction potential of modern reconstruction algorithms with ASIR-V and DLIR, and the use of VMI at 40–140 keV versus conventional polychromatic 120-kVp images.

Materials and methods

The measured object was a tissue-equivalent phantom [21, 22], composed of beeswax (density 0.95 g/cm³, fat-equivalent), salt-water solution (density 1.04 g/cm³, muscle-equivalent) and a bovine femoral bone (density 1.36 g/cm³). The phantom (length 19.40 cm, width 17.20–17.40 cm and height 15.20 cm) was placed into a separate phantom container (length 39.00 cm, width 32.00 cm and height 30.50 cm) and fastened.

Temperature regulation

A temperature-regulating water circuit system (Fig. 2) was developed for the systematic evaluation of the correlation of temperature and CT number. More than 35 L of water was heated in a water bath (Thermoboy c20, mgw Lauda) and subsequently pumped through a silicone hose into the phantom housing container. The phantom container features an overflow port located above the top edge of the phantom. Consequently, it is guaranteed that the phantom is completely immersed by water. When it reaches the overflow port, the water flows back into the water bath basin. This measurement setup constitutes a closed system; thus, the temperature can be kept constant over a long period of time. The temperature of the system was monitored by four T1 fiber-optic temperature probes (Rugged Fiber Optic Temperature Monitoring System). In the top center of the phantom is a 7.50-cm-deep blind hole that accommodates a temperature probe (probe 1) for monitoring the phantom’s core temperature. Another temperature probe (probe 2) is located on the side at a depth of 2 cm in the phantom. To avoid the risk of water infiltrating and thus distorting the temperature measurement, both holes were sealed with beeswax. Two additional temperature probes were located directly in the water in the phantom housing tank (probe 3) and in the water bath (probe 4).

CT measurements and evaluation

The measurements were performed in a temperature range from 18 up to 50°C (step size approx. 5°C) on a Revolution CT (GE Healthcare) by means of the GSI technology (voltage 80 kV and 140 kV, current 235–360 mA SmartmA, helical scan type, pitch 0.992). Complementary measurements were performed without GSI (voltage 120 kV, current 130–545 mA, helical scan type, pitch 0.992). Measurements were repeated (n = 14) for each temperature and CT protocol. To increase the reliability, the measurements were performed in thermodynamic equilibrium of the phantom. This required a waiting time of 4 to 10 h between measurements.

Fig. 1 Scheme of data acquisition and image reconstruction in the CT.

The raw data were reconstructed by FBP (= ASIR-V 0%), ASIR-V (50%, 100%) and DLIR-L (lowest filtering), DLIR-M (medium filtering) and DLIR-H (highest filtering) with the standard kernel and an image matrix of 512 × 512 pixels. For the evaluation of the GSI raw data, VMI at 40 keV, 50 keV, 60 keV, 70 keV, 80 keV, 100 keV, 120 keV and 140 keV were used. Additionally, the kernels bone, boneplus, soft and detail were examined for polychromatic images and GSI (VMI at 40 keV, 70 keV and 120 keV) with ASIR-V 100%. The standard kernel had to be used for DLIR.

Fig. 2 Left virtual monoenergetic image (VMI) at 70 keV of the tissueequivalent phantom with rectangle regions of interest (ROIs) for measurement of fat-equivalent (beeswax, 1), muscle-equivalent (salty water solution, 2) and bovine bone (3). Right temperature-regulating water circuit system with (a) phantom housing container with fixture for the (b) tissueequivalent phantom and (c) temperature probes. Water is pumped into the phantom housing container through (d) a silicone hose, (e) an overflow pipe returning the water back into the (f) water bath with heater (0–100 °C) and pump. During a series of measurements, the phantom container and the water bath were covered by lids so that the system constituted a closed circuit.

For each measurement, the average HU of rectangle regions of interest (ROIs) with mean areas of 8.24 cm² for fat (25 × 96 pixels) and muscle (40 × 60 pixels) and 2.58 cm² for bone (5 × 150 pixels) was analyzed (Fig. 2 left). The ROIs were set automatically; thus, the position of the ROIs was exactly at the same position for all reconstructions. The dependency of the average CT numbers as a function of temperature was analyzed using linear regression analysis. The software SPSS Statistics, version 26 (IBM), was used for statistical evaluation. The 95% and 50% confidence intervals (CI) were calculated for the temperature sensitivity.

Results

The measurement setup enabled reproducible measurements. During the measuring process, no movement artefacts caused by the circulating water system were observed (Fig. 3). The bone, in parts, showed beam hardening artefacts. The measurement uncertainty due to temperature gradients between the margin of the phantom and the core center was 0.59 ± 0.37 °C (temperature probes 1 and 2). The difference between water temperatures in the water bath and phantom housing container was 0.37 ± 0.28 °C (temperature probes 3 and 4).

The regression analysis showed a significant inverse linear dependency between temperature and average CT number (Fig. 4 left). Depending on the monoenergy, the VMI show a higher temperature sensitivity compared to polychromatic images (Fig. 4 right). For fat, the temperature sensitivities were 1.25 HU/°C (95% CI 1.18–1.33) for 120 kVp, 1.19 HU/°C (95% CI 1.08–1.31) for VMI at 40 keV and 1.35 HU/°C (95% CI 1.27–1.42) for VMI at 140 keV. For muscle, the temperature sensitivities were 0.38 HU/°C (95% CI 0.34–0.42) for 120 kVp, 0.47 HU/ °C (95% CI 0.34–0.59) for VMI at 40 keV and 0.37 HU/ °C (95% CI 0.33–0.42) for VMI at 140 keV. For bone, the temperature sensitivities were 1.15 HU/°C (95% CI 0.87–1.43) for 120 kVp, 3.58 HU/°C (95% CI 3.07– 4.09) for VMI at 50 keV and 0.35 HU/°C (95% CI 0.20–0.49) for VMI at 140 keV.

The application of ASIR-V and DLIR had no significant impact on the temperature sensitivity (Fig. 5 left). The reconstruction kernel has no effect on conventional polychromatic 120-kVp measurements (max difference 0.01 HU/°C); however, it can affect the temperature sensitivity for VMI (Fig. 5 right). For fat, the temperature sensitivity varies between 0.01 HU/°C (VMI at 70 keV) and 0.08 HU/°C (VMI at 40 keV: standard 1.19 HU/°C, bone 1.15 HU/°C, boneplus 1.15 HU/°C, soft 1.17 HU/ °C, detail 1.23 HU/°C; VMI at 120 keV: standard 1.34 HU/°C, bone 1.36 HU/°C, boneplus 1.36 HU/°C, soft 1.35 HU/°C, detail 1.33 HU/°C).

Fig. 3 Polychromatic 120-kVp images and virtual monoenergetic images (VMI) with three types of reconstruction (FBP, ASIR-V 100%, DLIR-H) are shown.

Fig. 4 Left dependency of the CT number on the temperature with linear regression analysis for virtual monoenergetic images (VMI) and polychromatic 120-kVp images. Right results of temperature sensitivity with 95% (gray) and 50% (black) confidence interval (CI) and coefficient of determination for fat, muscle and bone.

For muscle, compared to the standard kernel, only the bone kernel has a negative influence on the temperature sensitivity (VMI at 40 keV: standard 0.47 HU/°C, bone 0.39 HU/°C, boneplus 0.38 HU/°C; VMI at 70 keV: standard 0.40 HU/°C, bone 0.37 HU/°C, boneplus 0.37 HU/°C; VMI at 120 keV: standard 0.38 HU/°C, bone 0.37 HU/°C, boneplus 0.37 HU/°C). In contrast, the bone kernel may have advantages for representation of bones (VMI at 40 keV: standard 0.31 HU/°C, bone 0.51 HU/°C, boneplus 0.94 HU/°C, soft 0.30 HU/°C, detail 0.37 HU/°C; VMI at 70 keV: standard 1.72 HU/°C, bone 1.60 HU/°C, boneplus 1.67 HU/°C, soft 1.74 HU/°C, detail 1.70 HU/°C; VMI at 120 keV: standard 0.51 HU/°C, bone 0.55 HU/°C, boneplus 0.53 HU/°C, soft 0.50 HU/°C, detail 0.53 HU/°C).

Discussion

For evaluating the sensitivity of main human tissues in CT based thermometry, our study is one of the first integrating VMI. VMI proved suitable for CT-based thermometry and enable a higher temperature sensitivity (fat 8%, muscle 24%,  bone 211%) compared to conventional polychromatic 120kVp measurements. In addition, with conventional polychromatic measurements, there are no parameters to influence the temperature sensitivity.Tan et al [12] in 2019 did not  find any statistically significant difference in CT numbers with varied detection parameters for the same temperature setting, so that the scan parameters such as, e.g., voltage, tube current, pitch, gantry rotation time and slice thickness had no influence on the result of thermometry. With VMI, we found new parameters to modulate the temperature sensitivity. An interesting point is that, depending on the type of tissue, a certain monoenergy should be preferred (fat VMI at 140keV, muscle VMI at 40 keV and bone VMI at 50 keV).

Fig. 5 Results of temperature sensitivity with 95% (gray) and 50% (black) confidence interval (CI) for left different reconstruction types (FBP, ASIR-V 50%, ASIR-V 100%, DLIR-L, DLIR-M, DLIR-H) and right different reconstruction kernels (standard, bone, boneplus, soft, detail).

Furthermore, additional reconstructions can be obtained by dual-energy measurements. This leads to more detailed information about the material composition, e.g. material density images and data files (see Fig.1 options). However, VMI partly shows larger fluctuation of CT numbers, due to a slightly larger CI for the temperature sensitivity (Fig. 4 right). One reason could be beam hardening artefacts of bone, which also affect the measurements of fat and muscle. Separate calibration curves, depending on the tissue and monoenergy of the VMI, are required.

For clinical use, CT-based thermometry benefits from modern reconstruction algorithms such as ASIR-V and DLIR, which can drastically reduce the radiation dose [14, 16]. These kinds of reconstruction algorithms and the modification of other protocol parameters have no negative influence on CT numbers [12, 14, 17]. However, ASIR-V and DLIR reconstruction methods can reduce beam hardening artefacts (e.g. of bones) and thus improve the result for CT-based thermometry. Furthermore, the image noise can be drastically reduced, allowing more stable CT numbers for smaller ROIs.

The DLIR setup currently just works with the standard kernel. In future, more kernels will be added for DLIR from manufacturer, which should provide clear statements to kernel influences on CT-based thermometry in further investigations. A major problem of CT-based thermometry is the reproducibility of the results [13]. For this reason, we obtained multiple measurements for our investigations. The numerous repetitions lead to very reliable and stable correlations between temperature and CT number. In clinical application of CT-based thermometry, a slice should be measured multiple times to average the CT numbers.

The measured temperature sensitivity for the muscle imitating salt-water solution corresponds to literature sources [4, 5]. Fat density varies naturally between 0.90 and 0.97 g/cm³ [23]. Therefore, the homogeneous beeswax (density 0.95 g/cm³) is used as a homogeneous tissue equivalent to avoid influences due to density changes. We measure CT numbers between −53 HU and −97 HU for a temperature range from 18 to 50 °C. Porcine fat tissue samples, as comparable biotic tissue, show CT numbers between −105 and −135 HU for a temperature range from 20 to 45 °C [24]. The measurements with beeswax are limited to 50°C. Hot temperatures lead to rapid melting and faster changes in CT numbers than predicted by the linear regression (Fig. 4 left). Excluding the critical 50°C measurement, the temperature sensitivity is 1.02 HU/°C (120 kVp) and 1.13 HU/°C (VMI at 140 keV) for fat. Thus, VMI perform 11% better compared to 120 kVp.

There are some limitations to the present study. One of the major limitations is that the experiments were performed on a tissue-equivalent phantom. For example, physiological processes (e.g. blood perfusion) may cause local alterations in data points due to thermal conductivity and heat redistribution. In addition, we could not heat the phantom above 50 °C, due to the melting point of beeswax. However, it is possible to use the linear regression shown in Fig. 4 to estimate the CT numbers at higher temperatures.

In conclusion, with the application of dual-energy and VMI, CT-based thermometry has parameters to influence and improve temperature sensitivity. Using ASIR-V and DLIR has no negative influence on the results, opening up the potential of drastic dose reductions. The temperatureregulating water circuit system enables reproducible measurements of CT numbers as a function of temperature, e.g. for bone, fat and muscle tissue equivalents. The developed setup is a groundwork for further research concerning complex test objects regarding CT-based thermometry.

Funding: Open Access funding enabled and organized by Projekt DEAL.

Declarations

Guarantor: The scientific guarantor of this publication is Andreas Heinrich.

Conflict of Interest: The authors declare no competing interests.

Statistics and Biometry: At least one of the authors (Andreas Heinrich) has significant statistical expertise.

Informed Consent: Written informed consent was not required (phantom measurement).

Ethical Approval: Institutional Review Board approval was not required (phantom measurement).

Methodology

• prospective
• experimental
• performed at one institution

Open Access: This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Originally published at: European Radiology https://doi.org/10.1007/s00330-021-08206-z

References

1. Fani F, Schena E, Saccomandi P, Silvestri S (2014) CT-based thermometry: an overview. Int J Hyperthermia 30:219–227

2. Schena E, Giurazza F, Massaroni C, Fong Y, Park JJ, Saccomandi P (2017) Thermometry based on computed tomography images during microwave ablation: trials on ex vivo porcine liver. In: 2017 IEEE International Instrumentation and Measurement Technology Conference (I2MTC), pp 1–6

3. Paul J, Vogl TJ, Chacko A (2015) Dual energy computed tomography thermometry during hepatic microwave ablation in an exvivo porcine model. Phys Med 31:683–691

4. Jenne JW, Bahner M, Spoo J et al (1997) CT on-line monitoring of HIFU therapy. In: IEEE Ultrasonics Symposium Proceedings An International Symposium (Cat No97CH36118), vol 1372, pp 1377–1380

5. Fallone BG, Moran PR, Podgorsak EB (1982) Noninvasive thermometry with a clinical x-ray CT scanner. Med Phys 9:715–721

6. Bruners P, Levit E, Penzkofer T et al (2010) Multi-slice computed tomography: a tool for non-invasive temperature measurement? Int J Hyperthermia 26:359–365

7. Pandeya GD, Klaessens JHGM, Greuter MJW et al (2011) Feasibility of computed tomography based thermometry during interstitial laser heating in bovine liver. Eur Radiol 21:1733–1738

8. Bruners P, Pandeya GD, Levit E et al (2012) CT-based temperature monitoring during hepatic RF ablation: feasibility in an animal model. Int J Hyperthermia 28:55–61

9. Schena E, Massaroni C, Giurazza F et al (2016) Feasibility assessment and analysis of thermal sensitivity of CT-thermometry during microwave ablation of ex vivo porcine kidneys. In: 2016 Nanotechnology for Instrumentation and Measurement (NANOfIM), pp 5–16

10. Ruiter SJS, Heerink WJ, de Jong KP (2019) Liver microwave ablation: a systematic review of various FDA-approved systems. Eur Radiol 29:4026–4035

11. Kokuryo D, Kumamoto E, Kuroda K (2020) Recent technological advancements in thermometry. Adv Drug Deliv Rev 163-164:19– 39

12. Tan D, Mohamad NA, Wong YH et al (2019) Experimental assessment on feasibility of computed tomography-based thermometry for radiofrequency ablation on tissue equivalent polyacrylamide phantom. Int J Hyperthermia 36:553–560

13. Pohlan J, Kress W, Hermann KG et al (2020) Computed tomography thermography for ablation zone prediction in microwave ablation and cryoablation: advantages and challenges in an ex vivo porcine liver model. J Comput Assist Tomogr 44:744–749

14. Fan J, Yue M, Melnyk R (2014) Benefits of ASiR-V reconstruction for reducing patient radiation dose and preserving diagnostic quality in CT exams. White paper, GE Healthcare

15. Kim HG, Lee H-J, Lee S-K, Kim HJ, Kim M-J (2017) Head CT: image quality improvement with ASIR-V using a reduced radiation dose protocol for children. Eur Radiol 27:3609–3617

16. Greffier J, Hamard A, Pereira F et al (2020) Image quality and dose reduction opportunity of deep learning image reconstruction algorithm for CT: a phantom study. Eur Radiol 30:3951–3959. https://doi.org/10.1007/s00330-020-06724-w

17. JHsieh J, Liu E, Nett B, Tang J, Thibault JB, Sahney S (2019) A new era of image reconstruction: TrueFidelity™. Technical white paper on deep learning image reconstruction. GE Healthcare

18. Akagi M, Nakamura Y, Higaki T et al (2019) Deep learning reconstruction improves image quality of abdominal ultra-highresolution CT. Eur Radiol 29:6163–6171

19. Li J, Dong D, Fang M et al (2020) Dual-energy CT–based deep learning radiomics can improve lymph node metastasis risk prediction for gastric cancer. Eur Radiol 30:2324–2333

20. Rassouli N, Chalian H, Rajiah P, Dhanantwari A, Landeras L (2017) Assessment of 70-keV virtual monoenergetic spectral images in abdominal CT imaging: a comparison study to conventional polychromatic 120-kVp images. Abdom Radiol (NY) 42:2579– 2586

21. Hubig M, Schenkl S, Muggenthaler H et al (2018) Fully automatic CT-histogram-based fat estimation in dead bodies. Int J Legal Med 132:563–577. https://doi.org/10.1007/s00414-017-1757-5

22. Schenkl S, Muggenthaler H, Hubig M et al (2017) Automatic CTbased finite element model generation for temperature-based death time estimation: feasibility study and sensitivity analysis. Int J Leg Med 131:699–712

23. Lochner S, Graw M (2015) Physikalische Eigenschaften menschlichen Weichgewebes. Vergleichende Untersuchung an Erwachsenen und Kindern. 10. gemeinsames Symposium der DGVP und DGVM am 5. und 6. September 2014 in München: S.83.

24. Weiss N, Sosna J, Goldberg SN, Azhari H (2014) Non-invasive temperature monitoring and hyperthermic injury onset detection using X-ray CT during HIFU thermal treatment in ex vivo fatty tissue. Int J Hyperthermia 30:119–125