Faculty Bibliography

PURPOSE: Present a novel method for rapid prediction of temperature in vivo for a series of pulse sequences with differing levels and distributions of specific energy absorption rate (SAR). THEORY AND METHODS: After the temperature response to a brief period of heating is characterized, a rapid estimate of temperature during a series of periods at different heating levels is made using a linear heat equation and impulse-response (IR) concepts. Here the initial characterization and long-term prediction for a complete spine exam are made with the Pennes' bioheat equation where, at first, core body temperature is allowed to increase and local perfusion is not. Then corrections through time allowing variation in local perfusion are introduced. RESULTS: The fast IR-based method predicted maximum temperature increase within 1% of that with a full finite difference simulation, but required less than 3.5% of the computation time. Even higher accelerations are possible depending on the time step size chosen, with loss in temporal resolution. Correction for temperature-dependent perfusion requires negligible additional time and can be adjusted to be more or less conservative than the corresponding finite difference simulation. CONCLUSION: With appropriate methods, it is possible to rapidly predict temperature increase throughout the body for actual MR examinations. Magn Reson Med, 2015. (c) 2015 Wiley Periodicals, Inc.

Electromagnetic field simulations are increasingly used to assure RF safety of patients during MRI exams. In practice, however, tissue property distribution of the patient being imaged is not known, but may be represented with a pre-existing model. Repeatedly, agreement in transmit magnetic (B1 +) field distributions between two geometries has been used to suggest agreement in heating distributions. Here we examine relative effects of anatomical differences on B1 + distribution, Specific Absorption Rate (SAR) and temperature change (DeltaT). Numerical simulations were performed for a single surface coil positioned adjacent a homogeneous phantom and bovine phantom, each with slight geometric variations, and adjacent two different human body models. Experimental demonstration was performed on a bovine phantom using MR thermometry and B1 + mapping. Simulations and experiments demonstrate that B1 + distributions in different samples can be well correlated, while notable difference in maximum SAR and DeltaT occur. This work illustrates challenges associated with utilizing simulations or experiments for RF safety assurance purposes. Reliance on B1 + distributions alone for validation of simulations and/or experiments with a sample or subject for assurance of safety in another should be performed with caution.

In high field MRI, the spatial distribution of the radiofrequency magnetic ( B1) field is usually affected by the presence of the sample. For hardware design and to aid interpretation of experimental results, it is important both to anticipate and to accurately simulate the behavior of these fields. Fields generated by a radiofrequency surface coil were simulated using dyadic Green's functions, or experimentally measured over a range of frequencies inside an object whose electrical properties were varied to illustrate a variety of transmit [Formula: see text] and receive [Formula: see text] field patterns. In this work, we examine how changes in polarization of the field and interference of propagating waves in an object can affect the B1 spatial distribution. Results are explained conceptually using Maxwell's equations and intuitive illustrations. We demonstrate that the electrical conductivity alters the spatial distribution of distinct polarized components of the field, causing "twisted" transmit and receive field patterns, and asymmetries between [Formula: see text] and [Formula: see text]. Additionally, interference patterns due to wavelength effects are observed at high field in samples with high relative permittivity and near-zero conductivity, but are not present in lossy samples due to the attenuation of propagating EM fields. This work provides a conceptual framework for understanding B1 spatial distributions for surface coils and can provide guidance for RF engineers.

The efficiencies of the designs of the different coils used to produce the static, gradient, and radiofrequency magnetic fields, as well as assessing the various interactions between these fields and the subject being imaged, are fundamental to the performance and safety of MRI. Here we review the basic field equations for these types of interactions and how they can be used to simulate MRI performance and safety

PURPOSE: All methods presented to date to map both conductivity and permittivity rely on multiple acquisitions to compute quantitatively the magnitude of radiofrequency transmit fields, B1+. In this work, we propose a method to compute both conductivity and permittivity based solely on relative receive coil sensitivities ( B1-) that can be obtained in one single measurement without the need to neither explicitly perform transmit/receive phase separation nor make assumptions regarding those phases. THEORY AND METHODS: To demonstrate the validity and the noise sensitivity of our method we used electromagnetic finite differences simulations of a 16-channel transceiver array. To experimentally validate our methodology at 7 Tesla, multi compartment phantom data was acquired using a standard 32-channel receive coil system and two-dimensional (2D) and 3D gradient echo acquisition. The reconstructed electric properties were correlated to those measured using dielectric probes. RESULTS: The method was demonstrated both in simulations and in phantom data with correlations to both the modeled and bench measurements being close to identity. The noise properties were modeled and understood. CONCLUSION: The proposed methodology allows to quantitatively determine the electrical properties of a sample using any MR contrast, with the only constraint being the need to have 4 or more receive coils and high SNR. Magn Reson Med, 2014. (c) 2014 Wiley Periodicals, Inc.

PURPOSE: To predict signal-to-noise ratio (SNR) trends and absorbed energy in magnetic resonance imaging (MRI) of the brain up to 14T. MATERIALS AND METHODS: A human head in an eight-channel transmit/receive coil was simulated with Maxwell and Bloch equations to determine excitation homogeneity with radiofrequency (RF) shimming, image homogeneity, SNR, and absorbed energy in MRI from 1.5 to 14T considering realistic field distributions and relaxation properties. RESULTS: RF shimming alone achieved a standard deviation in excitation flip angle less than 10 degrees in mid-brain up to 14T, but produced a small region with low excitation on a lower slice. Current reconstruction methods may produce shading artifacts at 14T. SNR increases with a greater-than-linear rate for gradient recalled echo (GRE) sequences having short (2 msec) echo time (TE) and long relaxation time (TR) ( approximately 2.3-fold increase from 7T to 14T), but a less-than-linear rate if TE is 10 msec ( approximately 1.6-fold increase from 7T to 14T). Depending on the sequence, SNR per square root of imaging time may produce a less-than-linear increase with B0 . Whole-head absorbed energy shows a less-than-quadratic increase with B0 (1.7-fold increase from 7T to 14T). CONCLUSION: Numerical simulations indicate that with proper preparation and precautions, imaging of the human brain at up to 14T could be performed safely, with advantages in SNR. J. Magn. Reson. Imaging 2015;41:1432-1439. (c) 2014 Wiley Periodicals, Inc.

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. Magn Reson Med 73:1420-1431, 2015. (c) 2014 Wiley Periodicals, Inc.

Due to the strong dependence of tissue electrical properties on temperature, it is important to consider the potential effects of intense tissue heating on the RF electromagnetic fields during MRI, as can occur in MR-guided focused ultrasound surgery. In principle, changes of the RF electromagnetic fields could affect both efficacy of RF pulses, and the MRI-induced RF heating (SAR) pattern. In this study, the equilibrium temperature distribution in a whole-body model with 2 mm resolution before and during intense tissue heating up to 60 degrees C at the target region was calculated. Temperature-dependent electric properties of tissues were assigned to the model to establish a temperature-dependent electromagnetic whole-body model in a 3T MRI system. The results showed maximum changes in conductivity, permittivity, [Formula: see text] and SAR of about 25%, 6%, 2%, and 20%, respectively. Though the B1 field and SAR distributions are both temperature-dependent, the potential harm to patients due to higher SARs is expected to be minimal and the effects on the B1 field distribution should have minimal effect on images from basic MRI sequences.

Magnetic resonance imaging (MRI) with T2-weighting is routinely performed to assess intervertebral disc degeneration. Standard clinical evaluations of MR images are qualitative, however, and do not focus on region-specific alterations in the disc. Utilizing a rabbit needle puncture model, T2 mapping was performed on injured discs to develop a quantitative description of the degenerative process following puncture. To do so, an 18G needle was inserted into four discs per rabbit (L3/L4 to L6/L7) and T2 maps were generated pre- and 4 weeks post-injury. Individual T2 maps were normalized to a disc-specific coordinate system and then averaged for pre- and post-injury population composite T2 maps. We also developed a method to automatically segment the nucleus pulposus by fitting the NP region of the T2 maps with modified 2-D and 3-D Gaussian distribution functions. Puncture injury produced alterations in MR signal intensity in a region-specific manner mirroring human degeneration. Population average T2 maps provided a quantitative representation of the injury response, and identified deviations of individual degenerate discs from the pre-injury population. We found that the response to standardized injury was modest at lower lumbar levels, likely as a result of increased disc dimensions. These tools will be valuable for the quantitative characterization of disc degeneration in future clinical and pre-clinical studies. (c) 2014 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 33:140-148, 2015.