Nonlinear Inversion from Partial EIT Data: Computational Experiments

Electrical impedance tomography (EIT) is a non-invasive imaging method in which an unknown physical body is probed with electric currents applied on the boundary, and the internal conductivity distribution is recovered from the measured boundary voltage data. The reconstruction task is a nonlinear and ill-posed inverse problem, whose solution calls for special regularized algorithms, such as D-bar methods which are based on complex geometrical optics solutions (CGOs). In many applications of EIT, such as monitoring the heart and lungs of unconscious intensive care patients or locating the focus of an epileptic seizure, data acquisition on the entire boundary of the body is impractical, restricting the boundary area available for EIT measurements. An extension of the D-bar method to the case when data is collected only on a subset of the boundary is studied by computational simulation. The approach is based on solving a boundary integral equation for the traces of the CGOs using localized basis functions (Haar wavelets). The numerical evidence suggests that the D-bar method can be applied to partial-boundary data in dimension two and that the traces of the partial data CGOs approximate the full data CGO solutions on the available portion of the boundary, for the necessary small $k$ frequencies.Comment: 24 pages, 12 figure

Spatiospectral concentration in the Cartesian plane

We pose and solve the analogue of Slepian's time-frequency concentration problem in the two-dimensional plane, for applications in the natural sciences. We determine an orthogonal family of strictly bandlimited functions that are optimally concentrated within a closed region of the plane, or, alternatively, of strictly spacelimited functions that are optimally concentrated in the Fourier domain. The Cartesian Slepian functions can be found by solving a Fredholm integral equation whose associated eigenvalues are a measure of the spatiospectral concentration. Both the spatial and spectral regions of concentration can, in principle, have arbitrary geometry. However, for practical applications of signal representation or spectral analysis such as exist in geophysics or astronomy, in physical space irregular shapes, and in spectral space symmetric domains will usually be preferred. When the concentration domains are circularly symmetric in both spaces, the Slepian functions are also eigenfunctions of a Sturm-Liouville operator, leading to special algorithms for this case, as is well known. Much like their one-dimensional and spherical counterparts with which we discuss them in a common framework, a basis of functions that are simultaneously spatially and spectrally localized on arbitrary Cartesian domains will be of great utility in many scientific disciplines, but especially in the geosciences.Comment: 34 pages, 7 figures. In the press, International Journal on Geomathematics, April 14th, 201