EIT Forward Problem Parallel Simulation Environment with Anisotropic Tissue and Realistic Electrode Models

Electrical impedance tomography (EIT) is an imaging technology based on impedance measurements. To retrieve meaningful insights from these measurements, EIT relies on detailed knowledge of the underlying electrical properties of the body. This is obtained from numerical models of current flows therein. The nonhomogeneous and anisotropic electric properties of human tissues make accurate modeling and simulation very challenging, leading to a tradeoff between physical accuracy and technical feasibility, which at present severely limits the capabilities of EIT. This work presents a complete algorithmic flow for an accurate EIT modeling environment featuring high anatomical fidelity with a spatial resolution equal to that provided by an MRI and a novel realistic complete electrode model implementation. At the same time, we demonstrate that current graphics processing unit (GPU)-based platforms provide enough computational power that a domain discretized with five million voxels can be numerically modeled in about 30 s.

[1]  R. Freund,et al.  QMR: a quasi-minimal residual method for non-Hermitian linear systems , 1991 .

[2]  C Gabriel,et al.  The dielectric properties of biological tissues: I. Literature survey. , 1996, Physics in medicine and biology.

[3]  R. W. Lau,et al.  The dielectric properties of biological tissues: II. Measurements in the frequency range 10 Hz to 20 GHz. , 1996, Physics in medicine and biology.

[4]  R. W. Lau,et al.  The dielectric properties of biological tissues: III. Parametric models for the dielectric spectrum of tissues. , 1996, Physics in medicine and biology.

[5]  William R B Lionheart,et al.  A Matlab toolkit for three-dimensional electrical impedance tomography: a contribution to the Electrical Impedance and Diffuse Optical Reconstruction Software project , 2002 .

[6]  Richard M. Leahy,et al.  BrainSuite: An Automated Cortical Surface Identification Tool , 2000, MICCAI.

[7]  David Isaacson,et al.  Electrical Impedance Tomography , 2002, IEEE Trans. Medical Imaging.

[8]  R. R. Harrison,et al.  A low-power low-noise CMOS amplifier for neural recording applications , 2003, IEEE J. Solid State Circuits.

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

[10]  Andreas Schulze-Bonhage,et al.  How large must an epileptic focus be to cause an electrographic status epilepticus—a case report , 2004, Clinical Neurophysiology.

[11]  Daniel Rueckert,et al.  Image Processing Workflows Using A Grid Enabled Image Database , 2004 .

[12]  Lior Horesh,et al.  Stroke type differentiation by multi-frequency electrical impedance tomography—a feasibility study , 2005 .

[13]  J. Haueisen,et al.  Role of Soft Bone, CSF and Gray Matter in EEG Simulations , 2003, Brain Topography.

[14]  D. Djajaputra Electrical Impedance Tomography: Methods, History and Applications , 2005 .

[15]  P. Fillard,et al.  MEDINRIA : DT-MRI PROCESSING AND VISUALIZATION SOFTWARE , 2006 .

[16]  William R B Lionheart,et al.  Uses and abuses of EIDORS: an extensible software base for EIT , 2006, Physiological measurement.

[17]  Peter Andersen,et al.  9.4T human MRI: Preliminary results , 2006, Magnetic resonance in medicine.

[18]  Xavier Tricoche,et al.  Influence of tissue conductivity anisotropy on EEG/MEG field and return current computation in a realistic head model: A simulation and visualization study using high-resolution finite element modeling , 2006, NeuroImage.

[19]  Z.J. Koles,et al.  A High-Resolution Anisotropic Finite-Volume Head Model for EEG Source Analysis , 2006, 2006 International Conference of the IEEE Engineering in Medicine and Biology Society.

[20]  R H Bayford,et al.  Bioimpedance tomography (electrical impedance tomography). , 2006, Annual review of biomedical engineering.

[21]  Jeff H. Duyn,et al.  High-field MRI of brain cortical substructure based on signal phase , 2007, Proceedings of the National Academy of Sciences.

[22]  J. L. Davidson,et al.  Predicted EIT current densities in the brain using a 3D anatomically realistic model of the head , 2007 .

[23]  Bong Seok Kim,et al.  The Complete Electrode Model For Imaging and Electrode Contact Compensation in Electrical Impedance Tomography , 2007, 2007 29th Annual International Conference of the IEEE Engineering in Medicine and Biology Society.

[24]  Cícero R. de Lima,et al.  Electrical impedance tomography through constrained sequential linear programming: a topology optimization approach , 2007 .

[25]  R. Sadleir,et al.  Modeling Skull Electrical Properties , 2007, Annals of Biomedical Engineering.

[26]  Á. Pascual-Leone,et al.  Technology Insight: noninvasive brain stimulation in neurology—perspectives on the therapeutic potential of rTMS and tDCS , 2007, Nature Clinical Practice Neurology.

[27]  V. Walsh,et al.  Diffusion tensor MRI-based estimation of the influence of brain tissue anisotropy on the effects of transcranial magnetic stimulation , 2007, NeuroImage.

[28]  David Atkinson,et al.  Use of anisotropic modelling in electrical impedance tomography; Description of method and preliminary assessment of utility in imaging brain function in the adult human head , 2008, NeuroImage.

[29]  Erik Lindholm,et al.  NVIDIA Tesla: A Unified Graphics and Computing Architecture , 2008, IEEE Micro.

[30]  William R B Lionheart,et al.  Simple FEMs aren’t as good as we thought: experiences developing EIDORS v3.3 , 2008 .

[31]  Keith D. Paulsen,et al.  3D electric impedance tomography reconstruction on multi-core computing platforms , 2008, 2008 30th Annual International Conference of the IEEE Engineering in Medicine and Biology Society.

[32]  Thomas Jost,et al.  An efficient multi-algorithms sparse linear solver for GPUs , 2009, PARCO.

[33]  Roberto Guerrieri,et al.  Triangular matrix inversion on Graphics Processing Unit , 2009, Proceedings of the Conference on High Performance Computing Networking, Storage and Analysis.

[34]  Paulius Micikevicius,et al.  3D finite difference computation on GPUs using CUDA , 2009, GPGPU-2.

[35]  Wolfgang Straßer,et al.  A Parallel Preconditioned Conjugate Gradient Solver for the Poisson Problem on a Multi-GPU Platform , 2010, 2010 18th Euromicro Conference on Parallel, Distributed and Network-based Processing.

[36]  Michal M. Okoniewski,et al.  Adaptive Mesh Refinement Techniques for 3-D Skin Electrode Modeling , 2010, IEEE Transactions on Biomedical Engineering.

[37]  Allen D. Malony,et al.  A fast BiCG solver for the isotropic Poisson equation in the forward EIT problem in cylinder phantoms , 2010 .