University of Cambridge > Talks.cam > Isaac Newton Institute Seminar Series > Bayesian preconditioning for truncated Krylov subspace regularization with an application to Magnetoencephalography (MEG)

Bayesian preconditioning for truncated Krylov subspace regularization with an application to Magnetoencephalography (MEG)

Add to your list(s) Download to your calendar using vCal

If you have a question about this talk, please contact Mustapha Amrani.

Inverse Problems

Co-authors: Daniela Calvetti (Case Western Reserve University), Laura Homa (Case Western Reserve University)

We consider the computational problem arising in magnetoencephalography (MEG), where the goal is to estimate the electric activity within the brain non-invasively from extra-cranial measurements of the magnetic field components. The problem is severely ill-posed due to the intrinsic non-uniqueness of the solution, and suffer further from the challenges of starting from a weak data signal, its high dimensionality and complexity of the noise, part of which is due to the brain itself. We propose a new algorithm that is based on truncated conjugate gradient algorithm for least squares (CGLS) with statistically inspired left and right preconditioners. We demonstrate that by carefully accounting for the spatiotemporal statistical structure of the brain noise, and by adopting a suitable prior within the Bayesian framework, we can design a robust and efficient method for the numerical solution of the MEG inverse problem which can improve the spatial and temporal resolution of events of short duration.

This talk is part of the Isaac Newton Institute Seminar Series series.

Tell a friend about this talk:

This talk is included in these lists:

Note that ex-directory lists are not shown.

 

© 2006-2024 Talks.cam, University of Cambridge. Contact Us | Help and Documentation | Privacy and Publicity