3-D Contrast Source Inversion-Electrical Properties Tomography

Repeatability of brain phase-based magnetic resonance electric properties tomography methods and effect of compressed SENSE and RF shimming

Magnetic resonance electrical properties tomography (MREPT) is an emerging imaging modality to noninvasively measure tissue conductivity and permittivity. Implementation of MREPT in the clinic requires repeatable measurements at a short scan time and an appropriate protocol. The aim of this study was to investigate the repeatability of conductivity measurements using phase-based MREPT and the effects of compressed SENSE (CS), and RF shimming on the precision of conductivity measurements. Conductivity measurements using turbo spin echo (TSE) and three-dimensional balanced fast field echo (bFFE) with CS factors were repeatable. Conductivity measurement using bFFE phase showed smaller mean and variance that those measured by TSE. The conductivity measurements using bFFE showed minimal deviation with CS factors up to 8, with deviation increasing at CS factors > 8. Subcortical structures produced less consistent measurements than cortical parcellations at higher CS factors. RF shimming using full slice coverage 2D dual refocusing echo acquisition mode (DREAM) and full coverage 3D dual TR approaches further improved measurement precision. BFFE is a more optimal sequence than TSE for phase-based MREPT in brain. Depending on the area of the brain being measured, the scan can be safely accelerated with compressed SENSE without sacrifice of precision, offering the potential to employ MREPT in clinical research and applications. RF shimming with better field mapping further improves precision of the conductivity measures.

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 Method for Electrical Property Tomography Based on a Three-Dimensional Integral Representation of the Electric Field

Magnetic resonance electrical properties tomography (MREPT) noninvasively reconstructs high-resolution electrical property (EP) maps using MRI scanners and is useful for diagnosing cancerous tissues. However, conventional MREPT methods have limitations: sensitivity to noise in the numerical Laplacian operation, difficulty in reconstructing three-dimensional (3D) EPs and convergence not guaranteed in the iterative process. We propose a novel, iterative 3D reconstruction MREPT method without a numerical Laplacian operation. We derive an integral representation of the electric field using its Helmholtz decomposition with Maxwell's equations, under the assumption that the EPs are known on the boundary of the region of interest with the approximation that the unmeasurable magnetic field components are zero. Then, we solve the simultaneous equations composed of the integral representation and Ampere's law using a convex projection algorithm whose convergence is theoretically guaranteed. The efficacy of the proposed method was validated through numerical simulations and a phantom experiment. The results showed that this method is effective in reconstructing 3D EPs and is robust to noise. It was also shown that our proposed method with the unmeasurable component H^- enhances the accuracy of the EPs in a background and that with all the components of the magnetic field reduces the artifacts at the center of the slices except when all the components of the electric field are close to zero.

Combining deep learning and 3D contrast source inversion in MR-based electrical properties tomography

Magnetic resonance electrical properties tomography (MR-EPT) is a technique used to estimate the conductivity and permittivity of tissues from MR measurements of the transmit magnetic field. Different reconstruction methods are available; however, all these methods present several limitations, which hamper the clinical applicability. Standard Helmholtz-based MR-EPT methods are severely affected by noise. Iterative reconstruction methods such as contrast source inversion electrical properties tomography (CSI-EPT) are typically time-consuming and are dependent on their initialization. Deep learning (DL) based methods require a large amount of training data before sufficient generalization can be achieved. Here, we investigate the benefits achievable using a hybrid approach, that is, using MR-EPT or DL-EPT as initialization guesses for standard 3D CSI-EPT. Using realistic electromagnetic simulations at 3 and 7 T, the accuracy and precision of hybrid CSI reconstructions are compared with those of standard 3D CSI-EPT reconstructions. Our results indicate that a hybrid method consisting of an initial DL-EPT reconstruction followed by a 3D CSI-EPT reconstruction would be beneficial. DL-EPT combined with standard 3D CSI-EPT exploits the power of data-driven DL-based EPT reconstructions, while the subsequent CSI-EPT facilitates a better generalization by providing data consistency.

High-frequency conductivity at Larmor-frequency in human brain using moving local window multilayer perceptron neural network

Magnetic resonance electrical properties tomography (MREPT) aims to visualize the internal high-frequency conductivity distribution at Larmor frequency using the B1 transceive phase data. From the magnetic field perturbation by the electrical field associated with the radiofrequency (RF) magnetic field, the high-frequency conductivity and permittivity distributions inside the human brain have been reconstructed based on the Maxwell’s equation. Starting from the Maxwell’s equation, the complex permittivity can be described as a second order elliptic partial differential equation. The established reconstruction algorithms have focused on simplifying and/or regularizing the elliptic partial differential equation to reduce the noise artifact. Using the nonlinear relationship between the Maxwell’s equation, measured magnetic field, and conductivity distribution, we design a deep learning model to visualize the high-frequency conductivity in the brain, directly derived from measured magnetic flux density. The designed moving local window multi-layer perceptron (MLW-MLP) neural network by sliding local window consisting of neighboring voxels around each voxel predicts the high-frequency conductivity distribution in each local window. The designed MLW-MLP uses a family of multiple groups, consisting of the gradients and Laplacian of measured B1 phase data, as the input layer in a local window. The output layer of MLW-MLP returns the conductivity values in each local window. By taking a non-local mean filtering approach in the local window, we reconstruct a noise suppressed conductivity image while maintaining spatial resolution. To verify the proposed method, we used B1 phase datasets acquired from eight human subjects (five subjects for training procedure and three subjects for predicting the conductivity in the brain).

Improving phase-based conductivity reconstruction by means of deep learning-based denoising of B 1 + phase data for 3T MRI

PURPOSE

To denoise B 1 + phase using a deep learning method for phase-based in vivo electrical conductivity reconstruction in a 3T MR system.

METHODS

For B 1 + phase deep-learning denoising, a convolutional neural network (U-net) was chosen. Training was performed on data sets from 10 healthy volunteers. Input data were the real and imaginary components of single averaged spin-echo data (SNR = 45), which was used to approximate the B 1 + phase. For label data, multiple signal-averaged spin-echo data (SNR = 128) were used. Testing was performed on in silico and in vivo data. Reconstructed conductivity maps were derived using phase-based conductivity reconstructions. Additionally, we investigated the usability of the network to various SNR levels, imaging contrasts, and anatomical sites (ie, T₁ , T₂ , and proton density-weighted brain images and proton density-weighted breast images. In addition, conductivity reconstructions from deep learning-based denoised data were compared with conventional image filters, which were used for data denoising in electrical properties tomography (ie, the Gaussian filtering and the Savitzky-Golay filtering).

RESULTS

The proposed deep learning-based denoising approach showed improvement for B 1 + phase for both in silico and in vivo experiments with reduced quantitative error measures compared with other methods. Subsequently, this resulted in an improvement of reconstructed conductivity maps from the denoised B 1 + phase with deep learning.

CONCLUSION

The results suggest that the proposed approach can be used as an alternative preprocessing method to denoise B 1 + maps for phase-based conductivity reconstruction without relying on image filters or signal averaging.

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.

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.

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.

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.

Accuracy and precision of electrical permittivity mapping at 3T: the impact of three B 1 + mapping techniques

Purpose

To investigate the sequence‐specific impact of B1+ amplitude mapping on the accuracy and precision of permittivity reconstruction at 3T in the pelvic region.

Methods

B1+ maps obtained with actual flip angle imaging (AFI), Bloch–Siegert (BS), and dual refocusing echo acquisition mode (DREAM) sequences, set to a clinically feasible scan time of 5 minutes, were compared in terms of accuracy and precision with electromagnetic and Bloch simulations and MR measurements. Permittivity maps were reconstructed based on these B1+ maps with Helmholtz‐based electrical properties tomography. Accuracy and precision in permittivity were assessed. A 2‐compartment phantom with properties and size similar to the human pelvis was used for both simulations and measurements. Measurements were also performed on a female volunteer’s pelvis.

Results

Accuracy was evaluated with noiseless simulations on the phantom. The maximum B1+ bias relative to the true B1+ distribution was 1% for AFI and BS and 6% to 15% for DREAM. This caused an average permittivity bias relative to the true permittivity of 7% to 20% for AFI and BS and 12% to 35% for DREAM. Precision was assessed in MR experiments. The lowest standard deviation in permittivity, found in the phantom for BS, measured 22.4 relative units and corresponded to a standard deviation in B1+ of 0.2% of the B1+ average value. As regards B1+ precision, in vivo and phantom measurements were comparable.

Conclusions

Our simulation framework quantitatively predicts the different impact of B1+ mapping techniques on permittivity reconstruction and shows high sensitivity of permittivity reconstructions to sequence‐specific bias and noise perturbation in the B1+ map. These findings are supported by the experimental results.

Developments in Electrical-Property Tomography Based on the Contrast-Source Inversion Method

The main objective of electrical-property tomography (EPT) is to retrieve dielectric tissue parameters from B ^ 1 + data as measured by a magnetic-resonance (MR) scanner. This is a so-called hybrid inverse problem in which data are defined inside the reconstruction domain of interest. In this paper, we discuss recent and new developments in EPT based on the contrast-source inversion (CSI) method. After a short review of the basics of this method, two- and three-dimensional implementations of CSI-EPT are presented along with a very efficient variant of 2D CSI-EPT called first-order induced current EPT (foIC-EPT). Practical implementation issues that arise when applying the method to measured data are addressed as well, and the limitations of a two-dimensional approach are extensively discussed. Tissue-parameter reconstructions of an anatomically correct male head model illustrate the performance of two- and three-dimensional CSI-EPT. We show that 2D implementation only produces reliable reconstructions under very special circumstances, while accurate reconstructions can be obtained with 3D CSI-EPT.

3-D Microwave Tomography Using the Soft Prior Regularization Technique: Evaluation in Anatomically Realistic MRI-Derived Numerical Breast Phantoms

OBJECTIVE

Fusion of magnetic resonance imaging (MRI) breast images with microwave tomography is accomplished through a soft prior technique, which incorporates spatial information (from MRI), i.e., accurate boundary location of different regions of interest, into the regularization process of the microwave image reconstruction algorithm.

METHODS

Numerical experiments were completed on a set of three-dimensional (3-D) breast geometries derived from MR breast data with different parenchymal densities, as well as a simulated tumor to evaluate the performance over a range of breast shapes, sizes, and property distributions.

RESULTS

When the soft prior regularization technique was applied, both permittivity and conductivity relative root mean square error values decreased by more than 87% across all breast densities, except in two cases where the error decrease was only 55% and 78%. In addition, the incorporation of structural priors increased contrast between tumor and fibroglandular tissue by 59% in permittivity and 192% in conductivity.

CONCLUSION

This study confirmed that the soft prior algorithm is robust in 3-D and can function successfully across a range of complex geometries and tissue property distributions.

SIGNIFICANCE

This study demonstrates that our microwave tomography is capable of recovering accurate tissue property distributions when spatial information from MRI is incorporated through soft prior regularization.