Comparison of MR-thermography and planning calculations in phantoms

Adapting Temperature Predictions to MR Imaging in Treatment Position to Improve Simulation-Guided Hyperthermia for Cervical Cancer

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.

Assessment of the thermal tissue models for the head and neck hyperthermia treatment planning

PURPOSE

To compare different thermal tissue models for head and neck hyperthermia treatment planning, and to assess the results using predicted and measured applied power data from clinical treatments.

METHODS

Three commonly used temperature models from literature were analysed: "constant baseline", "constant thermal stress" and "temperature dependent". Power and phase data of 93 treatments of 20 head and neck patients treated with the HYPERcollar3D applicator were used. The impact on predicted median temperature T50 inside the target region was analysed with maximum allowed temperature of 44 °C in healthy tissue. The robustness of predicted T50 for the three models against the influence of blood perfusion, thermal conductivity and the assumed hotspot temperature level was analysed.

RESULTS

We found an average predicted T50 of 41.0 ± 1.3 °C (constant baseline model), 39.9 ± 1.1 °C (constant thermal stress model) and 41.7 ± 1.1 °C (temperature dependent model). The constant thermal stress model resulted in the best agreement between the predicted power (P = 132.7 ± 45.9 W) and the average power measured during the hyperthermia treatments (P = 129.1 ± 83.0 W).

CONCLUSION

The temperature dependent model predicts an unrealistically high T50. The power values for the constant thermal stress model, after scaling simulated maximum temperatures to 44 °C, matched best to the average measured powers. We consider this model to be the most appropriate for temperature predictions using the HYPERcollar3D applicator, however further studies are necessary for developing of robust temperature model for tissues during heat stress.

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.

Systematic review of pre-clinical and clinical devices for magnetic resonance-guided radiofrequency hyperthermia

Clinical trials have demonstrated the therapeutic benefits of adding radiofrequency (RF) hyperthermia (HT) as an adjuvant to radio- and chemotherapy. However, maximum utilization of these benefits is hampered by the current inability to maintain the temperature within the desired range. RF HT treatment quality is usually monitored by invasive temperature sensors, which provide limited data sampling and are prone to infection risks. Magnetic resonance (MR) temperature imaging has been developed to overcome these hurdles by allowing noninvasive 3D temperature monitoring in the target and normal tissues. To exploit this feature, several approaches for inserting the RF heating devices into the MR scanner have been proposed over the years. In this review, we summarize the status quo in MR-guided RF HT devices and analyze trends in these hybrid hardware configurations. In addition, we discuss the various approaches, extract best practices and identify gaps regarding the experimental validation procedures for MR - RF HT, aimed at converging to a common standard in this process.

Quantitative, Multi-institutional Evaluation of MR Thermometry Accuracy for Deep-Pelvic MR-Hyperthermia Systems Operating in Multi-vendor MR-systems Using a New Anthropomorphic Phantom

Clinical outcome of hyperthermia depends on the achieved target temperature, therefore target conformal heating is essential. Currently, invasive temperature probe measurements are the gold standard for temperature monitoring, however, they only provide limited sparse data. In contrast, magnetic resonance thermometry (MRT) provides unique capabilities to non-invasively measure the 3D-temperature. This study investigates MRT accuracy for MR-hyperthermia hybrid systems located at five European institutions while heating a centric or eccentric target in anthropomorphic phantoms with pelvic and spine structures. Scatter plots, root mean square error (RMSE) and Bland-Altman analysis were used to quantify accuracy of MRT compared to high resistance thermistor probe measurements. For all institutions, a linear relation between MRT and thermistor probes measurements was found with R² (mean ± standard deviation) of 0.97 ± 0.03 and 0.97 ± 0.02, respectively for centric and eccentric heating targets. The RMSE was found to be 0.52 ± 0.31 °C and 0.30 ± 0.20 °C, respectively. The Bland-Altman evaluation showed a mean difference of 0.46 ± 0.20 °C and 0.13 ± 0.08 °C, respectively. This first multi-institutional evaluation of MR-hyperthermia hybrid systems indicates comparable device performance and good agreement between MRT and thermistor probes measurements. This forms the basis to standardize treatments in multi-institution studies of MR-guided hyperthermia and to elucidate thermal dose-effect relations.

The potential of constrained SAR focusing for hyperthermia treatment planning: analysis for the head & neck region

Clinical trials have shown that hyperthermia is a potent adjuvant to conventional cancer treatments, but the temperatures currently achieved in the clinic are still suboptimal. Hyperthermia treatment planning simulations have potential to improve the heating profile of phased-array applicators. An important open challenge is the development of an effective optimization procedure that enables uniform heating of the target region while keeping temperature below a threshold in healthy tissues. In this work, we analyzed the effectiveness and efficiency of a recently proposed optimization approach, i.e. focusing via constrained power optimization (FOCO), using 3D simulations of twelve clinical patient specific models. FOCO performance was compared against a clinically used particle swarm based optimization approach. Evaluation metrics were target coverage at the 25% iso-SAR level, target hotspot quotient, median target temperature (T50) and computational requirements. Our results show that, on average, constrained power focusing performs slightly better than the clinical benchmark ([Formula: see text]T50 [Formula: see text] °C), but outperforms this clinical benchmark for large target volumes ([Formula: see text]40 cm[Formula: see text], [Formula: see text]T50 [Formula: see text] °C). In addition, the results are achieved in a shorter time ([Formula: see text]%) and are repeatable because the approach is formulated as a convex optimization problem.

Optimization of Single Voxel MR Spectroscopy Sequence Parameters and Data Analysis Methods for Thermometry in Deep Hyperthermia Treatments

OBJECTIVE

The difference in the resonance frequency of water and methylene moieties of lipids quantifies in magnetic resonance spectroscopy the absolute temperature using a predefined calibration curve. The purpose of this study was the investigation of peak evaluation methods and the magnetic resonance spectroscopy sequence (point-resolved spectroscopy) parameter optimization that enables thermometry during deep hyperthermia treatments.

MATERIALS AND METHODS

Different Lorentz peak-fitting methods and a peak finding method using singular value decomposition of a Hankel matrix were compared. Phantom measurements on organic substances (mayonnaise and pork) were performed inside the hyperthermia 1.5-T magnetic resonance imaging system for the parameter optimization study. Parameter settings such as voxel size, echo time, and flip angle were varied and investigated.

RESULTS

Usually all peak analyzing methods were applicable. Lorentz peak-fitting method in MATLAB proved to be the most stable regardless of the number of fitted peaks, yet the slowest method. The examinations yielded an optimal parameter combination of 8 cm³ voxel volume, 55 millisecond echo time, and a 90° excitation pulse flip angle.

CONCLUSION

The Lorentz peak-fitting method in MATLAB was the most reliable peak analyzing method. Measurements in homogeneous and heterogeneous phantoms resulted in optimized parameters for the magnetic resonance spectroscopy sequence for thermometry.

Simulation techniques in hyperthermia treatment planning

Abstract Clinical trials have shown that hyperthermia (HT), i.e. an increase of tissue temperature to 39-44 °C, significantly enhance radiotherapy and chemotherapy effectiveness [1]. Driven by the developments in computational techniques and computing power, personalised hyperthermia treatment planning (HTP) has matured and has become a powerful tool for optimising treatment quality. Electromagnetic, ultrasound, and thermal simulations using realistic clinical set-ups are now being performed to achieve patient-specific treatment optimisation. In addition, extensive studies aimed to properly implement novel HT tools and techniques, and to assess the quality of HT, are becoming more common. In this paper, we review the simulation tools and techniques developed for clinical hyperthermia, and evaluate their current status on the path from 'model' to 'clinic'. In addition, we illustrate the major techniques employed for validation and optimisation. HTP has become an essential tool for improvement, control, and assessment of HT treatment quality. As such, it plays a pivotal role in the quest to establish HT as an efficacious addition to multi-modality treatment of cancer.

Thermometry of red blood cell concentrate: magnetic resonance decoding warm up process

PURPOSE

Temperature is a key measure in human red blood cell concentrate (RBC) quality control. A precise description of transient temperature distributions in RBC units removed from steady storage exposed to ambient temperature is at present unknown. Magnetic resonance thermometry was employed to visualize and analyse RBC warm up processes, to describe time courses of RBC mean, surface and core temperatures by an analytical model, and to determine and investigate corresponding model parameters.

METHODS

Warm-up processes of 47 RBC units stored at 1-6°C and exposed to 21.25°C ambient temperature were investigated by proton resonance frequency thermometry. Temperature distributions were visualized and analysed with dedicated software allowing derivation of RBC mean, surface and core temperature-time courses during warm up. Time-dependence of mean temperature was assumed to fulfil a lumped capacitive model of heat transfer. Time courses of relative surface and core temperature changes to ambient temperature were similarly assumed to follow shifted exponential decays characterized by a time constant and a relative time shift, respectively.

RESULTS

The lumped capacitive model of heat transfer and shifted exponential decays described time-dependence of mean, surface and core temperatures close to perfect (mean R(2) were 0.999±0.001, 0.996±0.004 and 0.998±0.002, respectively). Mean time constants were τmean = 55.3±3.7 min, τsurface = 41.4±2.9 min and τcore = 76.8±7.1 min, mean relative time shifts were Δsurface = 0.07±0.02 and Δcore = 0.04±0.01. None of the constants correlated significantly with temperature differences between ambient and storage temperature.

CONCLUSION

Lumped capacitive model of heat transfer and shifted exponential decays represent simple analytical formulas to describe transient mean, surface and core temperatures of RBC during warm up, which might be a helpful tool in RBC temperature monitoring and quality control. Independence of constants on differences between ambient and storage temperature suggests validity of models for arbitrary storage and ambient temperatures.

Fat-referenced MR thermometry in the breast and prostate using IDEAL

PURPOSE

To demonstrate a three-echo fat-referenced MR thermometry technique that estimates and corrects for time-varying phase disturbances in heterogeneous tissues.

MATERIALS AND METHODS

Fat protons do not exhibit a temperature-dependent frequency shift. Fat-referenced thermometry methods exploit this insensitivity and use the signal from fat to measure and correct for magnetic field disturbances. In this study, we present a fat-referenced method that uses interpolation of the fat signal to correct for phase disturbances in fat free regions. Phantom and ex vivo tissue cool-down experiments were performed to evaluate the accuracy of this method in the absence of motion. Non-heated in vivo imaging of the breast and prostate was performed to demonstrate measurement robustness in the presence of systemic and motion-induced field disturbances. Measurement accuracy of the method was compared to conventional proton resonance frequency shift MR thermometry.

RESULTS

In the ex vivo porcine tissue experiment, maximum measurement error of the fat-referenced method was reduced 42% from 3.3 to 1.9°C when compared to conventional MR thermometry. In the breasts, measurement errors were reduced by up to 70% from 6.4 to 1.9°C.

CONCLUSION

Ex vivo and in vivo results show that the proposed method reduces measurement errors in the heterogeneous tissue experiments when compared to conventional MR thermometry.

A heterogeneous human tissue mimicking phantom for RF heating and MRI thermal monitoring verification

This paper describes a heterogeneous phantom that mimics a human thigh with a deep-seated tumor, for the purpose of studying the performance of radiofrequency (RF) heating equipment and non-invasive temperature monitoring with magnetic resonance imaging (MRI). The heterogeneous cylindrical phantom was constructed with an outer fat layer surrounding an inner core of phantom material mimicking muscle, tumor and marrow-filled bone. The component materials were formulated to have dielectric and thermal properties similar to human tissues. The dielectric properties of the tissue mimicking phantom materials were measured with a microwave vector network analyzer and impedance probe over the frequency range of 80-500 MHz and at temperatures of 24, 37 and 45 °C. The specific heat values of the component materials were measured using a differential scanning calorimeter over the temperature range of 15-55 °C. The thermal conductivity value was obtained from fitting the curves obtained from one-dimensional heat transfer measurement. The phantom was used to verify the operation of a cylindrical four-antenna annular phased array extremity applicator (140 MHz) by examining the proton resonance frequency shift (PRFS) thermal imaging patterns for various magnitude/phase settings (including settings to focus heating in tumors). For muscle and tumor materials, MRI was also used to measure T1/T2* values (1.5 T) and to obtain the slope of the PRFS phase change versus temperature change curve. The dielectric and thermal properties of the phantom materials were in close agreement to well-accepted published results for human tissues. The phantom was able to successfully demonstrate satisfactory operation of the tested heating equipment. The MRI-measured thermal distributions matched the expected patterns for various magnitude/phase settings of the applicator, allowing the phantom to be used as a quality assurance tool. Importantly, the material formulations for the various tissue types may be used to construct customized phantoms that are tailored for different anatomical sites.

Mathematical formulation and analysis of the nonlinear system reconstruction of the online image-guided adaptive control of hyperthermia

PURPOSE

A nonlinear system reconstruction can theoretically provide timely system reconstruction when designing a real-time image-guided adaptive control for multisource heating for hyperthermia. This clinical need motivates an analysis of the essential mathematical characteristics and constraints of such an approach.

METHODS

The implicit function theorem (IFT), the Karush-Kuhn-Tucker (KKT) necessary condition of optimality, and the Tikhonov-Phillips regularization (TPR) were used to analyze and determine the requirements of the optimal system reconstruction. Two mutually exclusive generic approaches were analyzed to reconstruct the physical system: The traditional full reconstruction and the recently suggested partial reconstruction. Rigorous mathematical analysis based on IFT, KKT, and TPR was provided for all four possible nonlinear reconstructions: (1) Nonlinear noiseless full reconstruction, (2) nonlinear noisy full reconstruction, (3) nonlinear noiseless partial reconstruction, and (4) nonlinear noisy partial reconstruction, when a class of nonlinear formulations of system reconstruction is employed.

RESULTS

Effective numerical algorithms for solving each of the aforementioned four nonlinear reconstructions were introduced and formal derivations and analyses were provided. The analyses revealed the necessity of adding regularization when partial reconstruction is used. Regularization provides the theoretical support for one to uniquely reconstruct the optimal system. It also helps alleviate the negative influences of unavoidable measurement noise. Both theoretical analysis and numerical examples showed the importance of having a good initial guess for accomplishing nonlinear system reconstruction.

CONCLUSIONS

Regularization is mandatory for partial reconstruction to make it well posed. The Tikhonov-Phillips regularized Gauss-Newton algorithm has nice theoretical performance for partial reconstruction of systems with and without noise. The Levenberg-Marquardt algorithm is a more robust algorithmic option compared to the Gauss-Newton algorithm for nonlinear full reconstruction. A severe limitation of nonlinear reconstruction is the time consuming calculations required for the derivatives of temperatures to unknowns. Developing a method of model reduction or implementing a parallel algorithm can resolve this. The results provided herein are applicable to hyperthermia with blood perfusion nonlinearly depending on temperature and in the presence of thermally significant blood vessels.

Accuracy of real time noninvasive temperature measurements using magnetic resonance thermal imaging in patients treated for high grade extremity soft tissue sarcomas

PURPOSE

To establish accuracy of real time noninvasive temperature measurements using magnetic resonance thermal imaging in patients treated for high grade extremity soft tissue sarcomas.

METHODS

Protocol patients with advanced extremity sarcomas were treated with external beam radiation therapy and hyperthermia. Invasive temperature measures were compared to noninvasive magnetic resonance thermal imaging (MRTI) at 1.5 T performed during hyperthermia. Volumetric temperature rise images were obtained using the proton resonance frequency shift (PRFS) technique during heating in a 140 MHz miniannular phased array applicator. MRTI temperature changes were compared to invasive measurements of temperature with a multisensor fiber optic probe inside a #15 g catheter in the tumor. Since the PRFS technique is sensitive to drifts in the primary imaging magnetic field, temperature change distributions were corrected automatically during treatment using temperature-stable reference materials to characterize field changes in 3D. The authors analyzed MRT images and compared, in evaluable treatments, MR-derived temperatures to invasive temperatures measured in extremity sarcomas. Small regions of interest (ROIs) were specified near each invasive sensor identified on MR images. Temperature changes in the interstitial sensors were compared to the corresponding ROI PRFS-based temperature changes over the entire treatment and over the steady-state period. Nonevaluable treatments (motion/imaging artifacts, noncorrectable drifts) were not included in the analysis.

RESULTS

The mean difference between MRTI and interstitial probe measurements was 0.91 degrees C for the entire heating time and 0.85 degrees C for the time at steady state. These values were obtained from both tumor and normal tissue ROIs. When the analysis is done on just the tumor ROIs, the mean difference for the whole power on time was 0.74 degrees C and during the period of steady state was 0.62 degrees C.

CONCLUSIONS

The data show that for evaluable treatments, excellent correlation (deltaT < 1 degrees C) of MRTI-ROI and invasive measurements can be achieved, but that motion and other artifacts are still serious challenges that must be overcome in future work.

A literature survey on indicators for characterisation and optimisation of SAR distributions in deep hyperthermia, a plea for standardisation

PURPOSE

To evaluate the predictive value of SAR indicators by assessing the correlation of a SAR indicator with the corresponding predicted temperature. Ultimately, this should lead to a number of verified SAR indicators for characterization and optimization of a predicted SAR distribution.

METHODS

A literature survey is followed by an evaluation of the SAR indicators on their functionality, using a set of heuristic classification criteria. To obtain an objective assessment of the predictive value for SAR characterisation, all SAR indicators are evaluated by correlating the value of the SAR indicator to the predicted target temperature when heated with the BSD2000 Sigma 60 applicator. Two methods were followed. First, the specificity of the SAR indicator to target temperature was assessed for each of the 36 patient-specific models, using 30 randomly chosen phase and amplitude settings. Secondly, each SAR indicator was used as a goal function to assess its suitability for optimisation purposes.

RESULTS

Only a selected number of SAR indicators correlate well with tumour/target-temperature. Hence, for target-related properties, an adequate set of SAR indicators is found in the literature. For hotspots, modifications are desirable. For optimisation purposes, improved objective functions have been defined.

CONCLUSIONS

From the correlation of the SAR indicators with tumour temperature, a preferred set of SAR indicators is derived: For target heating, 'average SAR ratio', 'Hotspot-target SAR ratio', and 'homogeneity coefficient' provide suitable objective criteria, while for hotspot reduction, 'Hotspot-target SAR ratio' is considered the most useful indicator. For optimisation procedures, 'Hotspot-target SAR ratio' is currently the most suitable objective function.

Real-time MRI-guided hyperthermia treatment using a fast adaptive algorithm

Magnetic resonance (MR) imaging is promising for monitoring and guiding hyperthermia treatments. The goal of this work is to investigate the stability of an algorithm for online MR thermal image guided steering and focusing of heat into the target volume. The control platform comprised a four-antenna mini-annular phased array (MAPA) applicator operating at 140 MHz (used for extremity sarcoma heating) and a GE Signa Excite 1.5 T MR system, both of which were driven by a control workstation. MR proton resonance frequency shift images acquired during heating were used to iteratively update a model of the heated object, starting with an initial finite element computed model estimate. At each iterative step, the current model was used to compute a focusing vector, which was then used to drive the next iteration, until convergence. Perturbation of the driving vector was used to prevent the process from stalling away from the desired focus. Experimental validation of the performance of the automatic treatment platform was conducted with two cylindrical phantom studies, one homogeneous and one muscle equivalent with tumor tissue (conductivity 50% higher) inserted, with initial focal spots being intentionally rotated 90 degrees and 50 degrees away from the desired focus, mimicking initial setup errors in applicator rotation. The integrated MR-HT treatment platform steered the focus of heating into the desired target volume in two quite different phantom tissue loads which model expected patient treatment configurations. For the homogeneous phantom test where the target was intentionally offset by 90 degrees rotation of the applicator, convergence to the proper phase focus in the target occurred after 16 iterations of the algorithm. For the more realistic test with a muscle equivalent phantom with tumor inserted with 50 degrees applicator displacement, only two iterations were necessary to steer the focus into the tumor target. Convergence improved the heating efficacy (the ratio of integral temperature in the tumor to integral temperature in normal tissue) by up to six-fold, compared to the first iteration. The integrated MR-HT treatment algorithm successfully steered the focus of heating into the desired target volume for both the simple homogeneous and the more challenging muscle equivalent phantom with tumor insert models of human extremity sarcomas after 16 and 2 iterations, correspondingly. The adaptive method for MR thermal image guided focal steering shows promise when tested in phantom experiments on a four-antenna phased array applicator.

Regional hyperthermia of the abdomen in conjunction with chemotherapy for peritoneal carcinomatosis: evaluation of two annular-phased-array applicators

BACKGROUND

Peritoneal carcinomatosis is a stage of gynecological and gastrointestinal malignancies with poor prognosis. Options for enhancing the effect of standard chemotherapy, such as aggressive surgery and intraperitoneal chemotherapy, have limitations. In this phase I/II study, we evaluated regional hyperthermia of the pelvis and abdomen using the annular-phased-array technique as an adjunct to chemotherapy.

METHODS

Forty-five patients with peritoneal carcinomatosis (with or without liver metastases) in colorectal cancer (CRC) (n = 16), ovarian cancer (OC) (n = 17), or gastric/pancreatic/biliary cancer (n = 12) underwent standard chemotherapy and regional hyperthermia. Most CRC patients received second-line chemotherapy. All OC patients were platinum resistant. Regional hyperthermia was applied using a SIGMA-60 applicator (OC), a SIGMA-Eye/MR applicator (CRC), or various ring applicators (gastric/pancreatic/biliary cancer).

RESULTS

Abdominal regional hyperthermia was well tolerated, with acceptable acute discomfort and no long-term morbidity. The SIGMA-Eye/MR applicator achieved higher systemic temperatures (associated with higher systemic stress) and more effective heating of the upper abdomen; the SIGMA-60 applicator achieved higher temperatures (and power densities) in the pelvis. Three-year overall survival was encouraging for patients with CRC (22%) and OC (29%) but not gastric/pancreatic/biliary cancer. For the SIGMA-60 applicator (patients with OC), higher measured temperatures at the vaginal stump correlated with better outcome. CONCLUSIONS. The SIGMA-60 and SIGMA-Eye/MR applicators are feasible for abdominal heating and have low toxicity. The SIGMA-60 applicator is specifically suitable for malignancies with high pelvic burden; the SIGMA-Eye/MR applicator better heats the upper abdomen, including the liver. Further randomized investigations are warranted.

Fast temperature optimization of multi-source hyperthermia applicators with reduced-order modeling of 'virtual sources'

The goal of this work is to build the foundation for facilitating real-time magnetic resonance image guided patient treatment for heating systems with a large number of physical sources (e.g. antennas). Achieving this goal requires knowledge of how the temperature distribution will be affected by changing each source individually, which requires time expenditure on the order of the square of the number of sources. To reduce computation time, we propose a model reduction approach that combines a smaller number of predefined source configurations (fewer than the number of actual sources) that are most likely to heat tumor. The source configurations consist of magnitude and phase source excitation values for each actual source and may be computed from a CT scan based plan or a simplified generic model of the corresponding patient anatomy. Each pre-calculated source configuration is considered a 'virtual source'. We assume that the actual best source settings can be represented effectively as weighted combinations of the virtual sources. In the context of optimization, each source configuration is treated equivalently to one physical source. This model reduction approach is tested on a patient upper-leg tumor model (with and without temperature-dependent perfusion), heated using a 140 MHz ten-antenna cylindrical mini-annular phased array. Numerical simulations demonstrate that using only a few pre-defined source configurations can achieve temperature distributions that are comparable to those from full optimizations using all physical sources. The method yields close to optimal temperature distributions when using source configurations determined from a simplified model of the tumor, even when tumor position is erroneously assumed to be approximately 2.0 cm away from the actual position as often happens in practical clinical application of pre-treatment planning. The method also appears to be robust under conditions of changing, nonlinear, temperature-dependent perfusion. The proposed approach of using virtual sources reduces the number of variables that must be optimized to achieve a tumor-focused temperature distribution, thereby reducing the calculation time required in real-time control applications to about 1/3 to 1/4 of that required for full optimization.

Online feedback focusing algorithm for hyperthermia cancer treatment

PURPOSE

Magnetic resonance (MR) imaging is increasingly being utilized to visualize the 3D temperature distribution in patients during treatment with hyperthermia or thermal ablation therapy. The goal of this work is to lay the foundation for improving the localization of heat in tumors with an online focusing algorithm that uses MR images as feedback to iteratively steer and focus heat into the target.

METHODS

The algorithm iteratively updates the model that quantifies the relationship between the source (antenna) settings and resulting tissue temperature distribution. At each step in the iterative process, optimal settings of power and relative phase of each antenna are computed to maximize averaged tumor temperature in the model. The MR-measured thermal distribution is then used to update/correct the model. This iterative procedure is repeated until convergence, i.e. until the model prediction and MR thermal image are in agreement. A human thigh tumor model heated in a 140 MHz four-antenna cylindrical mini-annular phased array is used for numerical validation of the proposed algorithm. Numerically simulated temperatures are used during the iterative process as surrogates for MR thermal images. Gaussian white noise with a standard deviation of 0.3 degrees C and zero mean is added to simulate MRI measurement uncertainty. The algorithm is validated for cases where the source settings for the first iteration are based on erroneous models: (1) tissue property variability, (2) patient position mismatch, (3) a simple idealized patient model built from CT-based actual geometry, and (4) antenna excitation uncertainty due to load dependent impedance mismatch and antenna cross-coupling. Choices of starting heating vector are also validated.

RESULTS

The algorithm successfully steers and focuses a tumor when there is no antenna excitation uncertainty. Temperature is raised to > or = 43 degrees C for more than about 90% of tumor volume, accompanied by less than about 20% of normal tissue volume being raised to a temperature > or = 41 degrees C. However, when there is antenna excitation uncertainty, about 40% to 80% of normal tissue volume is raised to a temperature > or = 41 degrees C. No significant tumor heating improvement is observed in all simulations after about 25 iteration steps.

CONCLUSIONS

A feedback control algorithm is presented and shown to be successful in iteratively improving the focus of tissue heating within a four-antenna cylindrical phased array hyperthermia applicator. This algorithm appears to be robust in the presence of errors in assumed tissue properties, including realistic deviations of tissue properties and patient position in applicator. Only moderate robustness was achieved in the presence of misaligned applicator/tumor positioning and antenna excitation errors resulting from load mismatch or antenna cross coupling.

Implications of clinical RF hyperthermia on protection limits in the RF range

The systemic temperature is meticulously regulated to 37-37.5 degrees C. Organ systems (skin, digestive system, muscles) have a considerable potential to regulate the perfusion for thermal regulation, physical activity, or digestion. While the regulation of the systemic temperature (37.5 degrees C) is quite strict, the tolerance and regulation potential with respect to local heat is more variable. Laboratory studies provided the relationship between thermal doses and cytotoxic effects. Tissue damage for short-term expositions (in the range of minutes) is only possible for temperatures above 50 degrees C. Radiofrequency radiation is utilized in cancer therapy, inducing local tissue temperatures in the range of 40-45 degrees C for 30-60 min. During local hyperthermia (with heated volumes <1 L) specific absorption rates (SARs) of 100-200 W kg, reactive perfusions of 20-40 mL/100 g/min, and tumor temperatures of 42-43 degrees C are achieved. Normally no side effects or damage in the normal tissue, such as muscle or skin, have been seen. During regional hyperthermia, SARs of 30-40 W kg are found in heated volumes of 10 L with temperatures of 41-42 degrees C in tumor-related measurement points. Then the reactive average perfusion is 6-9 mL/100 g/min (mean value 8 mL/100 g/min). Local temperatures even for higher SAR are regulated to values of not more than 40-42 degrees C. For these temperatures no damages in normal tissues have been found after regional hyperthermia in hundreds of patients. We conclude that the thermoregulatory potential for the whole body or large body regions is limited by the cardiac output, which can at least double the output from 5 to 10 L min. Even higher is the potential to compensate in smaller volumes. Here the perfusion in muscle can be increased from the basal value of 2-4 mL/100 g/min more than 5-10-fold.

Quantitative validation of the 3D SAR profile of hyperthermia applicators using the gamma method

For quality assurance of hyperthermia treatment planning systems, quantitative validation of the electromagnetic model of an applicator is essential. The objective of this study was to validate a finite-difference time-domain (FDTD) model implementation of the Lucite cone applicator (LCA) for superficial hyperthermia. The validation involved (i) the assessment of the match between the predicted and measured 3D specific absorption rate (SAR) distribution, and (ii) the assessment of the ratio between model power and real-world power. The 3D SAR distribution of seven LCAs was scanned in a phantom bath using the DASY4 dosimetric measurement system. The same set-up was modelled in SEMCAD X. The match between the predicted and the measured SAR distribution was quantified with the gamma method, which combines distance-to-agreement and dose difference criteria. Good quantitative agreement was observed: more than 95% of the measurement points met the acceptance criteria 2 mm/2% for all applicators. The ratio between measured and predicted power absorption ranged from 0.75 to 0.92 (mean 0.85). This study shows that quantitative validation of hyperthermia applicator models is feasible and is worth considering as a part of hyperthermia quality assurance procedures.