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

MR electrical properties mapping using vision transformers and canny edge detectors

PURPOSE

We developed a 3D vision transformer-based neural network to reconstruct electrical properties (EP) from magnetic resonance measurements.

THEORY AND METHODS

Our network uses the magnitude of the transmit magnetic field of a birdcage coil, the associated transceive phase, and a Canny edge mask that identifies the object boundaries as inputs to compute the EP maps. We trained our network on a dataset of 10 000 synthetic tissue-mimicking phantoms and fine-tuned it on a dataset of 11 000 realistic head models. We assessed performance in-distribution simulated data and out-of-distribution head models, with and without synthetic lesions. We further evaluated our network in experiments for an inhomogeneous phantom and a volunteer.

RESULTS

The conductivity and permittivity maps had an average peak normalized absolute error (PNAE) of 1.3% and 1.7% for the synthetic phantoms, respectively. For the realistic heads, the average PNAE for the conductivity and permittivity was 1.8% and 2.7%, respectively. The location of synthetic lesions was accurately identified, with reconstructed conductivity and permittivity values within 15% and 25% of the ground-truth, respectively. The conductivity and permittivity for the phantom experiment yielded 2.7% and 2.1% average PNAEs with respect to probe-measured values, respectively. The in vivo EP reconstruction truthfully preserved the subject's anatomy with average values over the entire head similar to the expected literature values.

CONCLUSION

We introduced a new learning-based approach for reconstructing EP from MR measurements obtained with a birdcage coil, marking an important step towards the development of clinically-usable in vivo EP reconstruction protocols.

A joint three-plane physics-constrained deep learning based polynomial fitting approach for MR electrical properties tomography

Magnetic resonance electrical properties tomography can extract the electrical properties of in-vivo tissue. To estimate tissue electrical properties, various reconstruction algorithms have been proposed. However, physics-based reconstructions are prone to various artifacts such as noise amplification and boundary artifact. Deep learning-based approaches are robust to these artifacts but need extensive training datasets and suffer from generalization to unseen data. To address these issues, we introduce a joint three-plane physics-constrained deep learning framework for polynomial fitting MR-EPT by merging physics-based weighted polynomial fitting with deep learning. Within this framework, deep learning is used to discern the optimal polynomial fitting weights for a physics based polynomial fitting reconstruction on the complex B1+ data. For the prediction of optimal fitting coefficients, three neural networks were separately trained on simulated heterogeneous brain models to predict optimal polynomial weighting parameters in three orthogonal planes. Then, the network weights were jointly optimized to estimate the polynomial weights in each plane for a combined conductivity reconstruction. Based on this physics-constrained deep learning approach, we achieved an improvement of conductivity estimation accuracy in comparison to a single plane estimation and a reduction of computational load. The results demonstrate that the proposed method based on 3D data exhibits superior performance in comparison to conventional polynomial fitting methods in terms of capturing anatomical detail and homogeneity. Crucially, in-vivo application of the proposed method showed that the method generalizes well to in-vivo data, without introducing significant errors or artifacts. This generalization makes the presented method a promising candidate for use in clinical applications.

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.

Unrolled Optimization via Physics-Assisted Convolutional Neural Network for MR-Based Electrical Properties Tomography: A Numerical Investigation

Magnetic Resonance imaging based Electrical Properties Tomography (MR-EPT) is a non-invasive technique that measures the electrical properties (EPs) of biological tissues. In this work, we present and numerically investigate the performance of an unrolled, physics-assisted method for 2D MR-EPT reconstructions, where a cascade of Convolutional Neural Networks is used to compute the contrast update. Each network takes in input the EPs and the gradient descent direction (encoding the physics underlying the adopted scattering model) and returns as output the updated contrast function. The network is trained and tested in silico using 2D slices of realistic brain models at 128 MHz. Results show the capability of the proposed procedure to reconstruct EPs maps with quality comparable to that of the popular Contrast Source Inversion-EPT, while significantly reducing the computational time.

A database for MR-based electrical properties tomography with in silico brain data-ADEPT

PURPOSE

Several reconstruction methods for MR-based electrical properties tomography (EPT) have been developed. However, the lack of common data makes it difficult to objectively compare their performances. This is, however, a necessary precursor for standardizing and introducing this technique in the clinical setting. To enable objective comparison of the performances of reconstruction methods and provide common data for their training and testing, we created ADEPT, a database of simulated data for brain MR-EPT reconstructions.

METHODS

ADEPT is a database containing in silico data for brain EPT reconstructions. This database was created from 25 different brain models, with and without tumors. Rigid geometric augmentations were applied, and different electrical properties were assigned to white matter, gray matter, CSF, and tumors to generate 120 different brain models. These models were used as input for finite-difference time-domain simulations in Sim4Life, used to compute the electromagnetic fields needed for MR-EPT reconstructions.

RESULTS

Electromagnetic fields from 84 healthy and 36 tumor brain models were simulated. The simulated fields relevant for MR-EPT reconstructions (transmit and receive RF fields and transceive phase) and their ground-truth electrical properties are made publicly available through ADEPT. Additionally, nonattainable fields such as the total magnetic field and the electric field are available upon request.

CONCLUSION

ADEPT will serve as reference database for objective comparisons of reconstruction methods and will be a first step toward standardization of MR-EPT reconstructions. Furthermore, it provides a large amount of data that can be exploited to train data-driven methods. It can be accessed from  https://doi.org/10.34894/V0HBJ8.

PIFON-EPT: MR-Based Electrical Property Tomography Using Physics-Informed Fourier Networks

We propose Physics-Informed Fourier Networks for Electrical Properties (EP) Tomography (PIFON-EPT), a novel deep learning-based method for EP reconstruction using noisy and/or incomplete magnetic resonance (MR) measurements. Our approach leverages the Helmholtz equation to constrain two networks, responsible for the denoising and completion of the transmit fields, and the estimation of the object's EP, respectively. We embed a random Fourier features mapping into our networks to enable efficient learning of high-frequency details encoded in the transmit fields. We demonstrated the efficacy of PIFON-EPT through several simulated experiments at 3 and 7 tesla(T) MR imaging, and showed that our method can reconstruct physically consistent EP and transmit fields. Specifically, when only 20% of the noisy measured fields were used as inputs, PIFON-EPT reconstructed the EP of a phantom with ≤ 5% error, and denoised and completed the measurements with ≤ 1% error. Additionally, we adapted PIFON-EPT to solve the generalized Helmholtz equation that accounts for gradients of EP between inhomogeneities. This yielded improved results at interfaces between different materials without explicit knowledge of boundary conditions. PIFON-EPT is the first method that can simultaneously reconstruct EP and transmit fields from incomplete noisy MR measurements, providing new opportunities for EPT research.

Data-driven electrical conductivity brain imaging using 3 T MRI

Magnetic resonance electrical properties tomography (MR-EPT) is a non-invasive measurement technique that derives the electrical properties (EPs, e.g., conductivity or permittivity) of tissues in the radiofrequency range (64 MHz for 1.5 T and 128 MHz for 3 T MR systems). Clinical studies have shown the potential of tissue conductivity as a biomarker. To date, model-based conductivity reconstructions rely on numerical assumptions and approximations, leading to inaccuracies in the reconstructed maps. To address such limitations, we propose an artificial neural network (ANN)-based non-linear conductivity estimator trained on simulated data for conductivity brain imaging. Network training was performed on 201 synthesized T2-weighted spin-echo (SE) data obtained from the finite-difference time-domain (FDTD) electromagnetic (EM) simulation. The dataset was composed of an approximated T2-w SE magnitude and transceive phase information. The proposed method was tested three in-silico and in-vivo on two volunteers and three patients' data. For comparison purposes, various conventional phase-based EPT reconstruction methods were used that ignoreB 1 + magnitude information, such as Savitzky-Golay kernel combined with Gaussian filter (S-G Kernel), phase-based convection-reaction EPT (cr-EPT), magnitude-weighted polynomial-fitting phase-based EPT (Poly-Fit), and integral-based phase-based EPT (Integral-based). From the in-silico experiments, quantitative analysis showed that the proposed method provides more accurate and improved quality (e.g., high structural preservation) conductivity maps compared to conventional reconstruction methods. Representatively, in the healthy brain in-silico phantom experiment, the proposed method yielded mean conductivity values of 1.97 ± 0.20 S/m for CSF, 0.33 ± 0.04 S/m for WM, and 0.52 ± 0.08 S/m for GM, which were closer to the ground-truth conductivity (2.00, 0.30, 0.50 S/m) than the integral-based method (2.56 ± 2.31, 0.39 ± 0.12, 0.68 ± 0.33 S/m). In-vivo ANN-based conductivity reconstructions were also of improved quality compared to conventional reconstructions and demonstrated network generalizability and robustness to in-vivo data and pathologies. The reported in-vivo brain conductivity values were in agreement with literatures. In addition, the proposed method was observed for various SNR levels (SNR levels = 10, 20, 40, and 58) and repeatability conditions (the eight acquisitions with the number of signal averages = 1). The preliminary investigations on brain tumor patient datasets suggest that the network trained on simulated dataset can generalize to unforeseen in-vivo pathologies, thus demonstrating its potential for clinical applications.

MARIE 2.0: A Perturbation Matrix Based Patient-Specific MRI Field Simulator

High static field MR scanners can produce human tissue images of astounding clarity, but rely on high frequency electromagnetic radiation that generates complicated in-tissue field patterns that are patient-specific and potentially harmful. Many such scanners use parallel transmitters to better control field patterns, but then adjust the transmitters based on general guidelines rather than optimizing for the specific patient, mostly because computing patient-specific fields was presumed far too slow. It was recently demonstrated that the combination of fast low-resolution tissue mapping and fast voxel-based field simulation can be used to perform a patient-specific MR safety check in minutes. However, the field simulation required several of those minutes, making it too slow to perform the dozens of simulations that would be needed for patient-specific optimization. In this paper we describe a compressed-perturbation-matrix technique that nearly eliminates the computational cost of including complex coils (or coils and shields) in voxel-based field simulation of tissue, thereby reducing simulation time from minutes to seconds. The approach is demonstrated on a wide variety of head+coil and head+coil+shield configurations, using the implementation in MARIE 2.0, the latest version of the open-source MR field simulator MARIE.

Repeatability and Reproducibility Uncertainty in Magnetic Resonance-Based Electric Properties Tomography of a Homogeneous Phantom

Uncertainty assessment is a fundamental step in quantitative magnetic resonance imaging because it makes comparable, in a strict metrological sense, the results of different scans, for example during a longitudinal study. Magnetic resonance-based electric properties tomography (EPT) is a quantitative imaging technique that retrieves, non-invasively, a map of the electric properties inside a human body. Although EPT has been used in some early clinical studies, a rigorous experimental assessment of the associated uncertainty has not yet been performed. This paper aims at evaluating the repeatability and reproducibility uncertainties in phase-based Helmholtz-EPT applied on homogeneous phantom data acquired with a clinical 3 T scanner. The law of propagation of uncertainty is used to evaluate the uncertainty in the estimated conductivity values starting from the uncertainty in the acquired scans, which is quantified through a robust James-Stein shrinkage estimator to deal with the dimensionality of the problem. Repeatable errors are detected in the estimated conductivity maps and are quantified for various values of the tunable parameters of the EPT implementation. The spatial dispersion of the estimated electric conductivity maps is found to be a good approximation of the reproducibility uncertainty, evaluated by changing the position of the phantom after each scan. The results underpin the use of the average conductivity (calculated by weighting the local conductivity values by their uncertainty and taking into account the spatial correlation) as an estimate of the conductivity of the homogeneous phantom.

A Hybrid Volume-Surface Integral Equation Method for Rapid Electromagnetic Simulations in MRI

OBJECTIVE

We developed a hybrid volume surface integral equation (VSIE) method based on domain decomposition to perform fast and accurate magnetic resonance imaging (MRI) simulations that include both remote and local conductive elements.

METHODS

We separated the conductive surfaces present in MRI setups into two domains and optimized electromagnetic (EM) modeling for each case. Specifically, interactions between the body and EM waves originating from local radiofrequency (RF) coils were modeled with the precorrected fast Fourier transform, whereas the interactions with remote conductive surfaces (RF shield, scanner bore) were modeled with a novel cross tensor train-based algorithm. We compared the hybrid-VSIE with other VSIE methods for realistic MRI simulation setups.

RESULTS

The hybrid-VSIE was the only practical method for simulation using 1 mm voxel isotropic resolution (VIR). For 2 mm VIR, our method could be solved at least 23 times faster and required 760 times lower memory than traditional VSIE methods.

CONCLUSION

The hybrid-VSIE demonstrated a marked improvement in terms of convergence times of the numerical EM simulation compared to traditional approaches in multiple realistic MRI scenarios.

SIGNIFICANCE

The efficiency of the novel hybrid-VSIE method could enable rapid simulations of complex and comprehensive MRI setups.

Novel Numerical Basis Sets for Electromagnetic Field Expansion in Arbitrary Inhomogeneous Objects

We investigated how to construct low-order subspace basis sets to accurately represent electromagnetic fields generated within inhomogeneous arbitrary objects by radio-frequency sources external to a Huygen's surface. The basis generation relies on the singular value decomposition of Green's functions integro-differential operators which makes it feasible to derive a reduced-order yet stable model. We present a detailed study of the theoretical and numerical requisites for generating such basis, and show how it can be used to calculate performance limits in magnetic resonance imaging applications. Finally, we propose a novel numerical framework for the computation of characteristic modes of arbitrary inhomogeneous objects. We validated accuracy and convergence properties of the numerical basis against a complete analytical basis in the case of a uniform spherical object. We showed that the discretization of the Huygens's surface has a minimal effect on the accuracy of the calculations, which mainly depended on the electromagnetic solver resolution and order of approximation.

A fast tool for the parametric analysis of human body exposed to LF electromagnetic fields in biomedical applications

A numerical procedure for analyzing electromagnetic (EM) fields interactions with biological tissues is presented. The proposed approach aims at drastically reducing the computational burden required by the repeated solution of large scale problems involving the interaction of the human body with EM fields, such as in the study of the time evolution of EM fields, uncertainty quantification, and inverse problems. The proposed volume integral equation (VIE), focused on low frequency applications, is a system of integral equations in terms of current density and scalar potential in the biological tissues excited by EM fields and/or electrodes connected to the human body. The proposed formulation requires the voxelization of the human body and takes advantage of the regularity of such discretization by speeding-up the computational procedure. Moreover, it exploits recent advancements in the solution of VIE by means of iterative preconditioned solvers and ad hoc parametric Model Order Reduction techniques. The efficiency of the proposed tool is demonstrated by applying it to a couple of realistic model problems: the assessment of the peripheral nerve stimulation, performed in terms of evaluation of the induced electric field, due to the gradient coils of a magnetic resonance imaging scanner during a clinical examination and the assessment of the exposure to environmental fields at 50 Hz of live-line workers with uncertain properties of the biological tissues. Thanks to the proposed method, uncertainty quantification analyses and time domain simulations are possible even for large scale problems and they can be performed on standard computers and reasonable computation time. Sample implementation of the method is made publicly available at https://github.com/UniPD-DII-ETCOMP/BioMOR.

Compression of volume-surface integral equation matrices via Tucker decomposition for magnetic resonance applications

In this work, we propose a method for the compression of the coupling matrix in volume-surface integral equation (VSIE) formulations. VSIE methods are used for electromagnetic analysis in magnetic resonance imaging (MRI) applications, for which the coupling matrix models the interactions between the coil and the body. We showed that these effects can be represented as independent interactions between remote elements in 3D tensor formats, and subsequently decomposed with the Tucker model. Our method can work in tandem with the adaptive cross approximation technique to provide fast solutions of VSIE problems. We demonstrated that our compression approaches can enable the use of VSIE matrices of prohibitive memory requirements, by allowing the effective use of modern graphical processing units (GPUs) to accelerate the arising matrix-vector products. This is critical to enable numerical MRI simulations at clinical voxel resolutions in a feasible computation time. In this paper, we demonstrate that the VSIE matrix-vector products needed to calculate the electromagnetic field produced by an MRI coil inside a numerical body model with 1 mm³ voxel resolution, could be performed in ~ 33 seconds in a GPU, after compressing the associated coupling matrix from ~ 80 TB to ~ 43 MB.

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.

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.

Deep Learning-Based Development of Personalized Human Head Model With Non-Uniform Conductivity for Brain Stimulation

Electromagnetic stimulation of the human brain is a key tool for neurophysiological characterization and the diagnosis of several neurological disorders. Transcranial magnetic stimulation (TMS) is a commonly used clinical procedure. However, personalized TMS requires a pipeline for individual head model generation to provide target-specific stimulation. This process includes intensive segmentation of several head tissues based on magnetic resonance imaging (MRI), which has significant potential for segmentation error, especially for low-contrast tissues. Additionally, a uniform electrical conductivity is assigned to each tissue in the model, which is an unrealistic assumption based on conventional volume conductor modeling. This study proposes a novel approach for fast and automatic estimation of the electric conductivity in the human head for volume conductor models without anatomical segmentation. A convolutional neural network is designed to estimate personalized electrical conductivity values based on anatomical information obtained from T1- and T2-weighted MRI scans. This approach can avoid the time-consuming process of tissue segmentation and maximize the advantages of position-dependent conductivity assignment based on the water content values estimated from MRI intensity values. The computational results of the proposed approach provide similar but smoother electric field distributions of the brain than that provided by conventional approaches.

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.

Deep learning-based reconstruction of in vivo pelvis conductivity with a 3D patch-based convolutional neural network trained on simulated MR data

PURPOSE

To demonstrate that mapping pelvis conductivity at 3T with deep learning (DL) is feasible.

METHODS

210 dielectric pelvic models were generated based on CT scans of 42 cervical cancer patients. For all dielectric models, electromagnetic and MR simulations with realistic accuracy and precision were performed to obtain B 1 + and transceive phase (ϕ^± ). Simulated B 1 + and ϕ^± served as input to a 3D patch-based convolutional neural network, which was trained in a supervised fashion to retrieve the conductivity. The same network architecture was retrained using only ϕ^± in input. Both network configurations were tested on simulated MR data and their conductivity reconstruction accuracy and precision were assessed. Furthermore, both network configurations were used to reconstruct conductivity maps from a healthy volunteer and two cervical cancer patients. DL-based conductivity was compared in vivo and in silico to Helmholtz-based (H-EPT) conductivity.

RESULTS

Conductivity maps obtained from both network configurations were comparable. Accuracy was assessed by mean error (ME) with respect to ground truth conductivity. On average, ME < 0.1 Sm^-1 for all tissues. Maximum MEs were 0.2 Sm^-1 for muscle and tumour, and 0.4 Sm^-1 for bladder. Precision was indicated with the difference between 90^th and 10^th conductivity percentiles, and was below 0.1 Sm^-1 for fat, bone and muscle, 0.2 Sm^-1 for tumour and 0.3 Sm^-1 for bladder. In vivo, DL-based conductivity had median values in agreement with H-EPT values, but a higher precision.

CONCLUSION

Anatomically detailed, noise-robust 3D conductivity maps with good sensitivity to tissue conductivity variations were reconstructed in the pelvis with DL.

Learning-based estimation of dielectric properties and tissue density in head models for personalized radio-frequency dosimetry

Radio-frequency dosimetry is an important process in assessments for human exposure safety and for compliance of related products. Recently, computational human models generated from medical images have often been used for such assessment, especially to consider the inter-subject variability. However, a common procedure to develop personalized models is time consuming because it involves excessive segmentation of several components that represent different biological tissues, which is a major obstacle in the inter-subject variability assessment of radiation safety. Deep learning methods have been shown to be a powerful approach for pattern recognition and signal analysis. Convolutional neural networks with deep architecture are proven robust for feature extraction and image mapping in several biomedical applications. In this study, we develop a learning-based approach for fast and accurate estimation of the dielectric properties and density of tissues directly from magnetic resonance images in a single shot. The smooth distribution of the dielectric properties in head models, which is realized using a process without tissue segmentation, improves the smoothness of the specific absorption rate (SAR) distribution compared with that in the commonly used procedure. The estimated SAR distributions, as well as that averaged over 10 g of tissue in a cubic shape, are found to be highly consistent with those computed using the conventional methods that employ segmentation.

Channel-combination method for phase-based |B1+| mapping techniques

PURPOSE

The aim of this study was to propose a channel combination method for |B₁⁺| mapping methods using phase difference to reconstruct |B₁⁺| map.

THEORY AND METHODS

Phase-based |B₁⁺| mapping methods commonly consider the phase difference of two scans to measure |B₁⁺|. Multiple receiver coils acquire a number of images and the phase difference at each channel is theoretically the same in the absence of noise. Affected by noise, phase difference is approximately governed by Gaussian distribution. Considering data from all channels as samples, estimation can be achieved by maximum likelihood method. With this method, all phase differences at each channel are combined into one. In this study, the proposed method is applied with Bloch-Siegert shift |B₁⁺| mapping method. Simulations are performed to illustrate the phase difference distribution and demonstrate the feasibility and facility of the proposed method. Phantom and vivo experiments are carried out at 1.5 T scanner equipped with 8-channel receiver coil. In all experiments, the proposed method is compared with weighted averaging (WA) method.

RESULTS

Simulations revealed appropriateness of approximating the distribution of phase difference to Gaussian distribution. Compared with WA method, the proposed method reduces errors of |B₁⁺| calculation. Phantom and vivo experiments provide further validation.

CONCLUSION

Considering phase noise distribution, the proposed method achieves channel combination by finding the estimation from data acquired by multiple receivers coil. The proposed method reduces |B₁⁺| reconstruction errors caused by noise.

Transceive phase mapping using the PLANET method and its application for conductivity mapping in the brain

Purpose

To demonstrate feasibility of transceive phase mapping with the PLANET method and its application for conductivity reconstruction in the brain.

Methods

Accuracy and precision of transceive phase (ϕ^±) estimation with PLANET, an ellipse fitting approach to phase‐cycled balanced steady state free precession (bSSFP) data, were assessed with simulations and measurements and compared to standard bSSFP. Measurements were conducted on a homogeneous phantom and in the brain of healthy volunteers at 3 tesla. Conductivity maps were reconstructed with Helmholtz‐based electrical properties tomography. In measurements, PLANET was also compared to a reference technique for transceive phase mapping, i.e., spin echo.

Results

Accuracy and precision of ϕ^± estimated with PLANET depended on the chosen flip angle and TR. PLANET‐based ϕ^± was less sensitive to perturbations induced by off‐resonance effects and partial volume (e.g., white matter + myelin) than bSSFP‐based ϕ^±. For flip angle = 25° and TR = 4.6 ms, PLANET showed an accuracy comparable to that of reference spin echo but a higher precision than bSSFP and spin echo (factor of 2 and 3, respectively). The acquisition time for PLANET was ~5 min; 2 min faster than spin echo and 8 times slower than bSSFP. However, PLANET simultaneously reconstructed T₁, T₂, B₀ maps besides mapping ϕ^±. In the phantom, PLANET‐based conductivity matched the true value and had the smallest spread of the three methods. In vivo, PLANET‐based conductivity was similar to spin echo‐based conductivity.

Conclusion

Provided that appropriate sequence parameters are used, PLANET delivers accurate and precise ϕ^± maps, which can be used to reconstruct brain tissue conductivity while simultaneously recovering T₁, T₂, and B₀ maps.