We have decided to discontinue the publication of preprints on our preprint server end of 2024. The publication culture within mathematics has changed so much due to the rise of repositories such as ArXiV (www.arxiv.org) that we are encouraging all institute members to make their preprints available there. An institute's repository in its previous form is, therefore, unnecessary. The preprints published to date will remain available here, but we will not add any new preprints here.
MiS Preprint
59/1998
Stable focal inverse source localization using combinatorial optimization techniques combined with regularization methods
Carsten H. Wolters, Rainer Beckmann, Adrian Rienäcker and Helmut Buchner
Abstract
The inverse problem arising from EEG and MEG is largely underdetermined. One strategy to alleviate this problem is the restriction to a limited number of pointlike sources, the focal source model. Although the singular value decomposition of the spatio-temporal data gives an estimate of the minimal number of dipoles contributing to the measurement, the exact number is unknown in advance and noise complicates the reconstruction. Classical nonlinear dipole fit algorithms do not give an estimate for the correct number because they are not stable with regard to an overestimation of this parameter. Too many sources may only describe noise, but can still attain a large magnitude during the inverse procedure and may be indiscernible from the true sources. This paper describes a nonlinear dipole fit reconstruction algorithm with a new regularization approach for the embedded linear problem, automatically controlled by the noise in the data and the condition of the occuring least square problems. The algorithm is stable with regard to source components which "nearly" lie in the kernel of the projection or lead field operator and it thus gives an estimate of the unknown number parameter. Simulation studies in a simulated sulcus structure are carried out for an instantaneous dipole model and spatial resolution in the sulcus and stability of the new method are compared with a classical reconstruction algorithm without regularization.