Fast Approximation of EEG Forward Problem and Application to Tissue Conductivity Estimation

Bioelectric source analysis in the human brain from scalp electroencephalography (EEG) signals is sensitive to the conductivity of the different head tissues. Conductivity values are subject dependent, so non-invasive methods for conductivity estimation are necessary to fine tune the EEG models. To do so, the EEG forward problem solution (so-called lead field matrix) must be computed for a large number of conductivity configurations. Computing one lead field requires a matrix inversion which is computationally intensive for realistic head models. Thus, the required time for computing a large number of lead fields can become impractical. In this work, we propose to approximate the lead field matrix for a set of conductivity configurations, using the exact solution only for a small set of basis points in the conductivity space. Our approach accelerates the computing time, while controlling the approximation error. Our method is tested for brain and skull conductivity estimation , with simulated and measured EEG data, corresponding to evoked somato-sensory potentials. This test demonstrates that the used approximation does not introduce any bias and runs significantly faster than if exact lead field were to be computed.Comment: Copyright (c) 2019 IEEE. Personal use of this material is permitted. However, permission to use this material for any other purposes must be obtained from the IEEE by sending a request to [email protected]

Inward continuation of the scalp potential distribution by means of the (vector) BEM

The vector Boundary Element Method (vBEM) is used for the calculation of a matrix that links the tangential components of the current density on the cortical and scalp surface. This so-called transfer matrix is compared to the transfer matrix that links the potential distribution on both surfaces. Forward and inverse calculations are performed to evaluate both types of transfer matrice

3D shape based reconstruction of experimental data in Diffuse Optical Tomography

Diffuse optical tomography (DOT) aims at recovering three-dimensional images of absorption and scattering parameters inside diffusive body based on small number of transmission measurements at the boundary of the body. This image reconstruction problem is known to be an ill-posed inverse problem, which requires use of prior information for successful reconstruction. We present a shape based method for DOT, where we assume a priori that the unknown body consist of disjoint subdomains with different optical properties. We utilize spherical harmonics expansion to parameterize the reconstruction problem with respect to the subdomain boundaries, and introduce a finite element (FEM) based algorithm that uses a novel 3D mesh subdivision technique to describe the mapping from spherical harmonics coefficients to the 3D absorption and scattering distributions inside a unstructured volumetric FEM mesh. We evaluate the shape based method by reconstructing experimental DOT data, from a cylindrical phantom with one inclusion with high absorption and one with high scattering. The reconstruction was monitored, and we found a 87% reduction in the Hausdorff measure between targets and reconstructed inclusions, 96% success in recovering the location of the centers of the inclusions and 87% success in average in the recovery for the volumes

EIT Reconstruction Algorithms: Pitfalls, Challenges and Recent Developments

We review developments, issues and challenges in Electrical Impedance Tomography (EIT), for the 4th Workshop on Biomedical Applications of EIT, Manchester 2003. We focus on the necessity for three dimensional data collection and reconstruction, efficient solution of the forward problem and present and future reconstruction algorithms. We also suggest common pitfalls or ``inverse crimes'' to avoid.Comment: A review paper for the 4th Workshop on Biomedical Applications of EIT, Manchester, UK, 200