Vector tomography for reconstructing electric elds with non-zero divergence in bounded domains
File(s) JCOMP-D-16-00257R1_accepted_version.pdf (749.95 KB)
Accepted version
Author(s)
Koulouri, A
Brookes, DM
Rimpiläinen, V
Type
Journal Article
Abstract
In vector tomography (VT), the aim is to reconstruct an unknown multi-dimensional vector field using line integral data. In the case of a 2-dimensional VT, two types of line integral data are usually required. These data correspond to integration of the parallel and perpendicular projection of the vector field along the integration lines and are called the longitudinal and transverse measurements, respectively. In most cases, however, the transverse measurements cannot be physically acquired. Therefore, the VT methods are typically used to reconstruct divergence-free (or source-free) velocity and flow fields that can be reconstructed solely from the longitudinal measurements. In this paper, we show how vector fields with non-zero divergence in a bounded domain can also be reconstructed from the longitudinal measurements without the need of explicitly evaluating the transverse measurements. To the best of our knowledge, VT has not previously been used for this purpose. In particular, we study low-frequency, time-harmonic electric fields generated by dipole sources in convex bounded domains which arise, for example, in electroencephalography (EEG) source imaging. We explain in detail the theoretical background, the derivation of the electric field inverse problem and the numerical approximation of the line integrals. We show that fields with non-zero divergence can be reconstructed from the longitudinal measurements with the help of two sparsity constraints that are constructed from the transverse measurements and the vector Laplace operator. As a comparison to EEG source imaging, we note that VT does not require mathematical modeling of the sources. By numerical simulations, we show that the pattern of the electric field can be correctly estimated using VT and the location of the source activity can be determined accurately from the reconstructed magnitudes of the field.
Date Issued
2016-10-19
Date Acceptance
2016-10-13
Citation
Journal of Computational Physics, 2016, 329, pp.73-90
ISSN
0021-9991
Publisher
Elsevier
Start Page
73
End Page
90
Journal / Book Title
Journal of Computational Physics
Volume
329
Copyright Statement
© 2016, Elsevier. Licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International http://creativecommons.org/licenses/by-nc-nd/4.0/
Subjects
Applied Mathematics
01 Mathematical Sciences
02 Physical Sciences
09 Engineering
Publication Status
Published
