Electrical Impedance Tomography for Piecewise Constant Domains Using Boundary Element Shape-Based Inverse Solutions

Shape-based solutions have recently received attention for certain ill-posed inverse problems. Their advantages include implicit imposition of relevant constraints and reduction in the number of unknowns, especially important for nonlinear ill-posed problems. We apply the shape-based approach to current-injection electrical impedance tomography (EIT) reconstructions. We employ a boundary element method (BEM) based solution for EIT. We introduce two shape models, one based on modified B-splines, and the other based on spherical harmonics, for BEM modeling of shapes. These methods allow us to parameterize the geometry of conductivity inhomogeneities in the interior of the volume. We assume the general shape of piecewise constant inhomogeneities is known but their conductivities and their exact location and shape is not. We also assume the internal conductivity profile is piecewise constant, meaning that each region has a constant conductivity. We propose and test three different regularization techniques to be used with either of the shape models. The performance of our methods is illustrated via both simulations in a digital torso model and phantom experiments when there is a single internal object. We observe that in the noisy environment, either simulated noise or real sources of noise in the experimental study, we get reasonable reconstructions. Since the signal-to-noise ratio (SNR) expected in modern EIT instruments is higher than that used in this study, these reconstruction methods may prove useful in practice

[1]  E. Somersalo,et al.  Statistical inverse problems: discretization, model reduction and inverse crimes , 2007 .

[2]  Masaru Ikehata,et al.  RECONSTRUCTION OF THE SHAPE OF THE INCLUSION BY BOUNDARY MEASUREMENTS , 1998 .

[3]  Dana H. Brooks,et al.  SCIRun/BioPSE: integrated problem solving environment for bioelectric field problems and visualization , 2004, 2004 2nd IEEE International Symposium on Biomedical Imaging: Nano to Macro (IEEE Cat No. 04EX821).

[4]  David Boas,et al.  Three-dimensional shape-based imaging of absorption perturbation for diffuse optical tomography. , 2003, Applied optics.

[5]  David Isaacson,et al.  Electric current computed tomography eigenvalues , 1990 .

[6]  G.J. Saulnier,et al.  ACT3: a high-speed, high-precision electrical impedance tomograph , 1991, IEEE Transactions on Biomedical Engineering.

[7]  Arto Voutilainen,et al.  Estimation of non-stationary region boundaries in EIT—state estimation approach , 2001 .

[8]  Rasmus R. Paulsen,et al.  Shape Modelling Using Markov Random Field Restoration of Point Correspondences , 2003, IPMI.

[9]  Simon R. Arridge,et al.  Three-dimensional reconstruction of shape and piecewise constant region values for optical tomography using spherical harmonic parametrization and a boundary element method , 2006 .

[10]  Eric L. Miller,et al.  Estimation-theoretic algorithms and bounds for three-dimensional polar shape-based imaging in diffuse optical tomography , 2006, 3rd IEEE International Symposium on Biomedical Imaging: Nano to Macro, 2006..

[11]  R. Kress,et al.  Inverse Acoustic and Electromagnetic Scattering Theory , 1992 .

[12]  A. Stroud Approximate calculation of multiple integrals , 1973 .

[13]  Ramani Duraiswami,et al.  Efficient 2D and 3D electrical impedance tomography using dual reciprocity boundary element techniques , 1998 .

[14]  J. Newell,et al.  Theory and performance of an adaptive current tomography system. , 1988, Clinical physics and physiological measurement : an official journal of the Hospital Physicists' Association, Deutsche Gesellschaft fur Medizinische Physik and the European Federation of Organisations for Medical Physics.

[15]  Armand Wirgin,et al.  The inverse crime , 2004, math-ph/0401050.

[16]  Sean S. B. Moore,et al.  FFTs for the 2-Sphere-Improvements and Variations , 1996 .

[17]  M. N. Morrow,et al.  The Effects of Thoracic Inhomogeneities on the Relationship Between Epicardial and Torso Potentials , 1986, IEEE Transactions on Biomedical Engineering.

[18]  Stephen P. Boyd,et al.  Control applications of nonlinear convex programming , 1998 .

[19]  Cristian Lorenz,et al.  Generation of Point-Based 3D Statistical Shape Models for Anatomical Objects , 2000, Comput. Vis. Image Underst..

[20]  Samuli Siltanen,et al.  Numerical method for finding the convex hull of an inclusion in conductivity from boundary measurements , 2000 .

[21]  Fetsje Bijma,et al.  In vivo measurement of the brain and skull resistivities using an EIT-based method and realistic models for the head , 2003, IEEE Transactions on Biomedical Engineering.

[22]  David P. Dobkin,et al.  The quickhull algorithm for convex hulls , 1996, TOMS.

[23]  Thomas F. Coleman,et al.  An Interior Trust Region Approach for Nonlinear Minimization Subject to Bounds , 1993, SIAM J. Optim..

[24]  Jan C. de Munck,et al.  The boundary element method in the forward and inverse problem of electrical impedance tomography , 2000, IEEE Transactions on Biomedical Engineering.

[25]  T. Oostendorp,et al.  Interpolation on a triangulated 3D surface , 1989 .

[26]  David Isaacson,et al.  Electrode boundary conditions and experimental validation for BEM-based EIT forward and inverse solutions , 2006, IEEE Transactions on Medical Imaging.

[27]  William R B Lionheart,et al.  A Matlab toolkit for three-dimensional electrical impedance tomography: a contribution to the Electrical Impedance and Diffuse Optical Reconstruction Software project , 2002 .

[28]  D. Isaacson,et al.  A 3-D boundary element solution to the forward problem of electrical impedance tomography , 2004, The 26th Annual International Conference of the IEEE Engineering in Medicine and Biology Society.

[29]  Carl de Boor,et al.  Efficient Computer Manipulation of Tensor Products , 1979, TOMS.

[30]  David Isaacson,et al.  Electrical Impedance Tomography , 1999, SIAM Rev..

[31]  S R Arridge,et al.  Simultaneous reconstruction of internal tissue region boundaries and coefficients in optical diffusion tomography , 2000, Physics in medicine and biology.

[32]  Martin Brühl,et al.  Explicit Characterization of Inclusions in Electrical Impedance Tomography , 2001, SIAM J. Math. Anal..

[33]  Gene H. Golub,et al.  Generalized cross-validation as a method for choosing a good ridge parameter , 1979, Milestones in Matrix Computation.

[34]  John A. Nelder,et al.  A Simplex Method for Function Minimization , 1965, Comput. J..

[35]  E. Miller,et al.  A shape reconstruction method for electromagnetic tomography using adjoint fields and level sets , 2000 .

[36]  Nadia Magnenat-Thalmann,et al.  Surface from Scattered Points: A Brief Survey of Recent Developments , 2005 .

[37]  Dianne P. O'Leary,et al.  The Use of the L-Curve in the Regularization of Discrete Ill-Posed Problems , 1993, SIAM J. Sci. Comput..

[38]  T. Savolainen,et al.  Estimating shapes and free surfaces with electrical impedance tomography , 2004 .

[39]  J. D. Munck,et al.  The electric resistivity of human tissues (100 Hz-10 MHz): a meta-analysis of review studies. , 1999, Physiological measurement.

[40]  David Isaacson,et al.  3-D Electrical Impedance Tomography for Piecewise Constant Domains With Known Internal Boundaries , 2007, IEEE Transactions on Biomedical Engineering.