Correction of motion-induced susceptibility artifacts and B0 drift during proton resonance frequency shift-based MR thermometry in the pelvis with background field removal methods

Magnetic resonance thermometry for hyperthermia in the oropharynx region

Magnetic resonance thermometry (MRT) can measure in-vivo 3D-temperature changes in real-time and noninvasively. However, for the oropharynx region and the entire head and neck, motion potentially introduces large artifacts. Considering long treatment times of 60-90 min, this study aims to evaluate whether MRT around the oropharynx is clinically feasible for hyperthermia treatments and quantify the effects of breathing and swallowing on MRT performance. A 3D-ME-FGRE sequence was used in a phantom cooling down and around the oropharynx of five volunteers over ∼75 min. The imaging protocol consisted of imaging with acceleration (ARC = 2), number of image averages (NEX = 1,2 and 3). For volunteers, the acquisitions included a breath-hold scan and scans with deliberate swallowing. MRT performance was quantified in neck muscle, spinal cord and masseter muscle, using mean average error (MAE), mean error (ME) and spatial standard deviation (SD). In phantom, an increase in NEX leads to a significant decrease in SD, but MAE and ME were unchanged. No significant difference was found in volunteers between the different scans. There was a significant difference between the regions evaluated: neck muscle had the best MAE (=1.96 °C) and SD (=0.82 °C), followed by spinal cord (MAE = 3.17 °C, SD = 0.92 °C) and masseter muscle (MAE = 4.53 °C, SD = 1.16 °C). Concerning the ME, spinal cord did best, then neck muscle and masseter muscle, with values of -0.64 °C, 1.15 °C and -3.05 °C respectively. Breathing, swallowing, and different ways of imaging (acceleration and NEX) do not significantly influence the MRT performance in the oropharynx region. The ROI selected however, leads to significant differences.

Technical advances in motion-robust MR thermometry

Proton resonance frequency shift (PRFS) MR thermometry is the most common method used in clinical thermal treatments because of its fast acquisition and high sensitivity to temperature. However, motion is the biggest obstacle in PRFS MR thermometry for monitoring thermal treatment in moving organs. This challenge arises because of the introduction of phase errors into the PRFS calculation through multiple methods, such as image misregistration, susceptibility changes in the magnetic field, and intraframe motion during MRI acquisition. Various approaches for motion correction have been developed for real-time, motion-robust, and volumetric MR thermometry. However, current technologies have inherent trade-offs among volume coverage, processing time, and temperature accuracy. These tradeoffs should be considered and chosen according to the thermal treatment application. In hyperthermia treatment, precise temperature measurements are of increased importance rather than the requirement for exceedingly high temporal resolution. In contrast, ablation procedures require robust temporal resolution to accurately capture a rapid temperature rise. This paper presents a comprehensive review of current cutting-edge MRI techniques for motion-robust MR thermometry, and recommends which techniques are better suited for each thermal treatment. We expect that this study will help discern the selection of motion-robust MR thermometry strategies and inspire the development of motion-robust volumetric MR thermometry for practical use in clinics.

MR Thermometry during Transcranial MR Imaging-Guided Focused Ultrasound Procedures: A Review

Interest in transcranial MR imaging-guided focused ultrasound procedures has recently grown. These incisionless procedures enable precise focal ablation of brain tissue using real-time monitoring by MR thermometry. This article will provide an updated review on clinically applicable technical underpinnings and considerations of proton resonance frequency MR thermometry, the most common clinically used MR thermometry sequence.

Straightforward Magnetic Resonance Temperature Measurements Combined with High Frame Rate and Magnetic Susceptibility Correction

Proton resonance frequency shift (PRFS) is an MRI-based simple temperature mapping method that exhibits higher spatial and temporal resolution than temperature mapping methods based on T1 relaxation time and diffusion. PRFS temperature measurements are validated against fiber-optic thermal sensors (FOSs). However, the use of FOSs may introduce temperature errors, leading to both underestimation and overestimation of PRFS measurements, primarily due to material susceptibility changes caused by the thermal sensors. In this study, we demonstrated susceptibility-corrected PRFS (scPRFS) with a high frame rate and accuracy for suitably distributed temperatures. A single-echo-based background removal technique was employed for phase variation correction, primarily owing to magnetic susceptibility, which enabled fast temperature mapping. The scPRFS was used to validate the temperature fidelity by comparing the temperatures of fiber-optic sensors and conventional PRFS through phantom-mimicked human and ex vivo experiments. This study demonstrates that scPRFS measurements in agar-gel are in good agreement with the thermal sensor readings, with a root mean square error (RMSE) of 0.33-0.36 °C in the phantom model and 0.12-0.16 °C in the ex vivo experiment. These results highlight the potential of scPRFS for precise thermal monitoring and ablation in both low- and high-temperature non-invasive therapies.

Image registration and mutual thresholding enable low interimage variability across dynamic MRI measurements of supraclavicular brown adipose tissue during mild cold exposure

PURPOSE

Activated brown adipose tissue (BAT) enhances lipid catabolism and improves cardiometabolic health. Quantitative MRI of the fat fraction (FF) of supraclavicular BAT (scBAT) is a promising noninvasive measure to assess BAT activity but suffers from high scan variability. We aimed to test the effects of coregistration and mutual thresholding on the scan variability in a fast (1 min) time-resolution MRI protocol for assessing scBAT FF changes during cold exposure.

METHODS

Ten volunteers (age 24.8 ± 3.0 years; body mass index 21.2 ± 2.1 kg/m² ) were scanned during thermoneutrality (32°C; 10 min) and mild cold exposure (18°C; 60 min) using a 12-point gradient-echo sequence (70 consecutive scans with breath-holds, 1.03 min per dynamic). Dynamics were coregistered to the first thermoneutral scan, which enabled drawing of single regions of interest in the scBAT depot. Voxel-wise FF changes were calculated at each time point and averaged across regions of interest. We applied mutual FF thresholding, in which voxels were included if their FF was greater than 30% FF in the reference scan and the registered dynamic. The efficacy of the coregistration was determined by using a moving average and comparing the mean squared error of residuals between registered and nonregistered data. Registered scBAT ΔFF was compared with single-scan thresholding using the moving average method.

RESULTS

Registered scBAT ΔFF had lower mean square error values than nonregistered data (0.07 ± 0.05% vs. 0.16 ± 0.14%; p < 0.05), and mutual thresholding reduced the scBAT ΔFF variability by 30%.

CONCLUSION

We demonstrate that coregistration and mutual thresholding improve stability of the data 2-fold, enabling assessment of small changes in FF following cold exposure.

Rectal MRI Interpretation After Neoadjuvant Therapy

In recent years, several key advances in the management of locally advanced rectal cancer have been made, including the implementation of total mesorectal excision as the standard surgical approach; use of neoadjuvant chemoradiotherapy in selected patients with a high risk of local recurrence, and finally, adoption of organ preservation strategies, through either local excision or nonoperative management in selected patients with clinical complete response following neoadjuvant chemoradiotherapy. This review aims to shed light on the role of rectal MRI in the assessment of treatment response after neoadjuvant therapy, which is especially important given the growing feasibility of nonoperative management. First, an overview of current neoadjuvant therapies and response assessment based on digital rectal examination, endoscopy, and MRI will be provided. Second, the use of a high-quality restaging rectal MRI protocol will be presented. Third, a step-by-step approach to assessing treatment response on restaging rectal MRI following neoadjuvant treatment will be outlined, acknowledging challenges faced by radiologists during MRI interpretation. Finally, research related to response assessment will be discussed. LEVEL OF EVIDENCE: 4 TECHNICAL EFFICACY: Stage 3.

A k-space-based method to measure and correct for temporal B0 field variations in MR temperature imaging

PURPOSE

Present a method to use change in phase in repeated Cartesian k-space measurements to monitor the change in magnetic field for dynamic MR temperature imaging.

METHODS

The method is applied to focused ultrasound heating experiments in a gelatin phantom and an ex vivo salt pork sample, without and with simulated respiratory motion.

RESULTS

In each experiment, phase variations due to B0 field drift and respiration were readily apparent in the measured phase difference. With correction, the SD of the temperature over time was reduced from 0.18°C to 0.14°C (no breathing) and from 0.81°C to 0.22°C (with breathing) for the gelatin phantom, and from 0.68°C to 0.13°C (no breathing) and from 1.06°C to 0.17°C (with breathing) for the pork sample. The accuracy in nonheated regions, assessed as the RMS error deviation from 0°C, improved from 1.70°C to 1.11°C (no breathing) and from 4.73°C to 1.47°C (with breathing) for the gelatin phantom, and from 5.95°C to 0.88°C (no breathing) and from 13.40°C to 1.73°C (with breathing) for the pork sample. The correction did not affect the temperature measurement accuracy in the heated regions.

CONCLUSION

This work demonstrates that phase changes resulting from variations in B0 due to drift and respiration, commonly seen in MR thermometry applications, can be measured directly from 3D Cartesian acquisition methods. The correction of temporal field variations using the presented technique improved temperature accuracy, reduced variability in nonheated regions, and did not reduce accuracy in heated regions.

Integrated thermal and magnetic susceptibility modeling for air-motion artifact correction in proton resonance frequency shift thermometry

PURPOSE

Hyperthermia treatments are successful adjuvants to conventional cancer therapies in which the tumor is sensitized by heating. To monitor and guide the hyperthermia treatment, measuring the tumor and healthy tissue temperature is important. The typical clinical practice heavily relies on intraluminal probe measurements that are uncomfortable for the patient and only provide spatially sparse temperature information. A solution may be offered through recent advances in magnetic resonance thermometry, which allows for three-dimensional internal temperature measurements. However, these measurements are not widely used in the pelvic region due to a low signal-to-noise ratio and presence of image artifacts.

METHODS

To advance the clinical integration of magnetic resonance-guided cancer treatments, we consider the problem of removing air-motion-induced image artifacts. Thereto, we propose a new combined thermal and magnetic susceptibility model-based temperature estimation scheme that uses temperature estimates to improve the removal of air-motion-induced image artifacts. The method is experimentally validated using a dedicated phantom that enables the controlled injection of air-motion artifacts and with in vivo thermometry from a clinical hyperthermia treatment.

RESULTS

We showed, using probe measurements in a heated phantom, that our method reduced the mean absolute error (MAE) by 58% compared to the state-of-the-art near a moving air volume. Moreover, with in vivo thermometry our method obtained a MAE reduction between 17% and 95% compared to the state-of-the-art.

CONCLUSION

We expect that the combined thermal and magnetic susceptibility modeling used in model-based temperature estimation can significantly improve the monitoring in hyperthermia treatments and enable feedback strategies to further improve MR-guided hyperthermia cancer treatments.

POD-Kalman filtering for improving noninvasive 3D temperature monitoring in MR-guided hyperthermia

Background

During resonance frequency (RF) hyperthermia treatment, the temperature of the tumor tissue is elevated to the range of 39–44°C. Accurate temperature monitoring is essential to guide treatments and ensure precise heat delivery and treatment quality. Magnetic resonance (MR) thermometry is currently the only clinical method to measure temperature noninvasively in a volume during treatment. However, several studies have shown that this approach is not always sufficiently accurate for thermal dosimetry in areas with motion, such as the pelvic region. Model‐based temperature estimation is a promising approach to correct and supplement 3D online temperature estimation in regions where MR thermometry is unreliable or cannot be measured. However, complete 3D temperature modeling of the pelvic region is too complex for online usage.

Purpose

This study aimed to evaluate the use of proper orthogonal decomposition (POD) model reduction combined with Kalman filtering to improve temperature estimation using MR thermometry. Furthermore, we assessed the benefit of this method using data from hyperthermia treatment where there were limited and unreliable MR thermometry measurements.

Methods

The performance of POD–Kalman filtering was evaluated in several heating experiments and for data from patients treated for locally advanced cervical cancer. For each method, we evaluated the mean absolute error (MAE) concerning the temperature measurements acquired by the thermal probes, and we assessed the reproducibility and consistency using the standard deviation of error (SDE). Furthermore, three patient groups were defined according to susceptibility artifacts caused by the level of intestinal gas motion to assess if the POD–Kalman filtering could compensate for missing and unreliable MR thermometry measurements.

Results

First, we showed that this method is beneficial and reproducible in phantom experiments. Second, we demonstrated that the combined method improved the match between temperature prediction and temperature acquired by intraluminal thermometry for patients treated for locally advanced cervical cancer. Considering all patients, the POD–Kalman filter improved MAE by 43% (filtered MR thermometry = 1.29°C, POD–Kalman filtered temperature = 0.74°C). Moreover, the SDE was improved by 47% (filtered MR thermometry = 1.16°C, POD–Kalman filtered temperature = 0.61°C). Specifically, the POD–Kalman filter reduced the MAE by approximately 60% in patients whose MR thermometry was unreliable because of the great amount of susceptibilities caused by the high level of intestinal gas motion.

Conclusions

We showed that the POD–Kalman filter significantly improved the accuracy of temperature monitoring compared to MR thermometry in heating experiments and hyperthermia treatments. The results demonstrated that POD–Kalman filtering can improve thermal dosimetry during RF hyperthermia treatment, especially when MR thermometry is inaccurate.

Phase-independent thermometry by Z-spectrum MR imaging

PURPOSE

Z-spectrum imaging, defined as the consecutive collection of images after saturating over a range of frequency offsets, has been recently proposed as a method to measure the fat-water fraction by the simultaneous detection of fat and water resonances. By incorporating a binomial pulse irradiated at each offset before the readout, the spectral selectivity of the sequence can be further amplified, making it possible to monitor the subtle proton resonance frequency shift that follows a change in temperature.

METHODS

We tested the hypothesis in aqueous and cream phantoms and in healthy mice, all under thermal challenge. The binomial module consisted of 2 sinc-shaped pulses of opposite phase separated by a delay. Such a delay served to spread out off-resonance spins, with the resulting excitation profile being a periodic function of the delay and the chemical shift.

RESULTS

During heating experiments, the water resonance shifted downfield, and by fitting the curve to a sine function it was possible to quantify the change in temperature. Results from Z-spectrum imaging correlated linearly with data from conventional MRI techniques like T₁ mapping and phase differences from spoiled GRE.

CONCLUSION

Because the measurement is performed solely on magnitude images, the technique is independent of phase artifacts and is therefore applicable in mixed tissues (e.g., fat). We showed that Z-spectrum imaging can deliver reliable temperature change measurement in both muscular and fatty tissues.

Susceptibility artifact correction in MR thermometry for monitoring of mild radiofrequency hyperthermia using total field inversion

PURPOSE

MR temperature monitoring of mild radiofrequency hyperthermia (RF-HT) of cancer exploits the linear resonance frequency shift of water with temperature. Motion-induced susceptibility distribution changes cause artifacts that we correct here using the total field inversion (TFI) approach.

METHODS

The performance of TFI was compared to two background field removal (BFR) methods: Laplacian boundary value (LBV) and projection onto dipole fields (PDF). Data sets with spatial susceptibility change and B 0 -drift were simulated, phantom heating experiments were performed, four volunteer data sets at thermoneutral conditions as well as data from one cervical cancer, two sarcoma, and one seroma patients undergoing mild RF-HT were corrected using the proposed methods.

RESULTS

Simulations and phantom heating experiments revealed that using BFR or TFI preserves temperature-induced phase change, while removing susceptibility artifacts and B 0 -drift. TFI resulted in the least cumulative error for all four volunteers. Temperature probe information from four patient data sets were best depicted by TFI-corrected data in terms of accuracy and precision. TFI also performed best in case of the sarcoma treatment without temperature probe.

CONCLUSION

TFI outperforms previously suggested BFR methods in terms of accuracy and robustness. While PDF consistently overestimates susceptibility contribution, and LBV removes valuable pixel information, TFI is more robust and leads to more accurate temperature estimations.

MR Thermometry Accuracy and Prospective Imaging-Based Patient Selection in MR-Guided Hyperthermia Treatment for Locally Advanced Cervical Cancer

Simple Summary

Monitoring and controlling the temperature distribution combined with precise energy delivery are key components for hyperthermia treatment success. Magnetic resonance (MR) imaging is used clinically to monitor the temperature of the treated volume non-invasively. However, there are no comprehensive systematic studies on MR thermometry accuracy during deep pelvic hyperthermia, and the few investigational studies suffer from a high probability of bias due to lacking objective criteria for data inclusion. This study presents the first systematic analysis and defines an imaging-based criterion for prospective patient selection to standardize clinical MR thermometry accuracy assessments.

Abstract

The efficacy of a hyperthermia treatment depends on the delivery of well-controlled heating; hence, accurate temperature monitoring is essential for ensuring effective treatment. For deep pelvic hyperthermia, there are no comprehensive and systematic reports on MR thermometry. Moreover, data inclusion generally lacks objective selection criteria leading to a high probability of bias when comparing results. Herein, we studied whether imaging-based data inclusion predicts accuracy and could serve as a tool for prospective patient selection. The accuracy of the MR thermometry in patients with locally advanced cervical cancer was benchmarked against intraluminal temperature. We found that gastrointestinal air motion at the start of the treatment, quantified by the Jaccard similarity coefficient, was a good predictor for MR thermometry accuracy. The results for the group that was selected for low gastrointestinal air motion improved compared to the results for all patients by 50% (accuracy), 26% (precision), and 80% (bias). We found an average MR thermometry accuracy of 2.0 °C when all patients were considered and 1.0 °C for the selected group. These results serve as the basis for comprehensive benchmarking of novel technologies. The Jaccard similarity coefficient also has good potential to prospectively determine in which patients the MR thermometry will be valuable.

Effect of gastrointestinal gas on the temperature distribution in pancreatic cancer hyperthermia treatment planning

PURPOSE

In pancreatic cancer treatment, hyperthermia can be added to increase efficacy of chemo- and/or radiotherapy. Gas in stomach, intestines and colon is often in close proximity to the target volume. We investigated the impact of variations in gastrointestinal gas (GG) on temperature distributions during simulated hyperthermia treatment (HT).

METHODS

We used sets of one CT and eight cone-beam CT (CBCT) scans obtained prior to/during fractionated image-guided radiotherapy in four pancreatic cancer patients. In Plan2Heat, we simulated locoregional heating by an ALBA-4D phased array radiofrequency system and calculated temperature distributions for (i) the segmented CT (sCT), (ii) sCT with GG replaced by muscle (sCT₀), (iii) sCT₀ with eight different GG distributions as visible on CBCT inserted (sCT_CBCT). We calculated cumulative temperature-volume histograms for the clinical target volume (CTV) for all ten temperature distributions for each patient and investigated the relationship between GG volume and change in ΔT₅₀ (temperature increase at 50% of CTV volume). We determined location and volume of normal tissue receiving a high thermal dose.

RESULTS

GG volume on CBCT varied greatly (9-991 cm³). ΔT₅₀ increased for increasing GG volume; maximum ΔT₅₀ difference per patient was 0.4-0.6 °C. The risk for GG-associated treatment-limiting hot spots appeared low. Normal tissue high-temperature regions mostly occurred anteriorly; their volume and maximum temperature showed moderate positive correlations with GG volume, while fat-muscle interfaces were associated with higher risks for hot spots.

CONCLUSIONS

Considerable changes in volume and position of gastrointestinal gas can occur and are associated with clinically relevant tumor temperature differences.

Clinical Performance and Future Potential of Magnetic Resonance Thermometry in Hyperthermia

Simple Summary

Hyperthermia is a treatment for cancer patients, which consists of heating the body to 43 °C. The temperature during treatment is usually measured by placing temperature probes intraluminal or invasively. The only clinically used option to measure temperature distributions non-invasively and in 3D is by MR thermometry (MRT). However, in order to be able to replace conventional temperature probes, MRT needs to become more reliable. In this review paper, we propose standardized performance thresholds for MRT, based on our experience of treating nearly 4000 patients. We then review the literature to assess to what extent these requirements are already being met in the clinic today and identify common problems. Lastly, using pre-clinical results in the literature, we assess where the biggest potential is to solve the problems identified. We hope that by standardizing MRT parameters as well as highlighting current and promising developments, progress in the field will be accelerated.

Abstract

Hyperthermia treatments in the clinic rely on accurate temperature measurements to guide treatments and evaluate clinical outcome. Currently, magnetic resonance thermometry (MRT) is the only clinical option to non-invasively measure 3D temperature distributions. In this review, we evaluate the status quo and emerging approaches in this evolving technology for replacing conventional dosimetry based on intraluminal or invasively placed probes. First, we define standardized MRT performance thresholds, aiming at facilitating transparency in this field when comparing MR temperature mapping performance for the various scenarios that hyperthermia is currently applied in the clinic. This is based upon our clinical experience of treating nearly 4000 patients with superficial and deep hyperthermia. Second, we perform a systematic literature review, assessing MRT performance in (I) clinical and (II) pre-clinical papers. From (I) we identify the current clinical status of MRT, including the problems faced and from (II) we extract promising new techniques with the potential to accelerate progress. From (I) we found that the basic requirements for MRT during hyperthermia in the clinic are largely met for regions without motion, for example extremities. In more challenging regions (abdomen and thorax), progress has been stagnating after the clinical introduction of MRT-guided hyperthermia over 20 years ago. One clear difficulty for advancement is that performance is not or not uniformly reported, but also that studies often omit important details regarding their approach. Motion was found to be the common main issue hindering accurate MRT. Based on (II), we reported and highlighted promising developments to tackle the issues resulting from motion (directly or indirectly), including new developments as well as optimization of already existing strategies. Combined, these may have the potential to facilitate improvement in MRT in the form of more stable and reliable measurements via better stability and accuracy.

Correction of motion-induced susceptibility artifacts and B0 drift during proton resonance frequency shift-based MR thermometry in the pelvis with background field removal methods

Purpose

The linear change of the water proton resonance frequency shift (PRFS) with temperature is used to monitor temperature change based on the temporal difference of image phase. Here, the effect of motion‐induced susceptibility artifacts on the phase difference was studied in the context of mild radio frequency hyperthermia in the pelvis.

Methods

First, the respiratory‐induced field variations were disentangled from digestive gas motion in the pelvis. The projection onto dipole fields (PDF) as well as the Laplacian boundary value (LBV) algorithm were applied on the phase difference data to eliminate motion‐induced susceptibility artifacts. Both background field removal (BFR) algorithms were studied using simulations of susceptibility artifacts, a phantom heating experiment, and volunteer and patient heating data.

Results

Respiratory‐induced field variations were negligible in the presence of the filled water bolus. Even though LBV and PDF showed comparable results for most data, LBV seemed more robust in our data sets. Some data sets suggested that PDF tends to overestimate the background field, thus removing phase attributed to temperature. The BFR methods even corrected for susceptibility variations induced by a subvoxel displacement of the phantom. The method yielded successful artifact correction in 2 out of 4 patient treatment data sets during the entire treatment duration of mild RF heating of cervical cancer. The heating pattern corresponded well with temperature probe data.

Conclusion

The application of background field removal methods in PRFS‐based MR thermometry has great potential in various heating applications and body regions to reduce motion‐induced susceptibility artifacts that originate outside the region of interest, while conserving temperature‐induced PRFS. In addition, BFR automatically removes up to a first‐order spatial B₀ drift.