# Localizing Brain Activity from Multiple Distinct Sources via EEG

Authors: George Dassios; Michael Doschoris; Konstantia Satrazemi
Journal: Journal of Applied Mathematics
 (2014)

---

## Abstract

An important question arousing in the framework of electroencephalography (EEG) is the possibility to recognize, by means of a recorded surface potential, the number of activated areas in the brain. In the present paper, employing a homogeneous spherical conductor serving as an approximation of the brain, we provide a criterion which determines whether the measured surface potential is evoked by a single or multiple localized neuronal excitations. We show that the uniqueness of the inverse problem for a single dipole is closely connected with attaining certain relations connecting the measured data. Further, we present the necessary and sufficient conditions which decide whether the collected data originates from a single dipole or from numerous dipoles. In the case where the EEG data arouses from multiple parallel dipoles, an isolation of the source is, in general, not possible.


---

## Body

## 1. Introduction

The electrochemically generated neuronal current in the brain and the secondary induction current sustained by the conductive cerebral tissue give rise to an electric potential, which can be recorded on the scalps surface. The dominant generators of the electric field measured on the surface of the head are graded synaptic potentials deriving mainly from the cerebral cortex [1], originating primarily from pyramidal cells uniformly oriented with apical dendrites perpendicular to the cortical surface [2]. The surface of the human cerebral cortex, comprising the outer covering layered structure of neuronal tissue of the brain, can be visualized as a highly convoluted sheet, strongly folded, consisting of fissures, sulci, and gyri. On the other hand, the basic fundamental model for the primary current distribution is a current dipole used as an equivalent source, summarizing the net effect of all microscopic currents located in a distinct region of the brain. This is a widely used approximation concept in the framework of neuroelectromagnetism [3, 4].Analysis of the registered surface electric potential provides the ability to examine cognitive processes via electroencephalography, a noninvasive, nonhazardous technique capable of following changes in neural activity with millisecond temporal resolution. An important aspect in the field of medical imaging and particular noninvasive brain imaging modalities such as EEG is the problem of identifying the position and moment of an equivalent dipole source. This is a crucial task in order to understand the mechanisms of brain response to various stimuli. However, there does not exist an exclusive source configuration for each set of electroencephalographic measurements. Hence, the corresponding inverse problem is nonunique.In addition, identifying the number of simultaneously active dipoles and their impact on the electric activity recorded on the surface of the head is a substantial nonunique task as well. Depending on the orientation of the dipoles the interaction of the related electric fields results in either an amplification of the registered EEG values or, on the other hand, cancellation effects occur [5, 6].On the computational level, detecting and localizing brain activity has been undertaken for over three decades [7–9] and various techniques and methods have been developed in order to obtain approximate solutions. These are mainly divided into two categories, namely, parametric and nonparametric methods, depending on if the number of dipole sources isa priori known or not. A review presenting those methods and techniques is given in [10, 11], whereas the key areas critically affecting the accuracy and precision of source localization are analyzed in [12].From the mathematical perspective, whenever the number of simultaneously active sources is inferred, analytic expressions determining the exact positions and moments of the corresponding dipoles can be derived for the spherical case [13] as well as for the mathematical complex but more realistic ellipsoidal geometry [14]. However, if the neuronal current exhibits two or more localized centers then a method to identify this possibility is essential.The present paper, employing a homogeneous spherical model for the brain, provides necessary and sufficient conditions for the identification of a single or multiple localized neuronal sources by reducing the problem to a set of simple algebraic equations. The simplicity of the demonstrated analysis is based on the manipulation of the exterior electric potential [15] delivering the recorded EEG data trivially as a limiting process. By linearity, the potential of every active dipole is collected and equating the matching coefficients, a set of algebraic relations are derived, establishing the basis of our approach. Finally, straightforward manipulations regarding the obtained linear system provide the conditions which have to be satisfied for the problem to be well-posed.The paper is organized as follows. In Section2 we reproduce the inversion algorithm for a single dipole providing uniqueness conditions interconnecting the EEG measurements. This development is then generalized to the case of several dipoles in Section 3. An explicit paradigm, investigating the conditions for the existence of two such dipoles, is also implemented. The results are collected in Section 4.


## 2. Inverse EEG for One Dipole

Assume a homogeneous spherical model of the brain of radiusa denoted by Ω
⊂
R
3 and let ∂
Ω be its boundary. The exterior to the brain region is denoted by Ω
c. Activation of a localized region inside the brain triggers a primary neuronal current J
p provoking an electric potential u
(
r
,
r
0
). In the case where the neuronal current is represented by a single equivalent dipole at the point r
0 with moment Q
0, we obtain J
p
=
Q
0
δ
(
r
-
r
0
), where δ denotes the dirac measure.Plonsey and Heppner [16] demonstrated that the electromagnetic activity of the brain is governed by the quasistatic theory of Maxwell’s equations; namely,
(1)
∇
×
E
=
0
,
(2)
∇
×
B
=
μ
0
(
J
p
+
σ
E
)
,
(3)
∇
·
B
=
0
,
where the magnetic permeability μ
0 is assumed to be constant everywhere in R
3.Equation (1) allows the introduction of an electric potential u such that
(4)
E
=
-
∇
u
.
Moreover, by taking the divergence of (2), we immediately conclude that theinterior electric potential u
- solves the following Neumann boundary-value problem in Ω:
(5)
σ
Δ
u
-
(
r
,
r
0
)
=
Q
0
·
∇
δ
(
r
-
r
0
)
,
r
∈
Ω
,
∂
r
u
-
(
r
,
r
0
)
=
0
,
r
∈
∂
Ω
,
where σ denotes the conductivity of the conducting medium occupying Ω. The solution regarding (5) is obtained by a straightforward expansion in terms of Legendre polynomials [15].Once the above problem is solved, knowledge of the solutionu
-
(
r
) leads to the exterior electric potential u
+
(
r
) satisfying the Dirichlet problem
(6)
Δ
u
+
(
r
,
r
0
)
=
0
,
r
∈
Ω
c
,
u
+
(
r
,
r
0
)
=
u
-
(
r
,
r
0
)
,
r
∈
∂
Ω
,
u
+
(
r
,
r
0
)
=
O
1
r
2
,
r
⟶
∞
,
provided as
(7)
u
+
(
r
,
r
0
)
=
1
σ
Q
0
·
∇
r
0
∑
n
=
1
∞
‍
∑
m
=
-
n
n
‍
1
n
r
0
n
r
n
+
1
Y
¯
n
m
(
r
^
0
)
Y
n
m
(
r
^
)
,
where Y
n
m
(
r
^
) denote the spherical harmonics, whereas the overline symbolizes complex conjugation.In order to pinpoint the positionr
0
=
(
x
01
,
x
02
,
x
03
) and moment Q
0
=
(
q
01
,
q
02
,
q
03
) of the dipole, six equations are required. As a result, the solution (7) is expanded for n
=
1,2 and the resulting relation is expressed in Cartesian coordinates, giving
(8)
u
+
r
,
r
0
=
1
4
π
σ
1
r
3
A
1
x
1
+
A
2
x
2
+
A
3
x
3
+
1
4
π
σ
1
r
5
B
1
x
1
2
+
B
2
x
2
2
+
B
3
x
3
2
+
B
4
x
1
x
2
B
1
x
1
2
+
B
2
x
2
2
+
B
3
x
3
2
+
B
4
x
1
x
2
+
B
5
x
2
x
3
+
B
6
x
3
x
1
+
O
1
r
4
,
where
(9)
(
A
1
,
A
2
,
A
3
)
=
3
(
q
01
,
q
02
,
q
03
)
,
(10)
B
1
=
15
2
q
01
x
01
-
5
2
Q
0
·
r
0
,
(11)
B
2
=
15
2
q
02
x
02
-
5
2
Q
0
·
r
0
,
(12)
B
3
=
15
2
q
03
x
03
-
5
2
Q
0
·
r
0
,
(13)
B
4
=
15
2
q
01
x
02
+
q
02
x
01
,
(14)
B
5
=
15
2
q
02
x
03
+
q
03
x
02
,
(15)
B
6
=
15
2
q
03
x
01
+
q
01
x
03
,
as well as the harmonicity condition
(16)
B
1
+
B
2
+
B
3
=
0
.Merging above equations into a single system yields(17)
A
r
0
=
B
given that
(18)
A
=
-
A
1
2
A
2
-
A
3
-
A
1
-
A
2
2
A
3
A
2
A
1
0
A
3
0
A
1
0
A
3
A
2
,
r
0
=
x
01
x
02
x
03
,
B
=
6
5
B
2
B
3
B
4
B
5
B
6
.A simple investigation of the latter shows that(19)
x
01
=
3
2
A
1
A
2
A
3
A
1
a
·
b
,
x
02
=
3
2
A
1
A
2
A
3
A
2
a
·
R
x
3
(
π
)
b
,
x
03
=
3
2
A
1
A
2
A
3
A
3
a
·
R
x
2
(
π
)
b
,
where a
=
(
A
1
,
A
2
,
A
3
)
⊤, b
=
(
-
B
4
,
B
5
,
B
6
)
⊤, and
(20)
R
x
3
(
α
)
=
cos
⁡
α
-
sin
α
0
sin
α
cos
⁡
α
0
0
0
1
,
R
x
2
(
α
)
=
cos
⁡
α
0
sin
α
0
1
0
-
sin
α
0
cos
⁡
α
are the rotation matrices about the x
3- and x
2-axis, respectively.Importantly, the above analysis reveals that the dipoles positionr
0 specified by relations (19) is uniqueonly if the recorded values for the coefficients A
l, l
=
1,2
,
3, and B
k, k
=
1,2
,
3,4
,
5,6, satisfy the following expressions:
(21)
2
A
1
A
2
A
3
(
B
3
-
B
2
)
+
3
A
1
A
2
2
-
A
3
2
B
5
+
3
A
2
2
+
A
3
2
A
3
B
4
-
A
2
B
6
=
0
,
(22)
-
2
A
1
A
2
A
3
(
B
2
+
2
B
3
)
+
3
A
1
2
+
A
3
2
A
1
B
5
-
A
3
B
4
+
3
A
2
A
3
2
-
A
1
2
B
6
=
0
.
Bearing in mind condition (16), relations (21) and (22) are easily obtained by solving any three equations of (10)–(15) with respect to r
0 and substituting the solution into the remaining two.Closing this section we note that a major drawback is that knowledge of the surface potentialu
(
a
r
^
,
r
0
) does not automatically imply knowledge of the coefficients present in (8) as well. A connection to the spherical case has to be made. This is achieved in view of the following formulas [15]:
(23)
u
(
a
r
^
,
r
0
)
=
∑
n
=
1
∞
‍
∑
m
=
-
n
n
‍
Γ
n
m
(
r
0
)
Y
n
m
(
r
^
)
,
Γ
n
m
(
r
0
)
=
1
σ
n
a
n
+
1
Q
0
·
∇
r
0
r
0
n
Y
¯
n
m
(
r
^
0
)
in conjunction with the orthogonality condition
(24)
Γ
n
m
r
0
=
∮
S
2
‍
u
a
r
^
,
r
0
Y
¯
n
m
r
^
d
Ω
r
^
,
where S
2 denotes the boundary of the unit sphere and d
Ω
(
r
^
) is the solid angle element.Relation (24) indicates the simple fact that measuring u
(
a
r
^
,
r
0
) provides information about the coefficients Γ
n
m
(
r
0
) alone. In order to obtain the values of A
l, l
=
1,2
,
3, as well as B
k, k
=
1,2
,
3,4
,
5,6, in the limit as r tends to a, relation (23) has to be expanded for n
=
1,2, providing the relations
(25)
A
j
=
3
4
π
a
2
∮
S
2
‍
u
(
a
r
^
,
r
0
)
x
j
d
Ω
(
r
^
)
,
j
=
1,2
,
3
,
B
1
=
5
8
π
a
4
∮
S
2
‍
u
(
a
r
^
,
r
0
)
2
x
1
2
-
x
2
2
-
x
3
2
d
Ω
(
r
^
)
,
B
2
=
5
8
π
a
4
∮
S
2
‍
u
(
a
r
^
,
r
0
)
2
x
2
2
-
x
3
2
-
x
1
2
d
Ω
(
r
^
)
,
B
3
=
5
8
π
a
4
∮
S
2
‍
u
(
a
r
^
,
r
0
)
2
x
3
2
-
x
1
2
-
x
2
2
d
Ω
(
r
^
)
,
B
4
=
15
4
π
a
4
∮
S
2
‍
u
(
a
r
^
,
r
0
)
x
1
x
2
d
Ω
(
r
^
)
,
B
5
=
15
4
π
a
4
∮
S
2
‍
u
(
a
r
^
,
r
0
)
x
2
x
3
d
Ω
(
r
^
)
,
B
6
=
15
4
π
a
4
∮
S
2
‍
u
(
a
r
^
,
r
0
)
x
3
x
1
d
Ω
(
r
^
)
.


## 3. EEG Inversion for Multiple Dipoles

In what follows, we consider a number of simultaneously active dipoles(
r
j
,
Q
j
), j
=
1,2
,
…
,
N. Due to linearity, the electric potential generated outside the sphere is
(26)
u
+
r
;
r
0
,
Q
0
=
∑
j
=
1
N
‍
u
+
r
;
r
j
,
Q
j
,
where the dipole (
r
0
,
Q
0
) represents the contribution of a fictitious dipole inside the brain, providing the same potential as the dipoles (
r
j
,
Q
j
), j
=
1,2
,
…
,
N, collectively.Equating coefficients in (26) given that r
j
=
(
x
j
1
,
x
j
2
,
x
j
3
) and Q
j
=
(
q
j
1
,
q
j
2
,
q
j
3
) leads to an identical system as (17) and (18),
(27)
(
A
1
,
A
2
,
A
3
)
=
3
∑
j
=
1
N
‍
q
j
1
,
∑
j
=
1
N
‍
q
j
2
,
∑
j
=
1
N
‍
q
j
3
,
(28)
B
1
=
∑
j
=
1
N
‍
15
2
q
j
1
x
j
1
-
5
2
Q
j
·
r
j
,
(29)
B
2
=
∑
j
=
1
N
‍
15
2
q
j
2
x
j
2
-
5
2
Q
j
·
r
j
,
(30)
B
3
=
∑
j
=
1
N
‍
15
2
q
j
3
x
j
3
-
5
2
Q
j
·
r
j
,
(31)
B
4
=
15
2
∑
j
=
1
N
‍
q
j
1
x
j
2
+
q
j
2
x
j
1
,
(32)
B
5
=
15
2
∑
j
=
1
N
‍
q
j
2
x
j
3
+
q
j
3
x
j
2
,
(33)
B
6
=
15
2
∑
j
=
1
N
‍
q
j
3
x
j
1
+
q
j
1
x
j
3
.
Regardless of the number of introduced dipole sources, the constraint (16) as well as the uniqueness conditions (21) and (22) are still binding.Let us investigate in the sequel the necessary circumstances regarding theN dipoles (
r
j
,
Q
j
) for which the uniqueness conditions (21) and (22) are satisfied. Substituting relations (28)–(33) into (21) and (22), respectively, and apprising the fact that the resulting equations are fulfilled if the coefficients of each variable vanishes lead to
(34)
q
02
2
+
q
03
2
q
03
q
j
2
-
q
02
q
j
3
=
0
,
(35)
q
01
q
02
q
02
q
j
3
-
q
03
q
j
2
+
q
02
q
03
q
02
q
j
1
-
q
01
q
j
2
+
q
03
2
q
03
q
j
1
-
q
01
q
j
3
=
0
,
(36)
q
01
q
03
q
02
q
j
3
-
q
03
q
j
2
+
q
02
q
03
q
01
q
j
3
-
q
03
q
j
1
+
q
02
2
q
01
q
j
2
-
q
02
q
j
1
=
0
,
(37)
q
01
q
03
q
02
q
j
1
-
q
01
q
j
2
+
q
01
q
02
q
03
q
j
1
-
q
01
q
j
3
+
q
03
2
q
02
q
j
3
-
q
03
q
j
2
=
0
,
(38)
q
01
2
+
q
03
2
q
01
q
j
3
-
q
03
q
j
1
=
0
,
(39)
q
01
q
03
q
03
q
j
2
-
q
02
q
j
3
+
q
02
q
03
q
03
q
j
1
-
q
01
q
j
3
+
q
01
2
q
01
q
j
2
-
q
02
q
j
1
=
0
,
for every j
=
1,2
,
…
,
N.Above equations are trivially satisfied in the case whereQ
0
=
0. However, this particular choice leads to u
+
(
r
;
r
0
,
Q
0
)
=
0 and turns the dipoles (
r
j
,
Q
j
), j
=
1,2
,
…
,
N, to silent sources, that is, brain activity impossible to be recorded via electroencephalographic measurements. On the other hand, whenever Q
0
≠
0  (34) and (35) are attained if
(40)
q
03
q
j
2
-
q
02
q
j
3
=
0
,
q
02
q
j
1
-
q
01
q
j
2
=
0
,
q
03
q
j
1
-
q
01
q
j
3
=
0
.
According to (9) and (27)(41)
q
0
s
=
∑
i
=
1
N
‍
q
i
s
,
s
=
1,2
,
3
and thus relations (40) read
(42)
q
j
2
∑
i
=
1
N
‍
q
i
3
-
q
j
3
∑
i
=
1
N
‍
q
i
2
=
0
,
j
=
1,2
,
…
,
N
,
q
j
1
∑
i
=
1
N
‍
q
i
2
-
q
j
2
∑
i
=
1
N
‍
q
i
1
=
0
,
j
=
1,2
,
…
,
N
,
q
j
1
∑
i
=
1
N
‍
q
i
3
-
q
j
3
∑
i
=
1
N
‍
q
i
1
=
0
,
j
=
1,2
,
…
,
N
.
It is not hard to show that aforementioned formulas are simultaneously satisfied provided that
(43)
x
^
s
·
Q
j
×
∑
i
=
1
i
≠
j
N
‍
Q
i
=
0
,
s
=
1,2
,
3
,
j
=
1,2
,
…
,
N
,
comprising N
-
1 equations for each direction x
^
s, s
=
1,2
,
3, or in other words, wherever the dipoles are parallel to each other; namely,
(44)
Q
i
×
Q
j
=
0
,
i
,
j
=
1,2
,
…
,
N
,
i
≠
j
.
Other occurrences, for example, q
01
=
0,  q
02
,
q
03
≠
0, and so forth, result in identical conclusions.In view of Section2, it is feasible to acquire the positions and moments of the N dipoles if at least 6
N of the corresponding coefficients regarding the RHS of (26), that is,
(45)
Γ
n
m
(
r
j
)
=
1
σ
n
a
n
+
1
∑
j
=
1
N
‍
Q
j
·
∇
r
j
r
j
n
Y
¯
n
m
(
r
^
j
)
are identified. In the interest of completeness, we note that a simple but strictly theoretical constraint which decides whether recorded EEG data could originate from multiple dipole sources is that the remaining coefficients (45) must be compatible with the given measurements [13].


### 3.1. Example: Inversion Algorithm for Two Dipoles

For the sake of demonstrating details of the presented analysis we establish to the full extent the case for two concurrently active dipoles.Equating coefficients in (26) leads to the following relations:
(46)
Q
0
=
Q
1
+
Q
2
,
(47)
3
q
01
x
01
-
Q
0
·
r
0
=
3
q
11
x
11
-
Q
1
·
r
1
+
3
q
21
x
21
-
Q
2
·
r
2
,
(48)
3
q
02
x
02
-
Q
0
·
r
0
=
3
q
12
x
12
-
Q
1
·
r
1
+
3
q
22
x
22
-
Q
2
·
r
2
,
(49)
3
q
03
x
03
-
Q
0
·
r
0
=
3
q
13
x
13
-
Q
1
·
r
1
+
3
q
23
x
23
-
Q
2
·
r
2
,
(50)
q
01
x
02
+
q
02
x
01
=
q
11
x
12
+
q
12
x
11
+
q
21
x
22
+
q
22
x
21
,
(51)
q
02
x
03
+
q
03
x
02
=
q
12
x
13
+
q
13
x
12
+
q
22
x
23
+
q
23
x
22
,
(52)
q
03
x
01
+
q
01
x
03
=
q
13
x
11
+
q
11
x
13
+
q
23
x
21
+
q
21
x
23
.Solving equations (50)–(52) with respect to r
0
=
(
x
01
,
x
02
,
x
03
) we arrive at
(53)
x
01
=
1
2
q
02
q
03
×
q
12
q
03
+
q
13
q
02
x
0
+
q
22
q
03
+
q
23
q
02
x
21
+
1
2
q
02
q
03
×
-
q
13
q
11
q
23
q
21
(
x
12
-
x
22
)
+
q
11
q
12
q
21
q
22
(
x
13
-
x
23
)
,
x
02
=
1
2
q
01
q
03
×
q
11
q
03
+
q
13
q
01
x
12
+
q
21
q
03
+
q
23
q
01
x
22
+
1
2
q
01
q
03
×
q
12
q
13
q
22
q
23
(
x
11
-
x
21
)
-
q
11
q
12
q
21
q
22
(
x
13
-
x
23
)
,
x
03
=
1
2
q
01
q
02
×
q
11
q
02
+
q
12
q
01
x
13
+
q
21
q
02
+
q
22
q
01
x
23
+
1
2
q
01
q
02
×
-
q
12
q
13
q
22
q
23
(
x
11
-
x
21
)
+
q
13
q
11
q
23
q
21
(
x
12
-
x
22
)
.On the other hand, subtracting (48) from (47) and combining the resulting equation with (50) to eliminate x
02 we obtain
(54)
x
~
01
=
1
q
01
2
+
q
02
2
×
q
11
q
01
+
q
12
q
02
x
11
+
q
21
q
01
+
q
22
q
02
x
21
+
1
q
01
2
+
q
02
2
q
11
q
12
q
21
q
22
(
x
12
-
x
22
)
.
Similarly, subtracting (49) from (48) and eliminating x
03 between the resulting equation and (51) we obtain
(55)
x
~
02
=
1
q
02
2
+
q
03
2
×
q
12
q
02
+
q
13
q
03
x
12
+
q
22
q
02
+
q
23
q
03
x
22
+
1
q
02
2
+
q
03
2
q
12
q
13
q
22
q
23
(
x
13
-
x
23
)
.
Finally, subtracting (47) from (49) and using the resulting equation to eliminate x
01 from (52) we arrive at
(56)
x
~
03
=
1
q
01
2
+
q
03
2
×
q
13
q
03
+
q
11
q
01
x
13
+
q
23
q
03
+
q
21
q
01
x
23
+
1
q
01
2
+
q
03
2
q
13
q
11
q
23
q
21
(
x
11
-
x
21
)
.
Obviously, the point r
0
=
(
x
01
,
x
02
,
x
03
) defined by relations (53) does not coincide with the point r
~
0
=
(
x
~
01
,
x
~
02
,
x
~
03
) defined from the solutions (54)–(56). This means that the system (47)–(52) is not compatible. That is, the solution we obtain depends on the choice of the equations we choose to calculate the x
0
j’s or, in other words, if we pick up any three equations among the six equations of the system, (with the exception of the first three which are linearly dependent) to calculate x
0
j, j
=
1,2
,
3, the solution we obtain does not necessarily satisfy the remaining equations. This observation tells us that the system (47)–(52) “knows” that it does not represent a single dipole and therefore it provides a criterion for deciding whether there is a single point excitation or it is an excitation that is due to more than one dipole. It is easily shown that if the moments of the two dipoles are parallel, that is, if
(57)
Q
2
=
γ
Q
1
,
γ
∈
R
,
then the representations for r
0 and r
~
0 coincide, and therefore, in this case, it is impossible to decide if the dipoles are one or two. In fact, under the condition (57), we obtain
(58)
Q
1
×
Q
2
=
0
,
which means that all the 2
×
2 determinants in expressions (53)–(56) vanish, and the common solution is given by
(59)
r
0
=
r
~
0
=
r
1
+
γ
r
2
1
+
γ
which is a point on the line segment that connects r
1 and r
2. Indeed,
(60)
r
0
-
r
1
×
r
0
-
r
2
=
0
,
which implies that the vectors r
0
-
r
1 and r
0
-
r
2 lie on the same line.From the physical point of view, the reason why parallel dipoles are indistinguishable is due the fact that the potential lines they generate are also similarly located and therefore in the exterior space the field has the pattern of a single current dipole. In any other case, this simple pattern of the field is destroyed and it is easy to identify a more complicated source. In order to see the discrepancy between the solutionr
0 and the solution r
~
0 we look at the difference
(61)
x
01
-
x
~
01
=
q
12
q
03
+
q
13
q
02
2
q
02
q
03
-
q
11
q
01
+
q
12
q
02
q
01
2
+
q
02
2
×
x
11
-
x
21
+
q
23
q
11
-
q
13
q
21
2
q
02
q
03
-
q
12
q
21
-
q
22
q
11
q
01
2
+
q
02
2
×
x
12
-
x
22
+
q
11
q
22
-
q
21
q
12
2
q
02
q
03
x
13
-
x
23
.
Since the moments Q
1 and Q
2 are independent of the positions r
1 and r
2 of the dipoles, it follows that x
01
=
x
~
01 if and only if the coefficients of x
11
-
x
21, x
12
-
x
22 and x
13
-
x
23 are equal to zero. After straightforward calculations we obtain Q
1
×
Q
2
=
0 if q
02
≠
0 and the two dipoles are parallel. The case where q
02
=
0 corresponds to Q
1
=
-
Q
2 and therefore to Q
0
=
0, namely, a significant cancellation effect takes place.The impact of individual dipole configurations (parallel or perpendicular) on the measured scalp potential is demonstrated in Figures1 and 2. Two dipoles of equal strength are situated in gyri in adjacent cortical regions (configurations Figures 1(b)–1(d)) and are compared with the corresponding single dipole excitation (configuration Figure 1(a)). We observe that, in the case of synchronous activation of two dipoles (Figure 1(b)) in close proximity, compared to a single excitation (Figure 1(a)), an amplification of the recorded scalp potential is detected. Locally, the increase of magnitude can be substantial (Figure 2(b)). On the contrary, instances (Figures 1(c) and 1(d)) result in a diminishing of magnitude (Figures 2(c) and 2(d)) as expected.Figure 1
Electroencephalographic recordings as a result of two simultaneously active dipoles in three distinct arrangements: (b) parallel, (c) perpendicular, and (d) antiparallel. The corresponding neuronal excitation occurs on the cortical surface in neighboring regions. The presented electric potential is plotted for a fixed azimuthal angle equal to3
π
/
8.The surface electric potential as a result of the activation of (a) a single dipole, (b) two parallel dipoles, (c) two perpendicular dipoles, and (d) two antiparallel dipoles in close proximity.
(a)
(b)
(c)
(d)


## 3.1. Example: Inversion Algorithm for Two Dipoles

For the sake of demonstrating details of the presented analysis we establish to the full extent the case for two concurrently active dipoles.Equating coefficients in (26) leads to the following relations:
(46)
Q
0
=
Q
1
+
Q
2
,
(47)
3
q
01
x
01
-
Q
0
·
r
0
=
3
q
11
x
11
-
Q
1
·
r
1
+
3
q
21
x
21
-
Q
2
·
r
2
,
(48)
3
q
02
x
02
-
Q
0
·
r
0
=
3
q
12
x
12
-
Q
1
·
r
1
+
3
q
22
x
22
-
Q
2
·
r
2
,
(49)
3
q
03
x
03
-
Q
0
·
r
0
=
3
q
13
x
13
-
Q
1
·
r
1
+
3
q
23
x
23
-
Q
2
·
r
2
,
(50)
q
01
x
02
+
q
02
x
01
=
q
11
x
12
+
q
12
x
11
+
q
21
x
22
+
q
22
x
21
,
(51)
q
02
x
03
+
q
03
x
02
=
q
12
x
13
+
q
13
x
12
+
q
22
x
23
+
q
23
x
22
,
(52)
q
03
x
01
+
q
01
x
03
=
q
13
x
11
+
q
11
x
13
+
q
23
x
21
+
q
21
x
23
.Solving equations (50)–(52) with respect to r
0
=
(
x
01
,
x
02
,
x
03
) we arrive at
(53)
x
01
=
1
2
q
02
q
03
×
q
12
q
03
+
q
13
q
02
x
0
+
q
22
q
03
+
q
23
q
02
x
21
+
1
2
q
02
q
03
×
-
q
13
q
11
q
23
q
21
(
x
12
-
x
22
)
+
q
11
q
12
q
21
q
22
(
x
13
-
x
23
)
,
x
02
=
1
2
q
01
q
03
×
q
11
q
03
+
q
13
q
01
x
12
+
q
21
q
03
+
q
23
q
01
x
22
+
1
2
q
01
q
03
×
q
12
q
13
q
22
q
23
(
x
11
-
x
21
)
-
q
11
q
12
q
21
q
22
(
x
13
-
x
23
)
,
x
03
=
1
2
q
01
q
02
×
q
11
q
02
+
q
12
q
01
x
13
+
q
21
q
02
+
q
22
q
01
x
23
+
1
2
q
01
q
02
×
-
q
12
q
13
q
22
q
23
(
x
11
-
x
21
)
+
q
13
q
11
q
23
q
21
(
x
12
-
x
22
)
.On the other hand, subtracting (48) from (47) and combining the resulting equation with (50) to eliminate x
02 we obtain
(54)
x
~
01
=
1
q
01
2
+
q
02
2
×
q
11
q
01
+
q
12
q
02
x
11
+
q
21
q
01
+
q
22
q
02
x
21
+
1
q
01
2
+
q
02
2
q
11
q
12
q
21
q
22
(
x
12
-
x
22
)
.
Similarly, subtracting (49) from (48) and eliminating x
03 between the resulting equation and (51) we obtain
(55)
x
~
02
=
1
q
02
2
+
q
03
2
×
q
12
q
02
+
q
13
q
03
x
12
+
q
22
q
02
+
q
23
q
03
x
22
+
1
q
02
2
+
q
03
2
q
12
q
13
q
22
q
23
(
x
13
-
x
23
)
.
Finally, subtracting (47) from (49) and using the resulting equation to eliminate x
01 from (52) we arrive at
(56)
x
~
03
=
1
q
01
2
+
q
03
2
×
q
13
q
03
+
q
11
q
01
x
13
+
q
23
q
03
+
q
21
q
01
x
23
+
1
q
01
2
+
q
03
2
q
13
q
11
q
23
q
21
(
x
11
-
x
21
)
.
Obviously, the point r
0
=
(
x
01
,
x
02
,
x
03
) defined by relations (53) does not coincide with the point r
~
0
=
(
x
~
01
,
x
~
02
,
x
~
03
) defined from the solutions (54)–(56). This means that the system (47)–(52) is not compatible. That is, the solution we obtain depends on the choice of the equations we choose to calculate the x
0
j’s or, in other words, if we pick up any three equations among the six equations of the system, (with the exception of the first three which are linearly dependent) to calculate x
0
j, j
=
1,2
,
3, the solution we obtain does not necessarily satisfy the remaining equations. This observation tells us that the system (47)–(52) “knows” that it does not represent a single dipole and therefore it provides a criterion for deciding whether there is a single point excitation or it is an excitation that is due to more than one dipole. It is easily shown that if the moments of the two dipoles are parallel, that is, if
(57)
Q
2
=
γ
Q
1
,
γ
∈
R
,
then the representations for r
0 and r
~
0 coincide, and therefore, in this case, it is impossible to decide if the dipoles are one or two. In fact, under the condition (57), we obtain
(58)
Q
1
×
Q
2
=
0
,
which means that all the 2
×
2 determinants in expressions (53)–(56) vanish, and the common solution is given by
(59)
r
0
=
r
~
0
=
r
1
+
γ
r
2
1
+
γ
which is a point on the line segment that connects r
1 and r
2. Indeed,
(60)
r
0
-
r
1
×
r
0
-
r
2
=
0
,
which implies that the vectors r
0
-
r
1 and r
0
-
r
2 lie on the same line.From the physical point of view, the reason why parallel dipoles are indistinguishable is due the fact that the potential lines they generate are also similarly located and therefore in the exterior space the field has the pattern of a single current dipole. In any other case, this simple pattern of the field is destroyed and it is easy to identify a more complicated source. In order to see the discrepancy between the solutionr
0 and the solution r
~
0 we look at the difference
(61)
x
01
-
x
~
01
=
q
12
q
03
+
q
13
q
02
2
q
02
q
03
-
q
11
q
01
+
q
12
q
02
q
01
2
+
q
02
2
×
x
11
-
x
21
+
q
23
q
11
-
q
13
q
21
2
q
02
q
03
-
q
12
q
21
-
q
22
q
11
q
01
2
+
q
02
2
×
x
12
-
x
22
+
q
11
q
22
-
q
21
q
12
2
q
02
q
03
x
13
-
x
23
.
Since the moments Q
1 and Q
2 are independent of the positions r
1 and r
2 of the dipoles, it follows that x
01
=
x
~
01 if and only if the coefficients of x
11
-
x
21, x
12
-
x
22 and x
13
-
x
23 are equal to zero. After straightforward calculations we obtain Q
1
×
Q
2
=
0 if q
02
≠
0 and the two dipoles are parallel. The case where q
02
=
0 corresponds to Q
1
=
-
Q
2 and therefore to Q
0
=
0, namely, a significant cancellation effect takes place.The impact of individual dipole configurations (parallel or perpendicular) on the measured scalp potential is demonstrated in Figures1 and 2. Two dipoles of equal strength are situated in gyri in adjacent cortical regions (configurations Figures 1(b)–1(d)) and are compared with the corresponding single dipole excitation (configuration Figure 1(a)). We observe that, in the case of synchronous activation of two dipoles (Figure 1(b)) in close proximity, compared to a single excitation (Figure 1(a)), an amplification of the recorded scalp potential is detected. Locally, the increase of magnitude can be substantial (Figure 2(b)). On the contrary, instances (Figures 1(c) and 1(d)) result in a diminishing of magnitude (Figures 2(c) and 2(d)) as expected.Figure 1
Electroencephalographic recordings as a result of two simultaneously active dipoles in three distinct arrangements: (b) parallel, (c) perpendicular, and (d) antiparallel. The corresponding neuronal excitation occurs on the cortical surface in neighboring regions. The presented electric potential is plotted for a fixed azimuthal angle equal to3
π
/
8.The surface electric potential as a result of the activation of (a) a single dipole, (b) two parallel dipoles, (c) two perpendicular dipoles, and (d) two antiparallel dipoles in close proximity.
(a)
(b)
(c)
(d)


## 4. Conclusions

An analytic criterion deciding the existence of a single dipole or multiple dipoles is presented. The analytic algorithm, on which the criterion resides, is derived by expressing the generated electric surface potential in Cartesian coordinates. Then, by equating the coefficients of the Cartesian monomials with the corresponding known coefficients from the Cartesian expansion of the recorded potential, we obtain a set of algebraic equations which establish the basis of our analysis. This system, obviously overdetermined, exhibits a unique solution only if certain constraints, provided via (17) and (18), are met. In the case of a single dipole relations (17) and (18) are trivially satisfied. On the other hand, if two or more dipoles are simultaneously active relations (17) and (18), which connect the electroencephalographic recordings, are fulfilled only if the dipoles are parallel. In this case it is impossible to decide if the EEG measurements are evoked by a single dipole or by a finite number of dipoles. In any other instance, the precise number of active sources can be decided. Summarizing, the necessary and sufficient conditions which decide whether the collected EEG data originates from a single dipole or from numerous dipoles are presented. In the event where the data arouses from multiple parallel dipoles, an isolation of the source is, in general, not possible.



---

*Source: 232747-2014-11-23.xml*