Realistic volumetric-approach to simulate transcranial electric stimulation—ROAST—a fully automated open-source pipeline

Research in the area of transcranial electrical stimulation (TES) often relies on computational models of current flow in the brain. Models are built on magnetic resonance images (MRI) of the human head to capture detailed individual anatomy. To simulate current flow, MRIs have to be segmented, virtual electrodes have to be placed on these anatomical models, the volume is tessellated into a mesh, and the finite element model is solved numerically to estimate the current flow. Various software tools are available for each step, as well as processing pipelines that connect these tools for automated or semi-automated processing. The goal of the present tool – ROAST – is to provide an end-to-end pipeline that can automatically process individual heads with realistic volumetric anatomy leveraging open-source software (SPM8, iso2mesh and getDP) and custom scripts to improve segmentation and execute electrode placement. When we compare the results on a standard head with other major commercial software tools for finite element modeling (ScanIP, Abaqus), ROAST only leads to a small difference of 9% in the estimated electric field in the brain. We obtain a larger difference of 47% when comparing results with SimNIBS, an automated pipeline that is based on surface segmentation of the head. We release ROAST as a fully automated pipeline available online as a open-source tool for TES modeling.

[1]  Stephen M Smith,et al.  Fast robust automated brain extraction , 2002, Human brain mapping.

[2]  Karl J. Friston,et al.  Statistical parametric mapping , 2013 .

[3]  M. Bikson,et al.  Transcranial current stimulation focality using disc and ring electrode configurations: FEM analysis , 2008, Journal of neural engineering.

[4]  Lucas C. Parra,et al.  Targeted electrical stimulation of active brain sources using reciprocity , 2016 .

[5]  D. Reato,et al.  Gyri-precise head model of transcranial direct current stimulation: Improved spatial focality using a ring electrode versus conventional rectangular pad , 2009, Brain Stimulation.

[6]  Dana H. Brooks,et al.  A pipeline for the simulation of transcranial direct current stimulation for realistic human head models using SCIRun/BioMesh3D , 2012, 2012 Annual International Conference of the IEEE Engineering in Medicine and Biology Society.

[7]  Anthony M. Norcia,et al.  In‐vivo measurements of human brain tissue conductivity using focal electrical current injection through intracerebral multicontact electrodes , 2017, Human brain mapping.

[8]  G. Crile,et al.  THE ELECTRICAL CONDUCTIVITY OF ANIMAL TISSUES UNDER NORMAL AND PATHOLOGICAL CONDITIONS , 1922 .

[9]  Carsten H. Wolters,et al.  A Discontinuous Galerkin Method to Solve the EEG Forward Problem Using the Subtraction Approach , 2015, SIAM J. Sci. Comput..

[10]  L. Parra,et al.  Inter-Individual Variation during Transcranial Direct Current Stimulation and Normalization of Dose Using MRI-Derived Computational Models , 2012, Front. Psychiatry.

[11]  Alexander Opitz,et al.  Determinants of the electric field during transcranial direct current stimulation , 2015, NeuroImage.

[12]  Alessia Paglialonga,et al.  Transcranial Direct Current Stimulation: Estimation of the Electric Field and of the Current Density in an Anatomical Human Head Model , 2011, IEEE Transactions on Biomedical Engineering.

[13]  Geddes La,et al.  Optimal stimulus duration for extracranial cortical stimulation. , 1987 .

[14]  Carsten H. Wolters,et al.  The Unfitted Discontinuous Galerkin Method for Solving the EEG Forward Problem , 2016, IEEE Transactions on Biomedical Engineering.

[15]  Markus Zahn,et al.  Transcranial direct current stimulation: A computer-based human model study , 2007, NeuroImage.

[16]  L. Geddes Optimal stimulus duration for extracranial cortical stimulation. , 1987, Neurosurgery.

[17]  Karl J. Friston,et al.  Unified segmentation , 2005, NeuroImage.

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

[19]  Yu Huang,et al.  Automated MRI segmentation for individualized modeling of current flow in the human head , 2013, Journal of neural engineering.

[20]  L. R. Dice Measures of the Amount of Ecologic Association Between Species , 1945 .

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

[22]  Won Hee Lee,et al.  Regional electric field induced by electroconvulsive therapy in a realistic finite element head model: Influence of white matter anisotropic conductivity , 2012, NeuroImage.

[23]  J. B. Ranck,et al.  Specific impedance of rabbit cerebral cortex. , 1963, Experimental neurology.

[24]  Adam J. Woods,et al.  Transcranial direct current stimulation in pediatric brain: A computational modeling study , 2012, 2012 Annual International Conference of the IEEE Engineering in Medicine and Biology Society.

[25]  Scott Makeig,et al.  Neuroelectromagnetic Forward Head Modeling Toolbox , 2010, Journal of Neuroscience Methods.

[26]  Lucas C. Parra,et al.  Optimal use of EEG recordings to target active brain areas with transcranial electrical stimulation , 2016, NeuroImage.

[27]  D. Louis Collins,et al.  Symmetric Atlasing and Model Based Segmentation: An Application to the Hippocampus in Older Adults , 2006, MICCAI.

[28]  Mark W. Woolrich,et al.  Advances in functional and structural MR image analysis and implementation as FSL , 2004, NeuroImage.

[29]  Yu Huang,et al.  ROAST: An Open-Source, Fully-Automated, Realistic Volumetric-Approach-Based Simulator For TES , 2018, 2018 40th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC).

[30]  Arthur W. Toga,et al.  Online resource for validation of brain segmentation methods , 2009, NeuroImage.

[31]  David A. Boas,et al.  Tetrahedral mesh generation from volumetric binary and grayscale images , 2009, 2009 IEEE International Symposium on Biomedical Imaging: From Nano to Macro.

[32]  Arthur W. Toga,et al.  Construction of a 3D probabilistic atlas of human cortical structures , 2008, NeuroImage.

[33]  Axel Thielscher,et al.  SimNIBS: A versatile toolbox for simulating fields generated by transcranial brain stimulation , 2015 .

[34]  W. Landau,et al.  Some relations between resistivity and electrical activity in the cerebral cortex of the cat. , 1955, Journal of cellular and comparative physiology.

[35]  Febo Cincotti,et al.  2019 41st Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC) , 2011, EMBC 2011.

[36]  Christophe Geuzaine,et al.  Gmsh: A 3‐D finite element mesh generator with built‐in pre‐ and post‐processing facilities , 2009 .

[37]  Marom Bikson,et al.  The Pursuit of DLPFC: Non-neuronavigated Methods to Target the Left Dorsolateral Pre-frontal Cortex With Symmetric Bicephalic Transcranial Direct Current Stimulation (tDCS) , 2015, Brain Stimulation.

[38]  Don Tucker The Scirun Inverse Eeg Pipeline | a Modeling and Simulation System for Cortical Mapping and Source Localization , 1998 .

[39]  H. C. Burger,et al.  Measurements of the specific Resistance of the human Body to direct Current , 2009 .

[40]  Abhishek Datta,et al.  State-of-art neuroanatomical target analysis of high-definition and conventional tDCS montages used for migraine and pain control , 2015, Front. Neuroanat..

[41]  Felipe Fregni,et al.  Classification of methods in transcranial Electrical Stimulation (tES) and evolving strategy from historical approaches to contemporary innovations , 2013, Journal of Neuroscience Methods.

[42]  Markus Zahn,et al.  Three-dimensional head model Simulation of transcranial magnetic stimulation , 2004, IEEE Transactions on Biomedical Engineering.

[43]  Kay Jann,et al.  In-vivo Imaging of Magnetic Fields Induced by Transcranial Direct Current Stimulation (tDCS) in Human Brain using MRI , 2016, Scientific Reports.

[44]  Lucas C Parra,et al.  The point spread function of the human head and its implications for transcranial current stimulation , 2012, Physics in medicine and biology.

[45]  Karl J. Friston,et al.  CHAPTER 2 – Statistical parametric mapping , 2007 .

[46]  F. Sánchez,et al.  Correlation between low-frequency electric conductivity and permittivity in the diaphysis of bovine femoral bone , 1992 .

[47]  Abhishek Datta,et al.  Transcranial direct current stimulation in patients with skull defects and skull plates: High-resolution computational FEM study of factors altering cortical current flow , 2010, NeuroImage.

[48]  Stefan Haufe,et al.  Large-scale EEG/MEG source localization with spatial flexibility , 2011, NeuroImage.

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

[50]  C. Gabriel Compilation of the Dielectric Properties of Body Tissues at RF and Microwave Frequencies. , 1996 .

[51]  Julie M. Baker,et al.  Individualized model predicts brain current flow during transcranial direct-current stimulation treatment in responsive stroke patient , 2011, Brain Stimulation.

[52]  D. Reato,et al.  Gyri – precise head model of transcranial DC stimulation : Improved spatial focality using a ring electrode versus conventional rectangular pad , 2010 .

[53]  Hartwig R. Siebner,et al.  Automatic skull segmentation from MR images for realistic volume conductor models of the head: Assessment of the state-of-the-art , 2018, NeuroImage.

[54]  Klaus Scheffler,et al.  Human in-vivo brain magnetic resonance current density imaging (MRCDI) , 2018, NeuroImage.

[55]  Marom Bikson,et al.  High‐Resolution Multi‐Scale Computational Model for Non‐Invasive Cervical Vagus Nerve Stimulation , 2018, Neuromodulation : journal of the International Neuromodulation Society.

[56]  Alexander Opitz,et al.  Electric field calculations in brain stimulation based on finite elements: An optimized processing pipeline for the generation and usage of accurate individual head models , 2013, Human brain mapping.

[57]  Alexander Opitz,et al.  Spatiotemporal structure of intracranial electric fields induced by transcranial electric stimulation in humans and nonhuman primates , 2016, Scientific Reports.

[58]  A. Dale,et al.  Cortical Surface-Based Analysis II: Inflation, Flattening, and a Surface-Based Coordinate System , 1999, NeuroImage.

[59]  Christophe Geuzaine,et al.  A general environment for the treatment of discrete problems and its application to the finite element method , 1998 .

[60]  Rosalind J. Sadleir,et al.  Transcranial direct current stimulation (tDCS) in a realistic head model , 2010, NeuroImage.

[61]  F. X. Bostick,et al.  Potential and current density distributions of cranial electrotherapy stimulation (CES) in a four-concentric-spheres model , 1996, IEEE Transactions on Biomedical Engineering.

[62]  Mark M. Stecker,et al.  Transcranial electric stimulation of motor pathways: a theoretical analysis , 2005, Comput. Biol. Medicine.

[63]  Yu Huang,et al.  Fully Automated Whole-Head Segmentation with Improved Smoothness and Continuity, with Theory Reviewed , 2015, PloS one.

[64]  Anders M. Dale,et al.  Cortical Surface-Based Analysis I. Segmentation and Surface Reconstruction , 1999, NeuroImage.

[65]  M. Fuchs,et al.  An improved boundary element method for realistic volume-conductor modeling , 1998, IEEE Transactions on Biomedical Engineering.

[66]  Axel Thielscher,et al.  Field modeling for transcranial magnetic stimulation: A useful tool to understand the physiological effects of TMS? , 2015, 2015 37th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC).

[67]  M. Hallett,et al.  Modeling the current distribution during transcranial direct current stimulation , 2006, Clinical Neurophysiology.

[68]  L. Parra,et al.  Measurements and models of electric fields in the in vivo human brain during transcranial electric stimulation , 2017, Brain Stimulation.

[69]  D. Griffiths Introduction to Electrodynamics , 2017 .