Single acquisition electrical property mapping based on relative coil sensitivities: A proof-of-concept demonstration

Changes of in vivo electrical conductivity in the brain and torso related to age, fat fraction and sex using MRI

This work was inspired by the observation that a majority of MR-electrical properties tomography studies are based on direct comparisons with ex vivo measurements carried out on post-mortem samples in the 90's. As a result, the in vivo conductivity values obtained from MRI in the megahertz range in different types of tissues (brain, liver, tumors, muscles, etc.) found in the literature may not correspond to their ex vivo equivalent, which still serves as a reference for electromagnetic modelling. This study aims to pave the way for improving current databases since the definition of personalized electromagnetic models (e.g. for Specific Absorption Rate estimation) would benefit from better estimation. Seventeen healthy volunteers underwent MRI of both brain and thorax/abdomen using a three-dimensional ultrashort echo-time (UTE) sequence. We estimated conductivity (S/m) in several classes of macroscopic tissue using a customized reconstruction method from complex UTE images, and give general statistics for each of these regions (mean-median-standard deviation). These values are used to find possible correlations with biological parameters such as age, sex, body mass index and/or fat volume fraction, using linear regression analysis. In short, the collected in vivo values show significant deviations from the ex vivo values in conventional databases, and we show significant relationships with the latter parameters in certain organs for the first time, e.g. a decrease in brain conductivity with age.

Phantom evaluation of electrical conductivity mapping by MRI: Comparison to vector network analyzer measurements and spatial resolution assessment

PURPOSE

To evaluate the performance of various MR electrical properties tomography (MR-EPT) methods at 3 T in terms of absolute quantification and spatial resolution limit for electrical conductivity.

METHODS

Absolute quantification as well as spatial resolution performance were evaluated on homogeneous phantoms and a phantom with holes of different sizes, respectively. Ground-truth conductivities were measured with an open-ended coaxial probe connected to a vector network analyzer (VNA). Four widely used MR-EPT reconstruction methods were investigated: phase-based Helmholtz (PB), phase-based convection-reaction (PB-cr), image-based (IB), and generalized-image-based (GIB). These methods were compared using the same complex images from a 1 mm-isotropic UTE sequence. Alternative transceive phase acquisition sequences were also compared in PB and PB-cr.

RESULTS

In large homogeneous phantoms, all methods showed a strong correlation with ground truth conductivities (r > 0.99); however, GIB was the best in terms of accuracy, spatial uniformity, and robustness to boundary artifacts. In the resolution phantom, the normalized root-mean-squared error of all methods grew rapidly (>0.40) when the hole size was below 10 mm, with simplified methods (PB and IB), or below 5 mm, with generalized methods (PB-cr and GIB).

CONCLUSION

VNA measurements are essential to assess the accuracy of MR-EPT. In this study, all tested MR-EPT methods correlated strongly with the VNA measurements. The UTE sequence is recommended for MR-EPT, with the GIB method providing good accuracy for structures down to 5 mm. Structures below 5 mm may still be detected in the conductivity maps, but with significantly lower accuracy.

Calibrated Frequency-Division Distorted Born Iterative Tomography for Real-Life Head Imaging

The clinical use of microwave tomography (MT) requires addressing the significant mismatch between simulated environment, which is used in the forward solver, and real-life system. To alleviate this mismatch, a calibrated tomography, which uses two homogeneous calibration phantoms and a modified distorted Born iterative method (DBIM), is presented. The two phantoms are used to derive a linear model that matches the forward solver to real-life measurements. Moreover, experimental observations indicate that signal quality at different frequencies varies between different antennas due to inevitably inconsistent manufacturing tolerance and variances in radio-frequency chains. An optimum frequency, at which the simulated and measured signals of the antenna present maximum similarity when irradiating the calibrated phantoms, is thus calculated for each antenna. A frequency-division DBIM (FD-DBIM), in which different antennas in the array transmit their corresponding optimum frequencies, is subsequently developed. A clinical brain scanner is then used to assess performance of the algorithm in lab and healthy volunteers' tests. The linear calibration model is first used to calibrate the measured data. After that FD-DBIM is used to solve the problem and map the dielectric properties of the imaged domain. The simulated and experimental results confirm validity of the presented approach and its superiority to other tomographic method.

Conductivity Tensor Imaging of the Human Brain Using Water Mapping Techniques

Conductivity tensor imaging (CTI) has been recently proposed to map the conductivity tensor in 3D using magnetic resonance imaging (MRI) at the frequency range of the brain at rest, i.e., low-frequencies. Conventional CTI mapping methods process the trans-receiver phase of the MRI signal using the MR electric properties tomography (MR-EPT) technique, which in turn involves the application of the Laplace operator. This results in CTI maps with a low signal-to-noise ratio (SNR), artifacts at tissue boundaries and a limited spatial resolution. In order to improve on these aspects, a methodology independent from the MR-EPT method is proposed. This relies on the strong assumption for which electrical conductivity is univocally pre-determined by water concentration. In particular, CTI maps are calculated by combining high-frequency conductivity derived from water maps and multi b-value diffusion tensor imaging (DTI) data. Following the implementation of a pipeline to optimize the pre-processing of diffusion data and the fitting routine of a multi-compartment diffusivity model, reconstructed conductivity images were evaluated in terms of the achieved spatial resolution in five healthy subjects scanned at rest. We found that the pre-processing of diffusion data and the optimization of the fitting procedure improve the quality of conductivity maps. We achieve reproducible measurements across healthy participants and, in particular, we report conductivity values across subjects of 0.55 ± 0.01Sm, 0.3 ± 0.01Sm and 2.15 ± 0.02Sm for gray matter (GM), white matter (WM), and cerebrospinal fluid (CSF), respectively. By attaining an actual spatial resolution of the conductivity tensor close to 1 mm in-plane isotropic, partial volume effects are reduced leading to good discrimination of tissues with similar conductivity values, such as GM and WM. The application of the proposed framework may contribute to a better definition of the head tissue compartments in electroencephalograpy/magnetoencephalography (EEG/MEG) source imaging and be used as biomarker for assessing conductivity changes in pathological conditions, such as stroke and brain tumors.

Transceive phase corrected 2D contrast source inversion-electrical properties tomography

PURPOSE

To remove the necessity of the tranceive phase assumption for CSI-EPT and show electrical properties maps reconstructed from measured data obtained using a standard 3T birdcage body coil setup.

METHODS

The existing CSI-EPT algorithm is reformulated to use the transceive phase rather than relying on the transceive phase assumption. Furthermore, the radio frequency (RF)-shield is numerically implemented to accurately model the RF fields inside the MRI scanner. We verify that the reformulated two-dimensional (2D) CSI-EPT algorithm can reconstruct electrical properties maps given 2D electromagnetic simulations. Afterward, the algorithm is tested with three-dimensional (3D) FDTD simulations to investigate if the 2D CSI-EPT can retrieve the electrical properties for 3D RF fields. Finally, an MR experiment at 3T with a phantom is performed.

RESULTS

From the results of the 2D simulations, it is seen that CSI-EPT can reconstruct the electrical properties using MRI accessible quantities. For 3D simulations, it is observed that the electrical properties are underestimated, nonetheless, CSI-EPT has a lower standard deviation than the standard Helmholtz based methods. Finally, the first CSI-EPT reconstructions based on measured data are presented showing comparable accuracy and precision to reconstructions based on simulated data, and demonstrating the feasibility of CSI-EPT.

CONCLUSIONS

The CSI-EPT algorithm was rewritten to use MRI accessible quantities. This allows for CSI-EPT to fully exploit the benefits of the higher static magnetic field strengths with a standard quadrature birdcage coil setup.

Decomposition of high-frequency electrical conductivity into extracellular and intracellular compartments based on two-compartment model using low-to-high multi-b diffusion MRI

Background

As an object’s electrical passive property, the electrical conductivity is proportional to the mobility and concentration of charged carriers that reflect the brain micro-structures. The measured multi-b diffusion-weighted imaging (Mb-DWI) data by controlling the degree of applied diffusion weights can quantify the apparent mobility of water molecules within biological tissues. Without any external electrical stimulation, magnetic resonance electrical properties tomography (MREPT) techniques have successfully recovered the conductivity distribution at a Larmor-frequency.

Methods

This work provides a non-invasive method to decompose the high-frequency conductivity into the extracellular medium conductivity based on a two-compartment model using Mb-DWI. To separate the intra- and extracellular micro-structures from the recovered high-frequency conductivity, we include higher b-values DWI and apply the random decision forests to stably determine the micro-structural diffusion parameters.

Results

To demonstrate the proposed method, we conducted phantom and human experiments by comparing the results of reconstructed conductivity of extracellular medium and the conductivity in the intra-neurite and intra-cell body. The phantom and human experiments verify that the proposed method can recover the extracellular electrical properties from the high-frequency conductivity using a routine protocol sequence of MRI scan.

Conclusion

We have proposed a method to decompose the electrical properties in the extracellular, intra-neurite, and soma compartments from the high-frequency conductivity map, reconstructed by solving the electro-magnetic equation with measured B1 phase signals.

Electrical Properties Tomography: A Methodological Review

Electrical properties tomography (EPT) is an imaging method that uses a magnetic resonance (MR) system to non-invasively determine the spatial distribution of the conductivity and permittivity of the imaged object. This manuscript starts by providing clear definitions about the data required for, and acquired in, EPT, followed by comprehensively formulating the physical equations underlying a large number of analytical EPT techniques. This thorough mathematical overview of EPT harmonizes several EPT techniques in a single type of formulation and gives insight into how they act on the data and what their data requirements are. Furthermore, the review describes machine learning-based algorithms. Matlab code of several differential and iterative integral methods is available upon request.

MR electrical properties imaging using a generalized image-based method

PURPOSE

To develop a fast and easy-to-use electrical properties tomography (EPT) method based on a single MR scan, avoiding both the need of a B₁ -map and transceive phase assumption, and that is robust against noise.

THEORY

Derived from Maxwell's equations, conductivity, and permittivity are reconstructed from a new partial differential equation involving the product of the RF fields and its derivatives. This also allows us to clarify and revisit the relevance of common assumptions of MREPT.

METHODS

Our new governing equation is solved using a 3D finite-difference scheme and compared to previous frameworks. The benefits of our method over selected existing MREPT methods are demonstrated for different simulation models, as well as for both an inhomogeneous agar phantom gel and in vivo brain data at 3T.

RESULTS

Simulation and experimental results are illustrated to highlight the merits of the proposed method over existing methods. We show the validity of our algorithm in versatile configurations, with many transition regions notably. Complex admittivity maps are also provided as a complementary MR contrast.

CONCLUSION

Because it avoids time-consuming RF field mapping and generalizes the use of standard MR image for electrical properties reconstruction, this contribution is promising as a new step forward for clinical applications.

Magnetic-Resonance-Based Electrical Property Mapping Using Global Maxwell Tomography With an 8-Channel Head Coil at 7 Tesla: A Simulation Study

OBJECTIVE

Global Maxwell Tomography (GMT) is a recently introduced volumetric technique for noninvasive estimation of electrical properties (EP) from magnetic resonance measurements. Previous work evaluated GMT using ideal radiofrequency (RF) excitations. The aim of this simulation study was to assess GMT performance with a realistic RF coil.

METHODS

We designed a transmit-receive RF coil with 8 decoupled channels for 7T head imaging. We calculated the RF transmit field ( B₁⁺) inside heterogeneous head models for different RF shimming approaches, and used them as input for GMT to reconstruct EP for all voxels.

RESULTS

Coil tuning/decoupling remained relatively stable when the coil was loaded with different head models. Mean error in EP estimation changed from [Formula: see text] to [Formula: see text] and from [Formula: see text] to [Formula: see text] for relative permittivity and conductivity, respectively, when changing head model without re-tuning the coil. Results slightly improved when an SVD-based RF shimming algorithm was applied, in place of excitation with one coil at a time. Despite errors in EP, RF transmit field ( B₁⁺) and absorbed power could be predicted with less than [Formula: see text] error over the entire head. GMT could accurately detect a numerically inserted tumor.

CONCLUSION

This work demonstrates that GMT can reliably reconstruct EP in realistic simulated scenarios using a tailored 8-channel RF coil design at 7T. Future work will focus on construction of the coil and optimization of GMT's robustness to noise, to enable in-vivo GMT experiments.

SIGNIFICANCE

GMT could provide accurate estimations of tissue EP, which could be used as biomarkers and could enable patient-specific estimation of RF power deposition, which is an unsolved problem for ultra-high-field magnetic resonance imaging.

Integral MR-EPT With the Calculation of Coil Current Distributions

Many integral equation (IE)-based magnetic resonance electrical property tomography (MR-EPT) methods use unloaded incident radio-frequency (RF) fields from simulations that may not fully reflect the real situation and thus lead to reconstruction errors. To improve the accuracy of IE-based MR-EPT methods, a novel approach that enables the calculation of loaded coil current distributions and avoids the explicit use of incident RF fields is presented in this paper. In the proposed method, a hybrid source composed of the current source from the coil and the contrast source from the subject are introduced in the integral equations. Because the loaded coil current distributions can be extracted from the reconstructed hybrid source, the simulated incident RF fields are eliminated from the problem formulations. To improve the convergence performance, a modified conjugate gradient (CG) scheme was used where the gradients of the current source and contrast source were balanced through using different weighting parameters. The proposed method was verified through full-wave simulations at 9.4 and 7 T involving a heterogeneous ball and an anatomical head phantom. The numerical results indicated that by using the proposed method, an accurate coil current distributions and EPs profiles can be reconstructed and the desirable robustness against noise can also be achieved.

Opening a new window on MR-based Electrical Properties Tomography with deep learning

In the radiofrequency (RF) range, the electrical properties of tissues (EPs: conductivity and permittivity) are modulated by the ionic and water content, which change for pathological conditions. Information on tissues EPs can be used e.g. in oncology as a biomarker. The inability of MR-Electrical Properties Tomography techniques (MR-EPT) to accurately reconstruct tissue EPs by relating MR measurements of the transmit RF field to the EPs limits their clinical applicability. Instead of employing electromagnetic models posing strict requirements on the measured MRI quantities, we propose a data driven approach where the electrical properties reconstruction problem can be casted as a supervised deep learning task (DL-EPT). DL-EPT reconstructions for simulations and MR measurements at 3 Tesla on phantoms and human brains using a conditional generative adversarial network demonstrate high quality EPs reconstructions and greatly improved precision compared to conventional MR-EPT. The supervised learning approach leverages the strength of electromagnetic simulations, allowing circumvention of inaccessible MR electromagnetic quantities. Since DL-EPT is more noise-robust than MR-EPT, the requirements for MR acquisitions can be relaxed. This could be a major step forward to turn electrical properties tomography into a reliable biomarker where pathological conditions can be revealed and characterized by abnormalities in tissue electrical properties.

Noninvasive Estimation of Electrical Properties From Magnetic Resonance Measurements via Global Maxwell Tomography and Match Regularization

OBJECTIVE

In this paper, we introduce global Maxwell tomography (GMT), a novel volumetric technique that estimates electric conductivity and permittivity by solving an inverse scattering problem based on magnetic resonance measurements.

METHODS

GMT relies on a fast volume integral equation solver, MARIE, for the forward path, and a novel regularization method, match regularization, designed specifically for electrical property estimation from noisy measurements. We performed simulations with three different tissue-mimicking numerical phantoms of different complexity, using synthetic transmit sensitivity maps with realistic noise levels as the measurements. We performed an experiment at 7 T using an eight-channel coil and a uniform phantom.

RESULTS

We showed that GMT could estimate relative permittivity and conductivity from noisy magnetic resonance measurements with an average error as low as 0.3% and 0.2%, respectively, over the entire volume of the numerical phantom. Voxel resolution did not affect GMT performance and is currently limited only by the memory of the graphics processing unit. In the experiment, GMT could estimate electrical properties within 5% of the values measured with a dielectric probe.

CONCLUSION

This work demonstrated the feasibility of GMT with match regularization, suggesting that it could be effective for accurate in vivo electrical property estimation. GMT does not rely on any symmetry assumption for the electromagnetic field, and can be generalized to estimate also the spin magnetization, at the expense of increased computational complexity.

SIGNIFICANCE

GMT could provide insight into the distribution of electromagnetic fields inside the body, which represents one of the key ongoing challenges for various diagnostic and therapeutic applications.

CONtrast Conformed Electrical Properties Tomography (CONCEPT) Based on Multi- Channel Transmission and Alternating Direction Method of Multipliers

In magnetic resonance-based electrical properties tomography (EPT), circularly polarized magnetic field B1 from a transmit radiofrequency (RF) coil is measured and utilized to infer the electrical conductivity and permittivity of biological tissues. Compared with a quadrature RF coil, a multi-channel transmit coil provides a plurality of unique transmit B1 patterns that help to alleviate the under-determinedness of EPT reconstruction problem, and it also allows to circumvent the "transceive phase assumption" that fails at ultra-high-field MRI. Here, a new approach, contrast conformed electrical properties tomography or CONCEPT, is proposed based on the multi-channel transmission that retrieves electrical properties (EPs) by solving a linear partial differential equation with discriminated L1 and L2 norm regularization informed by intermediate EP gradient. The theory of CONCEPT and a fast reconstruction algorithm based on the alternating direction method of multipliers are described and evaluated using numerical simulations, phantom experiment, and analysis of in vivo human brain data at 7 T MRI. Compared with the multi-channel gradient-based EPT (gEPT) method, this new technology does not require receive- B1 sensitivity profiles and does not rely on symmetry assumption regarding RF coil design and imaged target. Moreover, it is not dependent on external prior information, such as integration seed point or anatomical MRI, which can be sources of bias in reconstructed EP values. By deriving EPs from transmit B1 profiles only, CONCEPT can be used with RF coils that include receive-only arrays with large channel count which can, in turn, offer substantial gains in signal-to-noise ratio. It also holds potentials to image unsymmetrical body organs and diseased brain. CONCEPT provides solutions for the practical problems during the implementation of gEPT, thus representing a more generalized framework in the context of multi-channel RF transmission.

Mapping electrical properties heterogeneity of tumor using boundary informed electrical properties tomography (BIEPT) at 7T

PURPOSES

To develop and evaluate a boundary informed electrical properties tomography (BIEPT) technique for high-resolution imaging of tumor electrical properties (EPs) heterogeneity on a rodent tumor xenograft model.

METHODS

Tumor EP distributions were inferred from a reference area external to the tumor, as well as internal EP spatial variations derived from a plurality of relative transmit B1 measurements at 7T. Edge sparsity constraint was enforced to enhance numerical stability. Phantom experiments were performed to determine the imaging accuracy and sensitivity for structures of various EP values, as well as geometrical sizes down to 1.5 mm. Numerical simulation of a realistic rodent model was used to quantify the algorithm performance in the presence of noise. Eleven athymic rats with human breast cancer xenograft were imaged in vivo, and representative pathological samples were acquired for comparison.

RESULTS

Reconstructed EPs of the phantoms correspond well to the ground truth acquired from dielectric probe measurements, with the smallest structure reliably detectable being 3 mm. EPs heterogeneity inside a tumor is successfully retrieved in both simulated and experimental cases. In vivo tumor imaging results demonstrate similar local features and spatial patterns to anatomical MRI and pathological slides. The imaged conductivity of necrotic tissue is higher than that of viable tissues, which agrees with our expectation.

CONCLUSION

BIEPT enables robust detection of tumor EPs heterogeneity with high accuracy and sensitivity to small structures. The retrieved quantitative EPs reflect tumor pathological features (e.g., necrosis). These results provide strong rationale to further expand BIEPT studies toward pathological conditions where EPs may yield valuable, non-invasive biomarkers.

Electrical properties tomography: Available contrast and reconstruction capabilities

MR-based electrical properties tomography converts the MRI transmit/receive RF field measurements to tissue electrical property maps through dedicated reconstruction algorithms. Recent reports showed that despite limitations, electrical properties tomography holds promise for generating additional contrast for tumor detection and patient-specific modeling of tissue-RF field interactions. This review summarizes the available tissue electrical property contrasts and compares them with the capabilities of the most commonly used electrical properties tomography reconstruction method. Future directions and prospects of clinical translation are discussed.

Dictionary-based electric properties tomography

PURPOSE

To develop and validate a new algorithm called "dictionary-based electric properties tomography" (dbEPT) for deriving tissue electric properties from measured B₁ maps.

METHODS

Inspired by Magnetic Resonance fingerprinting, dbEPT uses a dictionary of local patterns ("atoms") of B₁ maps and corresponding electric properties distributions, derived from electromagnetic field simulations. For reconstruction, a pattern from a measured B₁ map is compared with the B₁ atoms of the dictionary. The B₁ atom showing the best match with the measured B₁ pattern yields the optimum electric properties pattern that is chosen for reconstruction. Matching was performed through machine learning algorithms. Two dictionaries, using transmit and transceive phases, were evaluated. The spatial distribution of local matching distance between optimal atom and measured pattern yielded a reconstruction reliability map. The method was applied to reconstruct conductivity of 4 volunteers' brains. A conventional, Helmholtz-based Electric properties tomography (EPT) reconstruction was performed for reference. Noise performance was studied through phantom simulations.

RESULTS

Quantitative values of conductivity agree with literature values. Results of the 2 dictionaries exhibit only minor differences. Somewhat larger differences are visible between dbEPT and Helmholtz-based EPT. Quantified by the correlation between conductivity and anatomic images, dbEPT depicts brain details more clearly than Helmholtz-based EPT. Matching distance is minimal in homogeneous brain ventricles and increases with tissue heterogeneity. Central processing unit time was approximately 2 minutes per dictionary training and 3 minutes per brain conductivity reconstruction using standard hardware equipment.

CONCLUSION

A new, dictionary-based approach for reconstructing electric properties is presented. Its conductivity reconstruction is able to overcome the EPT transceive-phase problem.

3-D Contrast Source Inversion-Electrical Properties Tomography

Contrast source inversion-electrical properties tomography (CSI-EPT) is an iterative reconstruction method to retrieve the electrical properties (EPs) of tissues from magnetic resonance data. The method is based on integral representations of the electromagnetic field and has been shown to allow EP reconstructions of small structures as well as tissue boundaries with compelling accuracy. However, to date, the CSI-EPT has been implemented for 2-D configurations only, which limits its applicability. In this paper, a full 3-D extension of the CSI-EPT method is presented, to enable CSI-EPT to be applied to realistic 3-D scenarios. Here, we demonstrate a proof-of-principle of 3-D CSI-EPT and present the reconstructions of a 3-D abdominal body section and a 3-D head model using different settings of the transmit coil. Numerical results show that the full 3-D approach yields accurate reconstructions of the EPs, even at tissue boundaries and is most accurate in regions where the absolute value of the electric field is highest.

Error analysis of helmholtz-based MR-electrical properties tomography

PURPOSE

MR electrical properties tomography (MR-EPT) aims to measure tissue electrical properties by computing spatial derivatives of measured B1+ data. This computation is very sensitive to spatial fluctuations caused, for example, by noise and Gibbs ringing. In this work, the error arising from the computation of spatial derivatives using finite difference kernels (FD error) has been investigated. In relation to this FD error, it has also been investigated whether mitigation strategies such as Gibbs ringing correction and Gaussian apodization can be beneficial for conductivity reconstructions.

METHODS

Conductivity reconstructions were performed on a phantom (by means of simulations and MR measurements at 3T) and on a human brain model. The accuracy was evaluated as a function of image resolution, FD kernel size, k-space windowing, and signal-to-noise ratio. The impact of mitigation strategies was also investigated.

RESULTS

The adopted small FD kernel is highly sensitive to spatial fluctuations, whereas the large FD kernel is more noise-robust. However, large FD kernels lead to extended numerical boundary error propagation, which severely hampers the MR-EPT reconstruction accuracy for highly spatially convoluted tissue structures such as the human brain. Mitigation strategies slightly improve the accuracy of conductivity reconstructions. For the adopted derivative kernels and the investigated scenario, MR-EPT conductivity reconstructions show low accuracy: less than 37% of the voxels have a relative error lower than 30%.

CONCLUSION

The numerical error introduced by the computation of spatial derivatives using FD kernels is one of the major causes of limited accuracy in Helmholtz-based MR-EPT reconstructions. Magn Reson Med 80:90-100, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Electrical Properties Tomography Based on $B_{{1}}$ Maps in MRI: Principles, Applications, and Challenges

OBJECTIVE

The purpose is to provide a comprehensive review of the electrical properties tomography (EPT) technique, which was introduced to image the electrical properties (EPs) of tissue noninvasively by exploiting the measured field data of MRI.

METHODS

We reviewed the principle of EPT, reconstruction methods, biomedical applications such as tumor imaging, and existing challenges. As a key application of EPT, the estimation of specific absorption rate (SAR) due to MRI was discussed in the background of elevated risk of tissue heating at high field.

RESULTS AND CONCLUSION

Since the originally proposed local, homogeneous Helmholtz equation-based reconstruction algorithm, advanced EPT algorithms have emerged to address the challenges of EPT, including reconstruction error near tissue boundaries, noise sensitivity, inaccurate phase estimation, and elimination of the unmeasurable component, along with demonstrations of in vivo experiments. EPT techniques have been applied to investigate EPs of both healthy and pathological tissues in vivo and factors contributing to various EP value, including sodium, water content, etc. More studies are anticipated to consolidate the current findings. EPT-based subject-specific SAR estimation has led to in vivo demonstration of its feasibility and prediction of temperature increase of phantom during MRI scans merely using measured data.

SIGNIFICANCE

EPT has the advantage of high resolution and practical feasibility in a clinical setup for imaging the biomedically interesting EPs of tissue in the radiofrequency range. EPT-based SAR estimation is another promising topic for predicting tissue heating of individual subjects during a specific MRI scan.

Electric properties tomography: Biochemical, physical and technical background, evaluation and clinical applications

Electric properties tomography (EPT) derives the patient's electric properties, i.e. conductivity and permittivity, using standard magnetic resonance (MR) systems and standard MR sequences. Thus, EPT does not apply externally mounted electrodes, currents or radiofrequency (RF) probes, as is the case in competing techniques. EPT is quantitative MR, i.e. it yields absolute values of conductivity and permittivity. This review summarizes the physical equations underlying EPT, the corresponding basic and advanced reconstruction techniques and practical numerical aspects to realize these reconstruction techniques. MR sequences which map the field information required for EPT are outlined, and experiments to validate EPT in phantom and in vivo studies are described. Furthermore, the review describes the clinical findings which have been obtained with EPT so far, and attempts to understand the physiologic background of these findings.

Fast Spin Echo Imaging-Based Electric Property Tomography With K-Space Weighting via ${T}_{2}$ Relaxation (rEPT)

Magnetic resonance electrical property tomography (MREPT) is a technique used to extract the electrical properties of tissues (conductivity in particular) using a magnetic resonance imaging system. In this paper, we propose an improved data acquisition scheme for the electrical property tomography technique by utilizing T ₂ modulation in fast spin echo (FSE) imaging. This technique was motivated by a numerical analysis of conductivity reconstruction in the frequency domain; results reveal the spatial frequency-dependent noise texture of conventional methods. A data-acquisition scheme using the FSE sequence was formulated to concentrate the signal within a specific frequency range where notable noise amplification is observed in the conventional method. Through numerical studies, the performance of the proposed acquisition was investigated. Furthermore, a compensation scheme was applied to reduce quantification errors due to tissue-specific T ₂ modulation, which is inherent in FSE imaging. The technique was applied to phantom and in vivo experiments. Results showed improved conductivity contrasts in both experiments, as compared with conventional MREPT methods.

Parallel transmission pulse design with explicit control for the specific absorption rate in the presence of radiofrequency errors

PURPOSE

A new framework for the design of parallel transmit (pTx) pulses is presented introducing constraints for local and global specific absorption rate (SAR) in the presence of errors in the radiofrequency (RF) transmit chain.

METHODS

The first step is the design of a pTx RF pulse with explicit constraints for global and local SAR. Then, the worst possible SAR associated with that pulse due to RF transmission errors ("worst-case SAR") is calculated. Finally, this information is used to re-calculate the pulse with lower SAR constraints, iterating this procedure until its worst-case SAR is within safety limits.

RESULTS

Analysis of an actual pTx RF transmit chain revealed amplitude errors as high as 8% (20%) and phase errors above 3° (15°) for spokes (spiral) pulses. Simulations show that using the proposed framework, pulses can be designed with controlled "worst-case SAR" in the presence of errors of this magnitude at minor cost of the excitation profile quality.

CONCLUSION

Our worst-case SAR-constrained pTx design strategy yields pulses with local and global SAR within the safety limits even in the presence of RF transmission errors. This strategy is a natural way to incorporate SAR safety factors in the design of pTx pulses. Magn Reson Med 75:2493-2504, 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.

Electrical conductivity and permittivity maps of brain tissues derived from water content based on T1 -weighted acquisition

PURPOSE

To develop an electrical properties tomography (EPT) technique that can provide in vivo electrical conductivity and permittivity images of biological tissue without performing complex-valued radiofrequency field measurements.

THEORY AND METHODS

Electrical conductivity and permittivity images are modeled as a monotonic function of tissues' water content (W) under the principle of Maxwell's mixture theory. Water content maps are estimated from two spin-echo images having different repetition times (TRs). For the modeling functions, physically measured parameters (electrical properties, water content, and T₁ ) of brain cerebrospinal fluid (CSF), gray matter, and white matter are used as landmark literature references. The formulations are validated by a developed electrolyte-protein phantom and by human brain studies at 3 Tesla (T).

RESULTS

The electrical properties (EPs) of the phantom estimated by the proposed method match well with the values measured on the bench. The conductivity and permittivity maps from all experiments show uncompromised spatial resolution without boundary artifacts and higher contrast when compared with water content maps.

CONCLUSIONS

Human brain and phantom EP images suggest that water content is a dominating factor in determining the electrical properties of tissues. Despite possible literature inaccuracies, the proposed method offers EP maps that can provide complementary information to current approaches, to facilitate EPT scans in clinical applications. Magn Reson Med 77:1094-1103, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Gradient-based electrical conductivity imaging using MR phase

PURPOSE

To develop a fast, practically applicable, and boundary artifact free electrical conductivity imaging method that does not use transceive phase assumption, and that is more robust against the noise.

THEORY

Starting from the Maxwell's equations, a new electrical conductivity imaging method that is based solely on the MR transceive phase has been proposed. Different from the previous phase based electrical properties tomography (EPT) method, a new formulation was derived by including the gradients of the conductivity into the equations.

METHODS

The governing partial differential equation, which is in the form of a convection-reaction-diffusion equation, was solved using a three-dimensional finite-difference scheme. To evaluate the performance of the proposed method numerical simulations, phantom and in vivo human experiments have been conducted at 3T.

RESULTS

Simulation and experimental results of the proposed method and the conventional phase-based EPT method were illustrated to show the superiority of the proposed method over the conventional method, especially in the transition regions and under noisy data.

CONCLUSION

With the contributions of the proposed method to the phase-based EPT approach, a fast and reliable electrical conductivity imaging appears to be feasible, which is promising for clinical diagnoses and local SAR estimation. Magn Reson Med 77:137-150, 2017. © 2016 Wiley Periodicals, Inc.

A new approach for electrical properties estimation using a global integral equation and improvements using high permittivity materials

Electrical Properties Tomography (EPT) using MRI is a technique that has been developed to provide a new contrast mechanism for in vivo imaging. Currently the most common method relies on the solution of the homogeneous Helmholtz equation, which has limitations in accurate estimation at tissue interfaces. A new method proposed in this work combines a Maxwell's integral equation representation of the problem, and the use of high permittivity materials (HPM) to control the RF field, in order to reconstruct the electrical properties image. The magnetic field is represented by an integral equation considering each point as a contrast source. This equation can be solved in an inverse method. In this study we use a reference simulation or scout scan of a uniform phantom to provide an initial estimate for the inverse solution, which allows the estimation of the complex permittivity within a single iteration. Incorporating two setups with and without the HPM improves the reconstructed result, especially with respect to the very low electric field in the center of the sample. Electromagnetic simulations of the brain were performed at 3T to generate the B1(+) field maps and reconstruct the electric properties images. The standard deviations of the relative permittivity and conductivity were within 14% and 18%, respectively for a volume consisting of white matter, gray matter and cerebellum.

Theoretical Investigation of Random Noise-Limited Signal-to-Noise Ratio in MR-Based Electrical Properties Tomography

In magnetic resonance imaging-based electrical properties tomography (MREPT), tissue electrical properties (EPs) are derived from the spatial variation of the transmit RF field (B1(+)). Here we derive theoretically the relationship between the signal-to-noise ratio (SNR) of the electrical properties obtained by MREPT and the SNR of the input B1(+) data, under the assumption that the latter is much greater than unity, and the noise in B1(+) at different voxels is statistically independent. It is shown that for a given B1(+) data, the SNR of both electrical conductivity and relative permittivity is proportional to the square of the linear dimension of the region of interest (ROI) over which the EPs are determined, and to the square root of the number of voxels in the ROI. The relationship also shows how the SNR varies with the main magnetic field (B0) strength. The predicted SNR is verified through numerical simulations on a cylindrical phantom with an analytically calculated B1(+) map, and is found to provide explanation of certain aspects of previous experimental results in the literature. Our SNR formula can be used to estimate minimum input data SNR and ROI size required to obtain tissue EP maps of desired quality.

Tissue electrical property mapping from zero echo-time magnetic resonance imaging

The capability of magnetic resonance imaging (MRI) to produce spatially resolved estimation of tissue electrical properties (EPs) in vivo has been a subject of much recent interest. In this work we introduce a method to map tissue EPs from low-flip-angle, zero-echo-time (ZTE) imaging. It is based on a new theoretical formalism that allows calculation of EPs from the product of transmit and receive radio-frequency (RF) field maps. Compared to conventional methods requiring separation of the transmit RF field (B(1)(+)) from acquired MR images, the proposed method has such advantages as: 1) reduced theoretical error, 2) higher acquisition speed, and 3) flexibility in choice of different transmit and receive RF coils. The method is demonstrated in electrical conductivity and relative permittivity mapping in a salt water phantom, as well as in vivo measurement of brain conductivity in healthy volunteers. The phantom results show the validity and scan-time efficiency of the proposed method applied to a piece-wise homogeneous object. Quality of in vivo EP results was limited by reconstruction errors near tissue boundaries, which highlights need for image segmentation in EP mapping in a heterogeneous medium. Our results show the feasibility of rapid EP mapping from MRI without B(1)(+) mapping.