A Subspace Pursuit-based Iterative Greedy Hierarchical solution to the neuromagnetic inverse problem

Magnetoencephalography (MEG) is an important non-invasive method for studying activity within the human brain. Source localization methods can be used to estimate spatiotemporal activity from MEG measurements with high temporal resolution, but the spatial resolution of these estimates is poor due to...

Full description

Bibliographic Details
Main Authors: Babadi, Behtash, Obregon-Henao, Gabriel, Lamus, Camilo, Brown, Emery N., Hamalainen, Matti S., Purdon, Patrick Lee
Other Authors: Institute for Medical Engineering and Science
Format: Article
Language:en_US
Published: Elsevier 2016
Online Access:http://hdl.handle.net/1721.1/102345
https://orcid.org/0000-0001-5651-5060
https://orcid.org/0000-0003-2668-7819
https://orcid.org/0000-0001-6841-112X
https://orcid.org/0000-0002-6777-7979
_version_ 1811082644690042880
author Babadi, Behtash
Obregon-Henao, Gabriel
Lamus, Camilo
Brown, Emery N.
Hamalainen, Matti S.
Purdon, Patrick Lee
author2 Institute for Medical Engineering and Science
author_facet Institute for Medical Engineering and Science
Babadi, Behtash
Obregon-Henao, Gabriel
Lamus, Camilo
Brown, Emery N.
Hamalainen, Matti S.
Purdon, Patrick Lee
author_sort Babadi, Behtash
collection MIT
description Magnetoencephalography (MEG) is an important non-invasive method for studying activity within the human brain. Source localization methods can be used to estimate spatiotemporal activity from MEG measurements with high temporal resolution, but the spatial resolution of these estimates is poor due to the ill-posed nature of the MEG inverse problem. Recent developments in source localization methodology have emphasized temporal as well as spatial constraints to improve source localization accuracy, but these methods can be computationally intense. Solutions emphasizing spatial sparsity hold tremendous promise, since the underlying neurophysiological processes generating MEG signals are often sparse in nature, whether in the form of focal sources, or distributed sources representing large-scale functional networks. Recent developments in the theory of compressed sensing (CS) provide a rigorous framework to estimate signals with sparse structure. In particular, a class of CS algorithms referred to as greedy pursuit algorithms can provide both high recovery accuracy and low computational complexity. Greedy pursuit algorithms are difficult to apply directly to the MEG inverse problem because of the high-dimensional structure of the MEG source space and the high spatial correlation in MEG measurements. In this paper, we develop a novel greedy pursuit algorithm for sparse MEG source localization that overcomes these fundamental problems. This algorithm, which we refer to as the Subspace Pursuit-based Iterative Greedy Hierarchical (SPIGH) inverse solution, exhibits very low computational complexity while achieving very high localization accuracy. We evaluate the performance of the proposed algorithm using comprehensive simulations, as well as the analysis of human MEG data during spontaneous brain activity and somatosensory stimuli. These studies reveal substantial performance gains provided by the SPIGH algorithm in terms of computational complexity, localization accuracy, and robustness.
first_indexed 2024-09-23T12:06:41Z
format Article
id mit-1721.1/102345
institution Massachusetts Institute of Technology
language en_US
last_indexed 2024-09-23T12:06:41Z
publishDate 2016
publisher Elsevier
record_format dspace
spelling mit-1721.1/1023452022-09-28T00:14:27Z A Subspace Pursuit-based Iterative Greedy Hierarchical solution to the neuromagnetic inverse problem Babadi, Behtash Obregon-Henao, Gabriel Lamus, Camilo Brown, Emery N. Hamalainen, Matti S. Purdon, Patrick Lee Institute for Medical Engineering and Science Harvard University--MIT Division of Health Sciences and Technology Massachusetts Institute of Technology. Department of Brain and Cognitive Sciences McGovern Institute for Brain Research at MIT Babadi, Behtash Lamus, Camilo Hamalainen, Matti S. Brown, Emery N. Purdon, Patrick Lee Magnetoencephalography (MEG) is an important non-invasive method for studying activity within the human brain. Source localization methods can be used to estimate spatiotemporal activity from MEG measurements with high temporal resolution, but the spatial resolution of these estimates is poor due to the ill-posed nature of the MEG inverse problem. Recent developments in source localization methodology have emphasized temporal as well as spatial constraints to improve source localization accuracy, but these methods can be computationally intense. Solutions emphasizing spatial sparsity hold tremendous promise, since the underlying neurophysiological processes generating MEG signals are often sparse in nature, whether in the form of focal sources, or distributed sources representing large-scale functional networks. Recent developments in the theory of compressed sensing (CS) provide a rigorous framework to estimate signals with sparse structure. In particular, a class of CS algorithms referred to as greedy pursuit algorithms can provide both high recovery accuracy and low computational complexity. Greedy pursuit algorithms are difficult to apply directly to the MEG inverse problem because of the high-dimensional structure of the MEG source space and the high spatial correlation in MEG measurements. In this paper, we develop a novel greedy pursuit algorithm for sparse MEG source localization that overcomes these fundamental problems. This algorithm, which we refer to as the Subspace Pursuit-based Iterative Greedy Hierarchical (SPIGH) inverse solution, exhibits very low computational complexity while achieving very high localization accuracy. We evaluate the performance of the proposed algorithm using comprehensive simulations, as well as the analysis of human MEG data during spontaneous brain activity and somatosensory stimuli. These studies reveal substantial performance gains provided by the SPIGH algorithm in terms of computational complexity, localization accuracy, and robustness. National Institutes of Health (U.S.) (New Innovator Award DP2-OD006454) National Institutes of Health (U.S.) (R01-EB006385-01) National Institutes of Health (U.S.) (P41-RR014075-11) 2016-05-02T13:34:09Z 2016-05-02T13:34:09Z 2013-09 Article http://purl.org/eprint/type/JournalArticle 10538119 http://hdl.handle.net/1721.1/102345 Babadi, Behtash, Gabriel Obregon-Henao, Camilo Lamus, Matti S. Hamalainen, Emery N. Brown, and Patrick L. Purdon. “A Subspace Pursuit-Based Iterative Greedy Hierarchical Solution to the Neuromagnetic Inverse Problem.” NeuroImage 87 (February 2014): 427–443. https://orcid.org/0000-0001-5651-5060 https://orcid.org/0000-0003-2668-7819 https://orcid.org/0000-0001-6841-112X https://orcid.org/0000-0002-6777-7979 en_US http://dx.doi.org/10.1016/j.neuroimage.2013.09.008 NeuroImage Creative Commons Attribution-NonCommercial-NoDerivs License http://creativecommons.org/licenses/by-nc-nd/4.0/ application/pdf Elsevier PMC
spellingShingle Babadi, Behtash
Obregon-Henao, Gabriel
Lamus, Camilo
Brown, Emery N.
Hamalainen, Matti S.
Purdon, Patrick Lee
A Subspace Pursuit-based Iterative Greedy Hierarchical solution to the neuromagnetic inverse problem
title A Subspace Pursuit-based Iterative Greedy Hierarchical solution to the neuromagnetic inverse problem
title_full A Subspace Pursuit-based Iterative Greedy Hierarchical solution to the neuromagnetic inverse problem
title_fullStr A Subspace Pursuit-based Iterative Greedy Hierarchical solution to the neuromagnetic inverse problem
title_full_unstemmed A Subspace Pursuit-based Iterative Greedy Hierarchical solution to the neuromagnetic inverse problem
title_short A Subspace Pursuit-based Iterative Greedy Hierarchical solution to the neuromagnetic inverse problem
title_sort subspace pursuit based iterative greedy hierarchical solution to the neuromagnetic inverse problem
url http://hdl.handle.net/1721.1/102345
https://orcid.org/0000-0001-5651-5060
https://orcid.org/0000-0003-2668-7819
https://orcid.org/0000-0001-6841-112X
https://orcid.org/0000-0002-6777-7979
work_keys_str_mv AT babadibehtash asubspacepursuitbasediterativegreedyhierarchicalsolutiontotheneuromagneticinverseproblem
AT obregonhenaogabriel asubspacepursuitbasediterativegreedyhierarchicalsolutiontotheneuromagneticinverseproblem
AT lamuscamilo asubspacepursuitbasediterativegreedyhierarchicalsolutiontotheneuromagneticinverseproblem
AT brownemeryn asubspacepursuitbasediterativegreedyhierarchicalsolutiontotheneuromagneticinverseproblem
AT hamalainenmattis asubspacepursuitbasediterativegreedyhierarchicalsolutiontotheneuromagneticinverseproblem
AT purdonpatricklee asubspacepursuitbasediterativegreedyhierarchicalsolutiontotheneuromagneticinverseproblem
AT babadibehtash subspacepursuitbasediterativegreedyhierarchicalsolutiontotheneuromagneticinverseproblem
AT obregonhenaogabriel subspacepursuitbasediterativegreedyhierarchicalsolutiontotheneuromagneticinverseproblem
AT lamuscamilo subspacepursuitbasediterativegreedyhierarchicalsolutiontotheneuromagneticinverseproblem
AT brownemeryn subspacepursuitbasediterativegreedyhierarchicalsolutiontotheneuromagneticinverseproblem
AT hamalainenmattis subspacepursuitbasediterativegreedyhierarchicalsolutiontotheneuromagneticinverseproblem
AT purdonpatricklee subspacepursuitbasediterativegreedyhierarchicalsolutiontotheneuromagneticinverseproblem