The direct and the inverse magnetic encephalography problem
-
Ivan P. Pakhnenko
Abstract
Magnetic encephalography provides a unique opportunity for non-invasive study of neural processes occurring in the brain, but at the same time produces a large amount of data. Processing this data in order to reconstruct the signal sources with a given accuracy is an extremely difficult task. Our work is devoted to solving the direct and inverse problems of magnetoencephalography. The inverse problem o magnetoencephalography is ill-posed and difficult for both analytical and numerical solutions. Additional complications arise from the volume (passive) currents and the associated magnetic fields, which strongly depend on the brain geometry. In this paper, we find approximate analytical solutions for the forward and the inverse problems in the spherical and spheroid geometry. To localize the source of activity, it is first necessary to consider a direct problem, which in our case has an exact analytical solution. For this purpose, magnetic induction was calculated on a given surface of a complex system of dipoles, and a program was written to visualize magnetic induction on the surface of the head. In certain models of the human head (spherical and spheroid), the Biot-Savart law and IC analysis provide the basis for a detailed study in order to develop an algorithm for localization of the primary motor cortex. For the general case of the spheroid model, the inverse problem can be approximately solved by neglecting the volumetric magnetic field near the points of the maximum magnetic field and taking into account only the primary magnetic field. This paper presents a step-by-step algorithm for obtaining a solution to the inverse MEG problem under the assumptions of discreteness of signal sources originating from different functional areas of the brain and the surface location of signal sources. The new method is based on our obtained exact analytical solution of the inverse problem for the one-dipole model. Within the adopted constraints, the method is also applicable for a larger number of dipoles. The case of two dipoles has been investigated so far. Using the ICA to separate the magnetic field into its components and applying the analytical formula for the case of one dipole, we can obtain an approximate solution of the inverse MEG problem with high accuracy. The localization methods presented in this paper are undoubtedly important in real clinical practice. Thus, during neurosurgical interventions, various areas of the brain can be damaged, including irreparable ones. Since the location of functional zones in the human brain is individual, the doctor needs to be able to localize these areas in the preoperative period with high accuracy. The methods developed in this paper serve to solve such an important task.
Abstract
Magnetic encephalography provides a unique opportunity for non-invasive study of neural processes occurring in the brain, but at the same time produces a large amount of data. Processing this data in order to reconstruct the signal sources with a given accuracy is an extremely difficult task. Our work is devoted to solving the direct and inverse problems of magnetoencephalography. The inverse problem o magnetoencephalography is ill-posed and difficult for both analytical and numerical solutions. Additional complications arise from the volume (passive) currents and the associated magnetic fields, which strongly depend on the brain geometry. In this paper, we find approximate analytical solutions for the forward and the inverse problems in the spherical and spheroid geometry. To localize the source of activity, it is first necessary to consider a direct problem, which in our case has an exact analytical solution. For this purpose, magnetic induction was calculated on a given surface of a complex system of dipoles, and a program was written to visualize magnetic induction on the surface of the head. In certain models of the human head (spherical and spheroid), the Biot-Savart law and IC analysis provide the basis for a detailed study in order to develop an algorithm for localization of the primary motor cortex. For the general case of the spheroid model, the inverse problem can be approximately solved by neglecting the volumetric magnetic field near the points of the maximum magnetic field and taking into account only the primary magnetic field. This paper presents a step-by-step algorithm for obtaining a solution to the inverse MEG problem under the assumptions of discreteness of signal sources originating from different functional areas of the brain and the surface location of signal sources. The new method is based on our obtained exact analytical solution of the inverse problem for the one-dipole model. Within the adopted constraints, the method is also applicable for a larger number of dipoles. The case of two dipoles has been investigated so far. Using the ICA to separate the magnetic field into its components and applying the analytical formula for the case of one dipole, we can obtain an approximate solution of the inverse MEG problem with high accuracy. The localization methods presented in this paper are undoubtedly important in real clinical practice. Thus, during neurosurgical interventions, various areas of the brain can be damaged, including irreparable ones. Since the location of functional zones in the human brain is individual, the doctor needs to be able to localize these areas in the preoperative period with high accuracy. The methods developed in this paper serve to solve such an important task.
Chapters in this book
- Frontmatter I
- Prologue I VII
- Prologue II XI
- Prologue III XIII
- Preface XVII
- Overview XIX
- Contents XXXIII
-
Part I: Theories
-
Part I-A: Overarching theory
- Introduction 1
- Universal axioms in classical Chinese philosophy 5
- Category theory for structural characterization 15
- Axiomatic bipolar dynamics and their control 45
-
Part I-B: Systems theories
- Introduction 75
- Stochastic formalization of agent-oriented systems 79
- Simplification of high-dimensional multitempo dynamic models 109
- Ideas of symmetry as a biophysical basis of system biomedicine 123
- Disorder of multiscale control 149
-
Part II: Person’s life-sphere
-
Part II-A: Person’s biosphere
- Introduction 185
- Mutations as activators of biological evolutionary processes at population levels 189
- Immunometabolism of T-cells in COVID-19 209
-
Part II-A.2: Body’s vital functions
- Introduction 245
- Structural modeling of vascular networks 249
- Mathematical modeling of AI application for the diagnosis of blood flow disorders 283
- Modeling of glucose and insulin regulation within the framework of a self-consistent model of the cardiovascular system 303
- Hemodynamics in residual myocardial ischemia 319
- The quasi-one-dimensional model of the lymph flow in the human lymphatic system 335
- An integrate-and-fire mechanism for modeling rhythmicity in the neuroendocrine system 365
- Kinetic network modeling of the neuroendocrine hypothalamic-pituitary-adrenal axis dynamics with particular attention on the role of alcohol as a digestif 377
- Inflammation and immune response in atherosclerosis 393
-
Part II-A.3: Body’s motor functions
- Introduction 423
- A magnetic resonance spectroscopy approach to quantitatively measure GABA and phosphorus level changes in the primary motor cortex elicited by transcranial direct current stimulation 427
-
Part II-A.4: Body’s operational functions
- Introduction 441
- The fermionic mind hypothesis–a category-theoretic verification of consciousness 445
- Cross-task cognitive workload measurement based on the sample selection of the EEG data 459
-
Part II-B: Person’s eco-sphere exposures
- Introduction 475
- The spread of SARS-CoV-2 in Russia and the evolution of the properties of the pathogen 479
- Agent-based modeling of epidemic spread via kinetic Monte Carlo method 491
- Control of SARS-nCoV outbreaks in China 2020 513
-
Part II-B.2: Civilization
- Introduction 531
- Pesticide exposure: Toward holistic environmental modeling 535
-
Part II-C: Person’s sociosphere exposures
- Introduction 559
- Evolution of the health system in Shanghai, China, 2016–2020 563
-
Part III: Technologies
- Introduction 577
- Design-process automation using functional process blocks 581
- Slow/fast dynamic models with applications to engineering problems 601
-
Part III-B: Information sciences
- Introduction 613
- Numerical modeling of medical ultrasound using the grid-characteristic method 617
- The direct and the inverse magnetic encephalography problem 635
-
Part III-C: Data-analytic sciences
- Introduction 653
- Assessing the bioequivalence of two different drugs with the same active ingredient 655
- Estimation of adjusted relative risks in log-binomial regression using the Bekhit–Schöpe–Wagenpfeil algorithm 665
-
Part IV: Clinical medicine
- Introduction 679
- Finding optimal two-stage combined treatment protocols for a blood cancer model 681
- Unraveling the mysteries: Mathematical perspectives on traditional Chinese medicine meridians 697
- Epilogue 721
- Index 723
Chapters in this book
- Frontmatter I
- Prologue I VII
- Prologue II XI
- Prologue III XIII
- Preface XVII
- Overview XIX
- Contents XXXIII
-
Part I: Theories
-
Part I-A: Overarching theory
- Introduction 1
- Universal axioms in classical Chinese philosophy 5
- Category theory for structural characterization 15
- Axiomatic bipolar dynamics and their control 45
-
Part I-B: Systems theories
- Introduction 75
- Stochastic formalization of agent-oriented systems 79
- Simplification of high-dimensional multitempo dynamic models 109
- Ideas of symmetry as a biophysical basis of system biomedicine 123
- Disorder of multiscale control 149
-
Part II: Person’s life-sphere
-
Part II-A: Person’s biosphere
- Introduction 185
- Mutations as activators of biological evolutionary processes at population levels 189
- Immunometabolism of T-cells in COVID-19 209
-
Part II-A.2: Body’s vital functions
- Introduction 245
- Structural modeling of vascular networks 249
- Mathematical modeling of AI application for the diagnosis of blood flow disorders 283
- Modeling of glucose and insulin regulation within the framework of a self-consistent model of the cardiovascular system 303
- Hemodynamics in residual myocardial ischemia 319
- The quasi-one-dimensional model of the lymph flow in the human lymphatic system 335
- An integrate-and-fire mechanism for modeling rhythmicity in the neuroendocrine system 365
- Kinetic network modeling of the neuroendocrine hypothalamic-pituitary-adrenal axis dynamics with particular attention on the role of alcohol as a digestif 377
- Inflammation and immune response in atherosclerosis 393
-
Part II-A.3: Body’s motor functions
- Introduction 423
- A magnetic resonance spectroscopy approach to quantitatively measure GABA and phosphorus level changes in the primary motor cortex elicited by transcranial direct current stimulation 427
-
Part II-A.4: Body’s operational functions
- Introduction 441
- The fermionic mind hypothesis–a category-theoretic verification of consciousness 445
- Cross-task cognitive workload measurement based on the sample selection of the EEG data 459
-
Part II-B: Person’s eco-sphere exposures
- Introduction 475
- The spread of SARS-CoV-2 in Russia and the evolution of the properties of the pathogen 479
- Agent-based modeling of epidemic spread via kinetic Monte Carlo method 491
- Control of SARS-nCoV outbreaks in China 2020 513
-
Part II-B.2: Civilization
- Introduction 531
- Pesticide exposure: Toward holistic environmental modeling 535
-
Part II-C: Person’s sociosphere exposures
- Introduction 559
- Evolution of the health system in Shanghai, China, 2016–2020 563
-
Part III: Technologies
- Introduction 577
- Design-process automation using functional process blocks 581
- Slow/fast dynamic models with applications to engineering problems 601
-
Part III-B: Information sciences
- Introduction 613
- Numerical modeling of medical ultrasound using the grid-characteristic method 617
- The direct and the inverse magnetic encephalography problem 635
-
Part III-C: Data-analytic sciences
- Introduction 653
- Assessing the bioequivalence of two different drugs with the same active ingredient 655
- Estimation of adjusted relative risks in log-binomial regression using the Bekhit–Schöpe–Wagenpfeil algorithm 665
-
Part IV: Clinical medicine
- Introduction 679
- Finding optimal two-stage combined treatment protocols for a blood cancer model 681
- Unraveling the mysteries: Mathematical perspectives on traditional Chinese medicine meridians 697
- Epilogue 721
- Index 723