Study of noise effects in electrical impedance tomography with resistor networks

We present a study of the numerical solution of the two dimensional electrical impedance tomography problem, with noisy measurements of the Dirichlet to Neumann map. The inversion uses parametrizations of the conductivity on optimal grids. The grids are optimal in the sense that finite volume discretizations on them give spectrally accurate approximations of the Dirichlet to Neumann map. The approximations are Dirichlet to Neumann maps of special resistor networks, that are uniquely recoverable from the measurements. Inversion on optimal grids has been proposed and analyzed recently, but the study of noise effects on the inversion has not been carried out. In this paper we present a numerical study of both the linearized and the nonlinear inverse problem. We take three different parametrizations of the unknown conductivity, with the same number of degrees of freedom. We obtain that the parametrization induced by the inversion on optimal grids is the most efficient of the three, because it gives the smallest standard deviation of the maximum a posteriori estimates of the conductivity, uniformly in the domain. For the nonlinear problem we compute the mean and variance of the maximum a posteriori estimates of the conductivity, on optimal grids. For small noise, we obtain that the estimates are unbiased and their variance is very close to the optimal one, given by the Cramer-Rao bound. For larger noise we use regularization and quantify the trade-off between reducing the variance and introducing bias in the solution. Both the full and partial measurement setups are considered.

[1]  V. Druskin,et al.  Resistor network approaches to electrical impedance tomography , 2011, 1107.0343.

[2]  Liliana Borcea,et al.  Pyramidal resistor networks for electrical impedance tomography with partial boundary measurements , 2010 .

[3]  Alexander V. Mamonov,et al.  Resistor networks and optimal grids for the numerical solution of electrical impedance tomography with partial boundary measurements , 2010 .

[4]  Matti Lassas,et al.  REGULARIZED D-BAR METHOD FOR THE INVERSE CONDUCTIVITY PROBLEM , 2009 .

[5]  Masahiro Yamamoto,et al.  Global uniqueness from partial Cauchy data in two dimensions , 2008, 0810.2286.

[6]  Liliana Borcea,et al.  Electrical impedance tomography with resistor networks , 2008 .

[7]  F. G. Vasquez On the parameterization of ill-posed inverse problems arising from elliptic partial differential equations , 2006 .

[8]  Faming Liang,et al.  Statistical and Computational Inverse Problems , 2006, Technometrics.

[9]  Liliana Borcea,et al.  On the continuum limit of a discrete inverse spectral problem on optimal finite difference grids , 2005 .

[10]  Sergio Vessella,et al.  Lipschitz stability for the inverse conductivity problem , 2005, Adv. Appl. Math..

[11]  Thomas A. Manteuffel,et al.  First-Order System Least Squares and Electrical Impedance Tomography , 2004, SIAM J. Numer. Anal..

[12]  M. Lassas,et al.  Calderóns' Inverse Problem for Anisotropic Conductivity in the Plane , 2004, math/0401410.

[13]  V. Druskin,et al.  Optimal finite difference grids for direct and inverse Sturm-Liouville problems , 2002 .

[14]  Jérôme Jaffré,et al.  Refinement and coarsening indicators for adaptive parametrization: application to the estimation of hydraulic transmissivities , 2002 .

[15]  Niculae Mandache,et al.  Exponential instability in an inverse problem for the Schrodinger equation , 2001 .

[16]  Juan Antonio Barceló,et al.  Stability of the Inverse Conductivity Problem in the Plane for Less Regular Conductivities , 2001 .

[17]  E. Somersalo,et al.  Statistical inversion and Monte Carlo sampling methods in electrical impedance tomography , 2000 .

[18]  Vladimir Druskin,et al.  Gaussian spectral rules for second order finite-difference schemes , 2000, Numerical Algorithms.

[19]  Vladimir Druskin,et al.  Optimal finite difference grids and rational approximations of the square root I. Elliptic problems , 2000 .

[20]  David V. Ingerman,et al.  Discrete and Continuous Dirichlet-to-Neumann Maps in the Layered Case , 2000, SIAM J. Math. Anal..

[21]  Vladimir Druskin,et al.  Application of the Difference Gaussian Rules to Solution of Hyperbolic Problems , 2000 .

[22]  Vladimir Druskin,et al.  Gaussian Spectral Rules for the Three-Point Second Differences: I. A Two-Point Positive Definite Problem in a Semi-Infinite Domain , 1999, SIAM J. Numer. Anal..

[23]  James A. Morrow,et al.  Circular planar graphs and resistor networks , 1998 .

[24]  H. Engl,et al.  Regularization of Inverse Problems , 1996 .

[25]  Y. C. Verdière,et al.  Reseaux électriques planaires II , 1994 .

[26]  C. Pommerenke Boundary Behaviour of Conformal Maps , 1992 .

[27]  E. Somersalo,et al.  Existence and uniqueness for electrode models for electric current computed tomography , 1992 .

[28]  B. Fitzpatrick Bayesian analysis in inverse problems , 1991 .

[29]  V. Druskin,et al.  Circular resistor networks for electrical impedance tomography with partial boundary measurements , 2010 .

[30]  H. Ben Ameur,et al.  Regularization of parameter estimation by adaptive discretization using refinement and coarsening indicators , 2002 .

[31]  Liliana Borcea,et al.  INVERSE PROBLEMS PII: S0266-5611(02)33630-X Optimal finite difference grids for direct and inverse Sturm–Liouville problems , 2002 .

[32]  David V. Ingerman,et al.  On a characterization of the kernel of the Dirichlet-to-Neumann map for a planar region , 1998 .

[33]  Gunther Uhlmann,et al.  Uniqueness in the inverse conductivity problem for nonsmooth conductivities in two dimensions , 1997 .

[34]  A. Nachman,et al.  Global uniqueness for a two-dimensional inverse boundary value problem , 1996 .

[35]  F. Y. Edgeworth,et al.  The theory of statistics , 1996 .

[36]  James A. Morrow,et al.  Finding the conductors in circular networks from boundary measurements , 1994 .

[37]  Y. C. Verdière,et al.  Réseaux électriques planaires I , 1994 .

[38]  Giovanni Alessandrini,et al.  Stable determination of conductivity by boundary measurements , 1988 .

[39]  D. Isaacson Distinguishability of Conductivities by Electric Current Computed Tomography , 1986, IEEE Transactions on Medical Imaging.

[40]  A. Seagar,et al.  Probing with low frequency electric currents. , 1983 .