Reducing computational costs in large scale 3D EIT by using a sparse Jacobian matrix with block-wise CGLS reconstruction

Electrical impedance tomography (EIT) is a fast and cost-effective technique to provide a tomographic conductivity image of a subject from boundary current–voltage data. This paper proposes a time and memory efficient method for solving a large scale 3D EIT inverse problem using a parallel conjugate gradient (CG) algorithm. The 3D EIT system with a large number of measurement data can produce a large size of Jacobian matrix; this could cause difficulties in computer storage and the inversion process. One of challenges in 3D EIT is to decrease the reconstruction time and memory usage, at the same time retaining the image quality. Firstly, a sparse matrix reduction technique is proposed using thresholding to set very small values of the Jacobian matrix to zero. By adjusting the Jacobian matrix into a sparse format, the element with zeros would be eliminated, which results in a saving of memory requirement. Secondly, a block-wise CG method for parallel reconstruction has been developed. The proposed method has been tested using simulated data as well as experimental test samples. Sparse Jacobian with a block-wise CG enables the large scale EIT problem to be solved efficiently. Image quality measures are presented to quantify the effect of sparse matrix reduction in reconstruction results.

[1]  Yue Xiu-li,et al.  Multifunctional magnetic nanoparticles for magnetic resonance image-guided photothermal therapy for cancer , 2014 .

[2]  Andy Adler,et al.  Addressing the computational cost of large EIT solutions , 2012, Physiological measurement.

[3]  William R B Lionheart,et al.  GREIT: a unified approach to 2D linear EIT reconstruction of lung images , 2009, Physiological measurement.

[4]  A. Adler,et al.  A measure of the information content of EIT data , 2008, Physiological Measurement.

[5]  Manuchehr Soleimani,et al.  A Narrow-Band Level Set Method Applied to EIT in Brain for Cryosurgery Monitoring , 2006, IEEE Transactions on Biomedical Engineering.

[6]  W. Lionheart.  EIT reconstruction algorithms: pitfalls, challenges and recent developments , 2003, Physiological measurement.

[7]  Yousef Saad,et al.  Iterative methods for sparse linear systems , 2003 .

[8]  Hugh McCann,et al.  Krylov subspace iterative techniques: on the detection of brain activity with electrical impedance tomography , 2002, IEEE Transactions on Medical Imaging.

[9]  Wei Wang,et al.  The number of electrodes and basis functions in EIT image reconstruction , 2002, Physiological measurement.

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

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

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

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

[14]  Brian H. Brown,et al.  Imaging spatial distributions of resistivity using applied potential tomography , 1983 .

[15]  Manuchehr Soleimani,et al.  Three-dimensional magnetic induction tomography imaging using a matrix free krylov subspace inversion algorithm , 2012 .

[16]  Manuchehr Soleimani,et al.  A trust region subproblem for 3D electrical impedance tomography inverse problem using experimental data , 2009 .