Partially parallel imaging with phase-sensitive data: Increased temporal resolution for magnetic resonance temperature imaging

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.

Magnetic resonance thermometry and its biological applications - Physical principles and practical considerations

Most parameters that influence the magnetic resonance imaging (MRI) signal experience a temperature dependence. The fact that MRI can be used for non-invasive measurements of temperature and temperature change deep inside the human body has been known for over 30 years. Today, MR temperature imaging is widely used to monitor and evaluate thermal therapies such as radio frequency, microwave, laser, and focused ultrasound therapy. In this paper we cover the physical principles underlying the biological applications of MR temperature imaging and discuss practical considerations and remaining challenges. For biological tissue, the MR signal of interest comes mostly from hydrogen protons of water molecules but also from protons in, e.g., adipose tissue and various metabolites. Most of the discussed methods, such as those using the proton resonance frequency (PRF) shift, T₁, T₂, and diffusion only measure temperature change, but measurements of absolute temperatures are also possible using spectroscopic imaging methods (taking advantage of various metabolite signals as internal references) or various types of contrast agents. Currently, the PRF method is the most used clinically due to good sensitivity, excellent linearity with temperature, and because it is largely independent of tissue type. Because the PRF method does not work in adipose tissues, T₁- and T₂-based methods have recently gained interest for monitoring temperature change in areas with high fat content such as the breast and abdomen. Absolute temperature measurement methods using spectroscopic imaging and contrast agents often offer too low spatial and temporal resolution for accurate monitoring of ablative thermal procedures, but have shown great promise in monitoring the slower and usually less spatially localized temperature change observed during hyperthermia procedures. Much of the current research effort for ablative procedures is aimed at providing faster measurements, larger field-of-view coverage, simultaneous monitoring in aqueous and adipose tissues, and more motion-insensitive acquisitions for better precision measurements in organs such as the heart, liver, and kidneys. For hyperthermia applications, larger coverage, motion insensitivity, and simultaneous aqueous and adipose monitoring are also important, but great effort is also aimed at solving the problem of long-term field drift which gets interpreted as temperature change when using the PRF method.

Fast MR thermometry using an echo-shifted sequence with simultaneous multi-slice imaging

PURPOSE

Real-time monitoring is important for the safety and effectiveness of high-intensity focused ultrasound (HIFU) therapy. Magnetic resonance imaging is the preferred imaging modality for HIFU monitoring, with its unique capability of temperature imaging. For real-time temperature imaging, higher temporal resolution and larger spatial coverage are needed. In this study, a sequence based on the echo-shifted RF-spoiled gradient echo (GRE) with simultaneous multi-slice (SMS) imaging was designed for fast temperature imaging.

METHODS

A phantom experiment was conducted to evaluate the accuracy of the echo-shifted sequence using a fluorescent fiber thermometer as reference. The temperature uncertainty of the echo-shifted sequence was compared with the traditional GRE sequence at room temperature through the ex vivo porcine muscle. Finally, the ex vivo porcine liver tissue experiment using HIFU heating was performed to demonstrate that the spatial coverage was increased without decreasing temporal resolution.

RESULTS

The echo-shifted sequence had a better temperature uncertainty performance compared with the traditional GRE sequence with the same temporal resolution. The ex vivo heating experiment confirmed that by combining the SMS technique and echo-shifted sequence, the spatial coverage was increased without decreasing the temporal resolution while maintaining high temperature measurement precision.

CONCLUSION

The proposed technique was validated as an effective real-time method for monitoring HIFU therapy.

MR techniques for guiding high-intensity focused ultrasound (HIFU) treatments

To make full use of the ability of magnetic resonance (MR) to guide high-intensity focused ultrasound (HIFU) treatment, effort has been made to improve techniques for thermometry, motion tracking, and sound beam visualization. For monitoring rapid temperature elevation with proton resonance frequency (PRF) shift, data acquisition and processing can be accelerated with parallel imaging and/or sparse sampling in conjunction with appropriate signal processing methods. Thermometry should be robust against tissue motion, motion-induced magnetic field variation, and susceptibility change. Thus, multibaseline, referenceless, or hybrid techniques have become important. In cases with adipose or bony tissues, for which PRF shift cannot be used, thermometry with relaxation times or signal intensity may be utilized. Motion tracking is crucial not only for thermometry but also for targeting the focus of an ultrasound in moving organs such as the liver, kidney, or heart. Various techniques for motion tracking, such as those based on an anatomical image atlas with optical-flow displacement detection, a navigator echo to seize the diaphragm position, and/or rapid imaging to track vessel positions, have been proposed. Techniques for avoiding the ribcage and near-field heating have also been examined. MR acoustic radiation force imaging (MR-ARFI) is an alternative to thermometry that can identify the location and shape of the focal spot and sound beam path. This technique could be useful for treating heterogeneous tissue regions or performing transcranial therapy. All of these developments, which will be discussed further in this review, expand the applicability of HIFU treatments to a variety of clinical targets while maintaining safety and precision.

LEVEL OF EVIDENCE

2 Technical Efficacy: Stage 4 J. Magn. Reson. Imaging 2018;47:316-331.

MR thermometry for focused ultrasound monitoring utilizing model predictive filtering and ultrasound beam modeling

Background

A major challenge in using magnetic resonance temperature imaging (MRTI) to monitor focused ultrasound (FUS) applications is achieving high spatio-temporal resolution over a large field of view (FOV). This is important to accurately monitor all ultrasound (US) power depositions. Magnetic resonance (MR) subsampling in conjunction with thermal model-based reconstruction of the MRTI utilizing Pennes bioheat transfer equation (PBTE) is one promising approach. The thermal properties used in the thermal model are often estimated from a pre-treatment, low-power sonication.

Methods

In this proof-of-concept study we investigate the use of US simulations computed using the hybrid angular spectrum (HAS) method to estimate the US power deposition density Q, thereby avoiding the pre-treatment sonication and any potential tissue damage. MRTI reconstructions are performed using a thermal model-based reconstruction method called model predictive filtering (MPF). Experiments are performed in a homogeneous gelatin phantom and in a gelatin phantom with embedded plastic skull. MPF reconstructions are compared to separate sonications imaged with fully sampled data over a smaller FOV. Temperature root-mean-square errors (RMSE) and focal spot positions and shapes are evaluated.

Results

HAS simulations accurately predict the location of the focal spot (to within 1 mm) in both phantoms. Accurate temperature maps (RMSE below 1 °C), where the location of the focal spot agrees well with fully sampled “truth” (to within 1 mm), are also achieved in both phantoms.

Conclusions

HAS simulations can be used to accurately predict the focal spot location in homogeneous media and when focusing through an aberrating plastic skull. The HAS simulated power deposition (Q) patterns can be used in the MPF thermal model-based reconstruction to obtain accurate temperature maps with high spatio-temporal resolution over large FOVs.

Model predictive filtering MR thermometry: Effects of model inaccuracies, k-space reduction factor, and temperature increase rate

PURPOSE

Evaluate effects of model parameter inaccuracies (thermal conductivity, k, and ultrasound power deposition density, Q), k-space reduction factor (R), and rate of temperature increase ( T˙) in a thermal model-based reconstruction for MR-thermometry during focused-ultrasound heating.

METHODS

Simulations and ex vivo experiments were performed to investigate the accuracy of the thermal model and the model predictive filtering (MPF) algorithm for varying R and T˙, and their sensitivity to errors in k and Q. Ex vivo data was acquired with a segmented EPI pulse sequence to achieve large field-of-view (192 × 162 × 96 mm) four-dimensional temperature maps with high spatiotemporal resolution (1.5 × 1.5 × 2.0 mm, 1.7 s).

RESULTS

In the simulations, 50% errors in k and Q resulted in maximum temperature root mean square errors (RMSE) of 6 °C for model only and 3 °C for MPF. Using recently developed methods, estimates of k and Q were accurate to within 3%. The RMSE between MPF and true temperature increased with R and T˙. In the ex vivo study the RMSE remained below 0.7 °C for R ranging from 4 to 12 and T˙ of 0.28-0.75 °C/s.

CONCLUSION

Errors in MPF temperatures occur due to errors in k and Q. These MPF temperature errors increase with increase in R and T˙, but are smaller than those obtained using the thermal model alone.

Accelerated MRI thermometry by direct estimation of temperature from undersampled k-space data

PURPOSE

Acceleration of magnetic resonance (MR) thermometry is desirable for several applications of MR-guided focused ultrasound, such as those requiring greater volume coverage, higher spatial resolution, or higher frame rates.

METHODS

We propose and validate a constrained reconstruction method that estimates focal temperature changes directly from k-space without spatial or temporal regularization. A model comprising fully-sampled baseline images is fit to undersampled k-space data, which removes aliased temperature maps from the solution space. Reconstructed temperature maps are compared to maps reconstructed using parallel imaging (iterative self-consistent parallel imaging reconstruction [SPIRiT]) and conventional hybrid thermometry, and temporally constrained reconstruction thermometry.

RESULTS

Temporal step response simulations demonstrate finer temporal resolution and lower error in 4×-undersampled radial k-space reconstructions compared to temporally constrained reconstruction. Simulations show that the k-space method can achieve higher accelerations with multiple receive coils. Phantom heating experiments further demonstrate the algorithm's advantage over reconstructions relying on parallel imaging alone to overcome undersampling artifacts. In vivo model error comparisons show the algorithm achieves low temperature error at higher acceleration factors (up to 32× with a radial trajectory) than compared reconstructions.

CONCLUSION

High acceleration factors can be achieved using the proposed temperature reconstruction algorithm, without sacrificing temporal resolution or accuracy.

Ultrasonic enhancement of drug penetration in solid tumors

Increasing the penetration of drugs within solid tumors can be accomplished through multiple ultrasound-mediated mechanisms. The application of ultrasound can directly change the structure or physiology of tissues or can induce changes in a drug or vehicle in order to enhance delivery and efficacy. With each ultrasonic pulse, a fraction of the energy in the propagating wave is absorbed by tissue and results in local heating. When ultrasound is applied to achieve mild hyperthermia, the thermal effects are associated with an increase in perfusion or the release of a drug from a temperature-sensitive vehicle. Higher ultrasound intensities locally ablate tissue and result in increased drug accumulation surrounding the ablated region of interest. Further, the mechanical displacement induced by the ultrasound pulse can result in the nucleation, growth and collapse of gas bubbles. As a result of such cavitation, the permeability of a vessel wall or cell membrane can be increased. Finally, the radiation pressure of the propagating pulse can translate particles or tissues. In this perspective, we will review recent progress in ultrasound-mediated tumor delivery and the opportunities for clinical translation.

Generalised polynomial chaos-based uncertainty quantification for planning MRgLITT procedures

PURPOSE

A generalised polynomial chaos (gPC) method is used to incorporate constitutive parameter uncertainties within the Pennes representation of bioheat transfer phenomena. The stochastic temperature predictions of the mathematical model are critically evaluated against MR thermometry data for planning MR-guided laser-induced thermal therapies (MRgLITT).

METHODS

The Pennes bioheat transfer model coupled with a diffusion theory approximation of laser tissue interaction was implemented as the underlying deterministic kernel. A probabilistic sensitivity study was used to identify parameters that provide the most variance in temperature output. Confidence intervals of the temperature predictions are compared to MR temperature imaging (MRTI) obtained during phantom and in vivo canine (n = 4) MRgLITT experiments. The gPC predictions were quantitatively compared to MRTI data using probabilistic linear and temporal profiles as well as 2-D 60 °C isotherms.

RESULTS

Optical parameters provided the highest variance in the model output (peak standard deviation: anisotropy 3.51 °C, absorption 2.94 °C, scattering 1.84 °C, conductivity 1.43 °C, and perfusion 0.94 °C). Further, within the statistical sense considered, a non-linear model of the temperature and damage-dependent perfusion, absorption, and scattering is captured within the confidence intervals of the linear gPC method. Multivariate stochastic model predictions using parameters with the dominant sensitivities show good agreement with experimental MRTI data.

CONCLUSIONS

Given parameter uncertainties and mathematical modelling approximations of the Pennes bioheat model, the statistical framework demonstrates conservative estimates of the therapeutic heating and has potential for use as a computational prediction tool for thermal therapy planning.

Magnetic resonance thermal imaging combined with SMASH navigators in the presence of motion

This study develops and tests an MR thermometry method combined with SMASH navigators in phantom experiments mimicking human liver motion with the purpose of detecting and correcting motion artifacts in thermal MR images. Experimental data were acquired on a 3T MRI scanner. Motion artifacts of mobile phantoms mimicking human liver motion were detected and corrected using the SMASH navigators and then MR temperature maps were obtained using a proton resonant frequency (PRF) shift method with complex image subtraction. Temperature acquired by MR thermal imaging was compared to that measured via thermocouples. MR thermal imaging combined with the SMASH navigator technique resulted in accurate temperature maps of the mobile phantoms compared to temperatures measured using the thermocouples. The differences between the obtained and measured temperatures varied from 8.2°C to 14.2°C and 2.2°C to 4.9°C without and with motion correction, respectively. Motion correction improved the temperature acquired by MR thermal imaging by >55%. The combination of the MR thermal imaging and SMASH navigator technique will enable monitoring and controlling heat distribution and temperature change in tissues during thermal therapies and will be a very important tool for cancer treatment in mobile organs.

PACS number: 87.57.‐s

Magnetic resonance thermometry at 7T for real-time monitoring and correction of ultrasound induced mild hyperthermia

While Magnetic Resonance Thermometry (MRT) has been extensively utilized for non-invasive temperature measurement, there is limited data on the use of high field (≥7T) scanners for this purpose. MR-guided Focused Ultrasound (MRgFUS) is a promising non-invasive method for localized hyperthermia and drug delivery. MRT based on the temperature sensitivity of the proton resonance frequency (PRF) has been implemented in both a tissue phantom and in vivo in a mouse Met-1 tumor model, using partial parallel imaging (PPI) to speed acquisition. An MRgFUS system capable of delivering a controlled 3D acoustic dose during real time MRT with proportional, integral, and derivative (PID) feedback control was developed and validated. Real-time MRT was validated in a tofu phantom with fluoroptic temperature measurements, and acoustic heating simulations were in good agreement with MR temperature maps. In an in vivo Met-1 mouse tumor, the real-time PID feedback control is capable of maintaining the desired temperature with high accuracy. We found that real time MR control of hyperthermia is feasible at high field, and k-space based PPI techniques may be implemented for increasing temporal resolution while maintaining temperature accuracy on the order of 1°C.

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.

Ultrafast 1D MR thermometry using phase or frequency mapping

OBJECT

To develop an ultrafast MRI-based temperature monitoring method for application during rapid ultrasound exposures in moving organs.

MATERIALS AND METHODS

A slice selective 90° - 180° pair of RF pulses was used to solicit an echo from a column, which was then sampled with a train of gradient echoes. In a gel phantom, phase changes of each echo were compared to standard gradient-echo thermometry, and temperature monitoring was tested during focused ultrasound sonications. Signal-to-noise ratio (SNR) performance was evaluated in vivo in a rabbit brain, and feasibility was tested in a human heart.

RESULTS

The correlation between each echo in the acquisition and MRI-based temperature measurements was good (R = 0.98 ± 0.03). A temperature sampling rate of 19 Hz was achieved at 3T in the gel phantom. It was possible to acquire the water frequency in the beating heart muscle with 5-Hz sampling rate during a breath hold.

CONCLUSION

Ultrafast thermometry via phase or frequency monitoring along single columns was demonstrated. With a temporal resolution around 50 ms, it may be possible to monitor focal heating produced by short ultrasound pulses.

Combining two-dimensional spatially selective RF excitation, parallel imaging, and UNFOLD for accelerated MR thermometry imaging

MR thermometry can be a very challenging application, as good resolution may be needed along spatial, temporal, and temperature axes. Given that the heated foci produced during thermal therapies are typically much smaller than the anatomy being imaged, much of the imaged field-of-view is not actually being heated and may not require temperature monitoring. In this work, many-fold improvements were obtained in terms of temporal resolution and/or 3D spatial coverage by sacrificing some of the in-plane spatial coverage. To do so, three fast-imaging approaches were jointly implemented with a spoiled gradient echo sequence: (1) two-dimensional spatially selective RF excitation, (2) unaliasing by Fourier encoding the overlaps using the temporal dimension (UNFOLD), and (3) parallel imaging. The sequence was tested during experiments with focused ultrasound heating in ex vivo tissue and a tissue-mimicking phantom. Temperature maps were estimated from phase-difference images based on the water proton resonance frequency shift. Results were compared to those obtained from a spoiled gradient echo sequence sequence, using a t-test. Temporal resolution was increased by 24-fold, with temperature uncertainty less than 1°C, while maintaining accurate temperature measurements (mean difference between measurements, as observed in gel = 0.1°C ± 0.6; R = 0.98; P > 0.05).

Theranostics with multifunctional magnetic gold nanoshells: photothermal therapy and t2* magnetic resonance imaging

OBJECTIVES

to investigate the multifunctional imaging and therapeutic capabilities of core-shell nanoparticles composed of a superparamagnetic iron oxide (SPIO) core and a gold shell (SPIO@AuNS).

MATERIALS AND METHODS

the magnetic/optical properties of SPIO@AuNS were examined both in an agar gel phantom and in vivo by evaluating contrast-enhanced magnetic resonance imaging (MRI) and by measuring near-infrared (NIR) light-induced temperature changes mediated by SPIO@AuNS. In addition, the biodistribution and pharmacokinetics of In-labeled SPIO@AuNS after intravenous injection in mice bearing A431 tumors were evaluated in the presence and absence of an external magnet.

RESULTS

: In agar phantoms containing SPIO@AuNS, a significant contrast enhancement in T2-weighted MRI was observed and a linear increase in temperature was observed with increasing concentration and laser output power when irradiated with NIR light centered at an 808 nm. In vivo, T2*-MRI delineated SPIO@AuNS and magnetic resonance temperature imaging of the same tumors revealed significant temperature elevations when intratumorally injected with SPIO@AuNS (1 × 10 particles/mouse) and irradiated with NIR light (65.70°C ± 0.69°C vs. 44.23°C ± 0.24°C for saline + laser). Biodistribution studies in mice intravenously injected with In-labeled-SPIO@AuNS(1 × 10 particles/mouse) had an approximately 2-fold increase in SPIO@AuNS delivered into tumors in the presence of an external magnet compared with tumors without the magnet.

CONCLUSIONS

owing to its ability to mediate efficient photothermal ablation of cancer cells under MRI guidance, as well as the ability to be directed to solid tumors with an external magnetic field gradient, multifunctional SPIO@AuNS is a promising theranostic nanoplatform.

Fast fat-suppressed reduced field-of-view temperature mapping using 2DRF excitation pulses

The purpose of this study is to develop a fast and accurate temperature mapping method capable of both fat suppression and reduced field-of-view (rFOV) imaging, using a two-dimensional spatially-selective RF (2DRF) pulse. Temperature measurement errors caused by fat signals were assessed, through simulations. An 11×1140μs echo-planar 2DRF pulse was developed and incorporated into a gradient-echo sequence. Temperature measurements were obtained during focused ultrasound (FUS) heating of a fat-water phantom. Experiments both with and without the use of a 2DRF pulse were performed at 3T, and the accuracy of the resulting temperature measurements were compared over a range of TE values. Significant inconsistencies in terms of measured temperature values were observed when using a regular slice-selective RF excitation pulse. In contrast, the proposed 2DRF excitation pulse suppressed fat signals by more than 90%, allowing good temperature consistency regardless of TE settings. Temporal resolution was also improved, from 12 frames per minute (fpm) with the regular pulse to 28 frames per minute with the rFOV excitation. This technique appears promising toward the MR monitoring of temperature in moving adipose organs, during thermal therapies.

Magnetic resonance temperature imaging validation of a bioheat transfer model for laser-induced thermal therapy

PURPOSE

Magnetic resonance-guided laser-induced thermal therapy (MRgLITT) is currently undergoing initial safety and feasibility clinical studies for the treatment of intracranial lesions in humans. As studies progress towards evaluation of treatment efficacy, predictive computational models may play an important role for prospective 3D treatment planning. The current work critically evaluates a computational model of laser induced bioheat transfer against retrospective multiplanar MR thermal imaging (MRTI) in a canine model of the MRgLITT procedure in the brain.

METHODS

A 3D finite element model of the bioheat transfer that couples Pennes equation to a diffusion theory approximation of light transport in tissue is used. The laser source is modelled conformal with the applicator geometry. Dirichlet boundary conditions are used to model the temperature of the actively cooled catheter. The MRgLITT procedure was performed on n = 4 canines using a 1-cm diffusing tip 15-W diode laser (980 nm). A weighted L₂norm is used as the metric of comparison between the spatiotemporal MR-derived temperature estimates and model prediction.

RESULTS

The normalised error history between the computational models and MRTI was within 1-4 standard deviations of MRTI noise. Active cooling models indicate that the applicator temperature has a strong effect on the maximum temperature reached, but does not significantly decrease the tissue temperature away from the active tip.

CONCLUSIONS

Results demonstrate the computational model of the bioheat transfer may provide a reasonable approximation of the laser-tissue interaction, which could be useful for treatment planning, but cannot readily replace MR temperature imaging in a complex environment such as the brain.

Quantitative comparison of thermal dose models in normal canine brain

PURPOSE

Minimally invasive thermal ablative therapies as alternatives to conventional surgical management of solid tumors and other pathologies is increasing owing to the potential benefits of performing these procedures in an outpatient setting with reduced complications and comorbidity. Magnetic resonance temperature imaging (MRTI) measurement allows existing thermal dose models to use the spatiotemporal temperature history to estimate the thermal damage to tissue. However, the various thermal dose models presented in the literature employ different parameters and thresholds, affecting the reliability of thermal dosimetry. In this study, the authors quantitatively compared three thermal dose models (Arrhenius rate process, CEM43, and threshold temperature) using the dice similarity coefficient (DSC).

METHODS

The DSC was used to compare the spatial overlap between the region of thermal damage as predicted by the models for in vivo normal canine brain during thermal therapy to the region of thermal damage as revealed by contrast-enhanced T1-weighted images acquired immediately after therapy (< 20 min). The outer edge of the hyperintense rim of the ablation region was used as the surrogate marker for the limits of thermal coagulation. The DSC was also used to investigate the impact of varying the thresholds on each models' ability to predict the zone of thermal necrosis.

RESULTS

At previously reported thresholds, the authors found that all three models showed good agreement (defined as DSC > 0.7) with post-treatment imaging. All three models examined across the range of commonly applied thresholds consistently showed highly accurate spatial overlap, low variability, and little dependence on temperature uncertainty. DSC values corresponding to cited thresholds were not significantly different from peak DSC values.

CONCLUSIONS

Thus, the authors conclude that the all three thermal dose models can be used as a reliable surrogate for postcontrast tissue damage verification imaging in rapid ablation procedures and can also be used to enhance the capability of MRTI to control thermal therapy in real time.

Autoregressive moving average modeling for spectral parameter estimation from a multigradient echo chemical shift acquisition

The authors investigated the performance of the iterative Steiglitz-McBride (SM) algorithm on an autoregressive moving average (ARMA) model of signals from a fast, sparsely sampled, multiecho, chemical shift imaging (CSI) acquisition using simulation, phantom, ex vivo, and in vivo experiments with a focus on its potential usage in magnetic resonance (MR)-guided interventions. The ARMA signal model facilitated a rapid calculation of the chemical shift, apparent spin-spin relaxation time (T2*), and complex amplitudes of a multipeak system from a limited number of echoes (< or equal 16). Numerical simulations of one- and two-peak systems were used to assess the accuracy and uncertainty in the calculated spectral parameters as a function of acquisition and tissue parameters. The measured uncertainties from simulation were compared to the theoretical Cramer-Rao lower bound (CRLB) for the acquisition. Measurements made in phantoms were used to validate the T2* estimates and to validate uncertainty estimates made from the CRLB. We demonstrated application to real-time MR-guided interventions ex vivo by using the technique to monitor a percutaneous ethanol injection into a bovine liver and in vivo to monitor a laser-induced thermal therapy treatment in a canine brain. Simulation results showed that the chemical shift and amplitude uncertainties reached their respective CRLB at a signal-to-noise ratio (SNR) > or =5 for echo train lengths (ETLs) > or =4 using a fixed echo spacing of 3.3 ms. T2* estimates from the signal model possessed higher uncertainties but reached the CRLB at larger SNRs and/or ETLs. Highly accurate estimates for the chemical shift (<0.01 ppm) and amplitude (<1.0%) were obtained with > or =4 echoes and for T2*(<1.0%) with > or =7 echoes. We conclude that, over a reasonable range of SNR, the SM algorithm is a robust estimator of spectral parameters from fast CSI acquisitions that acquire < or =16 echoes for one- and two-peak systems. Preliminary ex vivo and in vivo experiments corroborated the results from simulation experiments and further indicate the potential of this technique for MR-guided interventional procedures with high spatiotemporal resolution approximately 1.6 x 1.6 x 4 mm3 in < or =5 s.

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.

Dynamic chemical shift imaging for image-guided thermal therapy: analysis of feasibility and potential

A fast chemical shift imaging (CSI) technique based on a multiple gradient-recalled acquisition using a small number of echoes with intentional aliasing of the reference lipid peak is studied to determine its feasibility for temperature monitoring. Simulations were implemented to find parameters where the lipid and water peaks can be measured using a Fourier-based peak fitting approach as well as using an innovative autoregressive moving average technique. A phantom consisting of 50% mayonnaise/50% lemon juice was calibrated to temperature and compared to literature values. A porcine kidney was treated ex vivo with an external laser and imaged with the CSI technique with comparisons to temperature readings from a fluoroptic monitoring system and complex phase difference (CPD) calculations. To demonstrate the technique in vivo, a Balb/c mouse with a CT26 xenograft in the subcutaneous lower back was treated using gold-coated, silica-core nanoshells heated with an 808 nm interstitial laser. Compared to standard CPD techniques using a two-dimensional fast spoiled gradient recalled echo, this technique maintains spatiotemporal resolution, has high signal-to-noise ratio and accuracy over a wide range of T2* tissue values, can separate water and lipid signals, and additionally can use the lipid peak, when present, as an internal reference.

Magnetic resonance temperature imaging for focused ultrasound surgery: a review

Magnetic resonance temperature imaging (MRTI) is an enabling technology that has recently demonstrated the potential to bring the emerging minimally invasive image-guided thermal therapy procedures, such as radiofrequency, microwave, laser, ultrasound, and cryosurgery, into the clinical setting with a level of safety and efficacy not previously possible. By coupling the wealth of soft tissue contrast mechanisms available with magnetic resonance imaging with its intrinsic temperature sensitivity, magnetic resonance imaging is in a unique position to provide image-guided treatment planning and verification and quantitative or qualitative feedback during treatment delivery, heightening of the control the physician has over the method, and enhancement of the ability to deliver conformal treatments. The basic principles behind MRTI technology and its application to minimally invasive thermal therapy during ultrasound thermal therapy delivery are reviewed in this study.

Auto-calibrated dynamic parallel MRI with phase-sensitive data

A number of MRI applications rely on dynamic phase information embedded in the acquired images. Such applications often require multiple acquisitions, leading to possibly long scan time and low temporal resolution. Previously, SENSE method has been used for phase-sensitive data to shorten acquisition time. However, SENSE can be subject to artifacts due to inaccurate coil sensitivities and low SNR. In this paper, dynamic phase data are derived from self-calibrated parallel MRI and an optimal method is used to combine phase information from multiple receiver channels. Simulation results using 4-channel prostate imaging data show that it is possible to get a factor of 3 speedup and the new method is more accurate than the SENSE method in reconstructing the phase information, thus has potential to improve phase-sensitive MRI applications such as phase contrast velocity mapping, temperature mapping for thermal therapy, and Dixon water/fat imaging.