# Inverse conductivity recovery problem in a spherical geometry from EEG data: uniqueness, reconstruction and stability results

Maureen Clerc, Juliette Leblond, Jean-Paul Marmorat, Christos Papageorgakis

## ► To cite this version:

Maureen Clerc, Juliette Leblond, Jean-Paul Marmorat, Christos Papageorgakis. Inverse conductivity recovery problem in a spherical geometry from EEG data: uniqueness, reconstruction and stability results. TAMTAM 2017 - 8ème colloque, Tendances dans les Applications Mathématiques en Tunisie Algérie Maroc, May 2017, Hammamet, Tunisia. pp.1. hal-01526931

HAL Id: hal-01526931

https://inria.hal.science/hal-01526931

Submitted on 23 May 2017

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
INVERSE CONDUCTIVITY RECOVERY PROBLEM IN A SPHERICAL GEOMETRY
FROM EEG DATA: UNIQUENESS, RECONSTRUCTION AND STABILITY RESULTS

MAUREEN CLERC<sup>1,3</sup>, JULIETTE LEBLOND<sup>1,3</sup>, JEAN-PAUL MARMORAT<sup>2</sup>, CHRISTOS PAPAGEORGAKIS<sup>1,3</sup>.

<sup>1</sup> INRIA Sophia Antipolis Méditerranée, Équipes Apics et Athena, France.

<sup>2</sup> Ecole des Mines ParisTech, Centre de Mathématiques Appliquées, Sophia Antipolis, France.

<sup>3</sup> Université Côte d'Azur, France.

Abstract

Electroencephalography (EEG) is a non invasive imaging technique that measures the effect of the electric activity of active brain regions, called sources, through values of the electric potential furnished by a set of electrodes placed at the surface of the scalp. A fundamental problem there is the inverse problem of source localization which aims at locating the sources of the electric activity using the acquired EEG measurements [2]. The quality of the source estimation depends on the accuracy of the conductivity model used to solve the problem. Among the head tissues, the skull conductivity is the one that influences most the accuracy of EEG source localization [3]. Often, conductivity estimation is performed prior to the source estimation to determine the unknown conductivity using either supplementary EEG measurements or even measurements acquired by other imaging techniques. Such a technique is electrical impedance tomography (EIT) where current is injected through a pair of EEG electrodes while the unknown conductivities can be estimated by the resulting measurements at the rest electrode locations.

We examine the inverse skull conductivity estimation problem, which aims at recovering the electrical conductivity properties of the skull from measurements given at the surface of the head by EEG measurements. Our goal is to show uniqueness and a constructive scheme for the inverse skull conductivity estimation problem using partial boundary EEG data from a single experiment, in the preliminary case of an homogeneous skull conductivity. This is a version of the many inverse conductivity issues still under study nowadays [1].

The head is assumed to be an isotropic piecewise homogeneous medium and we examine a layered spherical head model made of three concentric nested spheres, each of them modelling scalp, skull and brain tissues (from the outermost to the innermost layer). Each of the three layers is supposed to have a constant conductivity. We also assume that the conductivities of the brain and the scalp are known, while the conductivity to be recovered is the one of the intermediate spherical layer (skull).

We solve the above conductivity estimation problem from the available EEG partial boundary data, expanded on the spherical harmonics basis, and transmitted over the spherical interfaces by transfer functions, while we consider that the source term is already estimated (through a number of coefficients of its spherical harmonics expansion). Linear algebra computations then allow us to find polynomials that possess a root which should coincide with the unknown skull conductivity, thus solving the estimation problem. We derive uniqueness properties and a reconstruction algorithm for the skull conductivity. A numerical study shows that the algorithm is able to accurately estimate the skull conductivity, with good robustness properties with respect to various levels of noise. The properties of the inverse conductivity estimation problem are also examined with various source configurations (partially known sources) and EIT measurements.

This work was supported by the Région Provence-Alpes-Côte d'Azur, France, and BESA GmbH, Germany.

REFERENCES

[1] G. ALESSANDRINI, *Generic uniqueness and size estimates in the inverse conductivity problem with one measurement*, Le Matematiche, 54 (1999), pp. 5-14.

[2] M. CLERC, J. LEBLOND, J.-P. MARMORAT AND T. PAPADOPOULO, *Source localization in EEG using rational approximation on plane sections*, Inverse Problems, 28, 055018 (2012).

[3] S. VALLAGHÉ AND M. CLERC, *A Global Sensitivity Analysis of Three-and Four-Layer EEG Conductivity Models*, IEEE Trans. Biomedical Engineering, 56(4) (2009), pp. 988-995.