Reconstruction of fully three-dimensional high spatial and temporal resolution MR temperature maps for retrospective applications

Advances in screening hyperthermic nanomedicines in 3D tumor models

Hyperthermic nanomedicines are particularly relevant for tackling human cancer, providing a valuable alternative to conventional therapeutics. The early-stage preclinical performance evaluation of such anti-cancer treatments is conventionally performed in flat 2D cell cultures that do not mimic the volumetric heat transfer occurring in human tumors. Recently, improvements in bioengineered 3D in vitro models have unlocked the opportunity to recapitulate major tumor microenvironment hallmarks and generate highly informative readouts that can contribute to accelerating the discovery and validation of efficient hyperthermic treatments. Leveraging on this, herein we aim to showcase the potential of engineered physiomimetic 3D tumor models for evaluating the preclinical efficacy of hyperthermic nanomedicines, featuring the main advantages and design considerations under diverse testing scenarios. The most recent applications of 3D tumor models for screening photo- and/or magnetic nanomedicines will be discussed, either as standalone systems or in combinatorial approaches with other anti-cancer therapeutics. We envision that breakthroughs toward developing multi-functional 3D platforms for hyperthermia onset and follow-up will contribute to a more expedited discovery of top-performing hyperthermic therapies in a preclinical setting before their in vivo screening.

Comparisons of healthy human brain temperature predicted from biophysical modeling and measured with whole brain MR thermometry

Brain temperature is an understudied parameter relevant to brain injury and ischemia. To advance our understanding of thermal dynamics in the human brain, combined with the challenges of routine experimental measurements, a biophysical modeling framework was developed to facilitate individualized brain temperature predictions. Model-predicted brain temperatures using our fully conserved model were compared with whole brain chemical shift thermometry acquired in 30 healthy human subjects (15 male and 15 female, age range 18-36 years old). Magnetic resonance (MR) thermometry, as well as structural imaging, angiography, and venography, were acquired prospectively on a Siemens Prisma whole body 3 T MR scanner. Bland-Altman plots demonstrate agreement between model-predicted and MR-measured brain temperatures at the voxel-level. Regional variations were similar between predicted and measured temperatures (< 0.55 °C for all 10 cortical and 12 subcortical regions of interest), and subcortical white matter temperatures were higher than cortical regions. We anticipate the advancement of brain temperature as a marker of health and injury will be facilitated by a well-validated computational model which can enable predictions when experiments are not feasible.

Progress in Understanding Radiofrequency Heating and Burn Injuries for Safer MR Imaging

RF electromagnetic wave exposure during MRI scans induces heat and occasionally causes burn injuries to patients. Among all the types of physical injuries that have occurred during MRI examinations, RF burn injuries are the most common ones. The number of RF burn injuries increases as the static magnetic field of MRI systems increases because higher RFs lead to higher heating. The commonly believed mechanisms of RF burn injuries are the formation of a conductive loop by the patient's posture or cables, such as an electrocardiogram lead; however, the mechanisms of RF burn injuries that occur at the contact points, such as the bore wall and the elbow, remain unclear. A comprehensive understanding of RF heating is needed to address effective countermeasures against all RF burn injuries for safe MRI examinations. In this review, we summarize the occurrence of RF burn injury cases by categorizing RF burn injuries reported worldwide in recent decades. Safety standards and regulations governing RF heating that occurs during MRI examinations are presented, along with their theoretical and physiological backgrounds. The experimental assessment techniques for RF heating are then reviewed, and the development of numerical simulation techniques is explained. In addition, a comprehensive theoretical interpretation of RF burn injuries is presented. By including the results of recent experimental and numerical simulation studies on RF heating, this review describes the progress achieved in understanding RF heating from the standpoint of MRI burn injury prevention.

A Breast-Specific MR Guided Focused Ultrasound Platform and Treatment Protocol: First-in-Human Technical Evaluation

OBJECTIVE

This paper presents and evaluates a breast-specific magnetic resonance guided focused ultrasound (MRgFUS) system. A first-in-human evaluation demonstrates the novel hardware, a sophisticated tumor targeting algorithm and a volumetric magnetic resonance imaging (MRI) protocol.

METHODS

At the time of submission, N = 10 patients with non-palpable T0 stage breast cancer have been treated with the breast MRgFUS system. The described tumor targeting algorithm is evaluated both with a phantom test and in vivo during the breast MRgFUS treatments. Treatments were planned and monitored using volumetric MR-acoustic radiation force imaging (MR-ARFI) and temperature imaging (MRTI).

RESULTS

Successful technical treatments were achieved in 80 % of the patients. All patients underwent the treatment with no sedation and 60 % of participants had analgesic support. The total MR treatment time ranged from 73 to 114 minutes. Mean error between desired and achieved targeting in a phantom was 2.9 ±1.8 mm while 6.2 ±1.9 mm was achieved in patient studies, assessed either with MRTI or MR-ARFI measurements. MRTI and MR-ARFI were successful in 60 % and 70 % of patients, respectively.

CONCLUSION

The targeting accuracy allows the accurate placement of the focal spot using electronic steering capabilities of the transducer. The use of both volumetric MRTI and MR-ARFI provides complementary treatment planning and monitoring information during the treatment, allowing the treatment of all breast anatomies, including homogeneously fatty breasts.

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.

A variable flip angle golden-angle-ordered 3D stack-of-radial MRI technique for simultaneous proton resonant frequency shift and T1 -based thermometry

PURPOSE

To develop and evaluate a variable-flip-angle golden-angle-ordered 3D stack-of-radial MRI technique for simultaneous proton resonance frequency shift (PRF) and T₁ -based thermometry in aqueous and adipose tissues, respectively.

METHODS

The proposed technique acquires multiecho radial k-space data in segments with alternating flip angles to measure 3D temperature maps dynamically on the basis of PRF and T₁ . A sliding-window k-space weighted image contrast filter is used to increase temporal resolution. PRF is measured in aqueous tissues and T₁ in adipose tissues using fat/water masks. The accuracy for T₁ quantification was evaluated in a reference T₁ /T₂ phantom. In vivo nonheating experiments were conducted in healthy subjects to evaluate the stability of PRF and T₁ in the brain, prostate, and breast. The proposed technique was used to monitor high-intensity focused ultrasound (HIFU) ablation in ex vivo porcine fat/muscle tissues and compared to temperature probe readings.

RESULTS

The proposed technique achieved 3D coverage with 1.1-mm to 1.3-mm in-plane resolution and 2-s to 5-s temporal resolution. During 20 to 30 min of nonheating in vivo scans, the temporal coefficient of variation for T₁ was <5% in the brain, prostate, and breast fatty tissues, while the standard deviation of relative PRF temperature change was within 3°C in aqueous tissues. During ex vivo HIFU ablation, the temperatures measured by PRF and T₁ were consistent with temperature probe readings, with an absolute mean difference within 2°C.

CONCLUSION

The proposed technique achieves simultaneous PRF and T₁ -based dynamic 3D MR temperature mapping in aqueous and adipose tissues. It may be used to improve MRI-guided thermal procedures.

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.

Simultaneous proton resonance frequency shift thermometry and T1 measurements using a single reference variable flip angle T1 method

PURPOSE

Implement simultaneous proton resonance frequency (PRF) shift and T₁ measurements with equivalent temporal resolution using a single reference variable flip angle method. This novel method allows for simultaneous thermometry in both aqueous and fatty tissue.

METHODS

This method acquires a single reference image at the lower flip angle and all dynamic images at the higher angle. T₁ is calculated using a single reference variable flip angle method, which accounts for the reference image temperature remaining constant. Monte Carlo simulations determined the optimal dynamic flip angle for combined PRF and T₁ measurements. This method was evaluated in MR-guided focused ultrasound heating experiments using a gelatin phantom and human cadaver breasts. In vivo measurement precision was demonstrated in healthy female volunteers under nonheating conditions.

RESULTS

Temperature rise during MR-guided focused ultrasound heating was measured in aqueous tissue with both PRF and T₁ . Both measures show good qualitative agreement in both space and time in aqueous tissue. The T₁ change due to temperature increase was measured in fat, demonstrating the expected temporal response. The dynamic flip angle that produces optimal SNR for PRF measurements is lower than the optimal angle for T₁ measurements, necessitating the selection of a compromise angle.

CONCLUSION

The single reference variable flip angle method provides a reliable way to simultaneously measure PRF temperature and T₁ change and overcomes PRF's inability to simultaneously monitor temperature in aqueous and adipose tissues. Future work will calibrate T₁ change to temperature, enabling real-time temperature in fat and increasing patient safety and treatment efficacy during thermal interventional treatments.

Predicting lesion size by accumulated thermal dose in MR-guided focused ultrasound for essential tremor

PURPOSE

To correlate the accumulated thermal dose (ATD) with lesion size in magnetic resonance (MR)-guided focused ultrasound (MRgFUS) thalamotomy to help guide future clinical treatments.

MATERIALS AND METHODS

Thirty-six patients with medication-refractory essential tremor were treated using a commercial MRgFUS brain system (ExAblate 4000, InSightec) in a 3T MR scanner (MR750, GE Healthcare). Intraoperative MR-thermometry was performed to measure the induced temperature and thermal dose distributions (thermal coefficient = -0.00909 ppm/°C). The ATD was calculated over multiple sonications with appropriate corrections for spatial-shifting artifacts. The ATD profile sizes obtained for dose values of 17, 40, 100, 200, and 240 cumulative equivalent minutes at 43°C (CEM) were correlated with the corresponding lesion sizes measured via axial T1- and T2-weighted MR images acquired 1 day post-treatment.

RESULTS

Of a total of 232 included sonications, 83 required corrections for off-resonance-induced spatial-shifting artifacts (correction range = [1.1,2.2] mm). The mean lesion sizes measured on T2-weighted MR images (6.2 ± 1.3 mm, mean ± SD) were 15% larger than those measured on corresponding T1-weighted MR images (5.3 ± 1.2 mm, mean ± SD). The ATD values that provided the best correlations with the measured lesion sizes on T2- and T1-weighted MR images were 100 and 200 CEM, respectively.

CONCLUSION

The ATD was correlated with lesion size measured 1 day following MRgFUS thalamotomy for essential tremor. These data provide useful information for predicting brain lesion size and determining treatment endpoints in future clinical MRgFUS procedures.

Effect of k-space-weighted image contrast and ultrasound focus size on the accuracy of proton resonance frequency thermometry

PURPOSE

To construct a predictive model that describes how the duration and symmetry of a k-space-weighted image contrast (KWIC) window affects the temporal resolution of differently sized ultrasound foci when using a pseudo-golden angle stack-of-stars acquisition.

METHODS

We performed a modulation analysis of proton resonance frequency temperature measurements to create the temporal modulation transfer function for KWIC windows of different symmetry and temporal duration. We reconstructed simulated ultrasound heating trajectories and stack-of-stars k-space data as well as experimental phantom data using the same trajectories. Images were reconstructed using symmetric and asymmetric KWIC windows of 3 different temporal durations. Simulated results were compared against temporal modulation transfer function predictions, experimental results, and the original simulated temperatures.

RESULTS

The temporal modulation transfer function shows that temporal resolution with KWIC reconstructions depend on the object size. The KWIC window duration affected SNR and severity of undersampling artifacts. Accuracy and response delay improved as the KWIC window duration decreased or the size of the heated region within the KWIC plane increased. Precision worsened as the window duration decreased. Using a symmetric window eliminated the response delay to heated region size but introduced a large reconstruction delay.

CONCLUSION

The accuracy and precision of proton resonance frequency temperature measurements from a stack-of-stars acquisition using a sliding KWIC window reconstruction are dependent on the size of the KWIC window and the size and shape of the heated region. The temporal modulation transfer function of KWIC reconstructions for any object size can predict the temporal response to changes in signal being acquired, such as temperature and contrast enhancement.

A spiral-based volumetric acquisition for MR temperature imaging

PURPOSE

To develop a rapid pulse sequence for volumetric MR thermometry.

METHODS

Simulations were carried out to assess temperature deviation, focal spot distortion/blurring, and focal spot shift across a range of readout durations and maximum temperatures for Cartesian, spiral-out, and retraced spiral-in/out (RIO) trajectories. The RIO trajectory was applied for stack-of-spirals 3D imaging on a real-time imaging platform and preliminary evaluation was carried out compared to a standard 2D sequence in vivo using a swine brain model, comparing maximum and mean temperatures measured between the two methods, as well as the temporal standard deviation measured by the two methods.

RESULTS

In simulations, low-bandwidth Cartesian trajectories showed substantial shift of the focal spot, whereas both spiral trajectories showed no shift while maintaining focal spot geometry. In vivo, the 3D sequence achieved real-time 4D monitoring of thermometry, with an update time of 2.9-3.3 s.

CONCLUSION

Spiral imaging, and RIO imaging in particular, is an effective way to speed up volumetric MR thermometry. Magn Reson Med 79:3122-3127, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Volumetric MRI thermometry using a three-dimensional stack-of-stars echo-planar imaging pulse sequence

PURPOSE

To measure temperature over a large brain volume with fine spatiotemporal resolution.

METHODS

A three-dimensional stack-of-stars echo-planar imaging sequence combining echo-planar imaging and radial sampling with golden angle spacing was implemented at 3T for proton resonance frequency-shift temperature imaging. The sequence acquires a 188x188x43 image matrix with 1.5x1.5x2.75 mm3 spatial resolution. Temperature maps were reconstructed using sensitivity encoding (SENSE) image reconstruction followed by the image domain hybrid method, and using the k-space hybrid method. In vivo temperature maps were acquired without heating to measure temperature precision in the brain, and in a phantom during high-intensity focused ultrasound sonication.

RESULTS

In vivo temperature standard deviation was less than 1°C at dynamic scan times down to 0.75 s. For a given frame rate, scanning at a minimum repetition time (TR) with minimum acceleration yielded the lowest standard deviation. With frame rates around 3 s, the scan was tolerant to a small number of receive coils, and temperature standard deviation was 48% higher than a standard two-dimensional Fourier transform temperature mapping scan, but provided whole-brain coverage. Phantom temperature maps with no visible aliasing were produced for dynamic scan times as short as 0.38 s. k-Space hybrid reconstructions were more tolerant to acceleration.

CONCLUSION

Three-dimensional stack-of-stars echo-planar imaging temperature mapping provides volumetric brain coverage and fine spatiotemporal resolution. Magn Reson Med 79:2003-2013, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Simultaneous MR thermometry and acoustic radiation force imaging using interleaved acquisition

PURPOSE

A novel and practical method for simultaneously performing MR acoustic radiation force imaging (ARFI) and proton resonance frequency (PRF)-shift thermometry has been developed and tested. This could be an important tool for evaluating the success of MR-guided focused ultrasound procedures for which MR-thermometry measures temperature and thermal dose and MR-ARFI detects changes in tissue mechanical properties.

METHODS

MR imaging was performed using a gradient recalled echo segmented echo-planar imaging pulse sequence with bipolar motion encoding gradients (MEG). Images with ultrasound pulses (ON) and without ultrasound pulses (OFF) during the MEG were interleaved at the repetition time (TR) level. ARFI displacements were calculated by complex subtraction of ON-OFF images, and PRF temperature maps were calculated by baseline subtraction. Evaluations in tissue-mimicking phantoms and ex vivo porcine brain tissue were performed. Constrained reconstruction improved the temporal resolution of dynamic measurements.

RESULTS

Simultaneous maps of displacement and temperature were acquired in 2D and 3D while keeping tissue heating < 1°C. Accuracy of the temperature maps was comparable to the standard PRF sequence. Using constrained reconstruction and subsampled k-space (R = 4.33), 3D simultaneous temperature and displacement maps can be acquired every 4.7 s.

CONCLUSION

This new sequence acquires simultaneous temperature and displacement maps with minimal tissue heating, and can be applied dynamically for monitoring tissue mechanical properties during ablation procedures. Magn Reson Med 79:1515-1524, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

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.

Accelerated magnetic resonance thermometry in the presence of uncertainties

A model-based information theoretic approach is presented to perform the task of magnetic resonance (MR) thermal image reconstruction from a limited number of observed samples on k-space. The key idea of the proposed approach is to optimally detect samples of k-space that are information-rich with respect to a model of the thermal data acquisition. These highly informative k-space samples can then be used to refine the mathematical model and efficiently reconstruct the image. The information theoretic reconstruction was demonstrated retrospectively in data acquired during MR-guided laser induced thermal therapy (MRgLITT) procedures. The approach demonstrates that locations with high-information content with respect to a model-based reconstruction of MR thermometry may be quantitatively identified. These information-rich k-space locations are demonstrated to be useful as a guide for k-space undersampling techniques. The effect of interactively increasing the predicted number of data points used in the subsampled model-based reconstruction was quantified using the L2-norm of the distance between the subsampled and fully sampled reconstruction. Performance of the proposed approach was also compared with uniform rectilinear subsampling and variable-density Poisson disk subsampling techniques. The proposed subsampling scheme resulted in accurate reconstructions using a small fraction of k-space points, suggesting that the reconstruction technique may be useful in improving the efficiency of thermometry data temporal resolution.

Specialized volumetric thermometry for improved guidance of MRgFUS in brain

PURPOSE

MR thermometry is critical for safe and effective transcranial focused ultrasound. The current single-slice MR thermometry sequence cannot achieve all desired treatment monitoring requirements. We propose an approach in which the imaging requirements of different aspects of treatment monitoring are met by optimizing multiple sequences.

METHODS

Imaging requirements were determined for three stages of MR-guided focused ultrasound brain treatment: 1) focal spot localization, 2) focal spot monitoring, and 3) background monitoring. Multiple-echo spiral thermometry sequences were optimized for each set of requirements and then validated with in vivo signal-to-noise ratio measurements and with phantom heating experiments.

RESULTS

Each of the proposed sequences obtained better precision than the current two-dimensional Fourier transform (2DFT) thermometry sequence. Five-slice focal spot localization achieved two-fold better resolution with 1.9-fold better precision but two-fold longer acquisition compared to 2DFT. Five-slice focal monitoring achieved 2.1-fold better precision with similar speed but 12% larger voxels than 2DFT. Full-brain background monitoring was demonstrated in both axial (7.1 s) and sagittal (11.4 s) orientations. Phantom heating time curves were consistent across all sequences after correcting for resolution.

CONCLUSION

Multiple-echo spiral imaging significantly improves MR thermometry efficiency, enabling multiple-slice monitoring. Optimizing multiple specialized sequences provides better performance than can be achieved by any single sequence. Magn Reson Med 78:508-517, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

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.

Prior image based temporally constrained reconstruction algorithm for magnetic resonance guided high intensity focused ultrasound

PURPOSE

A prior image based temporally constrained reconstruction (PITCR) algorithm was developed for obtaining accurate temperature maps having better volume coverage, and spatial, and temporal resolution than other algorithms for highly undersampled data in magnetic resonance (MR) thermometry.

METHODS

The proposed PITCR approach is an algorithm that gives weight to the prior image and performs accurate reconstruction in a dynamic imaging environment. The PITCR method is compared with the temporally constrained reconstruction (TCR) algorithm using pork muscle data.

RESULTS

The PITCR method provides superior performance compared to the TCR approach with highly undersampled data. The proposed approach is computationally expensive compared to the TCR approach, but this could be overcome by the advantage of reconstructing with fewer measurements. In the case of reconstruction of temperature maps from 16% of fully sampled data, the PITCR approach was 1.57× slower compared to the TCR approach, while the root mean square error using PITCR is 0.784 compared to 2.815 with the TCR scheme.

CONCLUSIONS

The PITCR approach is able to perform more accurate reconstructions of temperature maps compared to the TCR approach with highly undersampled data in MR guided high intensity focused ultrasound.

MRI-guided focused ultrasound surgery in musculoskeletal diseases: the hot topics

MRI-guided focused ultrasound surgery (MRgFUS) is a minimally invasive treatment guided by the most sophisticated imaging tool available in today's clinical practice. Both the imaging and therapeutic sides of the equipment are based on non-ionizing energy. This technique is a very promising option as potential treatment for several pathologies, including musculoskeletal (MSK) disorders. Apart from clinical applications, MRgFUS technology is the result of long, heavy and cumulative efforts exploring the effects of ultrasound on biological tissues and function, the generation of focused ultrasound and treatment monitoring by MRI. The aim of this article is to give an updated overview on a "new" interventional technique and on its applications for MSK and allied sciences.

Correcting heat-induced chemical shift distortions in proton resonance frequency-shift thermometry

PURPOSE

To reconstruct proton resonance frequency-shift temperature maps free of chemical shift distortions.

THEORY AND METHODS

Tissue heating created by thermal therapies such as focused ultrasound surgery results in a change in proton resonance frequency that causes geometric distortions in the image and calculated temperature maps, in the same manner as other chemical shift and off-resonance distortions if left uncorrected. We propose an online-compatible algorithm to correct these distortions in 2DFT and echo-planar imaging acquisitions, which is based on a k-space signal model that accounts for proton resonance frequency change-induced phase shifts both up to and during the readout. The method was evaluated with simulations, gel phantoms, and in vivo temperature maps from brain, soft tissue tumor, and uterine fibroid focused ultrasound surgery treatments.

RESULTS

Without chemical shift correction, peak temperature and thermal dose measurements were spatially offset by approximately 1 mm in vivo. Spatial shifts increased as readout bandwidth decreased, as shown by up to 4-fold greater temperature hot spot asymmetry in uncorrected temperature maps. In most cases, the computation times to correct maps at peak heat were less than 10 ms, without parallelization.

CONCLUSION

Heat-induced proton resonance frequency changes create chemical shift distortions in temperature maps resulting from MR-guided focused ultrasound surgery ablations, but the distortions can be corrected using an online-compatible algorithm. Magn Reson Med 76:172-182, 2016. © 2015 Wiley Periodicals, Inc.

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.

Complex difference constrained compressed sensing reconstruction for accelerated PRF thermometry with application to MRI-induced RF heating

PURPOSE

Introduce a novel compressed sensing reconstruction method to accelerate proton resonance frequency shift temperature imaging for MRI-induced radiofrequency heating evaluation.

METHODS

A compressed sensing approach that exploits sparsity of the complex difference between postheating and baseline images is proposed to accelerate proton resonance frequency temperature mapping. The method exploits the intra-image and inter-image correlations to promote sparsity and remove shared aliasing artifacts. Validations were performed on simulations and retrospectively undersampled data acquired in ex vivo and in vivo studies by comparing performance with previously published techniques.

RESULTS

The proposed complex difference constrained compressed sensing reconstruction method improved the reconstruction of smooth and local proton resonance frequency temperature change images compared to various available reconstruction methods in a simulation study, a retrospective study with heating of a human forearm in vivo, and a retrospective study with heating of a sample of beef ex vivo.

CONCLUSION

Complex difference based compressed sensing with utilization of a fully sampled baseline image improves the reconstruction accuracy for accelerated proton resonance frequency thermometry. It can be used to improve the volumetric coverage and temporal resolution in evaluation of radiofrequency heating due to MRI, and may help facilitate and validate temperature-based methods for safety assurance.

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.

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.

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.

Trans-cranial focused ultrasound without hair shaving: feasibility study in an ex vivo cadaver model

In preparing a patient for a trans-cranial magnetic resonance (MR)-guided focused ultrasound procedure, current practice is to shave the patient’s head on treatment day. Here we present an initial attempt to evaluate the feasibility of trans-cranial focused ultrasound in an unshaved, ex vivo human head model. A human skull filled with tissue-mimicking phantom and covered with a wig made of human hair was sonicated using 220- and 710-kHz head transducers to evaluate the feasibility of acoustic energy transfer. Heating at the focal point was measured by MR proton resonance shift thermometry. Results showed that the hair had a negligible effect on focal spot thermal rise at 220 kHz and a 17% drop in temperature elevation when using 710 kHz.

Toward real-time availability of 3D temperature maps created with temporally constrained reconstruction

PURPOSE

To extend the previously developed temporally constrained reconstruction (TCR) algorithm to allow for real-time availability of three-dimensional (3D) temperature maps capable of monitoring MR-guided high intensity focused ultrasound applications.

METHODS

A real-time TCR (RT-TCR) algorithm is developed that only uses current and previously acquired undersampled k-space data from a 3D segmented EPI pulse sequence, with the image reconstruction done in a graphics processing unit implementation to overcome computation burden. Simulated and experimental data sets of HIFU heating are used to evaluate the performance of the RT-TCR algorithm.

RESULTS

The simulation studies demonstrate that the RT-TCR algorithm has subsecond reconstruction time and can accurately measure HIFU-induced temperature rises of 20°C in 15 s for 3D volumes of 16 slices (RMSE = 0.1°C), 24 slices (RMSE = 0.2°C), and 32 slices (RMSE = 0.3°C). Experimental results in ex vivo porcine muscle demonstrate that the RT-TCR approach can reconstruct temperature maps with 192 × 162 × 66 mm 3D volume coverage, 1.5 × 1.5 × 3.0 mm resolution, and 1.2-s scan time with an accuracy of ±0.5°C.

CONCLUSION

The RT-TCR algorithm offers an approach to obtaining large coverage 3D temperature maps in real-time for monitoring MR-guided high intensity focused ultrasound treatments.

Extended Kalman filtering for continuous volumetric MR-temperature imaging

Real time magnetic resonance (MR) thermometry has evolved into the method of choice for the guidance of high-intensity focused ultrasound (HIFU) interventions. For this role, MR-thermometry should preferably have a high temporal and spatial resolution and allow observing the temperature over the entire targeted area and its vicinity with a high accuracy. In addition, the precision of real time MR-thermometry for therapy guidance is generally limited by the available signal-to-noise ratio (SNR) and the influence of physiological noise. MR-guided HIFU would benefit of the large coverage volumetric temperature maps, including characterization of volumetric heating trajectories as well as near- and far-field heating. In this paper, continuous volumetric MR-temperature monitoring was obtained as follows. The targeted area was continuously scanned during the heating process by a multi-slice sequence. Measured data and a priori knowledge of 3-D data derived from a forecast based on a physical model were combined using an extended Kalman filter (EKF). The proposed reconstruction improved the temperature measurement resolution and precision while maintaining guaranteed output accuracy. The method was evaluated experimentally ex vivo on a phantom, and in vivo on a porcine kidney, using HIFU heating. On the in vivo experiment, it allowed the reconstruction from a spatio-temporally under-sampled data set (with an update rate for each voxel of 1.143 s) to a 3-D dataset covering a field of view of 142.5×285×54 mm(3) with a voxel size of 3×3×6 mm(3) and a temporal resolution of 0.127 s. The method also provided noise reduction, while having a minimal impact on accuracy and latency.