An iterative Newton-Raphson method to solve the inverse admittivity problem

Regularization Solver Guided FISTA for Electrical Impedance Tomography

Electrical impedance tomography (EIT) is non-destructive monitoring technology that can visualize the conductivity distribution in the observed area. The inverse problem for imaging is characterized by a serious nonlinear and ill-posed nature, which leads to the low spatial resolution of the reconstructions. The iterative algorithm is an effective method to deal with the imaging inverse problem. However, the existing iterative imaging methods have some drawbacks, such as random and subjective initial parameter setting, very time consuming in vast iterations and shape blurring with less high-order information, etc. To solve these problems, this paper proposes a novel fast convergent iteration method for solving the inverse problem and designs an initial guess method based on an adaptive regularization parameter adjustment. This method is named the Regularization Solver Guided Fast Iterative Shrinkage Threshold Algorithm (RS-FISTA). The iterative solution process under the L1-norm regular constraint is derived in the LASSO problem. Meanwhile, the Nesterov accelerator is introduced to accelerate the gradient optimization race in the ISTA method. In order to make the initial guess contain more prior information and be independent of subjective factors such as human experience, a new adaptive regularization weight coefficient selection method is introduced into the initial conjecture of the FISTA iteration as it contains more accurate prior information of the conductivity distribution. The RS-FISTA method is compared with the methods of Landweber, CG, NOSER, Newton-Raphson, ISTA and FISTA, six different distributions with their optimal parameters. The SSIM, RMSE and PSNR of RS-FISTA methods are 0.7253, 3.44 and 37.55, respectively. In the performance test of convergence, the evaluation metrics of this method are relatively stable at 30 iterations. This shows that the proposed method not only has better visualization, but also has fast convergence. It is verified that the RS-FISTA algorithm is the better algorithm for EIT reconstruction from both simulation and physical experiments.

EIT Imaging of admittivities with a D-bar method and spatial prior: experimental results for absolute and difference imaging

Electrical impedance tomography (EIT) is an emerging imaging modality that uses harmless electrical measurements taken on electrodes at a body's surface to recover information about the internal electrical conductivity and or permittivity. The image reconstruction task of EIT is a highly nonlinear inverse problem that is sensitive to noise and modeling errors making the image reconstruction task challenging. D-bar methods solve the nonlinear problem directly, bypassing the need for detailed and time-intensive forward models, to provide absolute (static) as well as time-difference EIT images. Coupling the D-bar methodology with the inclusion of high confidence a priori data results in a noise-robust regularized image reconstruction method. In this work, the a priori D-bar method for complex admittivities is demonstrated effective on experimental tank data for absolute imaging for the first time. Additionally, the method is adjusted for, and tested on, time-difference imaging scenarios. The ability of the method to be used for conductivity, permittivity, absolute as well as time-difference imaging provides the user with great flexibility without a high computational cost.

Direct EIT reconstructions of complex admittivities on a chest-shaped domain in 2-D

Electrical impedance tomography (EIT) is a medical imaging technique in which current is applied on electrodes on the surface of the body, the resulting voltage is measured, and an inverse problem is solved to recover the conductivity and/or permittivity in the interior. Images are then formed from the reconstructed conductivity and permittivity distributions. In the 2-D geometry, EIT is clinically useful for chest imaging. In this work, an implementation of a D-bar method for complex admittivities on a general 2-D domain is presented. In particular, reconstructions are computed on a chest-shaped domain for several realistic phantoms including a simulated pneumothorax, hyperinflation, and pleural effusion. The method demonstrates robustness in the presence of noise. Reconstructions from trigonometric and pairwise current injection patterns are included.

A direct D-bar reconstruction algorithm for recovering a complex conductivity in 2-D

A direct reconstruction algorithm for complex conductivities in W^2,∞ (Ω), where Ω is a bounded, simply connected Lipschitz domain in ℝ², is presented. The framework is based on the uniqueness proof by Francini [Inverse Problems 20 2000], but equations relating the Dirichlet-to-Neumann to the scattering transform and the exponentially growing solutions are not present in that work, and are derived here. The algorithm constitutes the first D-bar method for the reconstruction of conductivities and permittivities in two dimensions. Reconstructions of numerically simulated chest phantoms with discontinuities at the organ boundaries are included.

Imbalances in Regional Lung Ventilation

Imbalances in regional lung ventilation, with gravity-dependent collapse and overdistention of nondependent zones, are likely associated to ventilator-induced lung injury. Electric impedance tomography is a new imaging technique that is potentially capable of monitoring those imbalances. The aim of this study was to validate electrical impedance tomography measurements of ventilation distribution, by comparison with dynamic computerized tomography in a heterogeneous population of critically ill patients under mechanical ventilation. Multiple scans with both devices were collected during slow-inflation breaths. Six repeated breaths were monitored by impedance tomography, showing acceptable reproducibility. We observed acceptable agreement between both technologies in detecting right-left ventilation imbalances (bias = 0% and limits of agreement = -10 to +10%). Relative distribution of ventilation into regions or layers representing one-fourth of the thoracic section could also be assessed with good precision. Depending on electrode positioning, impedance tomography slightly overestimated ventilation imbalances along gravitational axis. Ventilation was gravitationally dependent in all patients, with some transient blockages in dependent regions synchronously detected by both scanning techniques. Among variables derived from computerized tomography, changes in absolute air content best explained the integral of impedance changes inside regions of interest (r(2) > or = 0.92). Impedance tomography can reliably assess ventilation distribution during mechanical ventilation.

Elliptic Cylinder Geometry for Distinguishability Analysis in Impedance Tomography

Electrical impedance tomography (EIT) is a technique that computes the cross-sectional impedance distribution within the body by using current and voltage measurements made on the body surface. It has been reported that the image reconstruction is distorted considerably when the boundary shape is considered to be more elliptical than circular as a more realistic shape for the measurement boundary. This paper describes an alternative framework for determining the distinguishability region with a finite measurement precision for different conductivity distributions in a body modeled by elliptic cylinder geometry. The distinguishable regions are compared in terms of modeling error for predefined inhomogeneities with elliptical and circular approaches for a noncircular measurement boundary at the body surface. Since most objects investigated by EIT are noncircular in shape, the analytical solution for the forward problem for the elliptical cross section approach is shown to be useful in order to reach a better assessment of the distinguishability region defined in a noncircular boundary. This paper is concentrated on centered elliptic inhomogeneity in the elliptical boundary and an analytic solution for this type of forward problem. The distinguishability performance of elliptical cross section with cosine injected current patterns is examined for different parameters of elliptical geometry.

Electrical impedance tomography of human brain function using reconstruction algorithms based on the finite element method

Electrical impedance tomography (EIT) is a recently developed technique which enables the internal conductivity of an object to be imaged using rings of external electrodes. In a recent study, EIT during cortical evoked responses showed encouraging changes in the raw impedance measurements, but reconstructed images were noisy. A simplified reconstruction algorithm was used which modelled the head as a homogeneous sphere. In the current study, the development and validation of an improved reconstruction algorithm are described in which realistic geometry and conductivity distributions have been incorporated using the finite element method. Data from computer simulations and spherical or head-shaped saline-filled tank phantoms, in which the skull was represented by a concentric shell of plaster of Paris or a real human skull, have been reconstructed into images. There were significant improvements in image quality as a result of the incorporation of accurate geometry and extracerebral layers in the reconstruction algorithm. Image quality, assessed by blinded subjective expert observers, also improved significantly when data from the previous evoked response study were reanalysed with the new algorithm. In preliminary images collected during epileptic seizures, the new algorithm generated EIT conductivity changes which were consistent with the electrographic ictal activity. Incorporation of realistic geometry and conductivity into the reconstruction algorithm significantly improves the quality of EIT images and lends encouragement to the belief that EIT may provide a low-cost, portable functional neuroimaging system in the foreseeable future.

Direct estimation of Cole parameters in multifrequency EIT using a regularized Gauss-Newton method

A major drawback of electrical impedance tomography is the poor quality of the conductivity images, i.e., the low spatial resolution as well as large errors in the reconstructed conductivity values. The main reason is the necessity for regularization of the ill-conditioned inverse problem which results in excessive spatial low-pass filtering. A novel regularization method (SMORR (spectral modelling regularized reconstructor)) is proposed, which is based on the inclusion of spectral a priori information in the form of appropriate tissue models (e.g. Cole models). This approach reduces the ill-posedness of the inverse problem, when multifrequency data are available. An additional advantage is the direct reconstruction of the (physiological) tissue parameters of interest instead of the conductivities. SMORR was compared with posterior fitting of a Cole model to the conductivity spectra obtained with a classical iterative reconstruction scheme at various frequencies. SMORR performed significantly better than the reference method concerning robustness against noise in the data.

Estimation of anomaly location and size using electrical impedance tomography

We developed a new algorithm that estimates locations and sizes of anomalies in electrically conducting medium based on electrical impedance tomography (EIT) technique. When only the boundary current and voltage measurements are available, it is not practically feasible to reconstruct accurate high-resolution cross-sectional conductivity or resistivity images of a subject. In this paper, we focus our attention on the estimation of locations and sizes of anomalies with different conductivity values compared with the background tissues. We showed the performance of the algorithm from experimental results using a 32-channel EIT system and saline phantom. With about 1.73% measurement error in boundary current-voltage data, we found that the minimal size (area) of the detectable anomaly is about 0.72% of the size (area) of the phantom. Potential applications include the monitoring of impedance related physiological events and bubble detection in two-phase flow. Since this new algorithm requires neither any forward solver nor time-consuming minimization process, it is fast enough for various real-time applications in medicine and nondestructive testing.

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.

The number of electrodes and basis functions in EIT image reconstruction

In electrical impedance tomography, many factors affect the image reconstruction results. Among them are the number of electrodes (NOE) and the number of conductivity basis functions (NOCBF) for image reconstruction. The NOCBF generally reflects the density of the mesh with which images are reconstructed. How and to what extent do these factors affect the image reconstruction and corresponding images? In this area detailed analysis is still lacking. This study aims to address the above question. In this study, image reconstruction and its ill-posed condition were analysed by singular value decomposition (SVD) and spectral expansion theory with different configurations of NOE and NOCBF. The results in this study indicate that for a circular 2D plane object with electrodes evenly located around the boundary: (1) Under certain conditions, increasing the NOE enables us to improve the ill-posed condition in image reconstruction and hence improve the image quality. Generally more improvement is expected near the image periphery than in the image centre. (2) Increasing the NOCBF generally worsens the ill-posed condition. But it enables the solution to be sought in a finer subspace and may be able to improve the image quality on the periphery, while generally the result in image centre depends more on the prior information incorporated in the regularization.

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.

Preliminary results from an EIT breast imaging simulation system

An electrical impedance tomography (EIT) simulation system for breast imaging has been developed using impedance data from a previous study and a finite-element model (FEM). This system has the functionality to construct various models of the breast, image the boundary voltages developed from any injection schema and provide parametric Cole-Cole modelling. Simulations indicate that breast carcinoma can be imaged and multi-frequency Cole-Cole dispersion data can be extracted. This is intended as the first stage towards providing an artificial intelligence based system capable of producing clinically relevant analytical data.

Feasibility of biomedical applications of Hall effect imaging

Hall effect imaging is a new technique for mapping the electrical properties of a sample. Its principle has been demonstrated in two- and three-dimensional phantom images. Based on the experimental data and theoretical understanding of this technique developed over the past few years, this paper addresses the most relevant question for biomedical applications: whether Hall effect imaging is ultimately applicable to complex biological systems such as the human body. The arguments are given at the basic physics level, so that the conclusion is independent of current technology status. These arguments are corroborated with imaging data of an aorta sample. The conclusion is that Hall effect imaging is not suited for quantifying the electrical constants in complex biological samples. This technique is able to produce high-resolution volume images of samples in vitro that reflect their electrical heterogeneity. However, quantitative measurements of electrical constants are not practical for complex samples.

Electrical impedance tomography of complex conductivity distributions with noncircular boundary

Electrical impedance tomography (EIT) uses low-frequency current and voltage measurements made on the boundary of a body to compute the conductivity distribution within the body. Since the permittivity distribution inside the body also contributes significantly to the measured voltages, the present reconstruction algorithm images complex conductivity distributions. A finite element model (FEM) is used to solve the forward problem, using a 6017-node mesh for a piecewise-linear potential distribution. The finite element solution using this mesh is compared with the analytical solution for a homogeneous field and a maximum error of 0.05% is observed in the voltage distribution. The boundary element method (BEM) is also used to generate the voltage data for inhomogeneous conductivity distributions inside regions with noncircular boundaries. An iterative reconstruction algorithm is described for approximating both the conductivity and permittivity distributions from this data. The results for an off-centered inhomogeneity showed a 35% improvement in contrast from that seen with only one iteration, for both the conductivity and the permittivity values. It is also shown that a significant improvement in images results from accurately modeling a noncircular boundary. Both static and difference images are distorted by assuming a circular boundary and the amount of distortion increases significantly as the boundary shape becomes more elliptical. For a homogeneous field in an elliptical body with axis ratio of 0.73, an image reconstructed assuming the boundary to be circular has an artifact at the center of the image with an error of 20%. This error increased to 37% when the axis ratio was 0.64. A reconstruction algorithm which used a mesh with the same axis ratio as the elliptical boundary reduced the error in the conductivity values to within 0.5% of the actual values.