Very High Resolution NeuroElectric Brain Imaging Realized by Referee Consensus Processing

Electrical current flow within populations of neurons is a fundamental constituent of brain function. The resulting fluctuating magnetic fields may be sampled noninvasively with an array of magnetic field detectors positioned outside a patient’s head.  This is magnetoencephalography (MEG).  Each source may be characterized by 5-6 parameters, the xyz location and the xyz direction.  The magnetic field measurements are nonlinear in the location parameters; hence the source location is identifiable only via search of the brain volume.  When there is one or a very few sources, this may be practical; solutions for the general problem are weak. Referee consensus is a new method which enables identification of one source at a time regardless of the number and location of others. This “independence” enables solution of the general problem and insures suitability to grid computing.  The computation scales linearly with the number of nonlinear parameters. MEG recordings were obtained from volunteers while they performed a cognitive task The recordings were processed on the Open Science Grid (≈150 CPU hours/sec of data).  On average 500-1500 sources were active throughout.  Statistical analyses demonstrated < 2 mm resolving power [1] and very strong findings (p < 0.02 400 ) when testing for task specific information in the extracted virtual recordings from each individual. 3D maps of differential activation, neuroelectric tomography, provide a very high resolution functional imaging modality which compares favorably with functional MR imaging. Referee consensus is applicable widely to measurement systems including microwave telescope  imaging, seismic tomography, and financial market linkage identification. Applicability requires: (1) The measurements are linear in at least one parameter of each “source.” (2) Each source is detectable at multiple sensors. (3) A sequence of measurements in time is available.   [1] Linear dimensions are represented in this standard form.  Volume dimensions are represented throughout in terms of the length of a side, e.g. 8mm 3 instead of ½ cc, ½ cm 3 or 512 mm 3 .

[1]  Timothy Verstynen,et al.  In vivo mapping of microstructural somatotopies in the human corticospinal pathways. , 2011, Journal of neurophysiology.

[2]  Walter Schneider,et al.  Referee consensus: a platform technology for nonlinear optimization , 2013, XSEDE.

[3]  Anders M. Dale,et al.  Vector-based spatial–temporal minimum L1-norm solution for MEG , 2006, NeuroImage.

[4]  Gary H. Glover,et al.  Grand Challenges in Mapping the Human Brain: NSF Workshop Report , 2013, IEEE Transactions on Biomedical Engineering.

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

[6]  Wei Wang,et al.  Task-related MEG source localization via discriminant analysis , 2011, 2011 Annual International Conference of the IEEE Engineering in Medicine and Biology Society.

[7]  P. J. Charles,et al.  Real Time Signal Processing in the Clinical Setting , 1998, Annals of Biomedical Engineering.

[8]  D. Cheyne,et al.  Spatiotemporal mapping of cortical activity accompanying voluntary movements using an event‐related beamforming approach , 2006, Human brain mapping.

[9]  Raul Coimbra,et al.  Integrated imaging approach with MEG and DTI to detect mild traumatic brain injury in military and civilian patients. , 2009, Journal of neurotrauma.

[10]  Se Robinson,et al.  Functional neuroimaging by Synthetic Aperture Magnetometry (SAM) , 1999 .

[11]  J. Sarvas Basic mathematical and electromagnetic concepts of the biomagnetic inverse problem. , 1987, Physics in medicine and biology.

[12]  D. Lehmann,et al.  Low resolution electromagnetic tomography: a new method for localizing electrical activity in the brain. , 1994, International journal of psychophysiology : official journal of the International Organization of Psychophysiology.