William R. B. Lionheart

A Matlab toolkit for three-dimensional electrical impedance tomography: a contribution to the Electrical Impedance and Diffuse Optical Reconstruction Software project

The objective of the Electrical Impedance and Diffuse Optical Reconstruction Software project is to develop freely available software that can be used to reconstruct electrical or optical material properties from boundary measurements. Nonlinear and ill posed problems such as electrical impedance and optical tomography are typically approached using a finite element model for the forward calculations and a regularized nonlinear solver for obtaining a unique and stable inverse solution. Most of the commercially available finite element programs are unsuitable for solving these problems because of their conventional inefficient way of calculating the Jacobian, and their lack of accurate electrode modelling. A complete package for the two-dimensional EIT problem was officially released by Vauhkonen et al at the second half of 2000. However most industrial and medical electrical imaging problems are fundamentally three-dimensional. To assist the development we have developed and released a free toolkit of Matlab routines which can be employed to solve the forward and inverse EIT problems in three dimensions based on the complete electrode model along with some basic visualization utilities, in the hope that it will stimulate further development. We also include a derivation of the formula for the Jacobian (or sensitivity) matrix based on the complete electrode model.

Nonuniqueness in diffusion-based optical tomography

A condition on nonuniqueness in optical tomography is stated. The main result applies to steady-state (dc) diffusion-based optical tomography, wherein we demonstrate that simultaneous unique recovery of diffusion and absorption coefficients cannot be achieved. A specific example of two images that give identical dc data is presented. If the refractive index is considered an unknown, then nonuniqueness also occurs in frequency-domain and time-domain optical tomography, if the underlying model of the diffusion approximation is employed.

GREIT: A unified approach to 2D linear EIT reconstruction of lung images

Electrical impedance tomography (EIT) is an attractive method for clinically monitoring patients during mechanical ventilation, because it can provide a non-invasive continuous image of pulmonary impedance which indicates the distribution of ventilation. However, most clinical and physiological research in lung EIT is done using older and proprietary algorithms; this is an obstacle to interpretation of EIT images because the reconstructed images are not well characterized. To address this issue, we develop a consensus linear reconstruction algorithm for lung EIT, called GREIT (Graz consensus Reconstruction algorithm for EIT). This paper describes the unified approach to linear image reconstruction developed for GREIT. The framework for the linear reconstruction algorithm consists of (1) detailed finite element models of a representative adult and neonatal thorax, (2) consensus on the performance figures of merit for EIT image reconstruction and (3) a systematic approach to optimize a linear reconstruction matrix to desired performance measures. Consensus figures of merit, in order of importance, are (a) uniform amplitude response, (b) small and uniform position error, (c) small ringing artefacts, (d) uniform resolution, (e) limited shape deformation and (f) high resolution. Such figures of merit must be attained while maintaining small noise amplification and small sensitivity to electrode and boundary movement. This approach represents the consensus of a large and representative group of experts in EIT algorithm design and clinical applications for pulmonary monitoring. All software and data to implement and test the algorithm have been made available under an open source license which allows free research and commercial use.

Nonlinear image reconstruction for electrical capacitance tomography using experimental data

Electrical capacitance tomography (ECT) attempts to image the permittivity distribution of an object by measuring the electrical capacitances between sets of electrodes placed around its periphery. Image reconstruction in ECT is a nonlinear and ill-posed inverse problem. Although reconstruction techniques based on a linear approximation are fast, they are not adequate for all cases. In this paper, we study the nonlinearity of the inverse permittivity problem of ECT. A regularized Gauss-Newton scheme has been implemented for nonlinear image reconstruction. The forward problem has been solved at each iteration using the finite element method and the Jacobian matrix is recalculated using an efficient adjoint field method. Regularization techniques are required to overcome the ill-posedness: smooth generalized Tikhonov regularization for the smoothly varying case, and total variation (TV) regularization when there is a sharp transition of the permittivity have been used. The reconstruction results for experimental ECT data demonstrate the advantage of TV regularization for jump changes, and show improvement of the image quality by using nonlinear reconstruction methods.

A MATLAB package for the EIDORS project to reconstruct two-dimensional EIT images.

The EIDORS (electrical impedance and diffuse optical reconstruction software) project aims to produce a software system for reconstructing images from electrical or diffuse optical data. MATLAB is a software that is used in the EIDORS project for rapid prototyping, graphical user interface construction and image display. We have written a MATLAB package (http://venda.uku.fi/ vauhkon/) which can be used for two-dimensional mesh generation, solving the forward problem and reconstructing and displaying the reconstructed images (resistivity or admittivity). In this paper we briefly describe the mathematical theory on which the codes are based on and also give some examples of the capabilities of the package.

In Vivo Impedance Imaging With Total Variation Regularization

We show that electrical impedance tomography (EIT) image reconstruction algorithms with regularization based on the total variation (TV) functional are suitable for in vivo imaging of physiological data. This reconstruction approach helps to preserve discontinuities in reconstructed profiles, such as step changes in electrical properties at interorgan boundaries, which are typically smoothed by traditional reconstruction algorithms. The use of the TV functional for regularization leads to the minimization of a nondifferentiable objective function in the inverse formulation. This cannot be efficiently solved with traditional optimization techniques such as the Newton method. We explore two implementations methods for regularization with the TV functional: the lagged diffusivity method and the primal dual-interior point method (PD-IPM). First we clarify the implementation details of these algorithms for EIT reconstruction. Next, we analyze the performance of these algorithms on noisy simulated data. Finally, we show reconstructed EIT images of in vivo data for ventilation and gastric emptying studies. In comparison to traditional quadratic regularization, TV regulariza tion shows improved ability to reconstruct sharp contrasts.

Recovery of region boundaries of piecewise constant coefficients of an elliptic PDE from boundary data

In this study we consider the recovery of smooth region boundaries of piecewise constant coefficients of an elliptic PDE, - a+b = f, from data on the exterior boundary . The assumption made is that the values of the coefficients (a,b) are known a priori but the information about the geometry of the smooth region boundaries where a and b are discontinous is missing. For the full characterization of (a,b) it is then sufficient to find the region boundaries separating different values of the coefficients. This leads to a nonlinear ill-posed inverse problem. In this study we propose a numerical algorithm that is based on the finite-element method and subdivision of the discretization elements. We formulate the forward problem as a mapping from a set of coefficients representing boundary shapes to data on , and derive the Jacobian of this forward mapping. Then an iterative algorithm which seeks a boundary configuration minimizing the residual norm between measured and predicted data is implemented. The method is illustrated first for a general elliptic PDE and then applied to optical tomography where the goal is to find the diffusion and absorption coefficients of the object by transilluminating the object with visible or near-infrared light. Numerical test results for this specific application are given with synthetic data.

Uniqueness and reconstruction in magnetic resonance-electrical impedance tomography (MR-EIT)

Magnetic resonance-electrical impedance tomography (MR-EIT) was first proposed in 1992. Since then various reconstruction algorithms have been suggested and applied. These algorithms use peripheral voltage measurements and internal current density measurements in different combinations. In this study the problem of MR-EIT is treated as a hyperbolic system of first-order partial differential equations, and three numerical methods are proposed for its solution. This approach is not utilized in any of the algorithms proposed earlier. The numerical solution methods are integration along equipotential surfaces (method of characteristics), integration on a Cartesian grid, and inversion of a system matrix derived by a finite difference formulation. It is shown that if some uniqueness conditions are satisfied, then using at least two injected current patterns, resistivity can be reconstructed apart from a multiplicative constant. This constant can then be identified using a single voltage measurement. The methods proposed are direct, non-iterative, and valid and feasible for 3D reconstructions. They can also be used to easily obtain slice and field-of-view images from a 3D object. 2D simulations are made to illustrate the performance of the algorithms.

Absolute Conductivity Reconstruction in Magnetic Induction Tomography Using a Nonlinear Method

Magnetic induction tomography (MIT) attempts to image the electrical and magnetic characteristics of a target using impedance measurement data from pairs of excitation and detection coils. This inverse eddy current problem is nonlinear and also severely ill posed so regularization is required for a stable solution. A regularized Gauss-Newton algorithm has been implemented as a nonlinear, iterative inverse solver. In this algorithm, one needs to solve the forward problem and recalculate the Jacobian matrix for each iteration. The forward problem has been solved using an edge based finite element method for magnetic vector potential A and electrical scalar potential V, a so called A, A-V formulation. A theoretical study of the general inverse eddy current problem and a derivation, paying special attention to the boundary conditions, of an adjoint field formula for the Jacobian is given. This efficient formula calculates the change in measured induced voltage due to a small perturbation of the conductivity in a region. This has the advantage that it involves only the inner product of the electric fields when two different coils are excited, and these are convenient computationally. This paper also shows that the sensitivity maps change significantly when the conductivity distribution changes, demonstrating the necessity for a nonlinear reconstruction algorithm. The performance of the inverse solver has been examined and results presented from simulated data with added noise

A three-dimensional inverse finite-element method applied to experimental eddy-current imaging data

Eddy-current techniques can be used to create electrical conductivity mapping of an object. The eddy-current imaging system in this paper is a magnetic induction tomography (MIT) system. MIT images the electrical conductivity of the target based on impedance measurements from pairs of excitation and detection coils. The inverse problem here is ill-posed and nonlinear. Current state-of-the-art image reconstruction methods in MIT are generally based on linear algorithms. In this paper, a regularized Gauss-Newton scheme has been implemented based on an edge finite-element forward solver and an efficient formula for the Jacobian matrix. Applications of Tikhonov and total variation regularization have been studied. Results are presented from experimental data collected from a newly developed MIT system. The paper also presents further progress in using an MIT system for molten metal flow visualization in continuous casting by applying the proposed algorithm in a real experiment in a continuous casting pilot plant of Corus RD&T, Teesside Technology Centre.