EIT Reconstruction Algorithms: Pitfalls, Challenges and Recent Developments

We review developments, issues and challenges in Electrical Impedance Tomography (EIT), for the 4th Workshop on Biomedical Applications of EIT, Manchester 2003. We focus on the necessity for three dimensional data collection and reconstruction, efficient solution of the forward problem and present and future reconstruction algorithms. We also suggest common pitfalls or ``inverse crimes'' to avoid.Comment: A review paper for the 4th Workshop on Biomedical Applications of EIT, Manchester, UK, 200

High-order regularized regression in Electrical Impedance Tomography

We present a novel approach for the inverse problem in electrical impedance tomography based on regularized quadratic regression. Our contribution introduces a new formulation for the forward model in the form of a nonlinear integral transform, that maps changes in the electrical properties of a domain to their respective variations in boundary data. Using perturbation theory the transform is approximated to yield a high-order misfit unction which is then used to derive a regularized inverse problem. In particular, we consider the nonlinear problem to second-order accuracy, hence our approximation method improves upon the local linearization of the forward mapping. The inverse problem is approached using Newton's iterative algorithm and results from simulated experiments are presented. With a moderate increase in computational complexity, the method yields superior results compared to those of regularized linear regression and can be implemented to address the nonlinear inverse problem

Absolute electrical impedance tomography (aEIT) guided ventilation therapy in critical care patients: simulations and future trends

Thoracic electrical impedance tomography (EIT) is a noninvasive, radiation-free monitoring technique whose aim is to reconstruct a cross-sectional image of the internal spatial distribution of conductivity from electrical measurements made by injecting small alternating currents via an electrode array placed on the surface of the thorax. The purpose of this paper is to discuss the fundamentals of EIT and demonstrate the principles of mechanical ventilation, lung recruitment, and EIT imaging on a comprehensive physiological model, which combines a model of respiratory mechanics, a model of the human lung absolute resistivity as a function of air content, and a 2-D finite-element mesh of the thorax to simulate EIT image reconstruction during mechanical ventilation. The overall model gives a good understanding of respiratory physiology and EIT monitoring techniques in mechanically ventilated patients. The model proposed here was able to reproduce consistent images of ventilation distribution in simulated acutely injured and collapsed lung conditions. A new advisory system architecture integrating a previously developed data-driven physiological model for continuous and noninvasive predictions of blood gas parameters with the regional lung function data/information generated from absolute EIT (aEIT) is proposed for monitoring and ventilator therapy management of critical care patients

A Data-Driven Edge-Preserving D-bar Method for Electrical Impedance Tomography

In Electrical Impedance Tomography (EIT), the internal conductivity of a body is recovered via current and voltage measurements taken at its surface. The reconstruction task is a highly ill-posed nonlinear inverse problem, which is very sensitive to noise, and requires the use of regularized solution methods, of which D-bar is the only proven method. The resulting EIT images have low spatial resolution due to smoothing caused by low-pass filtered regularization. In many applications, such as medical imaging, it is known \emph{a priori} that the target contains sharp features such as organ boundaries, as well as approximate ranges for realistic conductivity values. In this paper, we use this information in a new edge-preserving EIT algorithm, based on the original D-bar method coupled with a deblurring flow stopped at a minimal data discrepancy. The method makes heavy use of a novel data fidelity term based on the so-called {\em CGO sinogram}. This nonlinear data step provides superior robustness over traditional EIT data formats such as current-to-voltage matrices or Dirichlet-to-Neumann operators, for commonly used current patterns.Comment: 24 pages, 11 figure

Comparison of linear and non-linear monotononicity-based shape reconstruction using exact matrix characterizations

Detecting inhomogeneities in the electrical conductivity is a special case of the inverse problem in electrical impedance tomography, that leads to fast direct reconstruction methods. One such method can, under reasonable assumptions, exactly characterize the inhomogeneities based on monotonicity properties of either the Neumann-to-Dirichlet map (non-linear) or its Fr\'echet derivative (linear). We give a comparison of the non-linear and linear approach in the presence of measurement noise, and show numerically that the two methods give essentially the same reconstruction in the unit disk domain. For a fair comparison, exact matrix characterizations are used when probing the monotonicity relations to avoid errors from numerical solution to PDEs and numerical integration. Using a special factorization of the Neumann-to-Dirichlet map also makes the non-linear method as fast as the linear method in the unit disk geometry.Comment: 18 pages, 5 figures, 1 tabl

The regularized monotonicity method: detecting irregular indefinite inclusions

In inclusion detection in electrical impedance tomography, the support of perturbations (inclusion) from a known background conductivity is typically reconstructed from idealized continuum data modelled by a Neumann-to-Dirichlet map. Only few reconstruction methods apply when detecting indefinite inclusions, where the conductivity distribution has both more and less conductive parts relative to the background conductivity; one such method is the monotonicity method of Harrach, Seo, and Ullrich. We formulate the method for irregular indefinite inclusions, meaning that we make no regularity assumptions on the conductivity perturbations nor on the inclusion boundaries. We show, provided that the perturbations are bounded away from zero, that the outer support of the positive and negative parts of the inclusions can be reconstructed independently. Moreover, we formulate a regularization scheme that applies to a class of approximative measurement models, including the Complete Electrode Model, hence making the method robust against modelling error and noise. In particular, we demonstrate that for a convergent family of approximative models there exists a sequence of regularization parameters such that the outer shape of the inclusions is asymptotically exactly characterized. Finally, a peeling-type reconstruction algorithm is presented and, for the first time in literature, numerical examples of monotonicity reconstructions for indefinite inclusions are presented.Comment: 28 pages, 7 figure