Noninvasive electrical conductivity measurement by MRI: a test of its validity and the electrical conductivity characteristics of glioma

Fast bias-corrected conductivity mapping using stimulated echoes

OBJECTIVE

To demonstrate the potential of a double angle stimulated echo (DA-STE) method for fast and accurate "full" homogeneous Helmholtz-based electrical properties tomography using a simultaneous B 1 + magnitude and transceive phase measurement.

METHODS

The combination of a spin and stimulated echo can be used to yield an estimate of both B 1 + magnitude and transceive phase and thus provides the means for "full" EPT reconstruction. An interleaved 2D acquisition scheme is used for rapid acquisition. The method was validated in a saline phantom and compared to a double angle method based on two single gradient echo acquisitions (GRE-DAM). The method was evaluated in the brain of a healthy volunteer.

RESULTS

The B 1 + magnitude obtained with DA-STE showed excellent agreement with the GRE-DAM method. Conductivity values based on the "full" EPT reconstruction also agreed well with the expectations in the saline phantom. In the brain, the method delivered conductivity values close to literature values.

DISCUSSION

The method allows the use of the "full" Helmholtz-based EPT reconstruction without the need for additional measurements. As a result, quantitative conductivity values are improved compared to phase-based EPT reconstructions. DA-STE is a fast complex- B 1 + mapping technique that could render EPT clinically relevant at 3 T.

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.

Bias correction for phase-based cr-MREPT using low resolution B1+ magnitude

ObjectiveFull-form Magnetic Resonance Electrical Properties Tomography (MREPT) requires bothB1+magnitude and phase information. SinceB1+phase can be obtained faster and with higher SNR compared toB1+magnitude, several phase-based methods have been developed for conductivity imaging. However, phase-based methods suffer from a concave bias due to the assumption that∇|B1+|is negligible in the ROI.ApproachIn this paper, we re-derive the central equation of phase-based cr-MREPT without assuming that∇|B1+|is negligible and thus propose a correction method directly integrated into the equation system.Main resultsProposed method successfully corrects the concave bias on both simulated and experimental data and significantly increases image quality.SignificanceThe proposed correction method depends on a very low-resolution|B1+|map, and therefore the imaging time does not increase significantly for obtainingB1+magnitude. Moreover, correction can be achieved using simulatedB1+magnitude, hence completely removing the additional imaging requirement.

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.

Feasibility of undersampled spiral trajectories in MREPT for fast conductivity imaging

PURPOSE

To investigate spiral-based imaging including trajectories with undersampling as a fast and robust alternative for phase-based magnetic resonance electrical properties tomography (MREPT) techniques.

METHODS

Spiral trajectories with various undersampling ratios were prescribed to acquire images from an experimental phantom and a healthy volunteer at 3T. The non-Cartesian acquisitions were reconstructed using SPIRiT, and conductivity maps were derived using phase-based cr-MREPT. The resulting maps were compared between different sampling trajectories. Additionally, a conductivity map was obtained using a Cartesian balanced SSFP acquisition from the volunteer to comparatively demonstrate the robustness of the proposed method.

RESULTS

The phantom and volunteer results illustrate the benefits of the spiral acquisitions. Specifically, undersampled spiral acquisitions display improved robustness against field inhomogeneity artifacts and lowered SD values with shortened readout times. Furthermore, average of conductivity values measured for the cerebrospinal fluid with the spiral acquisitions were 1.703 S/m, indicating a close agreement with the theoretical values of 1.794 S/m.

CONCLUSION

A spiral-based acquisition framework for conductivity imaging with and without undersampling is presented. Overall, spiral-based acquisitions improved robustness against field inhomogeneity artifacts, while achieving whole head coverage with multiple averages in less than a minute.

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.

An Adaptative Savitzky-Golay Kernel for Laplacian Estimation in Magnetic Resonance Electrical Property Tomography. ANNUAL INTERNATIONAL CONFERENCE OF THE IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. ANNUAL INTERNATIONAL CONFERENCE 2023;2023:1-4. [PMID: 38083553 DOI: 10.1109/embc40787.2023.10341200]

Magnetic Resonance electrical property tomography (MR-EPT) is a non-invasive imaging modality that reconstructs the living biological tissue's conductivity σ and εr permittivity using spatial derivatives of the measured RF field, also termed B1 data, in a magnetic resonance imaging system. The spatial derivative operator, particularly the Laplacian, amplifies the noise in the reconstructed electrical property (EP) maps, hence decreasing accuracy and increasing boundary artifacts. We propose a novel adaptative convolution kernel for generating numerical derivatives based on 3D Savitzky-Golay (SG) filters and local segmentation in a magnitude image. In comparison to typical SG kernel, the proposed kernel allows arbitrary shapes and sizes to vary with local tissue. It provides an automatic trade-off between noise and resolution, thereby significantly enhancing reconstruction accuracy and eliminating boundary artifacts.

REC-NN: A reconstruction error compensation neural network for Magnetic Resonance Electrical Property Tomography (MREPT)

Electrical properties (EPs) are expected as biomarkers for early cancer detection. Magnetic resonance electrical properties tomography (MREPT) is a technique to non-invasively estimate the EPs of tissues from MRI measurements. While noise sensitivity and artifact problems of MREPT are being solved progressively through recent efforts, the loss of tissue contrast emerges as an obstacle to the clinical applications of MREPT. To solve the problem, we propose a reconstruction error compensation neural network scheme (REC-NN) for a typical analytic MREPT method, Stab-EPT. Two NN structures: one with only ResNet blocks, and the other hybridizing ResNet blocks with an encoder-decoder structure. Results of experiments with digital brain phantoms show that, compared with Stab-EPT, and conventional NN based reconstruction, REC-NN improves both reconstruction accuracy and tissue contrast. It is found that, the encoder-decoder structure could improve the compensation accuracy of EPs in homogeneous region but showed worse reconstruction than only ResNet structure for tumorous tissues unseen in the training samples. Future research is required to address overcompensation problems, optimization of NN structure and application to clinical data.

Quantitative Assessment of Intervertebral Disc Composition by MRI: Sensitivity to Diurnal Variation

Whether diurnal variation exists in quantitative MRI indices such as the T1rho relaxation time (T1ρ) of the intervertebral disc (IVD) is yet to be explored. This prospective study aimed to evaluate the diurnal variation in T1ρ, apparent diffusion coefficient (ADC), and electrical conductivity (σ) of lumbar IVD and its relationship with other MRI or clinical indices. Lumbar spine MRI, including T1ρ imaging, diffusion-weighted imaging (DWI), and electric properties tomography (EPT), was conducted on 17 sedentary workers twice (morning and evening) on the same day. The T1ρ, ADC, and σ of IVD were compared between the time points. Their diurnal variation, if any, was tested for correlation with age, body mass index (BMI), IVD level, Pfirrmann grade, scan interval, and diurnal variation in IVD height index. The results showed a significant decrease in T1ρ and ADC and a significant increase in the σ of IVD in the evening. T1ρ variation had a weak correlation with age and scan interval, and ADC variation with scan interval. Diurnal variation exists for the T1ρ, ADC, and σ of lumbar IVD, which should be accounted for in image interpretation. This variation is thought to be due to diurnal variations in intradiscal water, proteoglycan, and sodium ion concentration.

Use of machine learning to improve the estimation of conductivity and permittivity based on longitudinal relaxation time T1 in magnetic resonance at 7 T

Electrical property tomography (EPT) is a noninvasive method that uses magnetic resonance imaging (MRI) to estimate the conductivity and permittivity of tissues, and hence, can be used as a biomarker. One branch of EPT is based on the correlation of water and relaxation time T1 with the conductivity and permittivity of tissues. This correlation was applied to a curve-fitting function to estimate electrical properties, it was found to have a high correlation between permittivity and T1 however the computation of conductivity based on T1 requires to estimate the water content. In this study, we developed multiple phantoms with several ingredients that modify the conductivity and permittivity and explored the use of machine learning algorithms to have a direct estimation of conductivity and permittivity based on MR images and the relaxation time T1. To train the algorithms, each phantom was measured using a dielectric measurement device to acquire the true conductivity and permittivity. MR images were taken for each phantom, and the T1 values were measured. Then, the acquired data were tested using curve fitting, regression learning, and neural fit models to estimate the conductivity and permittivity values based on the T1 values. In particular, the regression learning algorithm based on Gaussian process regression showed high accuracy with a coefficient of determination R² of 0.96 and 0.99 for permittivity and conductivity, respectively. The estimation of permittivity using regression learning demonstrated a lower mean error of 0.66% compared to the curve fitting method, which resulted in a mean error of 3.6%. The estimation of conductivity also showed that the regression learning approach had a lower mean error of 0.49%, whereas the curve fitting method resulted in a mean error of 6%. The findings suggest that utilizing regression learning models, specifically Gaussian process regression, can result in more accurate estimations for both permittivity and conductivity compared to other methods.

Technical note: Multi-receiver combination method for phase-based electrical property tomography of the breast

BACKGROUND

Phase-based electrical property tomography (EPT) is a technique that allows conductivity reconstruction with only phase of the B₁ field under the assumption that the magnitude of the B₁ fields are homogeneous. The more this assumption is violated, the less accurate the reconstructed conductivity. Thus, a method that ensures homogeneity of | B 1 - | $| {{\rm{B}}_1^ - } |$ field is important for breast image using multi-receiver coil.

PURPOSE

To develop a method for multi-receiver combination for phase-based EPT usable for breast EPT imaging in the clinic.

METHODS

Theory of the proposed method is presented. To validate the proposed method, the phantom and in-vivo experiments were conducted. Conductivity images were reconstructed using the transceive phase of the combined image and results were compared with another combination method.

RESULTS

The proposed method's conductivity results were more stable than those of the previous method when | B 1 + | $| {{\rm{B}}_1^ + } |$ was not homogeneous and when the homogeneous contrast region was small. The phantom and in-vivo results indicate that the proposed method produces improved conductivity images than the previous method. The proposed combination method also increased the conductivity contrast between benign and cancerous tissues.

CONCLUSION

The proposed method produced more stable multi-receiver combination for phase-based EPT of the breast in a clinical environment.

High frequency conductivity decomposition by solving physically constraint underdetermined inverse problem in human brain

The developed magnetic resonance electrical properties tomography (MREPT) can visualize the internal conductivity distribution at Larmor frequency by measuring the B1 transceive phase data from magnetic resonance imaging (MRI). The recovered high-frequency conductivity (HFC) value is highly complex and heterogeneous in a macroscopic imaging voxel. Using high and low b-value diffusion weighted imaging (DWI) data, the multi-compartment spherical mean technique (MC-SMT) characterizes the water molecule movement within and between intra- and extra-neurite compartments by analyzing the microstructures and underlying architectural organization of brain tissues. The proposed method decomposes the recovered HFC into the conductivity values in the intra- and extra-neurite compartments via the recovered intra-neurite volume fraction (IVF) and the diffusion patterns using DWI data. As a form of decomposition of intra- and extra-neurite compartments, the problem to determine the intra- and extra-neurite conductivity values from the HFC is still an underdetermined inverse problem. To solve the underdetermined problem, we use the compartmentalized IVF as a criterion to decompose the electrical properties because the ion-concentration and mobility have different characteristics in the intra- and extra-neurite compartments. The proposed method determines a representative apparent intra- and extra-neurite conductivity values by changing the underdetermined equation for a voxel into an over-determined minimization problem over a local window consisting of surrounding voxels. To suppress the noise amplification and estimate a feasible conductivity, we define a diffusion pattern distance to weight the over-determined system in the local window. To quantify the proposed method, we conducted a simulation experiment. The simulation experiments show the relationships between the noise reduction and the spatial resolution depending on the designed local window sizes and diffusion pattern distance. Human brain experiments (five young healthy volunteers and a patient with brain tumor) were conducted to evaluate and validate the reliability of the proposed method. To quantitatively compare the results with previously developed methods, we analyzed the errors for reconstructed extra-neurite conductivity using existing methods and indirectly verified the feasibility of the proposed method.

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.

Effects of Irradiation on Brain Tumors Using MR-Based Electrical Conductivity Imaging

Ionizing radiation delivers sufficient energy inside the human body to create ions, which kills cancerous tissues either by damaging the DNA directly or by creating charged particles that can damage the DNA. Recent magnetic resonance (MR)-based conductivity imaging shows higher sensitivity than other MR techniques for evaluating the responses of normal tissues immediately after irradiation. However, it is still necessary to verify the responses of cancer tissues to irradiation by conductivity imaging for it to become a reliable tool in evaluating therapeutic effects in clinical practice. In this study, we applied MR-based conductivity imaging to mouse brain tumors to evaluate the responses in irradiated and non-irradiated tissues during the peri-irradiation period. Absolute conductivities of brain tissues were measured to quantify the irradiation effects, and the percentage changes were determined to estimate the degree of response. The conductivity of brain tissues with irradiation was higher than that without irradiation for all tissue types. The percentage changes of tumor tissues with irradiation were clearly different than those without irradiation. The measured conductivity and percentage changes between tumor rims and cores to irradiation were clearly distinguished. The contrast of the conductivity images following irradiation may reflect the response to the changes in cellularity and the amounts of electrolytes in tumor tissues.

The Effect of Freezing and Thawing on Complex Permittivity of Bovine Tissues

The aim of this study was to investigate how the freezing and thawing of biological tissues affect their complex permittivity in the microwave frequency range from 0.5 MHz to 18 GHz. We measured the complex permittivity of ex vivo bovine tissues, including brain white and grey matter, liver, and muscle, using an open-ended coaxial probe. Bovine tissues were chosen for their availability and similarity to human tissue permittivity. The samples were measured at 25 °C, before they were frozen either in a commercial freezer below -18 °C or in liquid nitrogen, nominally at -196 °C. The measured permittivity before freezing was compared to the permittivity measured after freezing and thawing the tissues back to 25 °C. Statistical analysis of the results showed a statistically significant change in permittivity after freezing and thawing by both methods for all the measured tissues, at least in some parts of the measured frequency range. The largest difference was observed for the white matter, while the liver had the smallest percent change.

Physics Informed Neural Networks (PINN) for Low Snr Magnetic Resonance Electrical Properties Tomography (MREPT)

Electrical properties (EPs) of tissues facilitate early detection of cancerous tissues. Magnetic resonance electrical properties tomography (MREPT) is a technique to non-invasively probe the EPs of tissues from MRI measurements. Most MREPT methods rely on numerical differentiation (ND) to solve partial differential Equations (PDEs) to reconstruct the EPs. However, they are not practical for clinical data because ND is noise sensitive and the MRI measurements for MREPT are noisy in nature. Recently, Physics informed neural networks (PINNs) have been introduced to solve PDEs by substituting ND with automatic differentiation (AD). To the best of our knowledge, it has not been applied to MREPT due to the challenges in using PINN on MREPT as (i) a PINN requires part of ground-truth EPs as collocation points to optimize the network's AD, (ii) the noisy input data disrupts the optimization of PINNs despite the noise-filtering nature of NNs and additional denoising processes. In this work, we propose a PINN-MREPT model based on a canonical analytic MREPT model. A reference padding layer with known EPs was added to surround the region of interest for providing additive collocation points. Moreover, an optimizable diffusion coefficient was embedded in the analytic MREPT model used in the PINN-MREPT. The noise robustness of the proposed PINN-MREPT for single-sample reconstruction was tested by using numerical phantoms of human brain with extra tumor-like tissues at different noise levels. The results of numerical experiments show that PINN-MREPT outperforms two typical numerical MREPT methods in terms of reconstruction accuracy, sensitivity to the extra tissues, and the correlations of line profiles in the regions of interest. The advantage of the PINN-MREPT is shown by the results of an experiment on phantom measurement, too. Moreover, it is found that the diffusion term plays an important role to achieve a noise-robust PINN-MREPT. This is an important step moving forward to a clinical application of MREPT.

Application of High-Frequency Conductivity Map Using MRI to Evaluate It in the Brain of Alzheimer's Disease Patients

Background

The previous studies reported increased concentrations of metallic ions, imbalanced Na+ and K+ ions, and the increased mobility of protons by microstructural disruptions in Alzheimer's disease (AD).

Purpose

(1) to apply a high-frequency conductivity (HFC) mapping technique using a clinical 3T MRI system, (2) compare HFC values in the brains of participants with AD, amnestic mild cognitive impairment (MCI), and cognitively normal (CN) elderly people, (3) evaluate the relationship between HFC values and cognitive decline, and (4) explore usefulness of HFC values as an imaging biomarker to evaluate the differentiation of AD from CN.

Materials and Methods

This prospective study included 74 participants (23 AD patients, 27 amnestic MCI patients, and 24 CN elderly people) to explore the clinical application of HFC mapping in the brain from March 2019 to August 2021. We performed statistical analyses to compare HFC maps between the three participant groups, evaluate the association of HFC maps with Mini-Mental State Examination (MMSE) scores, and to evaluate the differentiation between the participant groups for HFC values for some brain areas.

Results

We obtained a good HFC map non-invasively. The HFC value was higher in the AD group than in the CN and MCI groups. MMSE scores were negatively associated with HFC values. Age was positively associated with HFC values. The HFC value in the insula has a high area under the receiver operating characteristic (ROC) curve (AUC) value to differentiate AD patients from the CN participants (Sensitivity [SE] = 82, Specificity [SP] =97, AUC = 0.902, p < 0.0001), better than gray matter volume (GMV) in hippocampus (SE = 79, SP = 83, AUC = 0.880, p < 0.0001). The classification for differentiating AD from CN was highest by adding the hippocampal GMV to the insular HFC value (SE = 87, SP = 87, AUC = 0.928, p < 0.0001).

Conclusion

High-frequency conductivity values were significantly increased in the AD group compared to the CN group and increased with age and disease severity. HFC values of the insula along with the GMV of the hippocampus can be used as an imaging biomarker to improve the differentiation of AD from CN.

Physics-Coupled Neural Network Magnetic Resonance Electrical Property Tomography (MREPT) for Conductivity Reconstruction

The electrical property (EP) of human tissues is a quantitative biomarker that facilitates early diagnosis of cancerous tissues. Magnetic resonance electrical properties tomography (MREPT) is an imaging modality that reconstructs EPs by the radio-frequency field in an MRI system. MREPT reconstructs EPs by solving analytic models numerically based on Maxwell's equations. Most MREPT methods suffer from artifacts caused by inaccuracy of the hypotheses behind the models, and/or numerical errors. These artifacts can be mitigated by adding coefficients to stabilize the models, however, the selection of such coefficient has been empirical, which limit its medical application. Alternatively, end-to-end Neural networks-based MREPT (NN-MREPT) learns to reconstruct the EPs from training samples, circumventing Maxwell's equations. However, due to its pattern-matching nature, it is difficult for NN-MREPT to produce accurate reconstructions for new samples. In this work, we proposed a physics-coupled NN for MREPT (PCNN-MREPT), in which an analytic model, cr-MREPT, works with diffusion and convection coefficients, learned by NNs from the difference between the reconstructed and ground-truth EPs to reduce artifacts. With two simulated datasets, three generalization experiments in which test samples deviate gradually from the training samples, and one noise-robustness experiment were conducted. The results show that the proposed PCNN-MREPT achieves higher accuracy than two representative analytic methods. Moreover, compared with an end-to-end NN-MREPT, the proposed method attained higher accuracy in two critical generalization tests. This is an important step to practical MREPT medical diagnoses.

In vivo electrical conductivity measurement of muscle, cartilage, and peripheral nerve around knee joint using MR-electrical properties tomography

This study aimed to investigate whether in vivo MR-electrical properties tomography (MR-EPT) is feasible in musculoskeletal tissues by evaluating the conductivity of muscle, cartilage, and peripheral nerve around the knee joint, and to explore whether these measurements change after exercise. This prospective study was approved by the institutional review board. On February 2020, ten healthy volunteers provided written informed consent and underwent MRI of the right knee using a three-dimensional balanced steady-state free precession (bSSFP) sequence. To test the effect of loading, the subjects performed 60 squatting exercises after baseline MRI, immediately followed by post-exercise MRI with the same sequences. After reconstruction of conductivity map based on the bSSFP sequence, conductivity of muscles, cartilages, and nerves were measured. Measurements between the baseline and post-exercise MRI were compared using the paired t-test. Test–retest reliability for baseline conductivity was evaluated using the intraclass correlation coefficient. The baseline and post-exercise conductivity values (mean ± standard deviation) [S/m] of muscles, cartilages, and nerves were 1.73 ± 0.40 and 1.82 ± 0.50 (p = 0.048), 2.29 ± 0.47 and 2.51 ± 0.37 (p = 0.006), and 2.35 ± 0.57 and 2.36 ± 0.57 (p = 0.927), respectively. Intraclass correlation coefficient for the baseline conductivity of muscles, cartilages, and nerves were 0.89, 0.67, and 0.89, respectively. In conclusion, in vivo conductivity measurement of musculoskeletal tissues is feasible using MR-EPT. Conductivity of muscles and cartilages significantly changed with an overall increase after exercise.

Brain Tissue Conductivity in Focal Cerebral Ischemia

BACKGROUND

Cerebral ischemia leads to oxygen depletion with rapid breakdown of transmembrane transporters and subsequent impaired electrolyte haemostasis. Electric properties tomography (EPT) is a new contrast in MRI which delivers information on tissue electrical conductivity. In the clinical realm it has been mostly used for tumour mapping. Ischemic cerebral stroke is another promising but neglected application. It might deliver additional information on tissue viability and possible response to therapy.

AIM

The aim of this study was to demonstrate tissue conductivity in a rodent model of stroke. Further, we aimed to compare electric conductivity in ischemic and non-ischemic cerebral tissue.

MATERIALS AND METHODS

Two male Wistar rats were used in this study and were subjected to permanent MCAO. The animals were scanned in a 3 Tesla system (Philips Achieva/Best, the Netherlands) using a dedicated solenoid animal coil (Philips/Hamburg, Germany). In addition to diffusion weighted imaging (DWI), EPT was performed using a steady-state free-precession (SSFP) sequence (repetition time/echo time = 4.5/2.3 ms, measured voxel size = 0.6 × 0.6 × 1.2 mm3, flip angle = 38°, number of excitations = 4). From the transceive phase ϕ of these SSFP scans, conductivity σ was estimated by the equation σ = Δϕ/(2μ0ω) with Δ the Laplacian operator, μ0 the magnetic permeability, and ω the Larmor frequency. Subsequently, a median filter was applied, which was locally restricted to voxels with comparable signal magnitude.

RESULTS

The animals exhibited an infarct as demonstrated on DWI. Conductivity within the infarcted region was 60-70 % of the conductivity of not affected contralateral tissue (0.39 ± 0.07 S/m and 0.31 ± 0.14 S/m vs. 0.64 ± 0.15 S/m and 0.66 ± 0.16 S/m, respectively).

DISCUSSION

Infarcted tissue exhibited decreased conductivity. Further in-vivo experiments with examination of the influence of reperfusion status and temporal evolution of the infarcted areas should be conducted. Depiction of the ischemic penumbra and possibly subclassification of the DWI lesion still seems to be a fruitful target for further studies.

Mapping electric bulk conductivity in the human heart

PURPOSE

To explore the technical feasibility of mapping the electric bulk conductivity in the human heart, and to determine quantitative conductivity values of myocardium and blood from a small group of volunteers.

METHODS

Using a 3T MR system, 6 healthy male volunteers were measured. For all volunteers, a time-resolved 2D sequence over the cardiac cycle was applied (electrocardiogram [ECG]-triggered SSFP acquired in breath-hold). From these data, a dedicated, so-called "2D conductivity" has been derived in the framework of electrical properties tomography (EPT). To validate the concept of 2D conductivity, a static 3D sequence (ECG-triggered and respiratory-gated SSFP 3D whole heart acquisition, allowing the full 3D reconstruction of conductivity) as well as a Q-flow sequence (for investigating the relation between flow and reconstruction errors of the conductivity) have been applied for one of the volunteers.

RESULTS

For both, blood and myocardium, quantitative values of obtained 2D conductivity were approximately two-thirds of the obtained 3D conductivity, as expected from Maxwell's equations. Furthermore, the quantitative conductivity values agreed with corresponding literature values. Conductivity of left-ventricular blood volume showed characteristic over- and under-shooting at specific time points during the cardiac cycle for all volunteers investigated. This over- and under-shooting correlated with the phase pattern caused by blood flow into/out of the ventricle.

CONCLUSION

The study demonstrated the technical feasibility of cardiac conductivity measurements using standard MR systems and standard MR sequences, and therefore, may open new options for MR-based cardiac diagnosis.

Comparison of Five Conductivity Tensor Models and Image Reconstruction Methods Using MRI

Imaging of the electrical conductivity distribution inside the human body has been investigated for numerous clinical applications. The conductivity tensors of biological tissue have been obtained from water diffusion tensors by applying several models, which may not cover the entire phenomenon. Recently, a new conductivity tensor imaging (CTI) method was developed through a combination of B1 mapping, and multi-b diffusion weighted imaging. In this study, we compared the most recent CTI method with the four existing models of conductivity tensors reconstruction. Two conductivity phantoms were designed to evaluate the accuracy of the models. Applied to five human brains, the conductivity tensors using the four existing models and CTI were imaged and compared with the values from the literature. The conductivity image of the phantoms by the CTI method showed relative errors between 1.10% and 5.26%. The images by the four models using DTI could not measure the effects of different ion concentrations subsequently due to prior information of the mean conductivity values. The conductivity tensor images obtained from five human brains through the CTI method were comparable to previously reported literature values. The images by the four methods using DTI were highly correlated with the diffusion tensor images, showing a coefficient of determination (R2) value of 0.65 to 1.00. However, the images by the CTI method were less correlated with the diffusion tensor images and exhibited an averaged R2 value of 0.51. The CTI method could handle the effects of different ion concentrations as well as mobilities and extracellular volume fractions by collecting and processing additional B1 map data. It is necessary to select an application-specific model taking into account the pros and cons of each model. Future studies are essential to confirm the usefulness of these conductivity tensor imaging methods in clinical applications, such as tumor characterization, EEG source imaging, and treatment planning for electrical stimulation.

High-frequency electrical properties tomography at 9.4T as a novel contrast mechanism for brain tumors

PURPOSE

To establish high-frequency magnetic resonance electrical properties tomography (MREPT) as a novel contrast mechanism for the assessment of glioblastomas using a rat brain tumor model.

METHODS

Six F98 intracranial tumor bearing rats were imaged longitudinally 8, 11 and 14 days after tumor cell inoculation. Conductivity and mean diffusivity maps were generated using MREPT and Diffusion Tensor Imaging. These maps were co-registered with T2 -weighted images and volumes of interests (VOIs) were segmented from the normal brain, ventricles, edema, viable tumor, tumor rim, and tumor core regions. Longitudinal changes in conductivity and mean diffusivity (MD) values were compared in these regions. A correlation analysis was also performed between conductivity and mean diffusivity values.

RESULTS

The conductivity of ventricles, edematous area and tumor regions (tumor rim, viable tumor, tumor core) was significantly higher (P < .01) compared to the contralateral cortex. The conductivity of the tumor increased over time while MD from the tumor did not change. A marginal positive correlation was noted between conductivity and MD values for tumor rim and viable tumor, whereas this correlation was negative for the tumor core.

CONCLUSION

We demonstrate a novel contrast mechanism based on ionic concentration and mobility, which may aid in providing complementary information to water diffusion in probing the microenvironment of brain tumors.

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.

Low conductivity on electrical properties tomography demonstrates unique tumor habitats indicating progression in glioblastoma

OBJECTIVES

Tissue conductivity measurements made with electrical properties tomography (EPT) can be used to define temporal changes in tissue habitats on longitudinal multiparametric MRI. We aimed to demonstrate the added insights for identifying tumor habitats obtained by including EPT with diffusion- and perfusion-weighted MRI, and to evaluate the use of these tumor habitats for determining tumor treatment response in post-treatment glioblastoma.

METHODS

Tumor habitats were developed from EPT, diffusion-weighted, and perfusion-weighted MRI in 60 patients with glioblastoma who underwent concurrent chemoradiotherapy. Voxels from EPT, apparent diffusion coefficient (ADC), and cerebral blood volume (CBV) maps were clustered into habitats, and each habitat was serially examined to assess its temporal change. The usefulness of temporal changes in tumor habitats for diagnosing tumor progression and treatment-related change was investigated using logistic regression. The performance of significant predictors was measured using the area under the curve (AUC) from receiver-operating-characteristics analysis with 1000-fold bootstrapping.

RESULTS

Five tumor habitats were identified, and of these, the hypervascular cellular habitat (odds ratio [OR] 5.45; 95% CI, 1.75-31.42; p = .02), hypovascular low conductivity habitat (OR 2.00; 95% CI, 1.45-3.05; p < .001), and hypovascular intermediate habitat (OR 1.57; 95% CI, 1.18-2.30; p = .006) were predictive of tumor progression. Low EPT and low CBV reflected a unique hypovascular low conductivity habitat that showed the highest diagnostic performance (AUC 0.86; 95% CI, 0.76-0.96). The combined habitats showed high performance (AUC 0.90; 95% CI, 0.82-0.98) in the differentiation of tumor progression from treatment-related change.

CONCLUSION

EPT reveals low conductivity habitats that can improve the diagnosis of tumor progression in post-treatment glioblastoma.

KEY POINTS

• Electrical properties tomography (EPT) demonstrated lower conductivity in tumor progression than in treatment-related change. • EPT allowed identification of a unique hypovascular low conductivity habitat when combined with cerebral blood volume mapping. • Tumor habitats with a hypovascular low conductivity habitat, hypervascular cellular habitat, and hypovascular intermediate habitat yielded high diagnostic performance for diagnosing tumor progression.

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.

Brain Tissue Conductivity Measurements with MR-Electrical Properties Tomography: An In Vivo Study

First in vivo brain conductivity reconstructions using Helmholtz MR-Electrical Properties Tomography (MR-EPT) have been published. However, a large variation in the reconstructed conductivity values is reported and these values differ from ex vivo conductivity measurements. Given this lack of agreement, we performed an in vivo study on eight healthy subjects to provide reference in vivo brain conductivity values. MR-EPT reconstructions were performed at 3 T for eight healthy subjects. Mean conductivity and standard deviation values in the white matter, gray matter and cerebrospinal fluid (σWM, σGM, and σCSF) were computed for each subject before and after erosion of regions at tissue boundaries, which are affected by typical MR-EPT reconstruction errors. The obtained values were compared to the reported ex vivo literature values. To benchmark the accuracy of in vivo conductivity reconstructions, the same pipeline was applied to simulated data, which allow knowledge of ground truth conductivity. Provided sufficient boundary erosion, the in vivo σWM and σGM values obtained in this study agree for the first time with literature values measured ex vivo. This could not be verified for the CSF due to its limited spatial extension. Conductivity reconstructions from simulated data verified conductivity reconstructions from in vivo data and demonstrated the importance of discarding voxels at tissue boundaries. The presented σWM and σGM values can therefore be used for comparison in future studies employing different MR-EPT techniques.

Impact of Peritumoral Edema During Tumor Treatment Field Therapy: A Computational Modelling Study

BACKGROUND

Tumor treatment fields (TTFie-lds) are an approved adjuvant therapy for glioblastoma (GBM). The magnitude of applied electrical field has been shown to be related to the anti-tumoral response. However, peritumoral edema may result in shunting of electrical current around the tumor, thereby reducing the intra-tumoral electric field. In this study, we systematically address this issue with computational simulations.

METHODS

Finite element models are created of a human head with varying amounts of peritumoral edema surrounding a virtual tumor. The electric field distribution was simulated using the standard TTFields electrode montage. Electric field magnitude was extracted from the tumor and related to edema thickness. Two patient specific models were created to confirm these results.

RESULTS

The inclusion of peritumoral edema decreased the average magnitude of the electric field within the tumor. In the model considering a frontal tumor and an anterior-posterior electrode configuration, ≥6 mm of peritumoral edema decreased the electric field by 52%. In the patient specific models, peritumoral edema decreased the electric field magnitude within the tumor by an average of 26%. The effect of peritumoral edema on the electric field distribution was spatially heterogenous, being most significant at the tissue interface between edema and tumor.

CONCLUSIONS

The inclusion of peritumoral edema during TTFields modelling may have a dramatic effect on the predicted electric field magnitude within the tumor. Given the importance of electric field magnitude for the anti-tumoral effects of TTFields, the presence of edema should be considered both in future modelling studies and when planning TTField therapy.

Low-frequency dominant electrical conductivity imaging of in vivo human brain using high-frequency conductivity at Larmor-frequency and spherical mean diffusivity without external injection current

Diffusion weighted imaging based on random Brownian motion of water molecules within a voxel provides information on the micro-structure of biological tissues through water molecule diffusivity. As the electrical conductivity is primarily determined by the concentration and mobility of ionic charge carriers, the macroscopic electrical conductivity of biological tissues is also related to the diffusion of electrical ions. This paper aims to investigate the low-frequency electrical conductivity by relying on a pre-defined biological model that separates the brain into the intracellular (restricted) and extracellular (hindered) compartments. The proposed method uses B1 mapping technique, which provides a high-frequency conductivity distribution at Larmor frequency, and the spherical mean technique, which directly estimates the microscopic tissue structure based on the water molecule diffusivity and neurite orientation distribution. The total extracellular ion concentration, which is separated from the high-frequency conductivity, is recovered using the estimated diffusivity parameters and volume fraction in each compartment. We propose a method to reconstruct the low-frequency dominant conductivity tensor by taking into consideration the extracted extracellular diffusion tensor map and the reconstructed electrical parameters. To demonstrate the reliability of the proposed method, we conducted two phantom experiments. The first one used a cylindrical acrylic cage filled with an agar in the background region and four anomalies for the effect of ion concentration on the electrical conductivity. The other experiment, in which the effect of ion mobility on the conductivity was verified, used cell-like materials with thin insulating membranes suspended in an electrolyte. Animal and human brain experiments were conducted to visualize the low-frequency dominant conductivity tensor images. The proposed method using a conventional MRI scanner can predict the internal current density map in the brain without directly injected external currents.

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.

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.

A robust electrical conductivity imaging method with total variation and wavelet regularization

PURPOSE

This study aims to develop and evaluate a robust conductivity imaging method that combines total variation and wavelet regularization to enhance the accuracy of conductivity maps.

THEORY AND METHODS

The proposed approach is based on a gradient-based method. The central equation is derived from Maxwell's equation and describes the relationship between conductivity and the transceive phase. A linear system equation is obtained via a finite-difference method and solved using a least-squares method. Total variation and wavelet transform regularization terms are added to the minimization problem and solved using the Split Bregman method to improve reconstruction stability. The proposed approach is compared with conventional and gradient-based methods. Numerical simulations are performed to validate the accuracy of the developed method, and the effects of noise are determined. Phantom and in vivo experiments are conducted at 3 T to verify the clinical applicability of the proposed method.

RESULTS

Numerical simulations show that the proposed method is more robust than other methods and can suppress the effects of noise. The quantitative conductivity value of the phantom experiment agrees with the measured value. The in vivo experiment results present a clear structure, and the conductivity value of the tumor region is significantly higher than that around healthy tissues.

CONCLUSION

The proposed electrical conductivity imaging method can improve the quality of conductivity reconstruction, and thus, has future clinical applications.

Magnetic resonance electrical property mapping at 21.1 T: a study of conductivity and permittivity in phantoms, ex vivo tissue and in vivo ischemia

Electrical properties (EP), namely conductivity and permittivity, can provide endogenous contrast for tissue characterization. Using electrical property tomography (EPT), maps of EP can be generated from conventional MRI data. This report investigates the feasibility and accuracy of EPT at 21.1 T for multiple RF coils and modes of operation using phantoms. Additionally, it demonstrates the EP of the in vivo rat brain with and without ischemia. Helmholtz-based EPT was implemented in its Full-form, which demands the complex [Formula: see text] field, and a simplified form requiring either just the [Formula: see text] field phase for conductivity or the [Formula: see text] field magnitude for permittivity. Experiments were conducted at 21.1 T using birdcage and saddle coils operated in linear or quadrature transceive mode, respectively. EPT approaches were evaluated using a phantom, ex and in vivo Sprague-Dawley rats under naïve conditions and ischemic stroke via transient middle cerebral artery occlusion. Different conductivity reconstruction approaches applied to the phantom displayed average errors of 12%-73% to the target acquired from dielectric probe measurements. Permittivity reconstructions showed higher agreement and an average 3%-8% error to the target depending on reconstruction approach. Conductivity and permittivity of ex and in vivo rodent brain were measured. Elevated EP in the ischemia region correlated with the increased sodium content and the influx of water intracellularly following ischemia in the lesion were detected. The Full-form technique generated from the linear birdcage provided the best accuracy for EP of the phantom. Phase-based conductivity and magnitude-based permittivity mapping provided reasonable estimates but also demonstrated the limitations of Helmholtz-based EPT at 21.1 T. Permittivity reconstruction was improved significantly over lower fields, suggesting a novel metric for in vivo brain studies. EPT applied to ischemic rat brain proved sensitivity to physiological changes, motivating the future application of more advanced reconstruction approaches.

In Vivo Measurement of Brain Tissue Response After Irradiation: Comparison of T2 Relaxation, Apparent Diffusion Coefficient, and Electrical Conductivity

Radiation therapy (RT) has been widely used as a powerful treatment tool to address cancerous tissues because of its ability to control cell growth. Its ionizing radiation damages the DNA of cancerous tissues, leading to cell death. Medical imaging, however, still has limitations regarding the reliability of its assessment of tissue response and in predicting the treatment effect because of its inability to provide contrast information on the gradual, minute tissue changes after RT. A recently developed magnetic resonance (MR)-based conductivity imaging method may provide direct, highly sensitive information on this tissue response because its contrast mechanism is based on the concentration and mobility of ions in intracellular and extracellular spaces. In this feasibility study, we applied T2-weighted, diffusion-weighted, and electrical conductivity imaging to mouse brain, thus, using the MR imaging to map the tissue response after radiation exposure. To evaluate the degree of response, we measured the T2 relaxation, apparent diffusion coefficient (ADC), and electrical conductivity of brain tissues before and after irradiation. The conductivity images, which showed significantly higher sensitivity than other MR imaging methods, indicated that the contrast is distinguishable in different ways at different areas of the brain. Future studies will focus on verifying these results and the long-term evaluation of conductivity changes using various irradiation methods for clinical applications.

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.

Automated gradient-based electrical properties tomography in the human brain using 7 Tesla MRI

Electrical properties of the brain tissues may yield useful biomarkers for neurological disorders and diseases, as well as contribute to safety assurance of ultra-high-field MRI. It has been reported that using B₁ maps from a multi-channel RF coil, the spatial variation of the electrical properties can be robustly retrieved. The absolute electrical property values can then be obtained by spatial integration, given that an integration seed point is assigned. In this study, we propose to exploit automatically detected seed points based on tissue piece-wise homogeneity (Helmholtz equation) for spatial integration. Numerical simulations of a numerical brain model and experiments involving 12 healthy volunteers were performed to demonstrate its feasibility and robustness in various noisy conditions and head positions. For in vivo imaging, we consistently observed higher conductivity and permittivity values in the white and gray matter compared to tabulated ex vivo probe measurement results found in the literature, a discrepancy that may be attributed to ex vivo experimental constraints. Our results suggest that the proposed technique produces consistent brain electrical properties in vivo that may contribute to improving diagnostic and therapeutic decisions.

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.

Variation in Reported Human Head Tissue Electrical Conductivity Values

Electromagnetic source characterisation requires accurate volume conductor models representing head geometry and the electrical conductivity field. Head tissue conductivity is often assumed from previous literature, however, despite extensive research, measurements are inconsistent. A meta-analysis of reported human head electrical conductivity values was therefore conducted to determine significant variation and subsequent influential factors. Of 3121 identified publications spanning three databases, 56 papers were included in data extraction. Conductivity values were categorised according to tissue type, and recorded alongside methodology, measurement condition, current frequency, tissue temperature, participant pathology and age. We found variation in electrical conductivity of the whole-skull, the spongiform layer of the skull, isotropic, perpendicularly- and parallelly-oriented white matter (WM) and the brain-to-skull-conductivity ratio (BSCR) could be significantly attributed to a combination of differences in methodology and demographics. This large variation should be acknowledged, and care should be taken when creating volume conductor models, ideally constructing them on an individual basis, rather than assuming them from the literature. When personalised models are unavailable, it is suggested weighted average means from the current meta-analysis are used. Assigning conductivity as: 0.41 S/m for the scalp, 0.02 S/m for the whole skull, or when better modelled as a three-layer skull 0.048 S/m for the spongiform layer, 0.007 S/m for the inner compact and 0.005 S/m for the outer compact, as well as 1.71 S/m for the CSF, 0.47 S/m for the grey matter, 0.22 S/m for WM and 50.4 for the BSCR.

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.