Matching pursuit for imaging high-contrast conductivity

We show that imaging an isotropic, high contrast, conductive medium is asymptotically equivalent to the identification of a unique resistor network, given measurements of currents and voltages at the boundary. We show that a matching pursuit approach can be used effectively towards the numerical solution of the high-contrast imaging problem, if the library of functions is constructed carefully and in accordance with the asymptotic theory. We also show how other libraries of functions that at first glance seem reasonable, in fact, do not work well. When the contrast in the conductivity is not so high, we show that wavelets can be used, especially non-orthogonal wavelet libraries. However, the library of functions that is based on the high-contrast asymptotic theory is more robust, even for intermediate contrasts, and especially so in the presence of noise.

[1]  Alan C. Tripp,et al.  INVERSE CONDUCTIVITY PROBLEM FOR INACCURATE MEASUREMENTS , 1996 .

[2]  Willis J. Tompkins,et al.  Comparing Reconstruction Algorithms for Electrical Impedance Tomography , 1987, IEEE Transactions on Biomedical Engineering.

[3]  Kohn,et al.  Variational constraints for electrical-impedance tomography. , 1990, Physical review letters.

[4]  H. Akaike.  A new look at the statistical model identification , 1974 .

[5]  Joseph B. Keller,et al.  Effective conductivity of periodic composites composed of two very unequal conductors , 1987 .

[6]  Fadil Santosa,et al.  Resolution and Stability Analysis of an Inverse Problem in Electrical Impedance Tomography: Dependence on the Input Current Patterns , 1994, SIAM J. Appl. Math..

[7]  Michael Vogelius,et al.  A backprojection algorithm for electrical impedance imaging , 1990 .

[8]  P. Zweifel.  Advanced Mathematical Methods for Scientists and Engineers , 1980 .

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

[10]  Randolph E. Bank,et al.  PLTMG - a software package for solving elliptic partial differential equations: users' guide 8.0 , 1998, Software, environments, tools.

[11]  J. Sylvester,et al.  A global uniqueness theorem for an inverse boundary value problem , 1987 .

[12]  V. Isakov.  Uniqueness and stability in multi-dimensional inverse problems , 1993 .

[13]  G. Uhlmann.  Inverse boundary value problems for partial differential equations. , 1998 .

[14]  Robert V. Kohn,et al.  Numerical implementation of a variational method for electrical impedance tomography , 1990 .

[15]  John Sylvester,et al.  The Dirichlet to Neumann map and applications , 1989 .

[16]  Joseph B. Keller,et al.  Conductivity of a Medium Containing a Dense Array of Perfectly Conducting Spheres or Cylinders or Nonconducting Cylinders , 1963 .

[17]  S. Kozlov,et al.  Geometric aspects of averaging , 1989 .

[18]  G. Papanicolaou,et al.  High-contrast impedance tomography , 1996 .

[19]  R. Kohn,et al.  Determining conductivity by boundary measurements II. Interior results , 1985 .

[20]  Liliana Borcea.  Direct and Inverse Problems for Transport in High Contrast Media , 1996 .

[21]  Hans D. Mittelmann,et al.  PLTMG: A Software Package for Solving Elliptic Partial Differential Equations (Randolph E. Bank) , 1994, SIAM Review.

[22]  Liliana Borcea,et al.  A hybrid numerical method for high contrast conductivity problems , 1997 .

[23]  Jorge P. Zubelli,et al.  Diffuse tomography: computational aspects of the isotropic case , 1992 .

[24]  M. Neuman,et al.  Impedance computed tomography algorithm and system. , 1985, Applied optics.

[25]  G. Schwarz.  Estimating the Dimension of a Model , 1978 .

[26]  Edward B Curtis,et al.  Determining the resistors in a network , 1990 .

[27]  Liliana Borcea,et al.  Network Approximation for Transport Properties of High Contrast Materials , 1998, SIAM J. Appl. Math..

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

[29]  David Isaacson,et al.  Layer stripping: a direct numerical method for impedance imaging , 1991 .

[30]  A velocity inversion problem involving an unknown source , 1990 .

[31]  R. Kohn,et al.  Relaxation of a variational method for impedance computed tomography , 1987 .

[32]  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..

[33]  Liliana Borcea,et al.  Asymptotic Analysis of Quasi-Static Transport in High Contrast Conductive Media , 1998, SIAM J. Appl. Math..

[34]  Per Christian Hansen,et al.  Analysis of Discrete Ill-Posed Problems by Means of the L-Curve , 1992, SIAM Rev..

[35]  M. Victor Wickerhauser,et al.  Wavelets: Algorithms and Applications (Yves Meyer) , 1994, SIAM Rev..

[36]  Robert V. Kohn,et al.  Determining conductivity by boundary measurements , 1984 .

[37]  Stéphane Mallat,et al.  Matching pursuits with time-frequency dictionaries , 1993, IEEE Trans. Signal Process..

[38]  J. Berryman.  Convexity properties of inverse problems with variational constraints , 1991 .

[39]  G. Batchelor,et al.  Thermal or electrical conduction through a granular material , 1977, Proceedings of the Royal Society of London. A. Mathematical and Physical Sciences.