Ohin Kwon

Magnetic resonance electrical impedance tomography (MREIT): simulation study of J-substitution algorithm

We developed a new image reconstruction algorithm for magnetic resonance electrical impedance tomography (MREIT). MREIT is a new EIT imaging technique integrated into magnetic resonance imaging (MRI) system. Based on the assumption that internal current density distribution is obtained using magnetic resonance imaging (MRI) technique, the new image reconstruction algorithm called J-substitution algorithm produces cross-sectional static images of resistivity (or conductivity) distributions. Computer simulations show that the spatial resolution of resistivity image is comparable to that of MRI. MREIT provides accurate high-resolution cross-sectional resistivity images making resistivity values of various human tissues available for many biomedical applications.

Reconstruction of conductivity and current density images using only one component of magnetic field measurements

Magnetic resonance current density imaging (MRCDI) is to provide current density images of a subject using a magnetic resonance imaging (MRI) scanner with a current injection apparatus. The injection current generates a magnetic field that we can measure from MR phase images. We obtain internal current density images from the measured magnetic flux densities via Amperes law. However, we must rotate the subject to acquire all of the three components of the induced magnetic flux density. This subject rotation is impractical in clinical MRI scanners when the subject is a human body. In this paper, we propose a way to eliminate the requirement of subject rotation by careful mathematical analysis of the MRCDI problem. In our new MRCDI technique, we need to measure only one component of the induced magnetic flux density and reconstruct both cross-sectional conductivity and current density images without any subject rotation.

Conductivity and current density image reconstruction using harmonic Bz algorithm in magnetic resonance electrical impedance tomography

Magnetic resonance electrical impedance tomography (MREIT) is to provide cross-sectional images of the conductivity distribution sigma of a subject. While injecting current into the subject, we measure one component Bz of the induced magnetic flux density B = (Bx, By, Bz) using an MRI scanner. Based on the relation between (inverted delta)2 Bz and inverted delta sigma, the harmonic Bz algorithm reconstructs an image of sigma using the measured Bz data from multiple imaging slices. After we obtain sigma, we can reconstruct images of current density distributions for any given current injection method. Following the description of the harmonic Bz algorithm, this paper presents reconstructed conductivity and current density images from computer simulations and phantom experiments using four recessed electrodes injecting six different currents of 26 mA. For experimental results, we used a three-dimensional saline phantom with two polyacrylamide objects inside. We used our 0.3 T (tesla) experimental MRI scanner to measure the induced Bz. Using the harmonic Bz algorithm, we could reconstruct conductivity and current density images with 82 x 82 pixels. The pixel size was 0.6 x 0.6 mm2. The relative L2 errors of the reconstructed images were between 13.8 and 21.5% when the signal-to-noise ratio (SNR) of the corresponding MR magnitude images was about 30. The results suggest that in vitro and in vivo experimental studies with animal subjects are feasible. Further studies are requested to reduce the amount of injection current down to less than 1 mA for human subjects.

J-substitution algorithm in magnetic resonance electrical impedance tomography (MREIT): phantom experiments for static resistivity images

Recently, a new static resistivity image reconstruction algorithm is proposed utilizing internal current density data obtained by magnetic resonance current density imaging technique. This new imaging method is called magnetic resonance electrical impedance tomography (MREIT). The derivation and performance of J-substitution algorithm in MREIT have been reported as a new accurate and high-resolution static impedance imaging technique via computer simulation methods. In this paper, we present experimental procedures, denoising techniques and image reconstructions using a 0.3-tesla (T) experimental MREIT system and saline phantoms. MREIT using J-substitution algorithm effectively utilizes the internal current density information resolving the problem inherent in a conventional EIT, that is, the low sensitivity of boundary measurements to any changes of internal tissue resistivity values. Resistivity images of saline phantoms show an accuracy of 6.8%-47.2% and spatial resolution of 64 /spl times/ 64. Both of them can be significantly improved by using an MRI system with a better signal-to-noise ratio.

Electrical conductivity imaging using gradient B/sub z/ 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.

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.

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.

Three-dimensional forward solver and its performance analysis for magnetic resonance electrical impedance tomography (MREIT) using recessed electrodes

In magnetic resonance electrical impedance tomography (MREIT), we try to reconstruct a cross-sectional resistivity (or conductivity) image of a subject. When we inject a current through surface electrodes, it generates a magnetic field. Using a magnetic resonance imaging (MRI) scanner, we can obtain the induced magnetic flux density from MR phase images of the subject. We use recessed electrodes to avoid undesirable artefacts near electrodes in measuring magnetic flux densities. An MREIT image reconstruction algorithm produces cross-sectional resistivity images utilizing the measured internal magnetic flux density in addition to boundary voltage data. In order to develop such an image reconstruction algorithm, we need a three-dimensional forward solver. Given injection currents as boundary conditions, the forward solver described in this paper computes voltage and current density distributions using the finite element method (FEM). Then, it calculates the magnetic flux density within the subject using the Biot-Savart law and FEM. The performance of the forward solver is analysed and found to be enough for use in MREIT for resistivity image reconstructions and also experimental designs and validations. The forward solver may find other applications where one needs to compute voltage, current density and magnetic flux density distributions all within a volume conductor.

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.

On a nonlinear partial differential equation arising in Magnetic Resonance Electrical Impedance Tomography

This paper considers the fundamental questions, such as existence and uniqueness, of a mathematical model arising in the MREIT system, which is an electrical impedance tomography technique integrated with magnetic resonance imaging. The mathematical model for MREIT is the Neumann problem of a nonlinear elliptic partial differential equation