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

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

• Somersalo, E (Case Western Reserve University)
• Tuesday 11 February 2014, 11:00-11:45
• Seminar Room 1, Newton Institute.

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.

This talk is included in these lists:

• All CMS events
• Featured lists
• INI info aggregator
• Isaac Newton Institute Seminar Series
• School of Physical Sciences
• Seminar Room 1, Newton Institute
• bld31

Note that ex-directory lists are not shown.