Comparison of algorithms for non-linear inverse 3D electrical tomography reconstruction.

Non-linear electrical impedance tomography reconstruction algorithms usually employ the Newton-Raphson iteration scheme to image the conductivity distribution inside the body. For complex 3D problems, the application of this method is not feasible any more due to the large matrices involved and their high storage requirements. In this paper we demonstrate the suitability of an alternative conjugate gradient reconstruction algorithm for 3D tomographic imaging incorporating adaptive mesh refinement and requiring less storage space than the Newton-Raphson scheme. We compare the reconstruction efficiency of both algorithms for a simple 3D head model. The results show that an increase in speed of about 30% is achievable with the conjugate gradient-based method without loss of accuracy.

[1]  J.P. Kaipio,et al.  Three-dimensional electrical impedance tomography based on the complete electrode model , 1999, IEEE Transactions on Biomedical Engineering.

[2]  R. Dutton,et al.  Delaunay triangulation and 3D adaptive mesh generation , 1997 .

[3]  Richard Bayford,et al.  Spherical and realistic head models in electrical impedance tomography , 2000 .

[4]  L M Heikkinen,et al.  A MATLAB package for the EIDORS project to reconstruct two-dimensional EIT images. , 2001, Physiological measurement.

[5]  D. C. Barber,et al.  Three-dimensional electrical impedance tomography , 1996, Nature.

[6]  Jürgen-Friedrich Hake.  Parallel algorithms for matrix operations and their performance on multiprocessor systems , 1992 .

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

[8]  Marc Molinari,et al.  Optimal imaging with adaptive mesh refinement in electrical impedance tomography. , 2002, Physiological measurement.

[9]  B H Blott,et al.  High fidelity imaging and high performance computing in nonlinear EIT. , 2000, Physiological measurement.

[10]  Ds Holder,et al.  Biomedical applications of EIT: a critical review , 1993 .

[11]  D. Isaacson,et al.  Electrical impedance tomography of complex conductivity distributions with noncircular boundary , 1997, IEEE Transactions on Biomedical Engineering.

[12]  M Molinari,et al.  Adaptive mesh refinement techniques for electrical impedance tomography. , 2001, Physiological measurement.

[13]  B H Blott,et al.  Nonlinear reconstruction constrained by image properties in electrical impedance tomography. , 1998, Physics in medicine and biology.

[14]  David S. Burnett,et al.  Finite Element Analysis: From Concepts to Applications , 1987 .

[15]  David P. Dobkin,et al.  The quickhull algorithm for convex hulls , 1996, TOMS.

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

[17]  S Arridge,et al.  A gradient-based optimisation scheme foroptical tomography. , 1998, Optics express.

[18]  Simon J. Cox,et al.  Efficient Non-Linear 3D Electrical Tomography Reconstruction , 2001 .

[19]  Gene H. Golub,et al.  Some History of the Conjugate Gradient and Lanczos Algorithms: 1948-1976 , 1989, SIAM Rev..

[20]  J. Shewchuk.  An Introduction to the Conjugate Gradient Method Without the Agonizing Pain , 1994 .

[21]  K. Foster,et al.  Dielectric properties of tissues and biological materials: a critical review. , 1989, Critical reviews in biomedical engineering.

[22]  A. J. Payne,et al.  Parallel Numerical Modelling of Ice Flow in Antarctica , 2000, PDPTA.