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.

Hyperthermia treatment consists of elevating the temperature of the tumor to increase the effectiveness of radiotherapy and chemotherapy. Hyperthermia treatment planning (HTP) is an important tool to optimize treatment quality using pre-treatment temperature predictions. The accuracy of these predictions depends on modeling uncertainties such as tissue properties and positioning. In this study, we evaluated if HTP accuracy improves when the patient is imaged inside the applicator at the start of treatment. Because perfusion is a major uncertainty source, the importance of accurate treatment position and anatomy was evaluated using different perfusion values. Volunteers were scanned using MR imaging without ('planning setup') and with the MR-compatible hyperthermia device ('treatment setup'). Temperature-based quality indicators were used to assess the differences between the standard, apparent and the optimized hyperthermia dose. We conclude that pre-treatment imaging can improve HTP predictions accuracy but also, that tissue perfusion modelling is crucial if temperature-based optimization is applied.

Purpose: We studied the differences between planning and treatment position, their impact on the accuracy of hyperthermia treatment planning (HTP) predictions, and the relevance of including true treatment anatomy and position in HTP based on magnetic resonance (MR) images. Materials and methods: All volunteers were scanned with an MR-compatible hyperthermia device, including a filled waterbolus, to replicate the treatment setup. In the planning setup, the volunteers were scanned without the device to reproduce the imaging in the current HTP. First, we used rigid registration to investigate the patient position displacements between the planning and treatment setup. Second, we performed HTP for the planning anatomy at both positions and the treatment mimicking anatomy to study the effects of positioning and anatomy on the quality of the simulated hyperthermia treatment. Treatment quality was evaluated using SAR-based parameters. Results: We found an average displacement of 2 cm between planning and treatment positions. These displacements caused average absolute differences of ∼12% for TC25 and 10.4%–15.9% in THQ. Furthermore, we found that including the accurate treatment position and anatomy in treatment planning led to an improvement of 2% in TC25 and 4.6%–10.6% in THQ. Conclusions: This study showed that precise patient position and anatomy are relevant since these affect the accuracy of HTP predictions. The major part of improved accuracy is related to implementing the correct position of the patient in the applicator. Hence, our study shows a clear incentive to accurately match the patient position in HTP with the actual treatment.

Purpose: MR thermometry (MRT) enables noninvasive temperature monitoring during hyperthermia treatments. MRT is already clinically applied for hyperthermia treatments in the abdomen and extremities, and devices for the head are under development. In order to optimally exploit MRT in all anatomical regions, the best sequence setup and post-processing must be selected, and the accuracy needs to be demonstrated. Methods: MRT performance of the traditionally used double-echo gradient-echo sequence (DE-GRE, 2 echoes, 2D) was compared to multi-echo sequences: a 2D fast gradient-echo (ME-FGRE, 11 echoes) and a 3D fast gradient-echo sequence (3D-ME-FGRE, 11 echoes). The different methods were assessed on a 1.5 T MR scanner (GE Healthcare) using a phantom cooling down from 59 °C to 34 °C and unheated brains of 10 volunteers. In-plane motion of volunteers was compensated by rigid body image registration. For the ME sequences, the off-resonance frequency was calculated using a multi-peak fitting tool. To correct for B0 drift, the internal body fat was selected automatically using water/fat density maps. Results: The accuracy of the best performing 3D-ME-FGRE sequence was 0.20 °C in phantom (in the clinical temperature range) and 0.75 °C in volunteers, compared to DE-GRE values of 0.37 °C and 1.96 °C, respectively. Conclusion: For hyperthermia applications, where accuracy is more important than resolution or scan-time, the 3D-ME-FGRE sequence is deemed the most promising candidate. Beyond its convincing MRT performance, the ME nature enables automatic selection of internal body fat for B0 drift correction, an important feature for clinical application.

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.

Temperature monitoring plays a central role in improving clinical effectiveness of adjuvant hyperthermia. The potential of magnetic resonance thermometry for treatment monitoring purposes led to several MR-guided hyperthermia approaches. However, the proposed solutions were sub-optimal due to technological and intrinsic limitations. These hamper achieving target conformal heating possibilities (applicator limitations) and accurate thermometry (inadequate signal-to-noise-ratio (SNR)). In this work, we studied proof of principle of a dual-function hyperthermia approach based on a coil array (64 MHz, 1.5 T) that is integrated in-between a phased array for heating (434 MHz) for maximum signal receive in order to improve thermometry accuracy. Hereto, we designed and fabricated a superficial hyperthermia mimicking planar array setup to study the most challenging interactions of generic phased-array setups in order to validate the integrated approach. Experiments demonstrated that the setup complies with the superficial hyperthermia guidelines for heating and is able to improve SNR at 2-4 cm depth by 17%, as compared to imaging using the body coil. Hence, the results showed the feasibility of our dual-function MR-guided hyperthermia approach as basis for the development of application specific setups.

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.

In hyperthermia, the general opinion is that pre-treatment optimization of treatment settings requires a patient-specific model. For deep pelvic hyperthermia treatment planning (HTP), tissue models comprising four tissue categories are currently discriminated. For head and neck HTP, we found that more tissues are required for increasing accuracy. In this work, we evaluated the impact of the number of segmented tissues on the predicted specific absorption rate (SAR) for the pelvic region. Highly detailed anatomical models of five healthy volunteers were selected from a virtual database. For each model, seven lists with varying levels of segmentation detail were defined and used as an input for a modeling study. SAR changes were quantified using the change in target-to-hotspot-quotient and maximum SAR relative differences, with respect to the most detailed patient model. The main finding of this study was that the inclusion of high water content tissues in the segmentation may result in a clinically relevant impact on the SAR distribution and on the predicted hyperthermia treatment quality when considering our pre-established thresholds. In general, our results underline the current clinical segmentation protocol and help to prioritize any improvements.

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 standard-ized 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.

Multichannel image registration is an important challenge in medical image analysis. Multichannel images result from modalities such as dual-energy CT or multispectral microscopy. Besides, multichannel feature images can be derived from acquired images, for instance, by applying multiscale feature banks to the original images to register. Multichannel registration techniques have been proposed, but most of them are applicable to only two multichannel images at a time. In the present study, we propose to formulate multichannel registration as a groupwise image registration problem. In this way, we derive a method that allows the registration of two or more multichannel images in a fully symmetric manner (i.e., all images play the same role in the registration procedure), and therefore, has transitive consistency by definition. The method that we introduce is applicable to any number of multichannel images, any number of channels per image, and it allows to take into account correlation between any pair of images and not just corresponding channels. In addition, it is fully modular in terms of dissimilarity measure, transformation model, regularisation method, and optimisation strategy. For two multimodal datasets, we computed feature images from the initially acquired images, and applied the proposed registration technique to the newly created sets of multichannel images. MIND descriptors were used as feature images, and we chose total correlation as groupwise dissimilarity measure. Results show that groupwise multichannel image registration is a competitive alternative to the pairwise multichannel scheme, in terms of registration accuracy and insensitivity towards registration reference spaces.

Image registration is an important task in medical image analysis. Whereas most methods are designed for the registration of two images (pairwise registration), there is an increasing interest in simultaneously aligning more than two images using groupwise registration. Multimodal registration in a groupwise setting remains difficult, due to the lack of generally applicable similarity metrics. In this work, a novel similarity metric for such groupwise registration problems is proposed. The metric calculates the sum of the conditional entropy between each image in the group and a representative template image constructed iteratively using principal component analysis. The proposed metric is validated in extensive experiments on synthetic and intrasubject clinical image data. These experiments showed equivalent or improved registration accuracy compared to other state-of-the-art (dis)similarity metrics and improved transformation consistency compared to pairwise mutual information.

In quantitative magnetic resonance imaging (qMRI), quantitative tissue properties can be estimated by fitting a signal model to the voxel intensities of a series of images acquired with different settings. To obtain reliable quantitative measures, it is necessary that the qMRI images are spatially aligned so that a given voxel corresponds in all images to the same anatomical location. The objective of the present study is to describe and evaluate a novel automatic groupwise registration technique using a dissimilarity metric based on an approximated form of total correlation. The proposed registration method is applied to five qMRI datasets of various anatomical locations, and the obtained registration performances are compared to these of a conventional pairwise registration based on mutual information. The results show that groupwise total correlation yields better registration performances than pairwise mutual information. This study also establishes that the formulation of approximated total correlation is quite analogous to two other groupwise metrics based on principal component analysis (PCA). Registration performances of total correlation and these two PCA-based techniques are therefore compared. The results show that total correlation yields performances that are analogous to these of the PCAbased techniques. However, compared to these PCA-based metrics, total correlation has two main advantages. Firstly, it is directly derived from a multivariate form of mutual information, while the PCA-based metrics were obtained empirically. Secondly, total correlation has the advantage of requiring no user-defined parameter.