Concurrency in electrical neuroinformatics: parallel computation for studying the volume conduction of brain electrical fields in human head tissues

Advances in human brain neuroimaging for high‐temporal and high‐spatial resolutions will depend on localization of electroencephalography (EEG) signals to their cortex sources. The source localization inverse problem is inherently ill‐posed and depends critically on the modeling of human head electromagnetics. We present a systematic methodology to analyze the main factors and parameters that affect the EEG source‐mapping accuracy. These factors are not independent, and their effect must be evaluated in a unified way. To do so requires significant computational capabilities to explore the problem landscape, quantify uncertainty effects, and evaluate alternative algorithms. Bringing high‐performance computing to this domain is necessary to open new avenues for neuroinformatics research. The head electromagnetics forward problem is the heart of the source localization inverse. We present two parallel algorithms to address tissue inhomogeneity and impedance anisotropy. Highly accurate head modeling environments will enable new research and clinical neuroimaging applications. Cortex‐localized dense‐array EEG analysis is the next‐step in neuroimaging domains such as early childhood reading, understanding of resting‐state brain networks, and models of full brain function. Therapeutic treatments based on neurostimulation will also depend significantly on high‐performance computing integration. Copyright © 2015 John Wiley & Sons, Ltd.

[1]  K. Foster,et al.  Dielectric Properties of Fluid-Saturated Bone - The Effect of Variation in Conductivity of Immersion Fluid , 1984, IEEE Transactions on Biomedical Engineering.

[2]  L. Parra,et al.  Optimized multi-electrode stimulation increases focality and intensity at target , 2011, Journal of neural engineering.

[3]  Richard Bowtell,et al.  Combining EEG and fMRI. , 2011, Methods in molecular biology.

[4]  H. Spekreijse,et al.  Mathematical dipoles are adequate to describe realistic generators of human brain activity , 1988, IEEE Transactions on Biomedical Engineering.

[5]  G. Huiskamp,et al.  The need for correct realistic geometry in the inverse EEG problem , 1999, IEEE Transactions on Biomedical Engineering.

[6]  P Berg,et al.  Multiple source analysis of interictal spikes: goals, requirements, and clinical value. , 1999, Journal of clinical neurophysiology : official publication of the American Electroencephalographic Society.

[7]  E. Frank.  Electric Potential Produced by Two Point Current Sources in a Homogeneous Conducting Sphere , 1952 .

[8]  R D Pascual-Marqui,et al.  Standardized low-resolution brain electromagnetic tomography (sLORETA): technical details. , 2002, Methods and findings in experimental and clinical pharmacology.

[9]  M. Scherg,et al.  Two bilateral sources of the late AEP as identified by a spatio-temporal dipole model. , 1985, Electroencephalography and clinical neurophysiology.

[10]  R. Pascual-Marqui.  Review of methods for solving the EEG inverse problem , 1999 .

[11]  J. D. Munck,et al.  A fast method to compute the potential in the multisphere model (EEG application) , 1993, IEEE Transactions on Biomedical Engineering.

[12]  Allen D. Malony,et al.  Next-generation human brain neuroimaging and the role of high-performance computing , 2013, 2013 International Conference on High Performance Computing & Simulation (HPCS).

[13]  Barry D. Van Veen,et al.  Statistical performance analysis of signal variance-based dipole models for MEG/EEG source localization and detection , 2003, IEEE Transactions on Biomedical Engineering.

[14]  S R Arridge,et al.  Recent advances in diffuse optical imaging , 2005, Physics in medicine and biology.

[15]  Kevin Whittingstall,et al.  Effects of dipole position, orientation and noise on the accuracy of EEG source localization , 2003, Biomedical engineering online.

[16]  D. Tucker,et al.  EEG coherency. I: Statistics, reference electrode, volume conduction, Laplacians, cortical imaging, and interpretation at multiple scales. , 1997, Electroencephalography and clinical neurophysiology.

[17]  D. Tucker.  Spatial sampling of head electrical fields: the geodesic sensor net. , 1993, Electroencephalography and clinical neurophysiology.

[18]  J. D. Munck.  The estimation of time varying dipoles on the basis of evoked potentials. , 1990 .

[19]  R Hoekema,et al.  Multigrid solution of the potential field in modeling electrical nerve stimulation. , 1998, Computers and biomedical research, an international journal.

[20]  P. Nicholson,et al.  Specific impedance of cerebral white matter. , 1965, Experimental neurology.

[21]  Á. Pascual-Leone,et al.  Noninvasive human brain stimulation. , 2007, Annual review of biomedical engineering.

[22]  L. Cohen,et al.  Transcranial direct current stimulation: State of the art 2008 , 2008, Brain Stimulation.

[23]  S. K. Law,et al.  Thickness and resistivity variations over the upper surface of the human skull , 2005, Brain Topography.

[24]  D. Tucker,et al.  Spatial sampling and filtering of EEG with spline Laplacians to estimate cortical potentials , 2005, Brain Topography.

[25]  Guang Cheng,et al.  Correlation Between Structure and Resistivity Variations of the Live Human Skull , 2008, IEEE Transactions on Biomedical Engineering.

[26]  Jing Li,et al.  Effects of holes on EEG forward solutions using a realistic geometry head model , 2007, Journal of neural engineering.

[27]  Jens Haueisen,et al.  Evaluation of the distortion of EEG signals caused by a hole in the skull mimicking the fontanel in the skull of human neonates , 2005, Clinical Neurophysiology.

[28]  D. A. Driscoll,et al.  EEG electrode sensitivity--an application of reciprocity. , 1969, IEEE transactions on bio-medical engineering.

[29]  C. D. Gelatt,et al.  Optimization by Simulated Annealing , 1983, Science.

[30]  R. Salmelin,et al.  Global optimization in the localization of neuromagnetic sources , 1998, IEEE Transactions on Biomedical Engineering.

[31]  J. J. Douglas.  Alternating direction methods for three space variables , 1962 .

[32]  Antoine Rémond,et al.  Methods of Analysis of Brain Electrical and Magnetic Signals , 1987 .

[33]  Chun-Chuan Chen,et al.  Accelerating Computation of DCM for ERP in MATLAB by External Function Calls to the GPU , 2013, PloS one.

[34]  Tadashi Yamazaki,et al.  Realtime cerebellum: A large-scale spiking network model of the cerebellum that runs in realtime using a graphics processing unit , 2013, Neural Networks.

[35]  R. Leahy,et al.  Mapping human brain function with MEG and EEG: methods and validation , 2004, NeuroImage.

[36]  W. Sutherling,et al.  Conductivities of Three-Layer Live Human Skull , 2004, Brain Topography.

[37]  J. D. Munck.  The potential distribution in a layered anisotropic spheroidal volume conductor , 1988 .

[38]  Alon Korngreen,et al.  Accelerating compartmental modeling on a graphical processing unit , 2013, Front. Neuroinform..

[39]  Anders Eklund,et al.  Medical image processing on the GPU - Past, present and future , 2013, Medical Image Anal..

[40]  E. R. Flynn,et al.  A Model for Frequency Dependence of Conductivities of the Live Human Skull , 2004, Brain Topography.

[41]  Chris R. Johnson,et al.  Volume Currents in Forward and Inverse Magnetoencephalographic Simulations Using Realistic Head Models , 2004, Annals of Biomedical Engineering.

[42]  Don M. Tucker,et al.  Regional head tissue conductivity estimation for improved EEG analysis , 2000, IEEE Transactions on Biomedical Engineering.

[43]  B. Vanrumste,et al.  Comparing iterative solvers for linear systems associated with the finite difference discretisation of the forward problem in electro-encephalographic source analysis , 2006, Medical and Biological Engineering and Computing.

[44]  Simon K. Warfield,et al.  Cortical Graph Smoothing: A Novel Method for Exploiting DWI-Derived Anatomical Brain Connectivity to Improve EEG Source Estimation , 2013, IEEE Transactions on Medical Imaging.

[45]  Allen D. Malony,et al.  Neuroanatomical segmentation in mri exploiting a priori knowledge , 2007 .

[46]  Daniel C. Javitt,et al.  Right hemisphere control of visuospatial attention: line-bisection judgments evaluated with high-density electrical mapping and source analysis☆ , 2003, NeuroImage.

[47]  J. C. DE MUNCK,et al.  A Parametric Method to Resolve the Ill‐Posed Nature of the EIT Reconstruction Problem: A Simulation Study , 1999 .

[48]  R. V. Uitert,et al.  Can a Spherical Model Substitute for a Realistic Head Model in Forward and Inverse MEG Simulations ? , 2002 .

[49]  M. Murray,et al.  EEG source imaging , 2004, Clinical Neurophysiology.

[50]  Y. D'Asseler,et al.  A finite difference method with reciprocity used to incorporate anisotropy in electroencephalogram dipole source localization , 2005 .

[51]  J. P. Ary,et al.  Location of Sources of Evoked Scalp Potentials: Corrections for Skull and Scalp Thicknesses , 1981, IEEE Transactions on Biomedical Engineering.

[52]  Hong Bao,et al.  GPGPU-Aided Ensemble Empirical-Mode Decomposition for EEG Analysis During Anesthesia , 2010, IEEE Transactions on Information Technology in Biomedicine.

[53]  Richard M. Leahy,et al.  Electromagnetic brain mapping , 2001, IEEE Signal Process. Mag..

[54]  Bart Vanrumste,et al.  Journal of Neuroengineering and Rehabilitation Open Access Review on Solving the Inverse Problem in Eeg Source Analysis , 2022 .

[55]  D. Lehmann,et al.  Handbook of electroencephalography and clinical neurophysiology , 1976 .

[56]  J. Pernier,et al.  A systematic evaluation of the spherical model accuracy in EEG dipole localization. , 1997, Electroencephalography and clinical neurophysiology.

[57]  Thomas R. Knösche,et al.  Influence of anisotropic conductivity on EEG source reconstruction: investigations in a rabbit model , 2006, IEEE Transactions on Biomedical Engineering.

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

[59]  E. Harth,et al.  Electric Fields of the Brain: The Neurophysics of Eeg , 2005 .

[60]  James Kozloski,et al.  Self-referential forces are sufficient to explain different dendritic morphologies , 2013, Front. Neuroinform..

[61]  D. A. Driscoll,et al.  Current Distribution in the Brain From Surface Electrodes , 1968, Anesthesia and analgesia.

[62]  J. Gotman,et al.  Non-uniform spatial sampling in EEG source analysis , 2001, 2001 Conference Proceedings of the 23rd Annual International Conference of the IEEE Engineering in Medicine and Biology Society.

[63]  M. Peters,et al.  Volume conduction effects in EEG and MEG. , 1998, Electroencephalography and clinical neurophysiology.

[64]  Mariano Sigman,et al.  CUDAICA: GPU Optimization of Infomax-ICA EEG Analysis , 2012, Comput. Intell. Neurosci..

[65]  Carlos H. Muravchik,et al.  Analysis of parametric estimation of head tissue conductivities using Electrical Impedance Tomography , 2013, Biomed. Signal Process. Control..

[66]  Á. Pascual-Leone,et al.  A Controlled Clinical Trial of Cathodal DC Polarization in Patients with Refractory Epilepsy , 2006, Epilepsia.

[67]  Christian Windischberger,et al.  Toward discovery science of human brain function , 2010, Proceedings of the National Academy of Sciences.

[68]  Allen D. Malony,et al.  A 3D Vector-Additive Iterative Solver for the Anisotropic Inhomogeneous Poisson Equation in the Forward EEG problem , 2009, ICCS.

[69]  M. Scherg,et al.  Evoked dipole source potentials of the human auditory cortex. , 1986, Electroencephalography and clinical neurophysiology.

[70]  Marisa O. Hollinshead,et al.  The organization of the human cerebral cortex estimated by intrinsic functional connectivity. , 2011, Journal of neurophysiology.

[71]  B. Gordon,et al.  Target Optimization in Transcranial Direct Current Stimulation , 2012, Front. Psychiatry.

[72]  Allen D. Malony,et al.  Conductivity Analysis for High-Resolution EEG , 2008, 2008 International Conference on BioMedical Engineering and Informatics.

[73]  P. Nunez,et al.  Electric fields of the brain , 1981 .

[74]  Lei Xing,et al.  GPU computing in medical physics: a review. , 2011, Medical physics.

[75]  Allen D. Malony,et al.  Use of Parallel Simulated Annealing for Computational Modeling of Human Head Conductivity , 2007, International Conference on Computational Science.

[76]  Zhigang Li,et al.  Development/global validation of a 6-month-old pediatric head finite element model and application in investigation of drop-induced infant head injury , 2013, Comput. Methods Programs Biomed..

[77]  Diego Fasoli,et al.  Three Applications of GPU Computing in Neuroscience , 2012, Computing in Science & Engineering.

[78]  Riaz A Khan,et al.  SOURCE LOCALIZATION OF BRAIN ELECTRICAL ACTIVITY VIA TIME-FREQUENCY LINEARLY CONSTRAINED MINIMUM VARIANCE METHOD , 2005 .