Electrical Impedance Tomography (EIT) is an imaging technique which aims to the reconstruction of the spatial electrical conductivity distribution of a section of the human body. In this paper, in order to solve the EIT forward and inverse problems, a finite difference approach to the solution of Maxwell’s equations, namely the admittance method, and the Newton-Raphson algorithm have been employed, respectively. The Jacobian matrix involved in the Newton-Raphson (N&R) algorithm has been computed by using both the admittance matrix of the forward problem and the socalled sensitivity matrix. The reconstruction problem has been solved using the standard Tikhonov regularization with various choices of the regularization matrix. The obtained results, for three thorax models of increasing complexity, approximating a section of the human chest, show that the lowest reconstruction errors are usually obtained by choosing as regularization matrix a discrete form of the Laplacian operator.
Implementation of the Newton-Raphson and admittance methods for EIT / D., Romano; Pisa, Stefano; Piuzzi, Emanuele. - In: INTERNATIONAL JOURNAL OF BIOELECTROMAGNETISM. - ISSN 1456-7857. - ELETTRONICO. - 12:(2010), pp. 12-20.
Implementation of the Newton-Raphson and admittance methods for EIT
PISA, Stefano;PIUZZI, Emanuele
2010
Abstract
Electrical Impedance Tomography (EIT) is an imaging technique which aims to the reconstruction of the spatial electrical conductivity distribution of a section of the human body. In this paper, in order to solve the EIT forward and inverse problems, a finite difference approach to the solution of Maxwell’s equations, namely the admittance method, and the Newton-Raphson algorithm have been employed, respectively. The Jacobian matrix involved in the Newton-Raphson (N&R) algorithm has been computed by using both the admittance matrix of the forward problem and the socalled sensitivity matrix. The reconstruction problem has been solved using the standard Tikhonov regularization with various choices of the regularization matrix. The obtained results, for three thorax models of increasing complexity, approximating a section of the human chest, show that the lowest reconstruction errors are usually obtained by choosing as regularization matrix a discrete form of the Laplacian operator.I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.