We present a subspace based, structure preserving solution method for the problem of Electrical Impedance Tomography, where the conductivity inside a simply connected 2-dimensional domain is sought from noisy and incomplete boundary data. Unlike conventional output-least squares algorithms that can be regarded as minimizing a certain error norm, solutions are recovered here as the minimizers of a closely related residual norm problem. An iterative solution scheme is shown to lead to a sequence of sparse matrix subproblems, with conditioning far more favourable than typically observed in output-least squares. We find that these sparse subproblems demonstrate a particular form of displacement structure that can be further elaborated to finally arrive upon an efficient computational implementation. In the second part of this contribution, we investigate in more detail the different sparsity patterns deriving from the problem formulation, and illustrate how these connect with the original concept of displacement structure.
Dimiccoli, L, Truyen, B & Cornelis, J 2006, A subspace based structure preserving solution method for the Electrical Impedance Tomography problem: Part II. in N Mastronardi, M Van Barel & R Vandebril (eds), Proceedings International Workshop Numerical Linear Algebra in Signals and Systems. pp. 21-21, International Workshop Numerical Linear Algebra in Signals and Systems, Monopoli, Italy, 11/09/06. <https://people.cs.kuleuven.be/~raf.vandebril/bari2006/abstracts.pdf>
Dimiccoli, L., Truyen, B., & Cornelis, J. (2006). A subspace based structure preserving solution method for the Electrical Impedance Tomography problem: Part II. In N. Mastronardi, M. Van Barel, & R. Vandebril (Eds.), Proceedings International Workshop Numerical Linear Algebra in Signals and Systems (pp. 21-21) https://people.cs.kuleuven.be/~raf.vandebril/bari2006/abstracts.pdf
@inproceedings{fc61593e0d074a54b1e1c2e4ca1a9270,
title = "A subspace based structure preserving solution method for the Electrical Impedance Tomography problem: Part II",
abstract = "We present a subspace based, structure preserving solution method for the problem of Electrical Impedance Tomography, where the conductivity inside a simply connected 2-dimensional domain is sought from noisy and incomplete boundary data. Unlike conventional output-least squares algorithms that can be regarded as minimizing a certain error norm, solutions are recovered here as the minimizers of a closely related residual norm problem. An iterative solution scheme is shown to lead to a sequence of sparse matrix subproblems, with conditioning far more favourable than typically observed in output-least squares. We find that these sparse subproblems demonstrate a particular form of displacement structure that can be further elaborated to finally arrive upon an efficient computational implementation. In the second part of this contribution, we investigate in more detail the different sparsity patterns deriving from the problem formulation, and illustrate how these connect with the original concept of displacement structure.",
keywords = "Electrical Impedance Tomography, sparsity patterns, displacement structure",
author = "Luca Dimiccoli and Bart Truyen and Jan Cornelis",
year = "2006",
month = sep,
day = "11",
language = "English",
pages = "21--21",
editor = "Nicola Mastronardi and {Van Barel}, Marc and Raf Vandebril",
booktitle = "Proceedings International Workshop Numerical Linear Algebra in Signals and Systems",
note = "International Workshop Numerical Linear Algebra in Signals and Systems ; Conference date: 11-09-2006 Through 15-09-2006",
}