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.

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).

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.

Angular resolution enhancement technique for diffusion-weighted imaging (DWI) using predicted diffusion gradient directions.

Anisotropic diffusion MRI techniques using single-shell or multi-shell acquisitions have been proposed as a means to overcome some limitations imposed by diffusion tensor imaging (DTI), especially in complex models of fibre orientation distribution in voxels. A long acquisition time for the angular resolution of diffusion MRI is a major obstacle to practical clinical implementations. In this paper, we propose a novel method to improve angular resolution of diffusion MRI acquisition using given diffusion gradient (DG) directions. First, we define a local diffusion pattern map of diffusion MR signals on a single shell in given DG directions. Using the local diffusion pattern map, we design a prediction scheme to determine the best DG direction to be synthesized within a nearest neighborhood DG directions group. Second, the local diffusion pattern map and the spherical distance on the shell are combined to determine a synthesized diffusion signal in the new DG direction. Using the synthesized and measured diffusion signals on a single sphere, we estimate a spin orientation distribution function (SDF) with human brain data. Although the proposed method is applied to SDF, a basic idea is to increase the angular resolution using the measured diffusion signals in various DG directions. The method can be applicable to different acquired multi-shell data or diffusion spectroscopic imaging (DSI) data. We validate the proposed method by comparing the recovered SDFs using the angular resolution enhanced diffusion signals with the recovered SDF using the measured diffusion data. The developed method provides an enhanced SDF resolution and improved multiple fiber structure by incorporating synthesized signals. The proposed method was also applied neurite orientation dispersion and density imaging (NODDI) using multi-shell acquisitions.

Microelectrode array analysis of hippocampal network dynamics following theta-burst stimulation via current source density reconstruction by Gaussian interpolation.

BACKGROUND: Multielectrode arrays (MEAs) have been used to understand electrophysiological network dynamics by recording real-time activity in groups of cells. The extent to which the collection of such data enables hypothesis testing on the level of circuits and networks depends largely on the sophistication of the analyses applied. NEW METHOD: We studied the systemic temporal variations of endogenous signaling within an organotypic hippocampal network following theta-burst stimulation (TBS) to the Schaffer collateral-commissural pathways. The recovered current source density (CSD) information from the raw grid of extracellular potentials by using a Gaussian interpolation method increases spatial resolution and avoids border artifacts by numerical differentials. RESULTS: We compared total sink and source currents in DG, CA3, and CA1; calculated accumulated correlation coefficients to compare pre- with post-stimulation CSD dynamics in each region; and reconstructed functional connectivity maps for regional cross-correlations with respect to temporal CSD variations. The functional connectivity maps for potential correlations pre- and post-TBS were compared to investigate the neural network as a whole, revealing differences post-TBS. COMPARISON WITH EXISTING METHOD(S): Previous MEA work on plasticity in hippocampal evoked potentials has focused on synchronicity across the hippocampus within isolated subregions. Such analyses ignore the complex relationships among diverse components of the hippocampal circuitry, thus failing to capture network-level behaviors integral to understanding hippocampal function. CONCLUSIONS: The proposed method of recovering current source density to examine whole-hippocampal function is sensitive to experimental manipulation and is worth further examination in the context of network-level analyses of neural signaling.

In vivo electrical conductivity imaging of a canine brain using a 3 T MREIT system.

Magnetic resonance electrical impedance tomography (MREIT) aims at producing high-resolution cross-sectional conductivity images of an electrically conducting object such as the human body. Following numerous phantom imaging experiments, the most recent study demonstrated successful conductivity image reconstructions of postmortem canine brains using a 3 T MREIT system with 40 mA imaging currents. Here, we report the results of in vivo animal imaging experiments using 5 mA imaging currents. To investigate any change of electrical conductivity due to brain ischemia, canine brains having a regional ischemic model were scanned along with separate scans of canine brains having no disease model. Reconstructed multi-slice conductivity images of in vivo canine brains with a pixel size of 1.4 mm showed a clear contrast between white and gray matter and also between normal and ischemic regions. We found that the conductivity value of an ischemic region decreased by about 10-14%. In a postmortem brain, conductivity values of white and gray matter decreased by about 4-8% compared to those in a live brain. Accumulating more experience of in vivo animal imaging experiments, we plan to move to human experiments. One of the important goals of our future work is the reduction of the imaging current to a level that a human subject can tolerate. The ability to acquire high-resolution conductivity images will find numerous clinical applications not supported by other medical imaging modalities. Potential applications in biology, chemistry and material science are also expected.

Measurement of induced magnetic flux density using injection current nonlinear encoding (ICNE) in MREIT.

Magnetic resonance electrical impedance tomography (MREIT) measures induced magnetic flux densities subject to externally injected currents in order to visualize conductivity distributions inside an electrically conducting object. Injection currents induce magnetic flux densities that appear in phase parts of acquired MR image data. In the conventional current injection method, we inject currents during the time segment between the end of the first RF pulse and the beginning of the reading gradient in order to ensure the gradient linearity. Noting that longer current injections can accumulate more phase changes, we propose a new pulse sequence called injection current nonlinear encoding (ICNE) where the duration of the injection current pulse is extended until the end of the reading gradient. Since the current injection during the reading gradient disturbs the gradient linearity, we first analyze the MR signal produced by the ICNE pulse sequence and suggest a novel algorithm to extract the induced magnetic flux density from the acquired MR signal. Numerical simulations and phantom experiments show that the new method is clearly advantageous in terms of the reduced noise level in measured magnetic flux density data. The amount of noise reduction depends on the choice of the data acquisition time and it was about 24% when we used a prolonged data acquisition time of 10.8 ms. The ICNE method will enhance the clinical applicability of the MREIT technique when it is combined with an appropriate phase artefact minimization method.

Optimization of current injection pulse width in MREIT.

In magnetic resonance electrical impedance tomography (MREIT), we inject electrical current into a volume conductor to induce a distribution of magnetic flux density. By measuring the internal magnetic flux density using an MR scanner, we reconstruct images of cross-sectional conductivity and current density distributions. One of the most important technical problems in MREIT is to reduce the noise level in the measured magnetic flux density data since it limits the quality of reconstructed images. The noise level is inversely proportional to the current injection pulse width and signal-to-noise ratio (SNR) of MR magnitude images. Knowing that we cannot simultaneously increase both factors for a chosen echo time, we show that there is an optimal current injection pulse width minimizing the noise level. Experimental results demonstrate that the optimal current injection pulse width and appropriately chosen data acquisition time considerably reduce the noise level. We suggest future works to reduce undesirable side effects due to an increased data acquisition time.

Phase artefact reduction in magnetic resonance electrical impedance tomography (MREIT).

Cross-sectional conductivity imaging in magnetic resonance electrical impedance tomography (MREIT) requires the measurement of internal magnetic flux density using an MRI scanner. Current injection MRI techniques have been used to induce magnetic flux density distributions that appear in phase parts of the obtained MR signals. Since any phase error, as well as noise, deteriorates the quality of reconstructed conductivity images, we must minimize them during the data acquisition process. In this paper, we describe a new method to correct unavoidable phase errors to reduce artefacts in reconstructed conductivity images. From numerical simulations and phantom experiments, we found that the zeroth- and first-order phase errors can be effectively minimized to produce better conductivity images. The promising results suggest that this technique should be employed together with improved MREIT pulse sequences in future studies of high-resolution conductivity imaging.

Conductivity image reconstruction from defective data in MREIT: numerical simulation and animal experiment.

Magnetic resonance electrical impedance tomography (MREIT) is designed to produce high resolution conductivity images of an electrically conducting subject by injecting current and measuring the longitudinal component, Bz, of the induced magnetic flux density B = (Bx, By, Bz). In MREIT, accurate measurements of Bz are essential in producing correct conductivity images. However, the measured Bz data may contain fundamental defects in local regions where MR magnitude image data are small. These defective Bz data result in completely wrong conductivity values there and also affect the overall accuracy of reconstructed conductivity images. Hence, these defects should be appropriately recovered in order to carry out any MREIT image reconstruction algorithm. This paper proposes a new method of recovering Bz data in defective regions based on its physical properties and neighboring information of Bz. The technique will be indispensable for conductivity imaging in MREIT from animal or human subjects including defective regions such as lungs, bones, and any gas-filled internal organs.

Basic setup for breast conductivity imaging using magnetic resonance electrical impedance tomography.

We present a new medical imaging technique for breast imaging, breast MREIT, in which magnetic resonance electrical impedance tomography (MREIT) is utilized to get high-resolution conductivity and current density images of the breast. In this work, we introduce the basic imaging setup of the breast MREIT technique with an investigation of four different imaging configurations of current-injection electrode positions and pathways through computer simulation studies. Utilizing the preliminary findings of a best breast MREIT configuration, additional numerical simulation studies have been carried out to validate breast MREIT at different levels of SNR. Finally, we have performed an experimental validation with a breast phantom on a 3.0 T MREIT system. The presented results strongly suggest that breast MREIT with careful imaging setups could be a potential imaging technique for human breast which may lead to early detection of breast cancer via improved differentiation of cancerous tissues in high-resolution conductivity images.

Harmonic decomposition in PDE-based denoising technique for magnetic resonance electrical impedance tomography.

Recent progress in magnetic resonance electrical impedance tomography (MREIT) research via simulation and biological tissue phantom studies have shown that conductivity images with higher spatial resolution and accuracy are achievable. In order to apply MREIT to human subjects, one of the important remaining problems to be solved is to reduce the amount of the injection current such that it meets the electrical safety regulations. However, by limiting the amount of the injection current according to the safety regulations, the measured MR data such as the z-component of magnetic flux density Bz in MREIT tend to have low SNR and get usually degraded in their accuracy due to the nonideal data acquisition system of an MR scanner. Furthermore, numerical differentiations of the measured Bz required by the conductivity image reconstruction algorithms tend to further deteriorate the quality and accuracy of the reconstructed conductivity images. In this paper, we propose a denoising technique that incorporates a harmonic decomposition. The harmonic decomposition is especially suitable for MREIT due to the physical characteristics of Bz. It effectively removes systematic and random noises, while preserving important key features in the MR measurements, so that improved conductivity images can be obtained. The simulation and experimental results demonstrate that the proposed denoising technique is effective for MREIT, producing significantly improved quality of conductivity images. The denoising technique will be a valuable tool in MREIT to reduce the amount of the injection current when it is combined with an improved MREIT pulse sequence.

Noise analysis in magnetic resonance electrical impedance tomography at 3 and 11 T field strengths.

In magnetic resonance electrical impedance tomography (MREIT), we measure the induced magnetic flux density inside an object subject to an externally injected current. This magnetic flux density is contaminated with noise, which ultimately limits the quality of reconstructed conductivity and current density images. By analysing and experimentally verifying the amount of noise in images gathered from two MREIT systems, we found that a carefully designed MREIT study will be able to reduce noise levels below 0.25 and 0.05 nT at main magnetic field strengths of 3 and 11 T, respectively, at a voxel size of 3 x 3 x 3 mm(3). Further noise level reductions can be achieved by optimizing MREIT pulse sequences and using signal averaging. We suggest two different methods to estimate magnetic flux noise levels, and the results are compared to validate the experimental setup of an MREIT system.

Identification of current density distribution in electrically conducting subject with anisotropic conductivity distribution.

Current density imaging (CDI) is able to visualize a three-dimensional current density distribution J inside an electrically conducting subject caused by an externally applied current. CDI may use a magnetic resonance imaging (MRI) scanner to measure the induced magnetic flux density B and compute J via the Ampere law [Formula: see text]. However, measuring all three components of B = (B(x), B(y), B(z)) has a technical difficulty due to the requirement of orthogonal rotations of the subject inside the MRI scanner. In this work, we propose a new method of reconstructing a current density image using only B(z) data so that we can avoid the subject rotation procedure. The method utilizes an auxiliary injection current to compensate the missing information of B(x) and B(y). The major advantage of the method is its applicability to a subject with an anisotropic conductivity distribution. Numerical experiments show the feasibility of the new technique.

Electrical conductivity images of biological tissue phantoms in MREIT.

We present cross-sectional conductivity images of two biological tissue phantoms. Each of the cylindrical phantoms with both diameter and height of 140 mm contained chunks of biological tissues such as bovine tongue and liver, porcine muscle and chicken breast within a conductive agar gelatin as the background medium. We attached four recessed electrodes on the sides of the phantom with equal spacing among them. Injecting current pulses of 480 or 120 mA ms into the phantom along two different directions, we measured the z-component Bz of the induced magnetic flux density B=(Bx, By, Bz) with a magnetic resonance electrical impedance tomography (MREIT) system based on a 3.0 T MRI scanner. Using the harmonic Bz algorithm, we reconstructed cross-sectional conductivity images from the measured Bz data. Reconstructed images clearly distinguish different tissues in terms of both their shapes and conductivity values. In this paper, we experimentally demonstrate the feasibility of the MREIT technique in producing conductivity images of different biological soft tissues with a high spatial resolution and accuracy when we use a sufficient amount of the injection current.

A mathematical model for breast cancer lesion estimation: electrical impedance technique using TS2000 commercial system.

We present a mathematical model to analyze transadmittance data for the detection of breast cancer using TransScan TS2000 commercial system. The model was constructed based on the assumption that a lesion exists near the surface of a breast region. The breast region that is considered as a background is assumed to be homogeneous at least near the surface where we attach a planar array of electrodes. Based on the model, we developed a lesion estimation algorithm utilizing single- or multifrequency transadmittance data. The approximate ratio of two conductivity values for the lesion and background needs to be known to estimate the size of the lesion even though the location estimate does not require this ratio. From the results of numerical simulations with added noise, we suggest better ways of interpreting TS2000 transadmittance images for the detection of breast cancer with improved accuracy. Since this study provides a rigorous mathematical modeling of TS2000 commercial system, it will be possible to apply the technique to lesion estimation problems based on more realistic models of breast regions in future studies.

Image reconstruction of anisotropic conductivity tensor distribution in MREIT: computer simulation study.

We describe a novel method of reconstructing images of an anisotropic conductivity tensor distribution inside an electrically conducting subject in magnetic resonance electrical impedance tomography (MREIT). MREIT is a recent medical imaging technique combining electrical impedance tomography (EIT) and magnetic resonance imaging (MRI) to produce conductivity images with improved spatial resolution and accuracy. In MREIT, we inject electrical current into the subject through surface electrodes and measure the z-component Bz of the induced magnetic flux density using an MRI scanner. Here, we assume that z is the direction of the main magnetic field of the MRI scanner. Considering the fact that most biological tissues are known to have anisotropic conductivity values, the primary goal of MREIT should be the imaging of an anisotropic conductivity tensor distribution. However, up to now, all MREIT techniques have assumed an isotropic conductivity distribution in the image reconstruction problem to simplify the underlying mathematical theory. In this paper, we firstly formulate a new image reconstruction method of an anisotropic conductivity tensor distribution. We use the relationship between multiple injection currents and the corresponding induced Bz data. Simulation results show that the algorithm can successfully reconstruct images of anisotropic conductivity tensor distributions. While the results show the feasibility of the method, they also suggest a more careful design of data collection methods and data processing techniques compared with isotropic conductivity imaging.

Magnetic resonance electrical impedance tomography at 3 Tesla field strength.

Magnetic resonance electrical impedance tomography (MREIT) is a recently developed imaging technique that combines MRI and electrical impedance tomography (EIT). In MREIT, cross-sectional electrical conductivity images are reconstructed from the internal magnetic field density data produced inside an electrically conducting object when an electrical current is injected into the object. In this work we present the results of electrical conductivity imaging experiments, and performance evaluations of MREIT in terms of noise characteristics and spatial resolution. The MREIT experiment was performed with a 3.0 Tesla MRI system on a phantom with an inhomogeneous conductivity distribution. We reconstructed the conductivity images in a 128 x 128 matrix format by applying the harmonic B(z) algorithm to the z-component of the internal magnetic field density data. Since the harmonic B(z) algorithm uses only a single component of the internal magnetic field data, it was not necessary to rotate the object in the MRI scan. The root mean squared (RMS) errors of the reconstructed images were between 11% and 35% when the injection current was 24 mA.

Static conductivity imaging using variational gradient Bz algorithm in magnetic resonance electrical impedance tomography.

A new image reconstruction algorithm is proposed to visualize static conductivity images of a subject in magnetic resonance electrical impedance tomography (MREIT). Injecting electrical current into the subject through surface electrodes, we can measure the induced internal magnetic flux density B = (Bx, By, Bz) using an MRI scanner. In this paper, we assume that only the z-component Bz is measurable due to a practical limitation of the measurement technique in MREIT. Under this circumstance, a constructive MREIT imaging technique called the harmonic Bz algorithm was recently developed to produce high-resolution conductivity images. The algorithm is based on the relation between inverted delta2Bz and the conductivity requiring the computation of inverted delta2Bz. Since twice differentiations of noisy Bz data tend to amplify the noise, the performance of the harmonic Bz algorithm is deteriorated when the signal-to-noise ratio in measured Bz data is not high enough. Therefore, it is highly desirable to develop a new algorithm reducing the number of differentiations. In this work, we propose the variational gradient Bz algorithm where Bz is differentiated only once. Numerical simulations with added random noise confirmed its ability to reconstruct static conductivity images in MREIT. We also found that it outperforms the harmonic Bz algorithm in terms of noise tolerance. From a careful analysis of the performance of the variational gradient Bz algorithm, we suggest several methods to further improve the image quality including a better choice of basis functions, regularization technique and multilevel approach. The proposed variational framework utilizing only Bz will lead to different versions of improved algorithms.

Electrical conductivity imaging using gradient B, decomposition algorithm in magnetic resonance electrical impedance tomography (MREIT).

In magnetic resonance electrical impedance tomography (MREIT), we try to visualize cross-sectional conductivity (or resistivity) images of a subject. We inject electrical currents into the subject through surface electrodes and measure the z component Bz of the induced internal magnetic flux density using an MRI scanner. Here, z is the direction of the main magnetic field of the MRI scanner. We formulate the conductivity image reconstruction problem in MREIT from a careful analysis of the relationship between the injection current and the induced magnetic flux density Bz. Based on the novel mathematical formulation, we propose the gradient Bz decomposition algorithm to reconstruct conductivity images. This new algorithm needs to differentiate Bz only once in contrast to the previously developed harmonic Bz algorithm where the numerical computation of (inverted delta)2Bz is required. The new algorithm, therefore, has the important advantage of much improved noise tolerance. Numerical simulations with added random noise of realistic amounts show the feasibility of the algorithm in practical applications and also its robustness against measurement noise.