Ramp-preserving denoising for conductivity image reconstruction in magnetic resonance electrical impedance tomography

In magnetic resonance electrical impedance tomography, among several conductivity image reconstruction algorithms, the harmonic B(z) algorithm has been successfully applied to B(z) data from phantoms and animals. The algorithm is, however, sensitive to measurement noise in B(z) data. Especially, in in vivo animal and human experiments where injection current amplitudes are limited within a few milliampere at most, measured B(z) data tend to have a low SNR. In addition, magnetic resonance (MR) signal void in outer layers of bones and gas-filled organs, for example, produces salt-pepper noise in the MR phase and, consequently, B(z) images. The B(z) images typically present areas of sloped transitions, which can be assimilated to ramps. Conductivity contrasts change ramp slopes in B(z) images and it is critical to preserve positions of those ramps to correctly recover edges in conductivity images. In this paper, we propose a ramp-preserving denoising method utilizing a structure tensor. Using an eigenvalue analysis, we identified local regions of salt-pepper noise. Outside the identified local regions, we applied an anisotropic smoothing to reduce noise while preserving their ramp structures. Inside the local regions of salt-pepper noise, we used an isotropic smoothing. After validating the proposed denoising method through numerical simulations, we applied it to in vivo animal imaging experiments. Both numerical simulation and experimental results show significant improvements in the quality of reconstructed conductivity images.     

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.     

A neural network image reconstruction technique for electrical impedance tomography

Reconstruction of images in electrical impedance tomography requires the solution of a nonlinear inverse problem on noisy data. This problem is typically ill-conditioned and requires either simplifying assumptions or regularization based on a priori knowledge. The authors present a reconstruction algorithm using neural network techniques which calculates a linear approximation of the inverse problem directly from finite element simulations of the forward problem. This inverse is adapted to the geometry of the medium and the signal-to-noise ratio (SNR) used during network training. Results show good conductivity reconstruction where measurement SNR is similar to the training conditions. The advantages of this method are its conceptual simplicity and ease of implementation, and the ability to control the compromise between the noise performance and resolution of the image reconstruction.     

Noninvasive imaging of bioimpedance distribution by means of current reconstruction magnetic resonance electrical impedance tomography

We have developed a novel magnetic resonance electrical impedance tomography (MREIT) algorithm-current reconstruction MREIT algorithm-for noninvasive imaging of electrical impedance distribution of a biological system using only one component of magnetic flux density. The newly proposed algorithm uses the inverse of Biot-Savart Law to reconstruct the current density distribution, and then, uses a modified J-substitution algorithm to reconstruct the conductivity image. A series of computer simulations has been conducted to evaluate the performance of the proposed current reconstruction MREIT algorithm with simulation settings for breast cancer imaging applications, with consideration of measurement noise, current injection strength, size of simulated tumors, spatial resolution, and position dependency. The present simulation results are highly promising, demonstrating the high spatial resolution, high accuracy in conductivity reconstruction, and robustness against noise of the proposed algorithm for imaging electrical impedance of a biological system. The present MREIT method may have potential applications to breast cancer imaging and imaging of other organs.     

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 x 64. Both of them can be significantly improved by using an MRI system with a better signal-to-noise ratio.     

Magnetic resonance driven electrical impedance tomography: a simulation study

Magnetic resonance electrical impedance tomography (MREIT) is a method for reconstructing a three-dimensional image of the conductivity distribution in a target volume using magnetic resonance (MR). In MREIT, currents are applied to the volume through surface electrodes and their effects on the MR induced magnetic fields are analyzed to produce the conductance image. However, current injection through surface electrodes poses technical problems such as the limitation on the safely applicable currents. In this paper, we present a new method called magnetic resonance driven electrical impedance tomography (MRDEIT), where the magnetic resonance in each voxel is used as the applied magnetic field source, and the resultant electromagnetic field is measured through surface electrodes or radio-frequency (RF) detectors placed near the surface. Because the applied magnetic field is at the RF frequency and eddy currents are the integral components in the method, a vector wave equation for the electric field is used as the basis of the analysis instead of a quasi-static approximation. Using computer simulations, it is shown that complex permittivity images can be reconstructed using MRDEIT, but that improvements in signal detection are necessary for detecting moderate complex permittivity changes.     

Comparison of PDE-based nonlinear diffusion approaches for image enhancement and denoising in optical coherence tomography

A comparison between two nonlinear diffusion methods for denoising OCT images is performed. Specifically, we compare and contrast the performance of the traditional nonlinear Perona-Malik filter with a complex diffusion filter that has been recently introduced by Gilboa et al.. The complex diffusion approach based on the generalization of the nonlinear scale space to the complex domain by combining the diffusion and the free Schridinger equation is evaluated on synthetic images and also on representative OCT images at various noise levels. The performance improvement over the traditional nonlinear Perona-Malik filter is quantified in terms of noise suppression, image structural preservation and visual quality. An average signal-to-noise ratio (SNR) improvement of about 2.5 times and an average contrast to noise ratio (CNR) improvement of 49% was obtained while mean structure similarity (MSSIM) was practically not degraded after denoising. The nonlinear complex diffusion filtering can be applied with success to many OCT imaging applications. In summary, the numerical values of the image quality metrics along with the qualitative analysis results indicated the good feature preservation performance of the complex diffusion process, as desired for better diagnosis in medical imaging processing.     

Magnetic resonance electrical impedance tomography for monitoring electric field distribution during tissue electroporation

Electroporation is a phenomenon caused by externally applied electric field of an adequate strength and duration to cells that results in the increase of cell membrane permeability to various molecules, which otherwise are deprived of transport mechanism. As accurate coverage of the tissue with a sufficiently large electric field presents one of the most important conditions for successful electroporation, applications based on electroporation would greatly benefit with a method of monitoring the electric field, especially if it could be done during the treatment. As the membrane electroporation is a consequence of an induced transmembrane potential which is directly proportional to the local electric field, we propose current density imaging (CDI) and magnetic resonance electrical impedance tomography (MREIT) techniques to measure the electric field distribution during electroporation. The experimental part of the study employs CDI with short high-voltage pulses, while the theoretical part of the study is based on numerical simulations of MREIT. A good agreement between experimental and numerical results was obtained, suggesting that CDI and MREIT can be used to determine the electric field during electric pulse delivery and that both of the methods can be of significant help in planning and monitoring of future electroporation based clinical applications.     

Weighted least-squares criteria for electrical impedance tomography

Methods are developed for the design of electrical impedance tomographic reconstruction algorithms with specified properties. Assuming a starting model with constant conductivity or some other specified background distribution, an algorithm with the following properties is found. (1) The optimum constant for the starting model is determined automatically. (2) The weighted least-squares error between the predicted and measured power dissipation data is as small as possible. (3) The variance of the reconstructed conductivity from the starting model is minimized. (4) Potential distributions with the largest volume integral of gradient squared have the least influence on the reconstructed conductivity, and therefore distributions most likely to be corrupted by contact impedance effects are deemphasized. (5) Cells that dissipate the most power during the current injection tests tend to deviate least from the background value. For a starting model with nonconstant conductivity, the reconstruction algorithm has analogous properties.     

Optimal experiments in electrical impedance tomography

Electrical impedance tomography (EIT) is a noninvasive imaging technique which aims to image the impedance within a body from electrical measurements made on the surface. The reconstruction of impedance images is a ill-posed problem which is both extremely sensitive to noise and highly computationally intensive. The authors define an experimental measurement in EIT and calculate optimal experiments which maximize the distinguishability between the region to be imaged and a best-estimate conductivity distribution. These optimal experiments can be derived from measurements made on the boundary. The analysis clarifies the properties of different voltage measurement schemes. A reconstruction algorithm based on the use of optimal experiments is derived. It is shown to be many times faster than standard Newton-based reconstruction algorithms, and results from synthetic data indicate that the images that it produces are comparable.     

New denoising scheme for magnetic resonance spectroscopy signals

A new scheme for denoising magnetic resonance spectroscopy (MRS) signals is presented. This scheme is based on projecting noisy MRS signals in different domains, consecutively, and performing noise filtering operations in these domains. The domains are chosen such that the noise portion, which is inseparable from the desired signal in one domain, is separable in the other. A set of stable, linear, time-frequency (SLTF) transforms with different resolutions was selected for these projections as an example. Scheme evaluation was performed using extensive MRS signals with various noise levels. Compared with one domain denoising, it was observed that the proposed scheme gives superior results that compensate for the excess computational requirements. The proposed scheme supersedes also the wavelet packet denoising schemes.     

A direct reconstruction algorithm for electrical impedance tomography

A direct (noniterative) reconstruction algorithm for electrical impedance tomography in the two-dimensional (2-D), cross-sectional geometry is reviewed. New results of a reconstruction of a numerically simulated phantom chest are presented. The algorithm is based on the mathematical uniqueness proof by A. I. Nachman [1996] for the 2-D inverse conductivity problem. In this geometry, several of the clinical applications include monitoring heart and lung function, diagnosis of pulmonary embolus, diagnosis of pulmonary edema, monitoring for internal bleeding, and the early detection of breast cancer.     

Reconstruction of three-dimensional data for electrical impedance tomography

A two-dimensional reconstruction algorithm based on a modified version of the method of sensitivity regions is used to reconstruct data obtained from a three-dimensional finite element model. By using data obtained from off-drive-plane measurements an improved image of changes in resistivity on the drive plane is obtained.<>     

On optimal current patterns for electrical impedance tomography

We develop a statistical criterion for optimal patterns in planar circular electrical impedance tomography. These patterns minimize the total variance of the estimation for the resistance or conductance matrix. It is shown that trigonometric patterns (Isaacson, 1986), originally derived from the concept of distinguishability, are a special case of our optimal statistical patterns. New optimal random patterns are introduced. Recovering the electrical properties of the measured body is greatly simplified when optimal patterns are used. The Neumann-to-Dirichlet map and the optimal patterns are derived for a homogeneous medium with an arbitrary distribution of the electrodes on the periphery. As a special case, optimal patterns are developed for a practical EIT system with a finite number of electrodes. For a general nonhomogeneous medium, with no a priori restriction, the optimal patterns for the resistance and conductance matrix are the same. However, for a homogeneous medium, the best current pattern is the worst voltage pattern and vice versa. We study the effect of the number and the width of the electrodes on the estimate of resistivity and conductivity in a homogeneous medium. We confirm experimentally that the optimal patterns produce minimum conductivity variance in a homogeneous medium. Our statistical model is able to discriminate between a homogenous agar phantom and one with a 2 mm air hole with error probability (p-value) 1/1000.     

Using compound electrodes in electrical impedance tomography

A compound electrode composed of a large outer electrode to inject current and a small inner electrode to sense voltage was developed and used to measure voltages from a physical phantom. The measured voltages were smaller in amplitude than those from conventional electrodes, demonstrating that the compound electrode can minimize contact impedance voltage drop from the measured data. A finite-element model was used for the compound electrode and incorporated into the regularized Newton-Raphson reconstruction algorithm. A sensitivity study showed that the reconstructed resistivity distributions are less dependent on the unknown contact resistance values for a compound electrode than a conventional electrode and that the use of a compound electrode results in improved images for the reconstruction algorithm     

Increasing image resolution in electrical impedance tomography

An effective approach to increase the image resolution in static electrical impedance tomography is proposed, in which the image with local high resolution is reconstructed by fine meshing only the impedance abnormal element in the finite element model based on a genetic algorithm. Experimental results from a laboratory phantom are presented     

An image reconstruction algorithm for three-dimensional electrical impedance tomography

Electrical impedance tomography (EIT) has been studied by many authors and in most of this work it has been considered to be a two-dimensional problem. Many groups are now turning their attention to the full three-dimensional case in which the computational demands become much greater. It is interesting to look for ways to reduce this demand and in this paper we describe an implementation of an algorithm that is able to achieve this by precomputing many of the quantities needed in the image reconstruction. The algorithm is based on a method called NOSER introduced some years ago by Cheney et al. [3]. In this paper we have significantly extended the method by introducing a more realistic electrode model into the analysis. We have given explicit formulae for the quantities involved so that the reader can reproduce our results.     

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.     

Measuring lung resistivity using electrical impedance tomography

We propose the use of electrical impedance tomography (EIT) imaging techniques in the measurement of lung resistivity for detection and monitoring of apnea and edema. In EIT, we inject currents into a subject using multiple electrodes and measure boundary voltages to reconstruct a cross-sectional image of internal resistivity distribution. We found that a simplified, therefore fast, version of the impedance imaging method can be used for detection and monitoring of apnea and edema. We have showed the feasibility of this method through computer simulations and human experiments. We speculate that the EIT imaging technique will be more reliable than the current impedance apnea monitoring method, since we are monitoring the change of internal lung resistivity. However, more study is required to verify that this method performs better in the presence of motion artifact than the conventional two-electrode impedance apnea monitoring method. Future work should include experiments which carefully simulate different kinds of motion artifacts.     

Estimation of optimal PDE-based denoising in the SNR sense.

This paper is concerned with finding the best partial differential equation-based denoising process, out of a set of possible ones. We focus on either finding the proper weight of the fidelity term in the energy minimization formulation or on determining the optimal stopping time of a nonlinear diffusion process. A necessary condition for achieving maximal SNR is stated, based on the covariance of the noise and the residual part. We provide two practical alternatives for estimating this condition by observing that the filtering of the image and the noise can be approximated by a decoupling technique, with respect to the weight or time parameters. Our automatic algorithm obtains quite accurate results on a variety of synthetic and natural images, including piecewise smooth and textured ones. We assume that the statistics of the noise were previously estimated. No a priori knowledge regarding the characteristics of the clean image is required. A theoretical analysis is carried out, where several SNR performance bounds are established for the optimal strategy and for a widely used method, wherein the variance of the residual part equals the variance of the noise.