Concentration maximization and local basis expansions (LBEX) for linear inverse problems

Linear inverse problems arise in biomedicine electro-encephalography and magnetoencephalography (EEG and MEG) and geophysics. The kernels relating sensors to the unknown sources are Green´s functions of some partial differential equation. This knowledge is obscured when treating the discretized kernels simply as matrices. Consequently, physical understanding of the fundamental resolution limits has been lacking. We relate the inverse problem to spatial Fourier analysis, and the resolution limits to uncertainty principles, providing conceptual links to underlying physics. Motivated by the spectral concentration problem and multitaper spectral analysis, our approach constructs local basis sets using maximally concentrated linear combinations of the measurement kernels