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.

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

Many areas of MR-guided thermal therapy research would benefit from temperature maps with high spatial and temporal resolution that cover a large three-dimensional volume. This article describes an approach to achieve these goals, which is suitable for research applications where retrospective reconstruction of the temperature maps is acceptable. The method acquires undersampled data from a modified three-dimensional segmented echo-planar imaging sequence and creates images using a temporally constrained reconstruction algorithm. The three-dimensional images can be zero-filled to arbitrarily small voxel spacing in all directions and then converted into temperature maps using the standard proton resonance frequency shift technique. During high intensity focused ultrasound heating experiments, the proposed method was used to obtain temperature maps with 1.5 mm × 1.5 mm × 3.0 mm resolution, 288 mm × 162 mm × 78 mm field of view, and 1.7 s temporal resolution. The approach is validated to demonstrate that it can accurately capture the spatial characteristics and time dynamics of rapidly changing high intensity focused ultrasound-induced temperature distributions. Example applications from MR-guided high intensity focused ultrasound research are shown to demonstrate the benefits of the large coverage fully three-dimensional temperature maps, including characterization of volumetric heating trajectories and near- and far-field heating.

The effects of spatial sampling choices on MR temperature measurements.

The purpose of this article is to quantify the effects that spatial sampling parameters have on the accuracy of magnetic resonance temperature measurements during high intensity focused ultrasound treatments. Spatial resolution and position of the sampling grid were considered using experimental and simulated data for two different types of high intensity focused ultrasound heating trajectories (a single point and a 4-mm circle) with maximum measured temperature and thermal dose volume as the metrics. It is demonstrated that measurement accuracy is related to the curvature of the temperature distribution, where regions with larger spatial second derivatives require higher resolution. The location of the sampling grid relative temperature distribution has a significant effect on the measured values. When imaging at 1.0 × 1.0 × 3.0 mm(3) resolution, the measured values for maximum temperature and volume dosed to 240 cumulative equivalent minutes (CEM) or greater varied by 17% and 33%, respectively, for the single-point heating case, and by 5% and 18%, respectively, for the 4-mm circle heating case. Accurate measurement of the maximum temperature required imaging at 1.0 × 1.0 × 3.0 mm(3) resolution for the single-point heating case and 2.0 × 2.0 × 5.0 mm(3) resolution for the 4-mm circle heating case.

A simulation technique for 3D MR-guided acoustic radiation force imaging.

PURPOSE: In magnetic resonance-guided focused ultrasound (MRgFUS) therapies, the in situ characterization of the focal spot location and quality is critical. MR acoustic radiation force imaging (MR-ARFI) is a technique that measures the tissue displacement caused by the radiation force exerted by the ultrasound beam. This work presents a new technique to model the displacements caused by the radiation force of an ultrasound beam in a homogeneous tissue model. METHODS: When a steady-state point-source force acts internally in an infinite homogeneous medium, the displacement of the material in all directions is given by the Somigliana elastostatic tensor. The radiation force field, which is caused by absorption and reflection of the incident ultrasound intensity pattern, will be spatially distributed, and the tensor formulation takes the form of a convolution of a 3D Green's function with the force field. The dynamic accumulation of MR phase during the ultrasound pulse can be theoretically accounted for through a time-of-arrival weighting of the Green's function. This theoretical model was evaluated experimentally in gelatin phantoms of varied stiffness (125-, 175-, and 250-bloom). The acoustic and mechanical properties of the phantoms used as parameters of the model were measured using independent techniques. Displacements at focal depths of 30- and 45-mm in the phantoms were measured by a 3D spin echo MR-ARFI segmented-EPI sequence. RESULTS: The simulated displacements agreed with the MR-ARFI measured displacements for all bloom values and focal depths with a normalized RMS difference of 0.055 (range 0.028-0.12). The displacement magnitude decreased and the displacement pattern broadened with increased bloom value for both focal depths, as predicted by the theory. CONCLUSIONS: A new technique that models the displacements caused by the radiation force of an ultrasound beam in a homogeneous tissue model theory has been rigorously validated through comparison with experimentally obtained 3D displacement data in homogeneous gelatin phantoms using a 3D MR-ARFI sequence. The agreement of the experimentally measured and simulated results demonstrates the potential to use MR-ARFI displacement data in MRgFUS therapies.