Simultaneous magnetic resonance phase and magnitude temperature maps in muscle

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.

An MRI-Guided Ring High-Intensity Focused Ultrasound System for Noninvasive Breast Ablation

High-intensity focused ultrasound (HIFU) has been used for noninvasive treatment of breast tumors, but the present magnetic resonance imaging (MRI)-guided HIFU (MRI-HIFU) systems encounter skin burn. In this study, a novel MRI-HIFU breast ablation system was developed to improve the above problem. The system consisted of the ring HIFU phased-array transducer, a commercial power amplifier, the mechanical positioner, and the graphical user interface control software. MRI thermometry was also established to monitor the temperature in the HIFU-treated tissue. Ablation of pork and the in vivo rabbit leg were carried out to validate the developed system. Results of fat-surrounding pork ablation showed that the ring HIFU system reached a safe margin of 3 mm without fat burn. Moreover, precision of the positioner moving the HIFU focal zone was within 6% error under MRI circumstances. The representative MRI temperature images show that the peak temperatures among the five ablations ranged between 66 °C and 91 °C, and their thermal doses were over 10000. The system could also ablate the biceps femoris of a rabbit without skin burn to form a lesion of 2.5 mm beneath the skin. With the HIFU dose of 315 W/10 s, the MRI temperature map revealed that the maximum temperature and the thermal dose were 60 °C and 3380, respectively. The MRI-guided ring HIFU system can ablate the target tissue near subcutaneous fat without fat burn. The system prototype is a promising tool for clinical implementation.

Recent technological advancements in thermometry

Thermometry is the key factor for achieving successful thermal therapy. Although invasive thermometry with a probe has been used for more than four decades, this method can only detect the local temperature within the probing volume. Noninvasive temperature imaging using a tomographic technique is ideal for monitoring hot-spot formation in the human body. Among various techniques, such as X-ray computed tomography, microwave tomography, echo sonography, and magnetic resonance (MR) imaging, the proton resonance frequency shift method of MR thermometry is the only method currently available for clinical practice because its temperature sensitivity is consistent in most aqueous tissues and can be easily observed using common clinical scanners. New techniques are being proposed to improve the robustness of this method against tissue motion. MR techniques for fat thermometry were also developed based on relaxation times. One of the latest non-MR techniques to attract attention is photoacoustic imaging.

Diffusion weighted magnetic resonance imaging for temperature measurements in catalyst supports with an axial gas flow

In situ thermometry of catalytic gas phase reactions allows to determine temperature profiles in catalyst beds. Diffusion weighted MRI is proposed as an alternative method for temperature measurements using capillaries filled with different liquids.

Correlation distance dependence of the resonance frequency of intermolecular zero quantum coherences and its implication for MR thermometry

PURPOSE

Because the resonance frequency of water-fat intermolecular zero-quantum coherences (iZQCs) reflects the water-fat frequency separation at the microscopic scale, these frequencies have been proposed and used as a mean to obtain more accurate temperature information. The purpose of this work was to investigate the dependence of the water-fat iZQC resonance frequency on sample microstructure and on the specific choice of the correlation distance.

METHODS

The effect of water-fat susceptibility gradients on the water-methylene iZQC resonance frequency was first computed and then measured for different water-fat emulsions and for a mixture of porcine muscle and fat. Similar measurements were also performed for mixed heteronuclear spin systems.

RESULTS

A strong dependence of the iZQC resonance frequency on the sample microstructure and on the specific choice of the correlation distance was found for spin systems like water and fat that do not mix, but not for spin systems that mix at the molecular level.

CONCLUSIONS

Because water and fat spins do not mix at the molecular level, the water-fat iZQC resonance frequency and its temperature coefficient are not only affected by sample microstructure but also by the specific choice of the correlation distance. Magn Reson Med 79:1429-1438, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Sampling strategies for subsampled segmented EPI PRF thermometry in MR guided high intensity focused ultrasound

PURPOSE

To investigate k-space subsampling strategies to achieve fast, large field-of-view (FOV) temperature monitoring using segmented echo planar imaging (EPI) proton resonance frequency shift thermometry for MR guided high intensity focused ultrasound (MRgHIFU) applications.

METHODS

Five different k-space sampling approaches were investigated, varying sample spacing (equally vs nonequally spaced within the echo train), sampling density (variable sampling density in zero, one, and two dimensions), and utilizing sequential or centric sampling. Three of the schemes utilized sequential sampling with the sampling density varied in zero, one, and two dimensions, to investigate sampling the k-space center more frequently. Two of the schemes utilized centric sampling to acquire the k-space center with a longer echo time for improved phase measurements, and vary the sampling density in zero and two dimensions, respectively. Phantom experiments and a theoretical point spread function analysis were performed to investigate their performance. Variable density sampling in zero and two dimensions was also implemented in a non-EPI GRE pulse sequence for comparison. All subsampled data were reconstructed with a previously described temporally constrained reconstruction (TCR) algorithm.

RESULTS

The accuracy of each sampling strategy in measuring the temperature rise in the HIFU focal spot was measured in terms of the root-mean-square-error (RMSE) compared to fully sampled "truth." For the schemes utilizing sequential sampling, the accuracy was found to improve with the dimensionality of the variable density sampling, giving values of 0.65 °C, 0.49 °C, and 0.35 °C for density variation in zero, one, and two dimensions, respectively. The schemes utilizing centric sampling were found to underestimate the temperature rise, with RMSE values of 1.05 °C and 1.31 °C, for variable density sampling in zero and two dimensions, respectively. Similar subsampling schemes with variable density sampling implemented in zero and two dimensions in a non-EPI GRE pulse sequence both resulted in accurate temperature measurements (RMSE of 0.70 °C and 0.63 °C, respectively). With sequential sampling in the described EPI implementation, temperature monitoring over a 192×144×135 mm3 FOV with a temporal resolution of 3.6 s was achieved, while keeping the RMSE compared to fully sampled "truth" below 0.35 °C.

CONCLUSIONS

When segmented EPI readouts are used in conjunction with k-space subsampling for MR thermometry applications, sampling schemes with sequential sampling, with or without variable density sampling, obtain accurate phase and temperature measurements when using a TCR reconstruction algorithm. Improved temperature measurement accuracy can be achieved with variable density sampling. Centric sampling leads to phase bias, resulting in temperature underestimations.

Toward volumetric MR thermometry with the MASTER sequence

MR temperature monitoring is an indispensable tool for high intensity focused ultrasound. In this paper, a new technique known as MASTER (multiple adjacent slice thermometry with excitation refocusing) is presented which improves the speed and accuracy of multiple-slice MR thermometry. Defocusing the magnetization after exciting a slice allows for multiple slices to be excited concurrently and stored in k-space. The magnetization from each excitation is then refocused and read in sequence. This approach increases TE for each slice, greatly improving temperature SNR as compared to conventional slice interleaving. Gradient sequence design optimization is required to minimize diffusion losses while maintaining high sequence efficiency. Flexibility in selecting position, update rate, accuracy, and voxel size for each slice independently allows for freedom in design to fit different application needs. Results are shown in phantom and in vivo validating the feasibility of the sequence, and comparing it to interleaved GRE. Sample design curves are presented that contrast the MASTER design space with that of interleaved GRE thermometry.

ARFI-prepared MRgHIFU in liver: simultaneous mapping of ARFI-displacement and temperature elevation, using a fast GRE-EPI sequence

MR acoustic radiation force imaging (ARFI) is an elegant adjunct to MR-guided high intensity focused ultrasound for treatment planning and optimization, permitting in situ assessment of the focusing and targeting quality. The thermal effect of high intensity focused ultrasound pulses associated with ARFI measurements is recommended to be monitored on line, in particular when the beam crosses highly absorbent structures or interfaces (e.g., bones or air-filled cavities). A dedicated MR sequence is proposed here, derived from a segmented gradient echo-echo planar imaging kernel by adding a bipolar motion encoding gradient with interleaved alternating polarities. Temporal resolution was reduced to 2.1 s, with in-plane spatial resolution of 1 mm. MR-ARFI measurements were executed during controlled animal breathing, with trans-costal successively steered foci, to investigate the spatial modulation of the focus intensity and the targeting offset. ARFI-induced tissue displacement measurements enabled the accurate localization, in vivo, of the high intensity focused ultrasound focal point in sheep liver, with simultaneous monitoring of the temperature elevation. ARFI-based precalibration of the focal point position was immediately followed by trans-costal MR-guided high intensity focused ultrasound ablation, monitored with a conventional proton resonance frequency shift MR thermometry sequence. The latter MR thermometry sequence had spatial resolution and geometrical distortion identical with the ARFI maps, hence no coregistration was required.

Multiresolution MRI temperature monitoring in a reduced field of view

PURPOSE

The purpose was to develop a new magnetic resonance imaging technique for fast temperature monitoring with extended volume coverage.

MATERIALS AND METHODS

The Multiple Resolutions Along Phase-Encode and Slice-Select Dimensions (MURPS) method was implemented in both a two-dimensional (2D) spoiled gradient echo (SPGR) sequence and a multishot echo-planar imaging (EPI) sequence. Both modified sequences were used to acquire image data from three slices with variable phase-encode resolution and slice thickness. In the SPGR sequence, a 2D resonant frequency pulse was also implemented to enable imaging within a reduced field of view, and this was used to monitor (at 1.5 T) the temperature changes in a live rabbit and in gel phantoms heated by focused ultrasound. A modified EPI sequence was tested during heating of a phantom undergoing motion.

RESULTS

The in vivo experiments demonstrated that temperature changes in unexpected locations away from the focal plane, such as near bone structures, could be detected due to the extra volume coverage afforded by the MURPS method. Temperature changes in a moving phantom were resolved using the MURPS EPI sequence with an acquisition rate of three slices every 300 ms.

CONCLUSION

The MURPS method enables temperature monitoring over multiple slices without loss of temporal resolution compared with single-slice imaging and, if combined with multishot EPI, enables volume temperature monitoring in moving organs.

Intradiscal temperature monitoring using double gradient-echo pulse sequences at 1.0T

PURPOSE

To validate an unspoiled gradient-recalled echo pulse sequence with dual echo acquisition as a means to increase temperature sensitivity while monitoring intradiscal laser ablation therapy.

MATERIALS AND METHODS

Phantom experiments as well as in vitro thermal ablation simulations were performed in an open 1.0T magnetic resonance (MR) scanner. Three methods of noninvasive MR-thermometry based on the signal void decrease caused by T1-relaxation time increase (T1), the temperature-dependent proton resonance frequency (PRF) shift, and a combination of both methods with complex differences (CD) were compared. Temperature accuracy and reliability of temperature distribution were the main assessment criteria.

RESULTS

The optimum temperature sensitivity was found using CD in phantom experiments. During in vitro experiments the PRF showed the smallest margin of error (T1: +/-1.64 degrees C, PRF: +/-1.23 degrees C, CD: +/-1.29 degrees C) and the best qualitative evaluation of temperature.

CONCLUSION

Intradiscal temperature monitoring with an unspoiled dual-echo sequence is most accurate with PRF-thermometry in combination with the long echo time. Magnitude images with an initial short echo time permit high image detail of the heat-induced lesion.

Noninvasive temperature mapping with MRI using chemical shift water-fat separation

Tissues containing both water and lipids, e.g., breast, confound standard MR proton reference frequency-shift methods for mapping temperatures due to the lack of temperature-induced frequency shift in lipid protons. Generalized Dixon chemical shift-based water-fat separation methods, such as GE's iterative decomposition of water and fat with echo asymmetry and least-squares estimation method, can result in complex water and fat images. Once separated, the phase change over time of the water signal can be used to map temperature. Phase change of the lipid signal can be used to correct for non-temperature-dependent phase changes, such as amplitude of static field drift. In this work, an image acquisition and postprocessing method, called water and fat thermal MRI, is demonstrated in phantoms containing 30:70, 50:50, and 70:30 water-to-fat by volume. Noninvasive heating was applied in an Off1-On-Off2 pattern over 50 min, using a miniannular phased radiofrequency array. Temperature changes were referenced to the first image acquisition. Four fiber optic temperature probes were placed inside the phantoms for temperature comparison. Region of interest (ROI) temperature values colocated with the probes showed excellent agreement (global mean +/- standard deviation: -0.09 +/- 0.34 degrees C) despite significant amplitude of static field drift during the experiments.

Real-time MR thermometry for monitoring HIFU ablations of the liver

A high-resolution and high-speed pulse sequence is presented for monitoring high-intensity focused ultrasound ablations in the liver in the presence of motion. The sequence utilizes polynomial-order phase saturation bands to perform outer volume suppression, followed by spatial-spectral excitation and three readout segmented echo-planar imaging interleaves. Images are processed with referenceless thermometry to create temperature-rise images every frame. The sequence and reconstruction were implemented in RTHawk and used to image stationary and moving sonications in a polyacrylamide gel phantom (62.4 acoustic W, 50 sec, 550 kHz). Temperature-rise images were compared between moving and stationary experiments. Heating spots and corresponding temperature-rise plots matched very well. The stationary sonication had a temperature standard deviation of 0.15 degrees C compared to values of 0.28 degrees C and 0.43 degrees C measured for two manually moved sonications at different velocities. Moving the phantom (while not heating) with respect to the transducer did not cause false temperature rises, despite susceptibility changes. The system was tested on nonheated livers of five normal volunteers. The mean temperature rise was -0.05 degrees C, with a standard deviation of 1.48 degrees C. This standard deviation is acceptable for monitoring high-intensity focused ultrasound ablations, suggesting real-time imaging of moving high-intensity focused ultrasound sonications can be clinically possible.

Hyperthermia MRI temperature measurement: evaluation of measurement stabilisation strategies for extremity and breast tumours

PURPOSE

MR thermometry using the proton resonance frequency shift (PRFS) method has been used to measure temperature changes during clinical hyperthermia treatment. However, frequency drift of the MRI system can add large errors to the measured temperature change. These drifts can be measured and corrected using oil references placed around the treatment region. In this study, the number and position of four or more oil references were investigated to obtain a practical approach to correct frequency drift during PRFS thermometry in phantoms and in vivo.

MATERIALS AND METHODS

Experiments were performed in a 140 MHz four antenna mini-annular phased array (MAPA) heat applicator (for treatment of extremity tumours) and an applicator for heating of the breast, with symmetric and asymmetric positioning of the oil references, respectively. Temperature change PRFS images were obtained during an hour or more of measurement with no application of heat. Afterwards, errors in calculating temperature change due to system drift were quantified with and without various oil reference correction arrangements.

RESULTS

Results showed good temperature correction in phantoms and in a human leg, with average errors of 0.28 degrees C and 0.94 degrees C respectively. There was further improvement in the leg when using eight or more oil references, reducing the average error to 0.44 degrees C, while the phantoms showed no significant improvement.

CONCLUSIONS

These results indicate that oil reference correction performs well in vivo, and that eight references can improve the correction by up to 0.5 degrees C compared to four references.

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.

Magnetic resonance imaging guidance for laser photothermal therapy

Temperature distribution is a crucial factor in determining the outcome of laser phototherapy in cancer treatment. Magnetic resonance imaging (MRI) is an ideal method for 3-D noninvasive temperature measurement. A 7.1-T MRI was used to determine laser-induced high thermal gradient temperature distribution of target tissue with high spatial resolution. Using a proton density phase shift method, thermal mapping is validated for in vivo thermal measurement with light-absorbing enhancement dye. Tissue-simulating phantom gels, biological tissues, and tumor-bearing animals were used in the experiments. An 805-nm laser was used to irradiate the samples, with laser power in the range of 1 to 3 W. A clear temperature distribution matrix within the target and surrounding tissue was obtained with a specially developed processing algorithm. The temperature mapping showed that the selective laser photothermal effect could result in temperature elevation in a range of 10 to 45 degrees C. The temperature resolution of the measurement was about 0.37 degrees C with 0.4-mm spatial resolution. The results of this study provide in vivo thermal information and future reference for optimizing laser dosage and dye concentration in cancer treatment.

MR thermometry

Minimally invasive thermal therapy as local treatment of benign and malignant diseases has received increasing interest in recent years. Safety and efficacy of the treatment require accurate temperature measurement throughout the thermal procedure. Noninvasive temperature monitoring is feasible with magnetic resonance (MR) imaging based on temperature-sensitive MR parameters such as the proton resonance frequency (PRF), the diffusion coefficient (D), T1 and T2 relaxation times, magnetization transfer, the proton density, as well as temperature-sensitive contrast agents. In this article the principles of temperature measurements with these methods are reviewed and their usefulness for monitoring in vivo procedures is discussed. Whereas most measurements give a temperature change relative to a baseline condition, temperature-sensitive contrast agents and spectroscopic imaging can provide absolute temperature measurements. The excellent linearity and temperature dependence of the PRF and its near independence of tissue type have made PRF-based phase mapping methods the preferred choice for many in vivo applications. Accelerated MRI imaging techniques for real-time monitoring with the PRF method are discussed. Special attention is paid to acquisition and reconstruction methods for reducing temperature measurement artifacts introduced by tissue motion, which is often unavoidable during in vivo applications.

From the RSNA refresher courses: MR imaging in hyperthermia

There is growing clinical evidence that the combination of radiation therapy and hyperthermia, when delivered at moderate temperatures (40 degrees-45 degrees C) for sustained times (30-90 minutes), is of benefit with regard to palliative relief of cancer, tumor response, local control, and survival. Adequate measurement of the temperature distribution achieved with the hyperthermia is a key element in successful therapy. Thermal dosimetry, even invasive dosimetry, is a complex topic when applied to the heterogeneous tissue of a tumor and associated organ systems. Imaging in hyperthermia therapy is performed primarily for estimation and control of temperature. Magnetic resonance (MR) imaging has unique parameter dependences that make it possible to monitor hyperthermia therapy by detection of proton resonant frequency changes or diffusion coefficient changes. In addition, MR imaging can be used to assess vascular parameters that not only allow selection of suitable patients for therapy but may also allow demonstration of response to therapy. Finally, as the use of thermally sensitive liposomes for delivery of chemotherapeutic agents is developed, MR imaging may allow determination of local drug dose.

Preliminary assessment of one-dimensional MR elastography for use in monitoring focused ultrasound therapy

The purpose of this work is to assess a fast technique that measures tissue stiffness and temperature during focused ultrasound thermal therapy (FUS). A one-dimensional (1D) MR elastography (MRE) pulse sequence was evaluated for the purpose of obtaining rapid measurements of thermally induced changes in tissue stiffness and temperature for monitoring FUS treatments. The accuracy of the 1D measurement was studied by comparing tissue displacements measured by 1D MRE with those measured by the well-established 2D MRE pulse sequence. The reproducibility of the 1D MRE measurement was assessed, in gel phantoms and ex vivo porcine tissue, for varied FUS intensity levels (31.5-199.9 W cm(-2)) and over a range of displacements at the focus (0.1-1 microm). Temperature elevations in agarose gel phantoms were measured using 1D MRE and calibrated using fiberoptic-thermometer-based measurements. The 1D MRE displacement measurements are highly correlated with those obtained with the 2D technique (R(2) = 0.88-0.93), indicating that 1D MRE can successfully measure tissue displacement. Ten repeated trials at each FUS power level yielded a minimum detectable displacement change of 0.2 microm in phantoms and 0.4 microm in tissue (at 95% confidence level). The 1D MRE temperature measurements correlated well with temperature changes measured simultaneously with fiberoptic thermometers (R(2) = 0.97). The 1D MRE technique is capable of detecting tissue displacements as low as 0.4 microm, which is an order of magnitude smaller than 5 microm displacements expected during FUS therapy (Le et al 2005 AIP Conf. Proc.: Ther. Ultrasound 829 186-90). Additionally, 1D MRE was shown to provide adequate measurements of temperature elevations in tissue. These findings indicate that 1D MRE may be an effective tool for monitoring FUS treatments.

Near-real-time feedback control system for liver thermal ablations based on self-referenced temperature imaging

Our challenge was to design and implement a dedicated temperature imaging feedback control system to guide and assist in a thermal liver ablation procedure in a double-donut 0.5T open MR scanner. This system has near-real-time feedback capability based on a newly developed "self-referenced" temperature imaging method using "moving-slab" and complex-field-fitting techniques. Two phantom validation studies and one ex vivo experiment were performed to compare the newly developed self-referenced method with the conventional subtraction method and evaluate the ability of the feedback control system in the same MR scanner. The near-real-time feedback system was achieved by integrating the following primary functions: (1) imaging of the moving organ temperature; (2) on-line needle tip tracking; (3) automatic turn-on/off the heating devices; (4) a Windows operating system-based novel user-interfaces. In the first part of the validation studies, microwave heating was applied in an agar phantom using a fast spoiled gradient recalled echo in a steady state sequence. In the second part of the validation and ex vivo study, target visualization, treatment planning and monitoring, and temperature and thermal dose visualization with the graphical user interface of the thermal ablation software were demonstrated. Furthermore, MR imaging with the "self-referenced" temperature imaging method has the ability to localize the hot spot in the heated region and measure temperature elevation during the experiment. In conclusion, we have demonstrated an interactively controllable feedback control system that offers a new method for the guidance of liver thermal ablation procedures, as well as improving the ability to assist ablation procedures in an open MR scanner.

Feasibility of simultaneous temperature and tissue stiffness detection by MRE

Temperature and tissue stiffness are two indices that can be used to monitor MRI-guided focused ultrasound thermal therapy. It would be beneficial to have both measures available to monitor treatment progression during thermal therapy. MR Elastography (MRE) has already been shown to provide tissue stiffness information; the purpose of this work is to demonstrate how temperature can be derived from the same MRE data acquisition. MRE data were acquired from 1.5% agarose phantoms and ex vivo porcine muscle tissue (from a grocery store) while they were heated slowly. The temperatures were measured using a fluorescent thermometer. The phase average from the MRE acquisition was used to calculate the phase shift induced by the proton resonance frequency shift associated with the temperature change. The results show that the phase shift due to temperature extracted from MRE data correlate well with the temperature change recorded by thermometer, yielding a temperature coefficient of -0.0096 ppm/ degrees C for the agarose phantom, and -0.0103 ppm/ degrees C for the ex vivo porcine tissue. These results indicate that it is possible to simultaneously measure both temperature and tissue shear stiffness using a new method of MRE data reconstruction.

Radiofrequency power deposition utilizing thermal imaging

Wavelength effects influence radiofrequency (RF) power deposition distributions and limit magnetic resonance (MR) medical applications at very high magnetic fields. The power depositions in spherical saline gel phantoms were deduced from proton resonance shift thermal maps at both 1.5 T and 3.0 T over a range of conductivities. Phase differences before and after RF heating were measured for both a quadrature head coil and a circular surface coil. A long echo time (TE) pulse sequence with a 3D phase unwrap algorithm provided increased thermal sensitivity. The measured thermal maps agreed with a model of eddy-current heating by circularly polarized oscillating RF fields in a conducting dielectric sphere. At 3.0 T, thermal maps were acquired with a <0.32 degrees C temperature rise at 4 W. Proton resonance shift thermal maps provided a measure of hot spots in very-high-field MR imaging (MRI), in which both the phase sensitivity and signal-to-noise ratio (SNR) were increased. The method provides a means of studying the heat distribution generated by RF coils excited by clinical pulse sequences.

High intensity focused ultrasound ablation of kidney guided by MRI

The effectiveness of magnetic resonance imaging (MRI) to monitor therapeutic protocols of high-intensity focused ultrasound (HIFU), in freshly excised pig kidney cortex is investigated. For high quality imaging, the pulse sequence fast spin echo (FSE) T1- and T2-weighted, and proton density were evaluated. For fast imaging, the pulse sequence T1-weighted fast spoiled gradient (FSPGR) was used. The main goal was to evaluate the MRI detection of large lesions (bigger than 1 cm x 1 cm x 1 cm) that is achieved by moving the transducer in a predetermined pattern. The contrast between lesion and kidney tissue is excellent with either T1-weighted or T2-weighted FSE. With T1-weighted FSE, the best contrast is observed for recovery time (TR) between 200 ms and 400 ms. With T2-weighted FSE best contrast can be achieved for echo time (TE) between 16 and 32 ms. T2-weighted FSE was proven as the best pulse sequence to detect cavitational activity. This advantage is attributed to the significant difference in signal intensity between air spaces and necrotic tissue. Air spaces appear brighter than thermal lesions. Therefore, for therapeutic protocols created using cavitational mode, T2-weighted FSE may be the optimum pulse sequence to use. The proton density pulse sequence does not provide any advantage over the T1- and T2-weighted pulse sequences. Using T1-weighted FSPGR, acquisition time as low as 5 s could be achieved. Good contrast and signal-to-noise ratio (SNR) are achieved with TR = 100 ms and flip angle between 75 to 90 degrees. The above techniques were very successful in detecting large lesion volumes.

Sliding window dual gradient echo (SW-dGRE): T1 and proton resonance frequency (PRF) calibration for temperature imaging in polyacrylamide gel

The aim of the work is to evaluate a magnetic resonance imaging (MRI) thermometry sequence suitable for targeting of focused ultrasound (FUS) when used in vascular occlusion studies. A sliding window dual gradient echo (SW-dGRE) sequence was used. This sequence has the capability of monitoring both T1 relaxation and phase changes, which vary with temperature. Preliminary work involved quantification of the changes in T1 relaxation time with temperature and obtaining the PRF shift coefficient in polyacrylamide gel as it underwent an exothermic reaction during polymerization (avoiding the use of an external heat source). Temperature changes were visualized using thermal maps acquired with the sequence. For FUS guidance a thermal imaging technique is required with a temporal resolution <5 s, a spatial resolution of approximately 1 mm and a temperature resolution of approximately 5 degrees C. The sequence was optimized to improve the CNR (contrast to noise ratio) and SNR (signal to noise ratio) in the phase and magnitude images respectively. The PRF coefficient obtained for the polyacrylamide gel was -9.98 +/- 0.24 ppb degrees C(-1), whilst deltaT1 and temperature change were related by a proportionality factor, the T1 temperature coefficient, of 102.3 +/- 2.9 ms degrees C(-1). The sequence produces an image at every 1.4 s interval. In both magnitude and phase data, the in-plane resolution is +/- 1.2 mm and the temperature resolution is approximately 2 degrees C. The advantage of this sequence is that the temperature obtained from the magnitude data can be confirmed independently using the phase data and vice versa. Thus the sequence can essentially be crosschecked.

MRI thermodosimetry in laser-induced interstitial thermotherapy

BACKGROUND AND OBJECTIVES

The aim of this study was to establish a correlation between a thermal measurement and a magnetic resonance imaging (MRI) signal during laser-induced interstitial thermotherapy (LITT) in liver.

STUDY DESIGN/MATERIALS AND METHODS

Pig liver was irradiated for 15 minutes with a diode laser at two different powers, 0.5 W (450 J) and 1.5 W (1,350 J). Tissue temperature was monitored every 20 seconds using thermocouples. Thermosensitive MRI sequences (T(1)-weighted Turbo-Flash) were acquired with the same irradiation parameters. Correlation between MRI signals (SI) and temperature measures was defined at two different distances from the fiber (5 and 10 mm).

RESULTS

At 0.5 W, temperatures rose progressively up to a maximum increase of 9.5 degrees C at 5 mm and 4 degrees C at 10 mm after 15 minutes. The corresponding MRI signal decreased progressively to -27.6 SI at 5 mm and -18.5 SI at 10 mm. At 1.5 W, temperatures rose dramatically at 5 mm, reaching a plateau. The temperature elevation measured at the end of the irradiation was of 30 degrees C whereas at 10 mm it was only 14.5 degrees C. The MRI signal varied accordingly, remaining inversely proportional to temperature (-76 SI at 5 mm and -35.5 SI at 10 mm).

CONCLUSIONS

An inversely proportional relationship was observed between MRI signal in sequential Turbo-Flash and temperature. MRI should allow to analyze heat diffusion in the liver, and thus to monitor real-time LITT treatments.

Feasibility of Internally Referenced Brain Temperature Imaging with a Metabolite Signal

The feasibility of using a metabolite signal as an internal reference for self-referenced temperature distribution measurement was examined. Line scan echo-planar spectroscopic imaging (LSEPSI) was applied to obtain quick multi-voxel spectroscopic measurements and to avoid possible spectral degradation from motion. Temperature distribution in a rabbit brain in vivo was successfully visualized by means of the chemical shift of water, which was measured by using naturally abundant (up to 10 mM) N-acetyl-aspartate (NAA) as the reference signal. Unlike the phase-mapping approach, this technique does not require a pixel-by-pixel subtraction. Therefore, in theory, it is more resistant to inter-scan motion or changes in susceptibility. The spatial and temporal resolutions of this technique are 1.5 cm3 and 4.5 min. A higher signal-to-noise ratio and optimization of the water and outer-volume suppression capabilities will be required to further enhance the temperature-mapping capabilities.

Fast MRI of RF heating via phase difference mapping

A method is presented for the rapid acquisition of temperature maps derived from phase difference maps. The temperature-dependent chemical shift coefficients (TDCSCs) of various concentrations of aqueous cobalt and dysprosium-based compounds were measured. The largest TDCSC calculated was for 100 mM DyEDTA, which had a TDCSC of -0.09 PPM/K; 160 mM CoCl2 had a TDCSC of -0.04 PPM/K. These temperature-dependent chemical shifts (TDCSs) result in phase changes in the MR signal with changing temperature. Agarose phantoms were constructed with each paramagnetic metal. A fast gradient-echo (FGRE) MR image was acquired to serve as the baseline image. A "test" MRI procedure was then performed on the phantom. Immediately afterwards, a second FGRE MR image was acquired, serving as the probing image. Proper image processing as a phase difference map between the probing image and the baseline image resulted in an image which quantitatively described the temperature increase of the phantom in response to a particular "test" imaging experiment. Applications of this technique in assessing the safety of pulse sequences and MR coils are discussed.

Magnetic resonance imaging-guided focused ultrasound thermal therapy in experimental animal models: correlation of ablation volumes with pathology in rabbit muscle and VX2 tumors

PURPOSE

To further investigate the use of magnetic resonance-guided focused ultrasound therapy (MRgFUS) as a noninvasive alternative to surgery in the local control of soft-tissue tumors by ablating prescribed volumes of VX2 rabbit tumors and comparing with ablation of normal tissue volumes.

MATERIALS AND METHODS

Small, ellipsoidal ablations at shallow depth were created using 5- to 15-second sonication pulses at radio frequency (RF) powers of 50-125 W using a spherical, air-backed transducer operating at 1.463 MHz under MR guidance in a 1.5-T clinical scanner.

RESULTS

Excellent correlation was observed between prescribed treatment volumes, MR thermal dosimetry, post-treatment verification MRI, and histopathology. Multifocal ablations of VX2 tumors in rabbits at depths of up to 2.5 cm resulted in complete ablation of the prescribed treatment volume.

CONCLUSION

MRgFUS is an effective technique for treating tumors in vivo. Techniques developed for treatments in homogeneous tissue volumes are applicable in the more complicated tumor environment if MR temperature feedback is available to modify treatment delivery parameters.

MR monitoring of tumour thermal therapy

Thermal therapy of tumour including hyperthermia and thermal ablation by heat or cold delivery requires on line monitoring. Due to its temperature sensitivity, Magnetic Resonance Imaging (MRI) allows thermal mapping at the time of the treatment. The different techniques of MR temperature monitoring based on water proton resonance frequency (PRF), longitudinal relaxation time T1, diffusion coefficient and MR Spectroscopic Imaging (MRSI) are reviewed and debated. The PRF method appears the most widely used and the most efficient at high magnetic field in spite of important drawbacks. The T1 method is the easiest method of visualisation of qualitative temperature distribution and quantitative measurement seems possible in the tissue surrounding the tumour up to a temperature of 45-65 degrees C. Despite its high temperature sensitivity, application of the diffusion method in vivo is restricted due to its high motion sensitivity. The recent MRSI technique seems very promising provided acquisition times can be reduced. Results from the literature indicate that MR temperature monitoring in vivo can be achieved in vivo with a precision of about 3 degrees C in 13 s for a voxel of 16 mm3 (1.5 x 1.5 x 7 mm) in 1.5 T scanners.

TmDOTA-: a sensitive probe for MR thermometry in vivo

The lanthanide complex, thulium 1,4,7,10-tetraazacyclododecane-1,4,7,10-tetraacetic acid (TmDOTA-), has been investigated as an agent for MR thermometry in vivo. The chemical shifts of the TmDOTA- protons were highly sensitive to temperature at a clinically relevant field strength, yet insensitive to pH and the presence of Ca2+. Given the excellent stability of lanthanide-DOTA complexes and high thermal sensitivity, TmDOTA- is expected to be a good candidate for MR thermometry in vivo.

Temperature quantification using the proton frequency shift technique: In vitro and in vivo validation in an open 0.5 tesla interventional MR scanner during RF ablation

Open magnetic resonance (MR) scanners allow MR-guided targeting of tumors, as well as temperature monitoring of radio frequency (RF) ablation. The proton frequency shift (PFS) technique, an accurate and fast imaging method for temperature quantification, was used to synthesize thermal maps after RF ablation in an open 0.5 T MR system under ex vivo and in vivo conditions. Calibration experiments with 1.5% agarose gel yielded a chemical shift factor of 0.011 +/- 0.001 ppm/ degrees C (r2 = 0.96). Three gradient echo (GRE) pulse sequences were tested for thermal mapping by comparison with fiberoptic thermometer (Luxtron Model 760) readings. Temperature uncertainty decreased from high to low bandwidths (BW): +/-5.9 degrees C at BW = 15.6 kHz, +/-1.4 degrees C at BW = 3.9 kHz, and +/-0.8 degrees C at BW = 2.5 kHz. In vitro experiments (N = 9) in the paraspinal muscle yielded a chemical shift factor of 0.008 +/- 0.001 ppm/ degrees C. Temperature uncertainty was determined as +/-2.7 degrees C (BW = 3.9 kHz, TE = 19.3 msec). The same experiments carried out in the paraspinal muscle (N = 9) of a fully anesthetized pig resulted in a temperature uncertainty of +/-4.3 degrees C (BW = 3.9 kHz, TE = 19.3 msec), which is higher than it is in vitro conditions (P < 0.15). Quantitative temperature monitoring of RF ablation is feasible in a 0.5 T open-configured MR scanner under ex vivo and in vivo conditions using the PFS technique.

Temperature monitoring with line scan echo planar spectroscopic imaging

A new magnetic resonance imaging method, line scan echo planar spectroscopic imaging (LSEPSI), is shown capable of providing rapid, internally referenced temperature monitoring from water and fat chemical shifts.

METHODS

Orthogonal 90 degrees and 180 degrees slice selective RF pulses inclined by 45 degrees from the image plane solicit a spin echo from a tissue column. The echo is read by asymmetric sampling of 32 gradient echoes spaced 1.4-1.8 ms apart. Sixty-four adjacent columns are sequentially sampled in 4.2-6.4 s with 4,096 voxels sampled with voxel volumes of 0.08-0.13 cm3. Mixed mayonnaise/water phantoms were used to correlate LSEPSI-derived chemical shifts and thermocouple-based temperature measurements from 23 to 60 degrees C with a 1.5 T scanner. Measurement artifacts unrelated to temperature were investigated with the phantom, as was the feasibility of applying the sequence in human breast in vivo.

RESULTS

The correlation between LSEPSI and thermocouple-based temperature measurements in the phantom was excellent (r2>0.99). Field drifts affecting the temperature measurements using the water peak alone were corrected by using the water/lipid peak difference. The sequence had an average temperature resolution of 1.4 degrees C in the phantom. The frequency difference measurement reduced the sensitivity to artifacts related to temperature. Both water and lipid peaks were detectable throughout many locations in the breast, suggesting the applicability of LSEPSI in this organ.

DISCUSSION

T1-saturation losses occur in conventional and echo-planar based 2D CSI sequences using phase encoding methods with short TR periods. These losses are eliminated when individual columns are sampled in snapshot fashion with LSEPSI since the effective TR becomes the time between scans rather than excitations. T1 saturation can make small spectral peaks difficult to detect at high temperatures and generally lowers the signal-to-noise ratio of the spectra. The rapid acquisition and insensitivity to T1 saturation effects make LSEPSI an attractive technique for monitoring thermal therapies in breast using the internally referenced fat/water frequency separation.

Image-guided therapy in urology

Technological breakthroughs have advanced the fields of urology, radiology, and minimally invasive surgery. Today, the various imaging modalities are increasingly applied to guiding therapy. Among the procedures now in use or under development are percutaneous cyst drainage or sclerotherapy; tissue ablation with high-intensity focused ultrasound, cold, heat, or photon radiation; and conformal radiation and brachytherapy. As current limitations are overcome, image-guided therapy will expand.

Intraprostatic temperature monitoring during transurethral microwave thermotherapy: status and future developments

Transurethral microwave thermotherapy is being applied as a minimally invasive treatment for alleviating the symptoms of urinary outlet obstruction associated with benign prostatic hyperplasia. Treatment progress has traditionally been guided in its effective power by rectally and urethrally measured temperatures, whereas intraprostatic temperatures would be preferred for feedback purposes. A critical evaluation is presented of intraprostatic thermometry techniques that have been suggested, the techniques currently being used and investigated, and the problems that remain to be solved. Techniques for noninvasive temperature measurement and detecting tissue response during thermal therapy are discussed in more detail. Results presented in the literature have shown magnetic resonance imaging and ultrasonic imaging to be adequate thermometry modalities. For treatment monitoring of transurethral microwave thermotherapy, ultrasonic imaging is especially promising. Future research will indicate whether the promise evolves into a sound clinical technique.

Magnetic resonance temperature imaging for guidance of thermotherapy

Continuous thermometry during a hyperthermic procedure may help to correct for local differences in heat conduction and energy absorption, and thus allow optimization of the thermal therapy. Noninvasive, three-dimensional mapping of temperature changes is feasible with magnetic resonance (MR) and may be based on the relaxation time T(1), the diffusion coefficient (D), or proton resonance frequency (PRF) of tissue water. The use of temperature-sensitive contrast agents and proton spectroscopic imaging can provide absolute temperature measurements. The principles and performance of these methods are reviewed in this paper. The excellent linearity and near-independence with respect to tissue type, together with good temperature sensitivity, make PRF-based temperature MRI the preferred choice for many applications at mid to high field strength (>/= 1 T). The PRF methods employ radiofrequency spoiled gradient-echo imaging methods. A standard deviation of less than 1 degrees C, for a temporal resolution below 1 second and a spatial resolution of about 2 mm, is feasible for a single slice for immobile tissues. Corrections should be made for temperature-induced susceptibility effects in the PRF method. If spin-echo methods are preferred, for example when field homogeneity is poor due to small ferromagnetic parts in the needle, the D- and T(1)-based methods may give better results. The sensitivity of the D method is higher that that of the T(1) methods provided that motion artifacts are avoided and the trace of D is evaluated. Fat suppression is necessary for most tissues when T(1), D, or PRF methods are employed. The latter three methods require excellent registration to correct for displacements between scans.

MR-based temperature monitoring for hot saline injection therapy

We applied magnetic resonance (MR) phase mapping methods to monitor the thermal frequency shift of water in order to study temperature changes from percutaneous hot saline injection therapy (PSIT) using in vitro swine livers and in vivo rabbit livers. The thermal coefficients calculated from the shifts of the water frequency with thermocouple based temperature measurements were -0.0085 +/- 0.0019 ppm/ degrees C for the in vitro studies and -0.0089 ppm/ degrees C for the in vivo studies. The error range was estimated to be +/- 3 degrees C and +/- 4.5 degrees C, respectively. Color-coded temperature maps were compared with macroscopic lesion sizes of the specimen. Regions defined using a 20 degrees C elevation in the initial images following hot saline injection (around 55 degrees C in absolute temperature) closely correlated with visible coagulation in size. We conclude that MR temperature monitoring of PSIT is quite feasible and may be helpful in expanding the clinical use of this thermal therapeutic tool for liver tumors.

Imaging temperature changes in an interventional 0.5 T magnet: in-vitro results

BACKGROUND AND OBJECTIVE

To evaluate the ability of monitoring laser induced temperature changes in an open, interventional 0.5 T magnet, adopting fast T1-weighted sequences.

MATERIALS AND METHODS

A fast gradient echo- (FGRE) and a fast spoiled gradient echo-sequence (FSPGR), both enabling an image update every 2.5 s, were investigated for their ability to visualize laser tissue effects at 5 Watt. Laser induced temperature was fluorooptically measured and correlated with signal intensity (SI) changes depicted by magnetic resonance imaging (MRI). MRI-lesions were compared with macroscopic findings.

RESULTS

SI changes on FGRE images appeared as early as 15 s following the onset of laser application and were significantly more pronounced than those seen on FSPGR images (p < .0001). A correlation of r = 0.94 (FGRE) and r = 0.92 (FSPGR) between temperature and SI loss was established. Owing to a steeper slope, the FGRE sequence was considered more sensitive to temperature changes. The areas of macroscopic tissue change correlated with those of SI loss, but lesion size was generally underestimated by MRI.

CONCLUSION

Laser monitoring is possible with rapid image updates in a midfield (0.5 T) interventional MRI environment using fast gradient echo sequence designs.

MR imaging of RF heating using a paramagnetic doped agarose phantom

In this paper, we present the first description of a technique to visualize and quantitate radiofrequency (RF) heating of a tissue phantom during a magnetic resonance imaging (MRI) procedure. We evaluated the heating patterns of four 10 cm diameter transmit/receive surface coils with differing degrees of distributed capacitance. The tissue phantom was a 6% agarose gel doped with 40 mM Na(4)HTm[DOTP], and possesses a conductivity intermediate to human muscle and fat. Heating was discerned via phase difference mapping using the large temperature dependent chemical shift coefficient for 23Na in Na(4)HTm[DOTP]. This coefficient is -0.5 ppm/ degrees C. Heating was highest where the phantom was closest to the surface coils, dropping off towards the center of the coil. No significant difference was observed in the heating patterns between the different surface coils. For the experimental setups used in this study, electric field 'hot spots' at the areas corresponding to the placement of the capacitor gaps were not observed.

Temperature monitoring in fat with MRI

The aim of the study was to test the hypothesis that fast spin echo T(1)-weighted images can be used to quantify the temperature in fat during thermal therapy in vivo. An MR compatible positioning device was used to manipulate focused ultrasound transducers in an MRI scanner. This system was used to sonicate fat tissue around the kidneys of 12 rabbits at various power levels for 10 to 20 sec. The scan parameters of T(1)-weighted fast spin echo (FSE) sequence were varied to optimize signal intensity characteristics while maintaining short scan times. An invasive optical probe was used to calibrate the temperature related signal intensity changes. For the T(1)-weighted FSE sequence, the signal intensity decreased with the temperature elevation at the rate of 0.97+/-0.02%/ degrees C. The single focused transducer produced a contrast-to-noise ratio more than 10 at power levels below the tissue damage threshold. The signal intensity was linearly dependent on the power, despite the measured temperatures being well above the coagulation threshold. This study demonstrates that T(1)-weighted FSE MRI sequences can be used to quantify the temperature elevation in fat in vivo during short focused ultrasound exposures. This can be very important for breast tumor surgery, fat ablation, and for treating deep seated tumors through superficial fat layers. Magn Reson Med 43:901-904, 2000.

Perturbations in hyperthermia temperature distributions associated with counter-current flow: numerical simulations and empirical verification

Two numerical techniques are used to calculate the effect of large vessel counter-current flow on hyperthermic temperature distributions. One is based on the Navier-Stokes equation for steady-state flow, and the second employs a convective-type boundary condition at the interface of the vessel walls. Steady-state temperature fields were calculated for two energy absorption rate distributions (ARD) in a cylindrical tissue model having two pairs of counter-current vessels (one pair with equal diameter vessels and another pair with unequal diameters). The first assumed a uniform ARD throughout cylinder; the second ARD was calculated for a tissue cylinder inside an existing four antenna Radiofrequency (RF) array. A tissue equivalent phantom was constructed to verify the numerical calculations. Temperatures induced with the RF array were measured using a noninvasive magnetic resonance imaging technique based on the chemical shift of water. Temperatures calculated using the two numerical techniques are in good agreement with the measured data. The results show: 1) the convective-type boundary condition technique reduces computation time by a factor of ten when compared to the fully conjugated method with little quantitative difference (approximately 0.3 degree C) in the numerical accuracy and 2) the use of noninvasive magnetic resonance imaging (thermal imaging) to quantitatively access the temperature perturbations near large vessels is feasible using the chemical shift technique.