An approach to rapid calculation of temperature change in tissue using spatial filters to approximate effects of thermal conduction

Numerical simulation in magnetic resonance imaging radiofrequency dosimetry

Magnetic Resonance Imaging (MRI) employs a radiofrequency electromagnetic field to create pictures on a computer. The prospective biological consequences of exposure to radiofrequency electromagnetic fields (RF EMFs) have not yet been demonstrated, and there is not enough evidence on biological hazards to offer a definite response concerning possible RF health dangers. Therefore, it is crucial to research the health concerns in reaction to RF EMFs, considering the entire exposure in terms of patients receiving MRI. Monitoring increases in temperaturein-vivothroughout MRI scan is extremely invasive and has resulted in a rise in the utilization of computational methods to estimate distributions of temperatures. The purpose of this study is to estimate the absorbed power of the brain exposed to RF in patients undergoing brain MRI scan. A three-dimensional Penne's bio-heat equation was modified to computationally analyze the temperature distributions and potential thermal effects within the brain during MRI scans in the 0.3 T to 1.5 T range (12.77 MHz to 63.87 MHz). The instantaneous temperature distributions of thein-vivotissue in the brain temperatures measured at a time, t = 20.62 s is 0.2 °C and t = 30.92 s is 0.4 °C, while the highest temperatures recorded at 1.03 min and 2.06 min were 0.4 °C and 0.6 °C accordingly. From the temperature distributions of thein-vivotissue in the brain temperatures measured, there is heat build-up in patients who are exposed to electromagnetic frequency ranges, and, consequently, temperature increases within patients are difficult to prevent. The study has, however, indicated that lengthier imaging duration appears to be related to increasing body temperature.

Uncertainty Quantification in SAR Induced by Ultra-High-Field MRI RF Coil via High-Dimensional Model Representation

As magnetic field strength in Magnetic Resonance Imaging (MRI) technology increases, maintaining the specific absorption rate (SAR) within safe limits across human head tissues becomes challenging due to the formation of standing waves at a shortened wavelength. Compounding this challenge is the uncertainty in the dielectric properties of head tissues, which notably affects the SAR induced by the radiofrequency (RF) coils in an ultra-high-field (UHF) MRI system. To this end, this study introduces a computational framework to quantify the impacts of uncertainties in head tissues' dielectric properties on the induced SAR. The framework employs a surrogate model-assisted Monte Carlo (MC) technique, efficiently generating surrogate models of MRI observables (electric fields and SAR) and utilizing them to compute SAR statistics. Particularly, the framework leverages a high-dimensional model representation technique, which constructs the surrogate models of the MRI observables via univariate and bivariate component functions, approximated through generalized polynomial chaos expansions. The numerical results demonstrate the efficiency of the proposed technique, requiring significantly fewer deterministic simulations compared with traditional MC methods and other surrogate model-assisted MC techniques utilizing machine learning algorithms, all while maintaining high accuracy in SAR statistics. Specifically, the proposed framework constructs surrogate models of a local SAR with an average relative error of 0.28% using 289 simulations, outperforming the machine learning-based surrogate modeling techniques considered in this study. Furthermore, the SAR statistics obtained by the proposed framework reveal fluctuations of up to 30% in SAR values within specific head regions. These findings highlight the critical importance of considering dielectric property uncertainties to ensure MRI safety, particularly in 7 T MRI systems.

Improved 1 H body imaging at 10.5 T: Validation and VOP-enabled imaging in vivo with a 16-channel transceiver dipole array

PURPOSE

To increase the RF coil performance and RF management for body imaging at 10.5 T by validating and evaluating a high-density 16-channel transceiver array, implementing virtual observation points (VOPs), and demonstrating specific absorption rate (SAR) constrained imaging in vivo.

METHODS

The inaccuracy of the electromagnetic model of the array was quantified based on B1 + and SAR data. Inter-subject variability was estimated using a new approach based on the relative SAR deviation of different RF shims between human body models. The pTx performance of the 16-channel array was assessed in simulation by comparison to a previously demonstrated 10-channel array. In vivo imaging of the prostate was performed demonstrating SAR-constrained static RF shimming and acquisition modes optimized for refocused echoes (AMORE).

RESULTS

The model inaccuracy of 29% and the inter-subject variability of 85% resulted in a total safety factor of 1.91 for pelvis studies. For renal and cardiac imaging, inter-subject variabilities of 121% and 141% lead to total safety factors of 2.25 and 2.45, respectively. The shorter wavelength at 10.5 T supported the increased element density of the 16-channel array which in turn outperformed the 10-channel version for all investigated metrics. Peak 10 g local SAR reduction of more than 25% without a loss of image quality was achieved in vivo, allowing a theoretical improvement in measurement efficiency of up to 66%.

CONCLUSIONS

By validating and characterizing a 16-channel dipole transceiver array, this work demonstrates, for the first time, a VOP-enabled RF coil for human torso imaging enabling increased pTx performance at 10.5 T.

Simulation-based evaluation of SAR and flip angle homogeneity for five transmit head arrays at 14 T

INTRODUCTION

Various research sites are pursuing 14 T MRI systems. However, both local SAR and RF transmit field inhomogeneity will increase. The aim of this simulation study is to investigate the trade-offs between peak local SAR and flip angle uniformity for five transmit coil array designs at 14 T in comparison to 7 T.

METHODS

Investigated coil array designs are: 8 dipole antennas (8D), 16 dipole antennas (16D), 8 loop coils (8D), 16 loop coils (16L), 8 dipoles/8 loop coils (8D8L) and for reference 8 dipoles at 7 T. Both RF shimming and k_T-points were investigated by plotting L-curves of peak SAR levels vs flip angle homogeneity.

RESULTS

For RF shimming, the 16L array performs best. For k_T-points, superior flip angle homogeneity is achieved at the expense of more power deposition, and the dipole arrays outperform the loop coil arrays.

DISCUSSION AND CONCLUSION

For most arrays and regular imaging, the constraint on head SAR is reached before constraints on peak local SAR are violated. Furthermore, the different drive vectors in k_T-points alleviate strong peaks in local SAR. Flip angle inhomogeneity can be alleviated by k_T-points at the expense of larger power deposition. For k_T-points, the dipole arrays seem to outperform loop coil arrays.

High-permittivity pads to enhance SNR and transmit efficiency in MRI of the heart at 7T: a simulation study

OBJECTIVE

High-permittivity pads have shown promising results in enhancing SNR and transmit efficiency when used for MRI of the brain, but fewer studies have been conducted to examine the performance of high-permittivity pads in other parts of the patient. In this work, we evaluate the impact on SNR and transmit efficiency distributions when high-permittivity pads with different thickness are positioned near the chest of the patient in combination with a transmit/receive array coil.

METHODS

The performance of the pads is evaluated through numerical simulations, and both the SNR distribution and the transmit efficiency maps are compared with those obtained when the pads are not present and the distance between the coils and the patient is minimal. The average improvement of SNR and transmit efficiency in the heart is also evaluated for different values of the permittivity of the pads.

RESULTS

In the scenario examined, high-permittivity pads can increase SNR and transmit efficiency in the heart volume by as much as 16% and 65%, respectively.

Bench to bore ramifications of inter-subject head differences on RF shimming and specific absorption rates at 7T

PURPOSE

This study shows how inter-subject variation over a dataset of 72 head models results in specific absorption rate (SAR) and B₁⁺ field homogeneity differences using common shim scenarios.

METHODS

MR-CT datasets were used to segment 71 head models into 10 tissue compartments. These head models were affixed to the shoulders and neck of the virtual family Duke model and placed within an 8 channel transmit surface-loop array to simulate the electromagnetic fields of a 7T imaging experiment. Radio frequency (RF) shimming using the Gerchberg-Saxton algorithm and Circularly Polarized shim weights over the entire brain and select slices of each model was simulated. Various SAR metrics and B₁⁺ maps were calculated to demonstrate the contribution of head variation to transmit inhomogeneity and SAR variability.

RESULTS

With varying head geometries the loading for each transmit loop changes as evidenced by changes in S-parameters. The varying shim conditions and head geometries are shown to affect excitation uniformity, spatial distributions of local SAR, and SAR averaging over different pulse sequences. The Gerchberg-Saxton RF shimming algorithm outperforms circularly polarized shimming for all head models. Peak local SAR within the coil most often occurs nearest the coil on the periphery of the body. Shim conditions vary the spatial distribution of SAR.

CONCLUSION

The work gives further support to the need for fast and more subject specific SAR calculations to maintain safety. Local SAR_10g is shown to vary spatially given shim conditions, subject geometry and composition, and position within the coil.

SAR and temperature distributions in a database of realistic human models for 7 T cardiac imaging

PURPOSE

To investigate inter-subject variability of B1+ , SAR and temperature rise in a database of human models using a local transmit array for 7 T cardiac imaging.

METHODS

Dixon images were acquired of 14 subjects and segmented in dielectric models with an eight-channel local transmit array positioned around the torso for cardiac imaging. EM simulations were done to calculate SAR distributions. Based on the SAR distributions, temperature simulations were performed for exposure times of 6 min and 30 min. Peak local SAR and temperature rise levels were calculated for different RF shim settings. A statistical analysis of the resulting peak local SAR and temperature rise levels was performed to arrive at safe power limits.

RESULTS

For RF shim vectors with random phase and uniformly distributed power, a safe average power limit of 35.7 W was determined (first level controlled mode). When RF amplitude and phase shimming was performed on the heart, a safe average power limit of 35.0 W was found. According to Pennes' model, our numerical study suggests a very low probability of exceeding the absolute local temperature limit of 40 °C for a total exposure time of 6 min and a peak local SAR of 20 W/kg. For a 30 min exposure time at 20 W/kg, it was shown that the absolute temperature limit can be exceeded in the case where perfusion does not change with temperature.

CONCLUSION

Safe power constraints were found for 7 T cardiac imaging with an eight-channel local transmit array, while considering the inter-subject variability of B1+ , SAR and temperature rise.

Temperature-based MRI safety simulations with a limited number of tissues

PURPOSE

Demonstrate ability to produce reasonable simulations of temperature using numerical models of the human body with a limited number of tissues.

METHODS

For both a male and female human body model, numerical simulations were used to calculate temperature distributions in three different models of the same human body: the original model with 35 tissues for the male model and 76 tissues for the female model, a simplified model having only three tissues (muscle, fat, and lung), and a simplified model having six tissues (muscle, fat, lung, bone, brain, and skin).

RESULTS

Although a three-tissue model gave reasonable specific absorption rate estimates in comparison to an original with many more tissues, because of tissue-specific thermal and physiological properties that do not affect specific absorption rate, such as rate of perfusion by blood, the three-tissue model did not provide temperature distributions similar to those of the original model. Inclusion of a few additional tissues, as in the six-tissue model, produced results in much better agreement with those from the original model.

CONCLUSION

Reasonable estimates of temperature can be simulated with a limited number of tissues, although this number is higher than the number of tissues required to produce reasonable simulations of specific absorption rate. For exposures primarily in the head and thorax, six tissues may be adequate for reasonable estimates of temperature.

RF heating measurement using MR thermometry and field monitoring: Methodological considerations and first in vivo results

PURPOSE

A MR thermometry (MRT) method with field monitoring is proposed to improve the measurement of small temperature variations induced in brain MRI exams.

METHODS

MR thermometry experiments were performed at 7 Tesla with concurrent field monitoring and RF heating. Images were reconstructed with nominal k-space trajectories and with first-order spherical harmonics correction. Experiments were performed in vitro with deliberate field disturbances and on an anesthetized macaque in 2 different specific absorption rate regimes, that is, at 50% and 100% of the maximal specific absorption rate level allowed in the International Electrotechnical Commission normal mode of operation. Repeatability was assessed by running a second separate session on the same animal.

RESULTS

Inclusion of magnetic field fluctuations in the reconstruction improved temperature measurement accuracy in vitro down to 0.02°C. Measurement precision in vivo was on the order of 0.15°C in areas little affected by motion. In the same region, temperature increase reached 0.5 to 0.8°C after 20 min of heating at 100% specific absorption rates and followed a rough factor of 2 with the 50% specific absorption rate scans. A horizontal temperature plateau, as predicted by Pennes bioheat model with thermal constants from the literature and constant blood temperature assumption, was not observed.

CONCLUSION

Inclusion of field fluctuations in image reconstruction was beneficial for the measurement of small temperature rises encountered in standard brain exams. More work is needed to correct for motion-induced field disturbances to extract reliable temperature maps.

Individualized SAR calculations using computer vision-based MR segmentation and a fast electromagnetic solver

PURPOSE

We propose a fast, patient-specific workflow for on-line specific absorption rate (SAR) supervision. An individualized electromagnetic model is created while the subject is on the table, followed by rapid SAR estimates for that individual. Our goal is an improved correspondence between the patient and model, reducing reliance on general anatomical body models.

METHODS

A 3D fat-water 3T acquisition (~2 minutes) is automatically segmented using a computer vision algorithm (~1 minute) into what we found to be the most important electromagnetic tissue classes: air, bone, fat, and soft tissues. We then compute the individual's EM field exposure and global and local SAR matrices using a fast electromagnetic integral equation solver. We assess the approach in 10 volunteers and compare to the SAR seen in a standard generic body model (Duke).

RESULTS

The on-the-table workflow averaged 7'44″. Simulation of the simplified Duke models confirmed that only air, bone, fat, and soft tissue classes are needed to estimate global and local SAR with an error of 6.7% and 2.7%, respectively, compared to the full model. In contrast, our volunteers showed a 16.0% and 20.3% population variability in global and local SAR, respectively, which was mostly underestimated by the Duke model.

CONCLUSION

Timely construction and deployment of a patient-specific model is computationally feasible. The benefit of resolving the population heterogeneity compared favorably to the modest modeling error incurred. This suggests that individualized SAR estimates can improve electromagnetic safety in MRI and possibly reduce conservative safety margins that account for patient-model mismatch, especially in non-standard patients.

Introduction of the snake antenna array: Geometry optimization of a sinusoidal dipole antenna for 10.5T body imaging with lower peak SAR

Purpose

To improve imaging performance for body MRI with a local transmit array at 10.5T, the geometry of a dipole antenna was optimized to achieve lower peak specific absorption rate (SAR) levels and a more uniform transmit profile.

Methods

Electromagnetic simulations on a phantom were used to evaluate the SAR and B1+‐performance of different dipole antenna geometries. The best performing antenna (the snake antenna) was simulated on human models in a 12‐channel array configuration for safety assessment and for comparison to a previous antenna design. This 12‐channel array was constructed after which electromagnetic simulations were validated by B1+‐maps and temperature measurements. After obtaining approval by the Food and Drug Administration to scan with the snake antenna array, in vivo imaging was performed on 2 volunteers.

Results

Simulation results on a phantom indicate a lower SAR and a higher transmit efficiency for the snake antenna compared to the fractionated dipole array. Similar results are found on a human body model: when comparing the trade‐off between uniformity and peak SAR, the snake antenna performs better for all imaging targets. Simulations and measurements are in good agreement. Preliminary imaging result were acquired in 2 volunteers with the 12‐channel snake antenna array.

Conclusion

By optimizing the geometry of a dipole antenna, peak SAR levels were lowered while achieving a more uniform transmit field as demonstrated in simulations on a phantom and a human body model. The array was constructed, validated, and successfully used to image 2 individuals at 10.5T.

Fast computation of soft tissue thermal response under deformation based on fast explicit dynamics finite element algorithm for surgical simulation

BACKGROUND AND OBJECTIVES

During thermal heating surgical procedures such as electrosurgery, thermal ablative treatment and hyperthermia, soft tissue deformation due to surgical tool-tissue interaction and patient movement can affect the distribution of thermal energy induced. Soft tissue temperature must be obtained from the deformed tissue for precise delivery of thermal energy. However, the classical Pennes bio-heat transfer model can handle only the static non-moving state of tissue. In addition, in order to enable a surgeon to visualise the simulated results immediately, the solution procedure must be suitable for real-time thermal applications.

METHODS

This paper presents a formulation of bio-heat transfer under the effect of soft tissue deformation for fast or near real-time tissue temperature prediction, based on fast explicit dynamics finite element algorithm (FED-FEM) for transient heat transfer. The proposed thermal analysis under deformation is achieved by transformation of the unknown deformed tissue state to the known initial static state via a mapping function. The appropriateness and effectiveness of the proposed formulation are evaluated on a realistic virtual human liver model with blood vessels to demonstrate a clinically relevant scenario of thermal ablation of hepatic cancer.

RESULTS

For numerical accuracy, the proposed formulation can achieve a typical 10^-3 level of normalised relative error at nodes and between 10^-4 and 10^-5 level of total errors for the simulation, by comparing solutions against the commercial finite element analysis package. For computation time, the proposed formulation under tissue deformation with anisotropic temperature-dependent properties consumes 2.518 × 10^-4 ms for one element thermal loads computation, compared to 2.237 × 10^-4 ms for the formulation without deformation which is 0.89 times of the former. Comparisons with three other formulations for isotropic and temperature-independent properties are also presented.

CONCLUSIONS

Compared to conventional methods focusing on numerical accuracy, convergence and stability, the proposed formulation focuses on computational performance for fast tissue thermal analysis. Compared to the classical Pennes model that handles only the static state of tissue, the proposed formulation can achieve fast thermal analysis on deformed states of tissue and can be applied in addition to tissue deformable models for non-linear heating analysis at even large deformation of soft tissue, leading to great translational potential in dynamic tissue temperature analysis and thermal dosimetry computation for computer-integrated medical education and personalised treatment.

Toward whole-cortex enhancement with an ultrahigh dielectric constant helmet at 3T

PURPOSE

To present a 3T brain imaging study using a conformal prototype helmet constructed with an ultra-high dielectric constant (uHDC; ε_r ~ 1000) materials that can be inserted into standard receive head-coils.

METHODS

A helmet conformal to a standard human head constructed with uHDC materials was characterized through electromagnetic simulations and experimental work. The signal-to-noise ratio (SNR), transmit efficiency, and power deposition with the uHDC helmet inserted within a 20-channel head coil were measured in vivo and compared with a 64-channel head coil and the 20-channel coil without the helmet. Seven healthy volunteers were analyzed.

RESULTS

Simulation and in vivo experimental results showed that transmit efficiency was improved by nearly 3 times within localized regions for a quadrature excitation, with a measured global increase of 58.21 ± 6.54% over 7 volunteers. The use of a parallel transmit spokes pulse compensated for severe degradation of B 1 + homogeneity, at the expense of higher global and local specific absorption rate levels. A SNR histogram analysis with statistical testing demonstrated that the uHDC helmet enhanced a 20-channel head coil to the level of the 64-channel head coil, with the improvements mainly within the cortical brain regions.

CONCLUSION

A prototype uHDC helmet enhanced the SNR of a standard head coil to the level of a high density 64-channel coil, although transmit homogeneity was compromised. Further improvements in SNR may be achievable with optimization of this technology, and could be a low-cost approach for future radiofrequency engineering work in the brain at 3T.

Optimization of the order and spacing of sequences in an MRI exam to reduce the maximum temperature and thermal dose

PURPOSE

Evaluate the possibility to reduce specific energy absorption rate (SAR)-induced maximum temperature and thermal dose by rearranging the order and spacing of sequences without increasing duration of the MRI examination.

METHODS

Using numerical simulations based on an actual SAR-intensive MRI examination, optimizations to reduce either maximum temperature or thermal dose were performed. For each permutation of groups of sequences having the same patient table position, temperature and thermal dose were computed very rapidly using recently published methods. Disposition of sequences was further adjusted by optimizing the spacing between each sequence without exceeding the original exam duration.

RESULTS

The maximum simulated temperature in the original exam was 42.38°C, and the maximum thermal dose was 3.23 cumulative effective minutes at 43°C (CEM43). After optimization to reduce maximum temperature, it was 41.77°C, and after optimization to minimize the thermal dose, it was 1.42 CEM43.

CONCLUSION

It is possible to reduce maximum temperature and thermal dose in the exam by changing the arrangement and spacing of the sequences without increasing the duration of the exam (by increasing TR or adding delays) or compromising image quality (by reducing flip angles).

Machine learning RF shimming: Prediction by iteratively projected ridge regression

PURPOSE

To obviate online slice-by-slice RF shim optimization and reduce B1+ mapping requirements for patient-specific RF shimming in high-field magnetic resonance imaging.

THEORY AND METHODS

RF Shim Prediction by Iteratively Projected Ridge Regression (PIPRR) predicts patient-specific, SAR-efficient RF shims with a machine learning approach that merges learning with training shim design. To evaluate it, a set of B1+ maps was simulated for 100 human heads for a 24-element coil at 7T. Features were derived from tissue masks and the DC Fourier coefficients of the coils' B1+ maps in each slice, which were used for kernelized ridge regression prediction of SAR-efficient RF shim weights. Predicted shims were compared to directly designed shims, circularly polarized mode, and nearest-neighbor shims predicted using the same features.

RESULTS

PIPRR predictions had 87% and 13% lower B1+ coefficients of variation compared to circularly polarized mode and nearest-neighbor shims, respectively, and achieved homogeneity and SAR similar to that of directly designed shims. Predictions were calculated in 4.92 ms on average.

CONCLUSION

PIPRR predicted uniform, SAR-efficient RF shims, and could save a large amount of B1+ mapping and computation time in RF-shimmed ultra-high field magnetic resonance imaging.

Parallel transmission RF pulse design with strict temperature constraints

RF safety in parallel transmission (pTx) is generally ensured by imposing specific absorption rate (SAR) limits during pTx RF pulse design. There is increasing interest in using temperature to ensure safety in MRI. In this work, we present a local temperature correlation matrix formalism and apply it to impose strict constraints on maximum absolute temperature in pTx RF pulse design for head and hip regions. Electromagnetic field simulations were performed on the head and hip of virtual body models. Temperature correlation matrices were calculated for four different exposure durations ranging between 6 and 24 min using simulated fields and body-specific constants. Parallel transmission RF pulses were designed using either SAR or temperature constraints, and compared with each other and unconstrained RF pulse design in terms of excitation fidelity and safety. The use of temperature correlation matrices resulted in better excitation fidelity compared with the use of SAR in parallel transmission RF pulse design (for the 6 min exposure period, 8.8% versus 21.0% for the head and 28.0% versus 32.2% for the hip region). As RF exposure duration increases (from 6 min to 24 min), the benefit of using temperature correlation matrices on RF pulse design diminishes. However, the safety of the subject is always guaranteed (the maximum temperature was equal to 39°C). This trend was observed in both head and hip regions, where the perfusion rates are very different.

Toward imaging the body at 10.5 tesla

PURPOSE

To explore the potential of performing body imaging at 10.5 Tesla (T) compared with 7.0T through evaluating the transmit/receive performance of similarly configured dipole antenna arrays.

METHODS

Fractionated dipole antenna elements for 10.5T body imaging were designed and evaluated using numerical simulations. Transmit performance of antenna arrays inside the prostate, kidneys and heart were investigated and compared with those at 7.0T using both phase-only radiofrequency (RF) shimming and multi-spoke pulses. Signal-to-noise ratio (SNR) comparisons were also performed. A 10-channel antenna array was constructed to image the abdomen of a swine at 10.5T. Numerical methods were validated with phantom studies at both field strengths.

RESULTS

Similar power efficiencies were observed inside target organs with phase-only shimming, but RF nonuniformity was significantly higher at 10.5T. Spokes RF pulses allowed similar transmit performance with accompanying local specific absorption rate increases of 25-90% compared with 7.0T. Relative SNR gains inside the target anatomies were calculated to be >two-fold higher at 10.5T, and 2.2-fold SNR gain was measured in a phantom. Gradient echo and fast spin echo imaging demonstrated the feasibility of body imaging at 10.5T with the designed array.

CONCLUSION

While comparable power efficiencies can be achieved using dipole antenna arrays with static shimming at 10.5T; increasing RF nonuniformities underscore the need for efficient, robust, and safe parallel transmission methods. Magn Reson Med 77:434-443, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Array-compressed parallel transmit pulse design

PURPOSE

To design array-compressed parallel transmit radiofrequency (RF) pulses and compare them to pulses designed with existing transmit array compression strategies.

THEORY AND METHODS

Array-compressed parallel RF pulse design is proposed as the joint optimization of a matrix of complex-valued compression weights that relate a full-channel physical array to a reduced-channel virtual array, along with a set of RF pulses for the virtual array. In this way, the physics of the RF pulse application determine the coil combination weights. Array-compressed pulse design algorithms are described for four parallel transmit applications: accelerated two-dimensional spiral excitation, multislice RF shimming, small-tip-angle kT -points excitation, and slice-selective spokes refocusing. Array-compressed designs are compared in simulations and an experiment to pulses designed using four existing array compression strategies.

RESULTS

In all cases, array-compressed pulses achieved the lowest root-mean-square excitation error among the array compression approaches. Low errors were generally achieved without increasing root-mean-square RF amplitudes or maximum local 10-gram specific absorption rate. Leave-one-out multisubject shimming simulations demonstrated that array-compressed RF shimming can identify useful fixed coil combination weights that perform well across a population.

CONCLUSION

Array-compressed pulse design jointly identifies the transmit coil array compression weights and RF pulses that perform best for a specific parallel excitation application. Magn Reson Med 76:1158-1169, 2016. © 2015 Wiley Periodicals, Inc.

B1-based SAR reconstruction using contrast source inversion-electric properties tomography (CSI-EPT)

Specific absorption rate (SAR) assessment is essential for safety purposes during MR acquisition. Online SAR assessment is not trivial and requires, in addition, knowledge of the electric tissue properties and the electric fields in the human anatomy. In this study, the potential of the recently developed CSI-EPT method to reconstruct SAR distributions is investigated. This method is based on integral representations for the electromagnetic field and attempts to reconstruct the tissue parameters and the electric field strength based on \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B_{1}^{ + }$$\end{document}B1+ field data only. Full three-dimensional FDTD simulations using a female pelvis model are used to validate two-dimensional CSI reconstruction results in the central transverse plane of a 3T body coil. Numerical experiments demonstrate that the reconstructed SAR distributions are in good agreement with the SAR distributions as determined via 3D FDTD simulations and show that these distributions can be computed very efficiently in the central transverse plane of a body coil with the two-dimensional approach of CSI-EPT.

A 16-channel combined loop-dipole transceiver array for 7 Tesla body MRI

PURPOSE

To develop a 16-channel transceive body imaging array at 7.0 T with improved transmit, receive, and specific absorption rate (SAR) performance by combining both loop and dipole elements and using their respective and complementary near and far field characteristics.

METHODS

A 16-channel radiofrequency (RF) coil array consisting of eight loop-dipole blocks (16LD) was designed and constructed. Transmit and receive performance was quantitatively investigated in phantom and human model simulations, and experiments on five healthy volunteers inside the prostate. Comparisons were made with 16-channel microstrip line (16ML) and 10-channel fractionated dipole antenna (10DA) arrays. The 16LD was used to acquire anatomic and functional images of the prostate, kidneys, and heart.

RESULTS

The 16LD provided > 14% improvements in the signal-to-noise ratio (SNR), peak B1+, B1+ transmit, and SAR efficiencies over the 16ML and 10DA in simulations inside the prostate. Experimentally, the 16LD had > 20% higher SNR and B1+ transmit efficiency compared with other arrays, and achieved up to 51.8% higher peak B1+ compared with 10DA.

CONCLUSION

Combining loop and dipole elements provided a body imaging array with high channel count and density while limiting inter-element coupling. The 16LD improved both near and far-field performance compared with existing 7.0T body arrays and provided high-quality MRI of the prostate kidneys and heart. Magn Reson Med 77:884-894, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Direct control of the temperature rise in parallel transmission by means of temperature virtual observation points: Simulations at 10.5 Tesla

PURPOSE

A method using parallel transmission to mitigate B1+ inhomogeneity while explicitly constraining the temperature rise is reported and compared with a more traditional SAR-constrained pulse design.

METHODS

Finite difference time domain simulations are performed on a numerical human head model and for a 16-channel coil at 10.5 Tesla. Based on a set of presimulations, a virtual observation point compression model for the temperature rise is derived. This compact representation is then used in a nonlinear programming algorithm for pulse design under explicit temperature rise constraints.

RESULTS

In the example of a time-of-flight sequence, radiofrequency pulse performance in some cases is increased by a factor of two compared with SAR-constrained pulses, while temperature rise is directly and efficiently controlled. Pulse performance can be gained by relaxing the SAR constraints, but at the expense of a loss of direct control on temperature.

CONCLUSION

Given the importance of accurate safety control at ultrahigh field and the lack of direct correspondence between SAR and temperature, this work motivates the need for thorough thermal studies in normal in vivo conditions. The tools presented here will possibly contribute to safer and more efficient MR exams.

Magnetic resonance temperature imaging-based quantification of blood flow-related energy losses

This study presents a new approach for evaluating bioheat transfer equation (BHTE) models used in treatment planning, control and evaluation of all thermal therapies. First, 3D magnetic resonance temperature imaging (MRTI) data are used to quantify blood flow-related energy losses, including the effects of perfusion and convection. Second, this information is used to calculate parameters of a BHTE model: in this paper the widely used Pennes BHTE. As a self-consistency check, the BHTE parameters are utilized to predict the temperatures from which they were initially derived. The approach is evaluated with finite-difference simulations and implemented experimentally with focused ultrasound heating of an ex vivo porcine kidney perfused at 0, 20 and 40 ml/min (n = 4 each). The simulation results demonstrate accurate quantification of blood flow-related energy losses, except in regions of sharp blood flow discontinuities, where the transitions are spatially smoothed. The smoothed transitions propagate into estimates of the Pennes perfusion parameter but have limited effect on the accuracy of temperature predictions using these estimates. Longer acquisition time periods mitigate the effects of MRTI noise, but worsen the effect of flow discontinuities. For the no-flow kidney experiments the estimates of a uniform, constant Pennes perfusion parameter are approximately zero, and at 20 and 40 ml/min the average estimates increase with flow rate to 3.0 and 4.2 kg/m(3) /s, respectively. When Pennes perfusion parameter values are allowed to vary spatially, but remain temporally constant, BHTE temperature predictions are more accurate than when using spatially uniform, constant Pennes perfusion values, with reductions in RMSE values of up to 79%. Locations with large estimated perfusion values correspond to high flow regions of the kidney observed in T1 -weighted MR images. This novel, MRTI-based technique holds promise for improving understanding of thermal therapy biophysics and for evaluating biothermal models.

Predicting long-term temperature increase for time-dependent SAR levels with a single short-term temperature response

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.

Diagnosis of breast tumor using thermal tomography q-r curve

Metabolic heat, the product following the metabolism of cells, is closely related to the pathological information of living organisms, which means there are strong connections between the heat distribution and the pathological state of the living organism. The mathematical function δ is introduced in the classical Pennes bioheat transfer equation as a point heat source, and by simplifying the boundary condition, a bioheat transfer model is established. Based on the temperature distribution of the human body surface, the q−r curve of heat intensity q varying with depth r is acquired while combining the fitting method of the Lorentz curve. According to 34,977 clinical confirmed cases and the corresponding classified statistics, diagnostic criteria (for breast diseases) for judging diseases by the q−r curve are proposed. The P -value of our statistics is <0.05 , which means our classified statistics are reliable. Six typical clinical examinations are performed, and the diagnosis results are very consistent with those of B-ultrasonic images, molybdenum target x-ray, and pathological examination, which suggests that the method of diagnosing diseases with a q−r curve has very good prospects for application. It is radiation free and noninvasive to the human body.

Q-r curve of thermal tomography and its clinical application on breast tumor diagnosis

Heat is the product following the metabolism of cells, and the metabolism is closely related with the pathological information of living organism. So, there are strong ties between the heat distribution and the pathological state in living organism. In this paper, the mathematical function δ is introduced in the classical Pennes bio-heat transfer equation as the point heat source. By simplifying the boundary conditions, a novel bio-heat transfer model is established and solved in a spherical coordinate system. Based on the temperature distribution of human body surface, the information of heat source is mined layer by layer, and the corresponding q-r curve of heat intensity varying with depth is acquired combining the fitting method of Lorentz curve. According to a large number of clinical confirmed cases and statistics, the diagnostic criteria judging diseases by q-r curve are proposed. Five typical clinical practices are performed and four of the diagnosis results are very consistent with those of molybdenum target (MT) X-ray, B-ultrasonic images and pathological examination, one gives the result of early stage malignant tumor that MT X-ray and B-ultrasonic can't check out. It is a radiation-free green method with noninvasive diagnostic procedure and accurate diagnosis result.

Combination of a multimode antenna and TIAMO for traveling-wave imaging at 9.4 Tesla

PURPOSE

To investigate the performance of a multimode antenna combined with time-interleaved acquisition of modes (TIAMO) for improved (1)H image homogeneity as compared to conventional traveling-wave imaging in the human brain at 9.4 Tesla (T).

METHODS

An adjustable three-port antenna was built to stimulate the propagation of three basic waveguide modes within a 9.4 T scanner bore. For TIAMO, two time-interleaved acquisitions using different linear combinations of these modes were optimized to achieve a homogeneous rooted sum-of-squares combination of their B1+ patterns ( B1,RSS+). The antenna's transmit and receive performance, as well as local specific absorption rate, were analyzed using experiments and numerical simulations.

RESULTS

The optimized TIAMO B1,RSS+ combination was superior to radiofrequency shimming. Across the entire brain, it improved the homogeneity of the excitation field by a factor of two and its maximum-to-minimum ratio by almost a factor of five as compared to the circularly polarized mode. The two-fold increase in "virtual" receive channels enhanced the parallel imaging performance and enabled the use of higher acceleration factors.

CONCLUSION

Despite the limited number of channels, a remote three-port antenna combined with TIAMO represents an easily implementable setup to achieve void-free (1)H images from the entire brain at 9.4 T, which can be used for anatomical localization and B0 shimming.

Numerical evaluation of image homogeneity, signal-to-noise ratio, and specific absorption rate for human brain imaging at 1.5, 3, 7, 10.5, and 14T in an 8-channel transmit/receive array

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° 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) (∼2.3-fold increase from 7T to 14T), but a less-than-linear rate if TE is 10 msec (∼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.

Measurement of SAR-induced temperature increase in a phantom and in vivo with comparison to numerical simulation

PURPOSE

To compare numerically simulated and experimentally measured temperature increase due to specific energy absorption rate from radiofrequency fields.

METHODS

Temperature increase induced in both a phantom and in the human forearm when driving an adjacent circular surface coil was mapped using the proton resonance frequency shift technique of magnetic resonance thermography. The phantom and forearm were also modeled from magnetic resonance image data, and both specific energy absorption rate and temperature change as induced by the same coil were simulated numerically.

RESULTS

The simulated and measured temperature increase distributions were generally in good agreement for the phantom. The relative distributions for the human forearm were very similar, with the simulations giving maximum temperature increase about 25% higher than measured.

CONCLUSION

Although a number of parameters and uncertainties are involved, it should be possible to use numerical simulations to produce reasonably accurate and conservative estimates of temperature distribution to ensure safety in magnetic resonance imaging. R01 EB006563