Inverse recovery in EEG

The inverse recovery in EEG is a Calderón-type inverse problem with the goal of recovering source terms and/or conductivity in layers of the human head from electroencephalographic measurements. Fundamentally, this inverse recovery seeks to solve the elliptic partial differential equation given by
or
where is the electric potential, is the conductivity, and represents primary current sources in the brain. Depending on the application, the inverse problem consists of either recovering from or recovering from .
Because the human head is highly inhomogeneous and composed of multiple layers with different conductivities, the inverse EEG problem is severely ill-posed and requires either analytical techniques or numerical approximations to obtain stable solutions.
Due to the form of the governing equation being a sort of "generalization" of the Laplace-Beltrami operator, this problem has deep connections to generalized analytic function theory, heat conduction, and broader electromagnetics. In fact, the problem can be seen as solving the Poisson equation for an inhomogeneous media, which is indeed how it is derived in the EEG problem.

Problem derivation

An overview on the physical foundations of the problem is given by Darbas and Lohrengel.
Consider the current density produced by neural activity,
where denotes primary current and the return current. Using the quasi-static approximation of Maxwell's equations,
The first of the above equation implies that the electric field is path-independent and thus we may write with the electric potential function. Then,
which gives
The primary current is modeled by pointwise sources located at some coordinate with dipolar moments . Using the Dirac delta distribution,
Thus, the source term is

Exact solutions

Constant conductivity

Many analytical studies assume that conductivity is piecewise constant in layers representing the brain, skull, and scalp. Let denote the head domain decomposed into three concentric spherical domains, where is a ball and are spherical shells. When is assumed piecewise constant, the governing equation reduces to a set of Laplace and Poisson equations:
Solutions can be obtained using Green's functions, spherical harmonics, and separation of variables.
A benchmark tool in the problem of piecewise constant conductivity is FindSource3D, a source localization tool which has shown perfect reconstruction when negligible error is assumed in source location. Such a tool has also been used for conductivity reconstruction, again with perfect accuracy under suitable error assumptions.

Non-constant conductivity

For non-constant or anisotropic conductivity, the inverse problem becomes substantially more difficult. Classical uniqueness and stability results for Calderón-type inverse problems apply under varying regularity assumptions on the boundary of ; for example, Alessandrini established uniqueness with Lipschitz boundaries, while Kohn and Vogelius proved stability under smooth boundary conditions.
Exact solutions for spatially varying are generally unavailable. Modern work therefore typically focuses on perturbative methods, integral equation formulations, or regularization techniques for ensuring stable recovery.

Numerical solutions

Because analytical solutions are available only in highly idealized geometries, numerical methods play a central role in practical EEG inverse recovery. Common approaches include typical Finite Element or Boundary Element Methods.
Numerical inverse recovery typically requires regularization methods such as Tikhonov regularization, sparsity constraints, or Bayesian inverse modeling to counteract ill-posedness. Practical EEG systems also incorporate noise modeling and electrode placement uncertainty.