Reconstructing cortical current density by exploring sparseness in the transform domain

In the present study, we have developed a novel electromagnetic source imaging approach to reconstruct extended cortical sources by means of cortical current density (CCD) modeling and a novel EEG imaging algorithm which explores sparseness in cortical source representations through the use of L1-norm in objective functions. The new sparse cortical current density (SCCD) imaging algorithm is unique since it reconstructs cortical sources by attaining sparseness in a transform domain (the variation map of cortical source distributions). While large variations are expected to occur along boundaries (sparseness) between active and inactive cortical regions, cortical sources can be reconstructed and their spatial extents can be estimated by locating these boundaries. We studied the SCCD algorithm using numerous simulations to investigate its capability in reconstructing cortical sources with different extents and in reconstructing multiple cortical sources with different extent contrasts. The SCCD algorithm was compared with two L2-norm solutions, i.e. weighted minimum norm estimate (wMNE) and cortical LORETA. Our simulation data from the comparison study show that the proposed sparse source imaging algorithm is able to accurately and efficiently recover extended cortical sources and is promising to provide high-accuracy estimation of cortical source extents.

[1]  C E Metz,et al.  Evaluation of receiver operating characteristic curve data in terms of information theory, with applications in radiography. , 1973, Radiology.

[2]  W. van Drongelen,et al.  EEG Source Imaging: Correlating Source Locations and Extents With Electrocorticography and Surgical Resections in Epilepsy Patients , 2007, Journal of clinical neurophysiology : official publication of the American Electroencephalographic Society.

[3]  T. Elbert,et al.  Phantom-limb pain as a perceptual correlate of cortical reorganization following arm amputation , 1995, Nature.

[4]  Barry D. Van Veen,et al.  Cortical patch basis model for spatially extended neural activity , 2006, IEEE Transactions on Biomedical Engineering.

[5]  Manbir Singh,et al.  An Evaluation of Methods for Neuromagnetic Image Reconstruction , 1987, IEEE Transactions on Biomedical Engineering.

[6]  R. Ilmoniemi,et al.  Magnetoencephalography-theory, instrumentation, and applications to noninvasive studies of the working human brain , 1993 .

[7]  K. Matsuura,et al.  Selective minimum-norm solution of the biomagnetic inverse problem , 1995, IEEE Transactions on Biomedical Engineering.

[8]  B He,et al.  Estimation of in vivo Human Brain-to-Skull Conductivity Ratio by means of Cortical Potential Imaging , 2005 .

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

[10]  W. Drongelen,et al.  Estimation of in vivo human brain-to-skull conductivity ratio from simultaneous extra- and intra-cranial electrical potential recordings , 2005, Clinical Neurophysiology.

[11]  Emmanuel J. Candès,et al.  Decoding by linear programming , 2005, IEEE Transactions on Information Theory.

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

[13]  V. Morozov.  On the solution of functional equations by the method of regularization , 1966 .

[14]  Jean Gotman,et al.  Evaluation of EEG localization methods using realistic simulations of interictal spikes , 2006, NeuroImage.

[15]  J S Ebersole,et al.  Noninvasive Localization of Epileptogenic Foci by EEG Source Modeling , 2000, Epilepsia.

[16]  Anders M. Dale,et al.  Improved Localization of Cortical Activity By Combining EEG and MEG with MRI Cortical Surface Reconstruction , 2002 .

[17]  Bin He,et al.  Modeling and Imaging of Bioelectrical Activity , 2005 .

[18]  Lei Ding,et al.  Sparse source imaging in electroencephalography with accurate field modeling , 2008, Human brain mapping.

[19]  Dmitry M. Malioutov,et al.  Optimal sparse representations in general overcomplete bases , 2004, 2004 IEEE International Conference on Acoustics, Speech, and Signal Processing.

[20]  Richard M. Leahy,et al.  MEG-based imaging of focal neuronal current sources , 1995, 1995 IEEE Nuclear Science Symposium and Medical Imaging Conference Record.

[21]  Eija Kalso,et al.  Phantom-limb pain , 1998, The Lancet.

[22]  R D Sidman,et al.  A method for localization of sources of human cerebral potentials evoked by sensory stimuli. , 1978, Sensory processes.

[23]  Stephen P. Boyd,et al.  Convex Optimization , 2004, Algorithms and Theory of Computation Handbook.

[24]  Toshimitsu Musha,et al.  Electric Dipole Tracing in the Brain by Means of the Boundary Element Method and Its Accuracy , 1987, IEEE Transactions on Biomedical Engineering.

[25]  M. Hämäläinen,et al.  Realistic conductivity geometry model of the human head for interpretation of neuromagnetic data , 1989, IEEE Transactions on Biomedical Engineering.

[26]  D. Donoho.  For most large underdetermined systems of linear equations the minimal 𝓁1‐norm solution is also the sparsest solution , 2006 .

[27]  Jos F. Sturm,et al.  A Matlab toolbox for optimization over symmetric cones , 1999 .

[28]  John W Belliveau,et al.  Monte Carlo simulation studies of EEG and MEG localization accuracy , 2002, Human brain mapping.

[29]  E. Somersalo,et al.  Visualization of Magnetoencephalographic Data Using Minimum Current Estimates , 1999, NeuroImage.

[30]  M. Fuchs,et al.  Smooth reconstruction of cortical sources from EEG or MEG recordings , 1996, NeuroImage.

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

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