- 1School of Mechanical Engineering, Xi’an Jiaotong University, Xi’an, China
- 2Shaanxi Key Laboratory of Intelligent Robot, Xi’an Jiaotong University, Xi’an, China
- 3School of Mechanical Engineering, Xinjiang University, Wulumuqi, China
Electroencephalogram (EEG) modeling in brain-computer interface (BCI) provides a theoretical foundation for its development. However, limited by the lack of guidelines in model parameter selection and the inability to obtain personal tissue information in practice, EEG modeling in BCI is mainly focused on the theoretical qualitative level which shows a gap between the theory and its application. Based on such problems, this work combined the surface EEG simulation with a converter based on the generative adversarial network (GAN), to establish the connection from simulated EEG to its application in BCI classification. For the scalp EEGs modeling, a mathematical model was built according to the physics of surface EEG, which consisted of the parallel 3-population neural mass model, the equivalent dipole, and the forward computation. For application, a converter based on the conditional GAN was designed, to transfer the simulated theoretical-only EEG to its practical version, in the lack of individual bio-information. To verify the feasibility, based on the latest microexpression-assisted BCI paradigm proposed by our group, the converted simulated EEGs were used in the training of BCI classifiers. The results indicated that, compared with training with insufficient real data, by adding the simulated EEGs, the overall performance showed a significant improvement (P = 0.04 < 0.05), and the test performance can be improved by 2.17% ± 4.23, in which the largest increase was up to 12.60% ± 1.81. Through this work, the link from theoretical EEG simulation to BCI classification has been initially established, providing an enhanced novel solution for the application of EEG modeling in BCI.
Introduction
Non-invasive brain-computer interface (BCI) as a branch in neuroscience, has continuously been a hotspot in recent decades. The focal point of the research has been shifted from the signal detection at first (Berger, 1929), to the paradigm proposal (Herrmann, 2001; Pfurtscheller and Neuper, 2001; Wolpaw et al., 2002; Polich, 2007), then the EEG decoding and classification (Koelstra et al., 2012; Chen et al., 2015; Muller-Putz et al., 2016). With the growing publication of related studies, attention began to be paid to the generation mechanism and modelling of EEG in BCI, attempting to provide more solid theoretical foundations.
Unlike clinical medicine (Sanz-Leon et al., 2015), EEG modeling launched by BCI usually adopted the neuronal population level-based modeling approach and has been more concerned with the theoretical qualitative simulation. The “neural population” had first been put forward in the 1970s, assuming that all nervous processes can be dependent upon the interaction of excitatory and inhibitory cells (Freeman, 1972, 1973, 1975; Wilson and Cowan, 1972, 1973). Later, it was adopted by Lopes da Silva in modeling the generation of rhythmic activity (Lopesdas et al., 1974; Lopes da Silva et al., 1976). In 1995, the neural mass model (NMM) was proposed by Jansen and Rit (1995) in the mathematical computation of visual evoked potential (Jansen et al., 1993). Since then, Jasen-Rit NMM established the connection between deep EEG modeling and the BCI paradigm, and became popular in terms of BCI-related EEG modeling. In 2006, Wendling NMM (Wendling et al., 2002), as an improved Jasen-Rit, was utilized in the simulation of cortical activity during motor tasks (Zavaglia et al., 2006). Afterward, NMM was applied in the EEG simulation based on the scene graph steady-state visual evoked potentials BCI paradigm (Li et al., 2018a) and the EEG driven by facial expression (Zhang et al., 2016).
Neural mass model-based EEG modeling does evidence the signal response under the BCI paradigm to some extent, but by far it has not formed guidance for the paradigm design and it is still a certain distance away from application in BCI. Firstly, according to the anatomy, the output of NMM is the postsynaptic membrane potential, which can be regarded as the deep source of surface EEG, but not the collected one. Certainly, under well-selected parameters, phenomena in surface EEG can be reproduced (i.e., significant frequency distribution (Jansen and Rit, 1995), rhythmic activity (Ursino et al., 2010), energy undulation (Zavaglia et al., 2006)); However, the selection of NMM parameters in BCI relied largely on manual adjustment and consumed large efforts. Another challenge is the forward computation from the inner source to the scalp signal, which especially requires one’s tissue conductivity and geometry (Cosandier-Rimele et al., 2010). Although the simplified three-shell concentric spherical head model can be adopted (Berg and Scherg, 1994), without extra modification the result still stayed at the theoretical level (meanwhile, it is unaffordable and impossible for investigators to obtain complex tissue data from every user in practice). Obstacles above (i.e., the selection of NMM parameters, and the lack of individual biometric data) limited the application of EEG simulation in BCI. However, without the application, modeling, and simulation of EEG, it lacked practical meaning, was stuck in the theoretical level, and felt within an inch of enough.
In the other aspect, building the link between EEG simulation and its application in BCI is also valuable for BCI. Simulated EEG has its natural advantage in data amount. If the simulated signal can be applied to BCI training, it can become another feasible solution to reduce the amount of pre-collected training data and alleviate the problem of data insufficiency. To build such a link, the core problem lies in how to realize the personal scalp EEG simulation without detailed individual tissue data and the fine tuning of the NMM parameter. Encouraged by the work of probabilistic forecasting (Koochali et al., 2019), data augmentation (Fan et al., 2020; Luo et al., 2020), and EEG feature generating (Krishna et al., 2021) with generative adversarial networks (GAN), in this work, we proposed an EEG simulation method with conditional GAN combined to compensate for the lack of tissue information, hence establishing the link between EEG modeling and its application.
In this paper, emphasis was laid on surface EEG modeling and realizing its application in BCI classification. In the methodology section, the generation physics of surface EEG were explained firstly thus deriving the mathematical model of scalp EEG; Then, focusing on the realistic dilemma, a conditional GAN converter was proposed to transform the simulated EEGs from the theoretical-only to its practical version, to overcome the lack of individual bio-information. In the data acquisition section, the experimental EEGs acquisition and the paradigm-related parameters were stated. In the results section, the feasibility of applying the EEG simulation to the BCI classification was estimated. As followed, the discussion section discussed the performance and the limitation of this work, and the conclusion was summarized in the conclusion section.
Methodology
To build the connection between EEG modeling and its application in BCI, the method started from the scalp EEG modeling, and then a converter based on the conditional GAN was combined to transfer the theoretical-only simulated-EEG to its practical version. The overall schematics diagram is illustrated in Figure 1.
The Mathematical Model of Surface Electroencephalogram
Electroencephalogram measures the large-scale simultaneous activation of the brain neurons through the surface electrodes placed on the scalp (Berger, 1929). It is generally believed that the extracellular potential field generated by postsynaptic potentials is the source of an EEG (Baillet et al., 2001; Hallez et al., 2007). Since only the regular arrangements of neurons can amplify their extracellular potential field to a measurable extent, the extracellular potential fields raised by the postsynaptic potential of pyramidal cells are generally admitted to be the generator of the EEG. The extracellular potential is caused by the migration of positively charged ions (Schaul, 1998), as illustrated in Figure 2. From the extracellular environment to the brain surface, the spread pathway contains the gray matter, the white matter, the cerebrospinal fluid (CSF), the skull, and the scalp. Only by passing through various tissues, can the electrical activity deep within the brain be finally detected by EEG electrodes. According to the physics of surface EEG, the modeling process of scalp EEG can be divided into three nodes: the postsynaptic membrane potential, the extracellular electric field and the forward computation from deep EEG to scalp EEG. Following are the details.
Modeling of Postsynaptic Membrane Potential
According to the physics, the modeling starts from the postsynaptic membrane potential. For being generated by large-scaled simultaneous activation, instead of isolated neurons, the cortical activity is simulated with the population-level based model: NMM.
In NMM, one population represents neurons that are lumped together with the same membrane potential. A single population includes the interaction among four neural subgroups (Wendling et al., 2002), pyramidal cells, excitatory interneurons, inhibitory interneurons with slow synaptic kinetics, and inhibitory interneurons with faster synaptic kinetics, as illustrated in Figure 3.
In Figure 3, a static nonlinearity sigmoid relationship is used to convert the average postsynaptic membrane potentials to an average spike density. Then a second-order transfer function h(t) transfers the presynaptic spike density to the postsynaptic membrane potential. This second-order representation of synaptic effect was firstly derived by Rotterdam (van Rotterdam et al., 1982) and is still in use. After simplification, the mathematical equations of Wendling’s NMM (Wendling et al., 2002) are
where A, B, & G are the average gain in synaptic effect, and a, b, & g are the time constant; C1∼C7 represent the connectivity constants between neuron groups; P(t) is the sum of exogenous contributions. The sigmoid function is
where e0 determines the maximum firing rate of the neural population, s0 is the mean firing threshold, and r is the slope.
Since a single population can only produce a unimodal spectrum, multiple coupled neural populations with different ways of connecting were studied (Wendling et al., 2000) to describe the neural activity in a wide bandwidth. By separating typical EEG into three bands as the low (4-12 Hz), the medium (12-30 Hz), and the high (>30 Hz) (Zavaglia et al., 2006), parallel feedforward schema of three populations were arranged to mimic the overall complexity of EEG, as illustrated in Figure 4. The postsynaptic membrane potential Vout under multi-parallel-populations is
where superscript L, M, and H indicate the three populations in representative EEG bands, and wk(k = L, M, H) stands for the weights.
Equivalent Model of Extracellular Electric Field
The generator of the extracellular electric field can be viewed as a simplified electric model with two current monopoles, as demonstrated in Figure 5. At the apical dendrite side, a current sink is used to describe the influx of positive ions. At the cell body, a current source is placed to represent the injection of positively charged ions into the extracellular environment.
The frequency range and time-variation of measured EEGs indicate that there is no charge piled up in the conducting extracellular volume (Plonsey and Heppner, 1967). Thereby the simplified electric model can be regarded as a current dipole that consists of two equal anisotropic charges. Following the population level-based modeling approach in “Modeling of Postsynaptic Membrane Potential,” the simultaneous electrically active pyramidal cells in a small patch can be represented as one equivalent dipole (He et al., 2002).
The electric-field intensity of the dipole in a quasi-static field can be driven utilizing the gradient operator
where E is the vector electric-field intensity in V/m, ∇ is the gradient operator, and U is the scalar potential field. For a dipole, the scalar potential in a vacuum at any point a is
where k is the scale factor, r1 r2 r are distances apart from the positive, the negative, and the dipole center respectively, p is distance between two monopoles, I is the current injected or removed, and d is the dipole moment. By combining (4) and (5), the equipotential lines and electric-field intensity in vicinity is illustrated in Figure 6.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/727394/fnins-15-727394-HTML-r1/image_m/fnins-15-727394-g006.jpg)
Figure 6. The simplified extracellular electric field and its equipotential lines generated by dipole.
Forward Computation of Surface Electroencephalogram
From the extracellular field to the scalp, charges pass through the gray matter, the white matter, the cerebrospinal fluid (CSF), the skull, and the scalp. Tissues on this pathway have their unique conductivities, resulting in a smearing effect on the forward potential computation (Wolters et al., 2006). Among them, the conductivities for gray matter, scalp, and CSF are isotropic (Baumann et al., 1997); While the white matter and the skull have anisotropic conductivity (Rush, 1900; Nicholson, 1965).
By multiplying the electric field E and the conductivity, the current density is introduced to describe the total flow of charge per time over a cross-section of area. According to Ohm’s law, the current density is
where vector field J in A/m2 indicates the current density perpendicular to the cross-sectional area of the conductor, and κ is the electrolytic conductivity with units S/m. For isotropic conductivities κ is a scalar; While for anisotropic conductivities, κ ∈ ℝ3×3 is a position dependent conductivity tensor as
By applying the divergence operator to the current density, according to Poisson’s equation, relationships can be defined
where the ∇ ⋅ J is often called the current source density and symbolized with Im in Enderle (2012). To compute the potential field on different tissues, combing equation (4), (6), and (8), a general formed Poisson’s differential equation can be obtained
For a dipole, the current source density at any point can be written as a delta function with its distance from the positive and the negative.
where r1 r2 are the same as in (5).
In the Cartesian coordinate system, setting the positive monopole at position (x1, y1, z1) and the negative at (x2, y2, z2), at any point (x, y, z) (10) becomes
for anisotropic conductivities. For isotropic conductivities with a scalar κ, (10) can be written as
GAN Converter: From Simulated Signal to Practical Electroencephalograms
To accurately simulate the scalp EEGs, one’s unique tissue conductivity, geometry, and the location of the EEG generator are crucial in the forward computation. Such data have to be measured with specific equipment (e.g., Magnetic Resonance Imaging). But in practice, for BCI, it is impossible to get this biological information from every user. Without detailed personal data, however, the simulated scalp EEGs can only stop at theory without being used in realistic applications. To mimic one’s real EEGs in the lack of individual bio-information, the generative adversarial network (GAN) was adopted to convert the simulated signal into its practical version.
Conditional Generative Adversarial Network for Practical Electroencephalograms Imitation
Generative adversarial network was firstly proposed by Goodfellow et al. (2014), which provides the ability to counterfeit images that are statistically indistinguishable from real ones. A GAN consists of two networks: a Generator and a Discriminator. The Generator is trained to produce counterfeits that can deceive the Discriminator. As an adversary, the Discriminator aims at distinguishing the true from the fake. Inspired by GAN’s genius in image forgery, it was adopted in this work as a powerful converter to transform the simulated scalp EEGs from the theoretical-only to the practical.
To transform the simulated EEGs hence applying in the training of BCI classifier, conditional GAN was adopted to output signals with labels as the specified intention. For traditional GANs, the input of the Generator is a random vector, and several up-sampling computations are used to expand pixels. However, as a converter in this research, the simulated scalp EEGs were fed as input instead, without extra up-sampling.
According to our previous experimental setup in the BCI serial studies (Lu et al., 2018a, 2020; Zhang et al., 2021b), to realize the real-time decoding, continuous collected EEGs (4 s) were sliced into short-windowed (100 ms) to form the training-set. To reproduce the dataset forming procedure and augment the real EEGs samples, the output of the Generator was sliced to window-length size before entering the Discriminator. The conditional GAN in this work is illustrated in Figure 7.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/727394/fnins-15-727394-HTML-r1/image_m/fnins-15-727394-g007.jpg)
Figure 7. The conditional GAN for simulated EEGs conversion, where the number of trials, timestamps, channels, and labels, and the window length is in accordance with the experimental setup (latterly stated in RESULTS).
Generator and Discriminator
In general cases, the input for the Generator is one-dimensional noise. Thus, one-hot labels for conditional GAN can be concatenated just after the noise. But in this work, with inputs being multi-channeled simulated scalp EEGs, one-hot labels were fed and transformed as an extra channel concatenated with the simulated signal. The transform of the label input is demonstrated in Figure 8, of which the Generator and the Discriminator shared the same method but a different number of dense nodes.
The detailed design of the Generator and the Discriminator is listed in Tables 1, 2. In both networks, the selection of kernel size followed the principle from a rough glance to a fine adjustment. The optimizer and the cost function of GAN are shown in Table 3.
In the training process of GAN, the output of the Generator was unclear at the beginning. Then the awkward outputs mixed with the real collected EEGs were fed into the Discriminator to train a pair of discerning eyes. Subsequently, the Discriminator was frozen, and the Generator was trained to confuse the Discriminator according to the loss of the GAN. After iterations, with a comparable capability of the two networks, the Generator was able to successfully convert the simulated EEGs to the convincingly share the same characteristics as the real collected data.
Experimental Data
In the experiment, the latest EEG-based control paradigm assisted by microexpressions (ME-BCI) (Zhang et al., 2021b) proposed by our team was adopted, to verify the feasibility of scalp EEG modeling and its application in BCI training.
Experimental Electroencephalograms Acquisition
The commercial wireless EEG acquisition system (Neuracle Technology Co., Ltd.) with 30 EEG channels and 1000 Hz sampling rate is adopted in this study, as illustrated in Figure 9. The electrode placement is demonstrated also in Figure 9, in accordance with the international 10-20 location system, in which AFz is the ground and Cpz is the reference by default. Eight healthy subjects (25-38 years of age) participated in the experiment (World Medical Association, 2013).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/727394/fnins-15-727394-HTML-r1/image_m/fnins-15-727394-g009.jpg)
Figure 9. The Neuracle wireless EEG acquisition system and its electrode placement, in which AFz is the ground and Cpz is the reference by default.
Following our latest progress on ME-BCI (Zhang et al., 2021b), real EEGs were collected. Four microexpressions were selected: micro raise-brow (mRB), micro furrow-brow (mFB), micro left-smirk (mLS), and micro right-smirk (mRS), as illustrated in Figure 10. Each microexpression was conducted for 4 sessions (including 6 trials) per subject. In each trial, 3 s countdown, 4 s microexpression, and 2 s rest were contained in turn, as demonstrated in Figure 11. During data collection, subjects were asked to sit quietly and avoid extra body movements.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/727394/fnins-15-727394-HTML-r1/image_m/fnins-15-727394-g010.jpg)
Figure 10. Illustration of four selected microexpressions. (A) Micro raise-brow. (B) Micro furrow-brow. (C) Micro left-smirk. (D) Micro right-smirk.
In the ME-BCI related serial works (Lu et al., 2018a, 2020; Zhang et al., 2021b), for the real-time requirement, each result was decoded with 100 ms EEG input. Thus, in this work, the EEG sequences were still sliced into pieces with short-window-length. Real collected EEG data were divided into two parts: 1st session as a template for simulated EEGs conversion, 2nd ∼ 4th sessions as testing-set for feasibility verification.
The Head Models and the Dipole Position
Facing the lack of one’s unique tissue conductivity and geometry, the three-shell standard mesh of the human head in FieldTrip (Oostenveld et al., 2010) was adopted as a template. In this standard mesh, the head structure is simplified into several surfaces, including brain, skull, and scalp; On the scalp are the EEG electrodes. Figure 12 shows the surfaces of the standard head model, and the electrodes in accordance with the Neuracle Device.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/727394/fnins-15-727394-HTML-r1/image_m/fnins-15-727394-g012.jpg)
Figure 12. Triangulated surfaces of the standard head model, and the electrodes location of the Neuracle device. (A) Brain. (B) Skull. (C) Scalp. (D) Electrodes location according to the Neuracle device.
As for the location of the EEG generator, since the precise positioning cannot be achieved without the extra professional equipment, the dipole was set in accordance with the mechanism. In our serial research on ME-BCI (Zhang et al., 2016, 2021b; Li et al., 2018b; Lu et al., 2018a, b, 2020), data-driven brain connectivity analysis demonstrated the main involvement of the motor cortex (Lu et al., 2018a; Zhang et al., 2021b), which conformed to the contralateral control facts. Meanwhile, evidence showed the frontal lobe and limbic system also participate in facial-expression processing (Price and Drevets, 2010; Li et al., 2018b; Lu et al., 2018b). Thereby, the dipole was set orthogonalization to the surface in the frontal lobe/motor cortex, separately for different microexpressions (Figure 13).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/727394/fnins-15-727394-HTML-r1/image_m/fnins-15-727394-g013.jpg)
Figure 13. Dipole position of scalp EEGs simulation in (A) frontal lobe, (B) left motor cortex, and (C) right motor cortex. For mRB and mFB, dipole was set in the frontal lobe; for mLS was the right motor cortex, and for mRS was the left motor cortex.
Results
Simulation of Multi-Channeled Scalp Electroencephalograms
Postsynaptic Membrane Potential
To simulate postsynaptic membrane potentials in three representative EEG bands as low (4-12 Hz), medium (12-30 Hz) and high (>30 Hz), different parameter values were taken in tuning the peaks. Part of parameters in (1) and (2) are common for each band whose values are determined according to anatomical facts (Table 4; Freeman, 1987; Jansen and Rit, 1995).
Using the 4th order Runge-Kutta algorithm, differential equations in (1) can obtain accurate numerical solutions. Under certain parameter combinations, the peak can be fine-tuned by adjusting only G (the average gain for fast inhibitory interneurons) (Zavaglia et al., 2006). Setting the exogenous input P(t) as a random noise (mean value = 60, variance = 100), Table 5 lists the system parameters to generate peaks approximately in three bands. Figure 14 illustrates the spectrum of each population.
To mimic the subject’s EEG, population weights are calculated via one’s real EEGs (1st session of experimental data). The uniformly scaled energies within the band’s range in Real EEGs were set as the population weights, as listed in Table 6. The mimic postsynaptic membrane potentials are shown in Figure 15.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/727394/fnins-15-727394-HTML-r1/image_m/fnins-15-727394-g015.jpg)
Figure 15. Simulated postsynaptic membrane potentials of subject S1 under micro-raise-brow, micro-furrow-brow, micro-left-smirk, and micro-right-smirk.
Multi-Channelled Scalp Electroencephalograms
With the standard mesh of the human head, the forward computation was calculated with the boundary element method (BEM) for its low computational needs (Fuchs et al., 2002). According to the dipole setup, the simulated multi-channeled scalp EEGs were demonstrated in Figure 16.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/727394/fnins-15-727394-HTML-r1/image_m/fnins-15-727394-g016.jpg)
Figure 16. Simulated micro right-smirk multi-channeled scalp EEGs (in aesthetics consideration, only parts of channels are shown).
However, due to the blankness of the individual bio-information, the simulated multi-channeled scalp EEGs lacked realism in the temporal domain. But literature showed its capability in mimicking the spectral characteristics (Li et al., 2018a). Table 7 lists the Pearson correlation coefficient between the simulated scalp EEGs and the real collected 1st session in spectrum distribution. Figure 17 demonstrated the spectrum comparison between simulated scalp EEGs and the real collected one.
Unlike visual evoked potentials, the frequency spectrum distribution of the ME-BCI paradigm showed no typical features (Zhang et al., 2016). As shown in Table 7, without a fine tuning in NMM parameters, dipole positions, and tissue configurations, the simulation ability of computation is limited for spectra without significant characteristics. The value in Table 7 indicated that the correlation existed, but was not high.
Conversion and Application of Simulated Scalp Electroencephalograms
Given the lack of biological information, the theoretical simulation of multi-channeled scalp EEGs can only reproduce some general phenomenon. To complete the conversion of simulated EEGs to its useful practical version in BCI, the simulated signals were fed into the conditional GAN. The training losses of the Discriminator and the GAN are illustrated in Figure 18.
With the random exogenous input P(t), the population weights, the dipole position, and the conditional label, endless practical simulated EEGs can be generated through the conversion. To verify the feasibility and capability of this method, the simulated EEGs mixed with the real collected 1st session were used to train a ME-BCI classifier, and the rest of the 2nd to 4th sessions were used to test its performance. The division of datasets is shown in Figure 19.
All EEGs, both the real and the simulated one, were detrended and filtered into [2 Hz, 55 Hz] with 4th-order Butterworth bandpass-filter. To form an intuitive comparison, same as previous research (Lu et al., 2018a, 2020), EEGs were sliced into short-windowed (100 ms) segments, and the same feature extraction method Common Spatial Pattern (CSP) was adopted to emphasize the performance between “training with only real EEGs” and “training with fake augmented EEGs.” The spatial filter of CSP was calculated via the training-set. The ‘one versus rest’ strategy was adopted in the CSP. Various commonly used classifiers were used to calculate the validation accuracy and test accuracy. The validation and test performances under ‘training with fake augmented EEGs’ are listed in Table 8.
As shown in Table 8, classifiers such as K-nearest neighbor (KNN) and some ensembles (i.e., bagged trees, and subspace KNN) achieved higher validation accuracies. As for test accuracies, due to the individual differences, the personal best varied from 73.7 to 98.1%. The mean test accuracy among classifiers is 84.93% ± 1.25. The overall performance in Table 8 indicates the feasibility of applying the simulated EEGs to the BCI classifier training, through a GAN conversion. For a detailed comparison, accuracies under the training-set with only the 1st session real EEGs were also estimated. The validation and test performances under ‘training with only the real 1st session EEGs’ are shown in Table 9.
Comparing Table 9 with Table 8, most accuracies have been improved after the data augmentation with simulated EEGs, both the validation and the test. The two-sample T-test shows that after adding the simulated EEGs, the overall performance (both validation and test) has a significant improvement (P = 0.04 < 0.05), compared with insufficient real data. By marking the improvement in test accuracy as positive, and the decline as negative, the overall performance under augmenting with fake EEGs is 2.17% ± 4.23. Among all subjects, test accuracy of S5 achieved the largest improvement as 12.60% ± 1.81 among 6 classifiers. The comparison between Tables 8, 9 confirms the feasibility and capability of scalp EEGs simulation as a new data augmentation method in BCI.
Figure 20 demonstrates the confusion matrix of six classifiers under training with “fake EEGs augmented” vs. “only real EEGs”. The blue high-lighted blocks in Figure 20B indicate that the performance has been improved after being augmented by the converted simulated EEGs. Here, the improvements mean a raise in true positive, or a decline in false positive. In all six classifiers, the distribution between “fake EEGs augmented” and “only real EEGs” basically remained the same, and showed no obvious imbalance among categories. The confusion result indicates that this training-set augmentation method with simulated EEGs shows no bias and maintains the same level of performance among categories.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/727394/fnins-15-727394-HTML-r1/image_m/fnins-15-727394-g020.jpg)
Figure 20. The confusion matrix of test accuracy under (A) training set with “fake EEGs augmented” versus (B) training set “only real EEGs.” The blue highlight represents that the performance has been improved with the simulation contained.
Discussion
In this work, a mathematical model of scalp EEG is presented, and a conditional GAN is put forward to convert the simulated scalp EEGs from theoretical-only to practical, hence applying in the training of BCI classifiers. Through this study, we have taken a step forward from the theoretical modeling of EEG, and initially established the connection between simulated EEGs and its practical applications in BCI. The pros and cons existing in this work are discussed as followed.
The Significance and the Performance
The Tuning of Neural Mass Model Parameters, and the Structure of Coupled Neural Mass Model
Neural mass model is represented as a group of differential equations with several parameters related to physiological mechanisms. Apart from the anatomy-related common parameters (C, e0, s0, and r), for a single population, Jasen-Rit NMM with two feedback interneurons has four adjustable parameters (Jansen and Rit, 1995), Wendling NMM with three feedbacks has six (Wendling et al., 2002), and the more in Ursiono NMM (Ursino et al., 2010). For a multi-population model, the adjustable parameters grow exponentially, and the coupling structure also affects the output (Wendling et al., 2000). However, the peak of NMM output shows a complex non-linear relationship to the coupling of parameters. To strictly reproduce the event within collected EEGs, large efforts have to be paid or an extra fitting algorithm has to be added to tune the NMM parameters. Another difficulty in parameter tuning is the chaos of NMM (Huang et al., 2011). Even with changeless parameters, the output of NMM remains fluctuated and shows instability to different initial values and random exogenous input.
In this work, with emphasis put on the link between EEG modeling and its application in BCI, the careful selection of the NMM parameters for each individual seems unaffordable and less practical. With a GAN converter added, the strict re-produce of EEG characteristics is no longer obsessed, thereby largely ease the NMM parameter tuning step.
Since single NMM can only generate a unimodal spectrum, for cortical responses, multi-coupled NMM was adopted in the modeling of postsynaptic membrane potential. In Jansen and Rit, 1995; Zhang et al., 2016; Li et al., 2018a, a double-column model was used to model the generation of neural activity potential. In this work, the triple-population coupled NMM structure was adopted; Compared to the double-column model, the simulated output of triple-population showed more details within the EEG band, and can reproduce the fast activity within EEG. As for the studies that involve triple coupled model, such as Zavaglia et al. (2006) and Zhang et al. (2021a), different detailed population weights were set to adjust the coupling proportion within multi-populations. Such population weights were calculated via collected personal EEG data, thus improving the realism and enhancing the personal characteristics of the cortical response.
Overcome the Lack of Personal Brain Tissue Data
From the postsynaptic membrane potential to the scalp EEGs, personal brain tissue data matters a lot. These data include the geometric shape of the brain surface, the tissue conductivity, and the position of the cortical sources. Even by using the BEM method, further combing a distributed source model (Dale et al., 2000), the realistic cortical shape is still crucial in the forward computation. This information by far has to be determined from Magnetic Resonance Imaging (MRI). An MRI system costs millions of dollars. For non-medical department students or even some colleges, there is little chance of accessing the MRI. Even if it can be detected, this personal and costly modeling method is not practical in the field of BCI. After all, in application, it is impossible for every user to undergo MRI detection.
To overcome these obstacles we faced during research, the GAN converter was used in this work to make up for the default in personal data. Focusing on the BCI field, not only simulation but application, with the method in this work, large numbers of applicable scalp EEGs can be generated. Meanwhile, it should be made feasible to apply the simulated EEGs to the training of BCI classifier.
Build the Connection Between Electroencephalogram Simulation and Application in BCI
Electroencephalogram modelling in BCI area is related to specific mental tasks or BCI paradigms, but it by far has not formed guidance for the paradigm design or built any connection with its application in BCI; In other words, the EEG modelling in BCI derives from BCI but exists independently of BCI, lacking integrity. Jansen and Rit (1995) simulated the generation of visual evoked potential in EEG, which outputted signals with approximate phase response and, in detail, exported the influence of NMM parameters. But, strictly, the output ended up as the postsynaptic membrane potential, and no more properties other than phase were discussed. Zavaglia et al. (2006) reproduced the power spectral density and the temporal changes of cortical EEG during finger movement; However, the use of large fitting methods (in optimizing NMM parameters) and the MRI measurement made it difficult to establish the connection with the application in BCI. Delphine’s work (Cosandier-Rimele et al., 2010) established the computational modeling of epileptic activity, proposed a solid study in scalp EEG simulation; Although this work was unrelated to BCI paradigms, it indicated that the scalp EEG simulation are very sensitive to the geometry and electrical properties of the different head tissues (mainly, the brain, skull, and scalp), discussing the gap between the simulated-EEG and the actual-EEG and its difficulty of application. Recently, Li et al. (2018a) simulated the EEG frequency response of the scene graph steady-state visual evoked potentials; Similarly, the output ended up at the postsynaptic membrane potential, discussed no other properties besides the visual response frequency, and formed no guidance for paradigm design or application. In the latest, Hanzhe (Zhang et al., 2021a) built the mathematical model of EEG for lower limb voluntary movement intention based on NMM; However, the output did not propagate to the scalp EEG, and made no comparison with the real collected signal.
The works above did make large amounts of progress in EEG modeling related to different BCI paradigms or mental tasks, however a large number of them stopped at the postsynaptic membrane potential, and focused mainly on the reproduction of little EEG response characteristics. Compared to them, this work focused on establishing the connection between the EEG modeling and its application, by proposing the scalp EEG model and applying the simulated scalp EEG into the training of BCI classifiers, so as to add to the modeling study; The 2.17% ± 4.23 improvement of test accuracy demonstrated the feasibility. Meanwhile, in order to establish such a connection, compared to those precise modeling methods that required individual tissue data, this work adopted the standard head mesh during the simulation and proposed an additional converter to overcome such a data default problem.
The Limitations and Further Work
The Overfit Problem in Classification
In the transformer from the awkward simulated EEGs to the applicable version, a conditional GAN was used. During the learning process of GAN, the signals in the real collected 1st session were used as the templates. The Discriminator would be deceived if the input showed similar features to the real templates. According to this logic, the Generator would finally output converted signals with the same features based on the judgment condition of the Discriminator. Therefore, the converted simulated signals would show distinct characteristics similar to templates.
It is widely recognized that surface EEG is a non-linear and unstable signal with a low signal-to-noise ratio, and is easily affected by one’s thinking or emotions. This attribute makes the EEG tend to show different feature distribution with time-shifting. In other words, the feature space or data distribution extracted under partial time has poor generalization. Thus, the BCI classifier would show large overfitting if trained with samples lacking diversity.
However, the simulated EEGs converted from GAN would have similar characteristics as the 1st session, so does the feature space. Since the focal point in this work is not emphasized on the EEG decoding, the feature extraction algorithm and the machine learning method are not discussed. Comparing Table 8 to Table 9, the accuracy in Table 8 shows a larger descent between the validation and the test than in Table 9. It can be found that, with the simulated EEGs added to the training, the test performance has been improved, but the over-fitting problem has been exacerbated at the same time. In future work, more attention should be paid to the diversity of the converted EEGs. The attentional mechanism can be involved (Lashgari et al., 2021), to focus more on the main characteristics so that to enable more randomness in other secondary characteristics.
The Lower Adaptability Among Classifiers
In Table 8, taking subject S4 as an example, among various classifiers, the validation accuracies range from 27.9 to 98.9%. For all subjects, classifiers such as “coarse tree,” “medium Gaussian SVM,” “coarse Gaussian SVM,” and “subspace discriminant” achieved low validation accuracies. While several KNNs performed better. The same result did not appear in the case of training with all real collected EEGs (Lu et al., 2020), although there are difference between real datasets.
As discussed above, the simulated EEGs converted by GAN shall appear the characteristic that accords with the discriminator’s judgment basis. According to Table 8’s result, such characteristics are more suitable with the computation logic of KNN. Further, it indicates that there are flaws in the design of the Discriminator in this article, perhaps the network structure or the cost function. Additionally, the symmetrical network structure tends to be adopted in the design of GAN, thus the same defects also exist in the Generator. In the next studies, more attention will be paid to the design of converter GAN, to ensure a better adaptability of simulated EEGs’ application among different classifiers; Additionally, more kinds of decoding algorithms will also be adopted to estimate the feasibility, such as methods with hybrid-domain features (Ieracitano et al., 2021), end-to-end deep learning methods, transfer learning methods and so on.
Conclusion
This work established the connection from the EEG simulation to its application in BCI, and proved its feasibility in improving the test accuracy compared to insufficient real data. In this paper, a mathematical model of surface EEG was presented, and a GAN converter was proposed to transfer the theoretical simulated EEGs to its applicable version in BCI training.
In accordance with the physics of surface EEG, the mathematical computation is firstly started with a triple-population coupled NMM. Then the dipole and the forward computation were followed to model the propagation. To overcome the lack of individual biological data and build a bridge between the simulation and the application, a GAN converter was established to transfer the simulated scalp signal to its applicable version. In application, the converted simulated EEGs were used in the training of BCI classifiers. With the simulated EEGs added into the training-set, compared with only insufficient real collected data, the overall performance improved significantly (P = 0.04 < 0.05), and the test performance showed an overall 2.17% ± 4.23 increase. Among all subjects, the largest increase is 12.60% ± 1.81. Through this work, we hope to provide a novel feasible solution for the application of surface EEG modeling in BCI.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics Statement
The studies involving human participants were reviewed and approved by the institutional review board of Xi’an Jiaotong University (No. 20211452), and all experiments were conducted in accordance with the Declaration of Helsinki. The patients/participants provided their written informed consent to participate in this study. Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.
Author Contributions
XZ supervised this work and revised the manuscript. ZL proposed and did the research, and wrote the manuscript. TZ, HL, and YW organized and carried out the experiments. QT revised the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This work was supported in part by the National Key Research and Development Program of China (No. 2017YFB1300303), and partially supported by grants from the Science and Technology Supporting Xinjiang Project (No. 2020E0259).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Baillet, S., Mosher, J. C., and Leahy, R. M. (2001). Electromagnetic brain mapping. IEEE Signal Process. Mag. 18, 14–30. doi: 10.1109/79.962275
Baumann, S. B., Wozny, D. R., Kelly, S. K., and Meno, F. M. (1997). The electrical conductivity of human cerebrospinal fluid at body temperature. IEEE Trans. Biomed. Eng. 44, 220–223. doi: 10.1109/10.554770
Berg, P., and Scherg, M. (1994). A fase method for forward computation of multiple-shell spherical head models. Electroencephal. Clin. Neurophysiol. 90, 58–64. doi: 10.1016/0013-4694(94)90113-9
Berger, H. (1929). Über das elektrenkephalogramm des menschen. Arch. Psychiatr. Nervenkr. 87, 527–570.
Chen, X. G., Wang, Y. J., Nakanishi, M., Gao, X. R., Jung, T. P., and Gao, S. K. (2015). High-speed spelling with a noninvasive brain-computer interface. Proc. Natl. Acad. Sci. U.S.A. 112, E6058–E6067. doi: 10.1073/pnas.1508080112
Cosandier-Rimele, D., Merlet, I., Bartolomei, F., Badier, J.-M., and Wendling, F. (2010). computational modeling of epileptic activity: from cortical sources to EEG signals. J. Clin. Neurophysiol. 27, 465–470. doi: 10.1097/WNP.0b013e3182005dcd
Dale, A. M., Liu, A. K., Fischl, B. R., Buckner, R. L., Belliveau, J. W., Lewine, J. D., et al. (2000). Dynamic statistical parametric mapping: combining fMRI and MEG for high-resolution imaging of cortical activity. Neuron 26, 55–67. doi: 10.1016/S0896-6273(00)81138-1
Enderle, J. D. (2012). “Chapter 12 - Bioelectric phenomena,” in Introduction to Biomedical Engineering, 2nd Edn, eds J. D. Enderle and J. D. Bronzino (Boston, MA: Academic Press), 747–815.
Fan, J. H., Sun, C., Chen, C., Jiang, X., Liu, X., Zhao, X., et al. (2020). EEG data augmentation: towards class imbalance problem in sleep staging tasks. J. Neural Eng. 17:056017. doi: 10.1088/1741-2552/abb5be
Freeman, W. J. (1973). “A model of the olfactory system,” in Neural Modeling, Brain Information Service Res. Report 1, eds M. A. B. Brazier, D. O. Walter, and D. Schneider (Los Angeles, CA: University of California).
Freeman, W. J. (1972). Linear analysis of dynamics of neural masses. Annu. Rev. Biophys. Bioeng. 1, 225–256. doi: 10.1146/annurev.bb.01.060172.001301
Freeman, W. J. (1975). “CHAPTER 7–signal processing by neural mass actions,” in Mass Action in the Nervous System, ed. W. J. Freeman (San Diego, CA: Academic Press), 402–461.
Freeman, W. J. (1987). Simulation of chaotic EEG patterns with a dynamic model of the olfactory system. Biol. Cybern. 56, 139–150. doi: 10.1007/BF00317988
Fuchs, M., Kastner, J., Wagner, M., Hawes, S., and Ebersole, J. S. (2002). A standardized boundary element method volume conductor model. Clin. Neurophysiol. 113, 702–712. doi: 10.1016/s1388-2457(02)00030-5
Goodfellow, I. J., Pouget-Abadie, J., Mirza, M., Xu, B., Warde-Farley, D., Ozair, S., et al. (2014). Generative adversarial networks. arXiv [Preprint]. arXiv:1406.2661.
Hallez, H., Vanrumste, B., Grech, R., Muscat, J., De Clercq, W., Vergult, A., et al. (2007). Review on solving the forward problem in EEG source analysis. J. Neuroeng. Rehabil. 4:46. doi: 10.1186/1743-0003-4-46
He, B., Yao, D., and Lian, J. (2002). High-resolution EEG: on the cortical equivalent dipole layer imaging. Clin. Neurophysiol. 113, 227–235. doi: 10.1016/S1388-2457(01)00740-4
Herrmann, C. S. (2001). Human EEG responses to 1-100 Hz flicker: resonance phenomena in visual cortex and their potential correlation to cognitive phenomena. Exp. Brain Res. 137, 346–353. doi: 10.1007/s002210100682
Huang, G., Zhang, D., Meng, J., and Zhu, X. (2011). Interactions between two neural populations: a mechanism of chaos and oscillation in neural mass model. Neurocomputing 74, 1026–1034. doi: 10.1016/j.neucom.2010.11.019
Ieracitano, C., Morabito, F. C., Hussain, A., and Mammone, N. (2021). A hybrid-domain deep learning-based BCI for discriminating hand motion planning from EEG sources. Int. J. Neural syst. 31:2150038. doi: 10.1142/s0129065721500386
Jansen, B. H., and Rit, V. G. (1995). Electroencephalogram and visual evoked potential generation in a mathematical model of coupled cortical columns. Biol. Cybern. 73, 357–366. doi: 10.1007/BF00199471
Jansen, B. H., Zouridakis, G., and Brandt, M. E. (1993). A neurophysiologically-based mathematical-model of flash visual evoked-potentials. Biol. Cybern. 68, 275–283. doi: 10.1007/bf00224863
Koelstra, S., Mühl, C., Soleymani, M., Lee, J., Yazdani, A., Ebrahimi, T., et al. (2012). DEAP: a database for emotion analysis using physiological signals. IEEE Trans. Affect. Comput. 3, 18–31. doi: 10.1109/t-affc.2011.15
Koochali, A., Schichtel, P., Dengel, A., and Ahmed, S. (2019). Probabilistic forecasting of sensory data with generative adversarial networks–ForGAN. IEEE Access 7, 63868–63880. doi: 10.1109/ACCESS.2019.2915544
Krishna, G., Co, T., Carnahan, M., Yan, H., and Tewfik, A. H. (2021). “Generating EEG features from acoustic features,” in Proceeding of the 2020 28th European Signal Processing Conference, Amsterdam, 1100–1104.
Lashgari, E., Ott, J., Connelly, A., Baldi, P., and Maoz, U. (2021). An end-to-end CNN with attentional mechanism applied to raw EEG in a BCI classification task. J. Neural Eng. 18:0460e3. doi: 10.1088/1741-2552/ac1ade
Li, R., Zhang, X. D., Li, H. Z., Zhang, L. M., Lu, Z. F., and Chen, J. C. (2018a). An approach for brain-controlled prostheses based on scene graph steady-state visual evoked potentials. Brain Res. 1692, 142–153. doi: 10.1016/j.brainres.2018.05.018
Li, R., Zhang, X., Lu, Z., Liu, C., Li, H., Sheng, W., et al. (2018b). An approach for brain-controlled prostheses based on a facial expression paradigm. Front. Neurosci. 12:943. doi: 10.3389/fnins.2018.00943
Lopes da Silva, F. H., van Rotterdam, A., Barts, P., van Heusden, E., and Burr, W. (1976). “Models of neuronal populations: the basic mechanisms of rhythmicity,” in Progress in Brain Research, eds M. A. Corner and D. F. Swaab (Amsterdam: Elsevier), 281–308.
Lopesdas, Fh, Hoeks, A., Smits, H., and Zetterberg, L. H. (1974). Model of brain rhythmic activity: alpha-rhythm of thalamus. Kybernetik 15, 27–37. doi: 10.1007/bf00270757
Lu, Z., Zhang, X., Li, H., Li, R., and Chen, J. (2018a). “A real-time brain control method based on facial expression for prosthesis operation,” in Proceedings of the 2018 IEEE International Conference on Robotics and Biomimetics, Kuala Lumpur.
Lu, Z., Zhang, X., Li, H., Zhang, T., and Zhang, Y. (2020). “A Semi-Asynchronous Real-Time Facial Expression Assisted Brain Control Method: an Extension,” in Proceedings of the 2020 The 10th IEEE International Conference on CYBER Technology in Automation, Control, and Intelligent Systems, Xi’an.
Lu, Z., Zhang, X., Li, R., and Guo, J. (2018b). A brain control method for prosthesises based on facial expression. China Mech. Eng. 29, 1454–1459; 1474.
Luo, Y., Zhu, L. Z., Wan, Z. Y., and Lu, B. L. (2020). Data augmentation for enhancing EEG-based emotion recognition with deep generative models. J. Neural Eng. 17:056021. doi: 10.1088/1741-2552/abb580
Muller-Putz, G. R., Schwarz, A., Pereira, J., and Ofner, P. (2016). “From classic motor imagery to complex movement intention decoding: the noninvasive Graz-BCI approach,” in Brain-Computer Interfaces: Lab Experiments To Real-World Applications, ed. D. Coyle (Amsterdam: Elsevier), 39–70.
Nicholson, P. W. (1965). Specific impedance of cerebral white matter. Exp. Neurol. 13, 386–401. doi: 10.1016/0014-4886(65)90126-3
Oostenveld, R., Fries, P., Maris, E., and Schoffelen, J.-M. (2010). FieldTrip: open source software for advanced analysis of MEG, EEG, and invasive electrophysiological data. Comput. Intell. Neurosci. 2011:156869. doi: 10.1155/2011/156869
Pfurtscheller, G., and Neuper, C. (2001). Motor imagery and direct brain-computer communication. Proc. IEEE 89, 1123–1134. doi: 10.1109/5.939829
Plonsey, R., and Heppner, D. B. (1967). Considerations of quasi-stationarity in electrophysiological systems. Bull. Math. Biophys. 29, 657–664. doi: 10.1007/BF02476917
Polich, J. (2007). Updating p300: an integrative theory of P3a and P3b. Clin. Neurophysiol. 118, 2128–2148.
Price, J. L., and Drevets, W. C. (2010). Neurocircuitry of mood disorders. Neuropsychopharmacology 35, 192–216. doi: 10.1038/npp.2009.104
Rush, S. (1900). Current distribution in the brain from surface electrodes. Anesth. Analg. 47, 717–723.
Sanz-Leon, P., Knock, S. A., Spiegler, A., and Jirsa, V. K. (2015). Mathematical framework for large-scale brain network modeling in The Virtual Brain. Neuroimage 111, 385–430. doi: 10.1016/j.neuroimage.2015.01.002
Schaul, N. (1998). The fundamental neural mechanisms of electroencephalography. Electroencephalogr. Clin. Neurophysiol. 106, 101–107. doi: 10.1016/S0013-4694(97)00111-9
Ursino, M., Cona, F., and Zavaglia, M. (2010). The generation of rhythms within a cortical region: analysis of a neural mass model. Neuroimage 52, 1080–1094. doi: 10.1016/j.neuroimage.2009.12.084
van Rotterdam, A., Lopes da Silva, F. H., van den Ende, J., Viergever, M. A., and Hermans, A. J. (1982). A model of the spatial-temporal characteristics of the alpha rhythm. Bull. Math. Biol. 44, 283–305. doi: 10.1016/S0092-8240(82)80070-0
Wendling, F., Bartolomei, F., Bellanger, J.-J., and Chauvel, P. (2002). Epileptic fast activity can be explained by a model of impaired GABAergic dendritic inhibition. Eur. J. Neurosci. 15, 1499–1508. doi: 10.1046/j.1460-9568.2002.01985.x
Wendling, F., Bellanger, J. J., Bartolomei, F., and Chauvel, P. (2000). Relevance of nonlinear lumped-parameter models in the analysis of depth-EEG epileptic signals. Biol. Cybern. 83, 367–378. doi: 10.1007/s004220000160
Wilson, H. R., and Cowan, J. D. (1972). Excitatory and inhibitory interactions in localized populations of model neurons. Biophys. J. 12, 1–24. doi: 10.1016/s0006-3495(72)86068-5
Wilson, H. R., and Cowan, J. D. (1973). A mathematical theory of functional dynamics of cortical and thalamic nervous-tissue. Kybernetik 13, 55–80. doi: 10.1007/bf00288786
Wolpaw, J. R., Birbaumer, N., McFarland, D. J., Pfurtscheller, G., and Vaughan, T. M. (2002). Brain-computer interfaces for communication and control. Clin. Neurophysiol. 113, 767–791.
Wolters, C. H., Anwander, A., Tricoche, X., Weinstein, D., Koch, M. A., and MacLeod, R. S. (2006). Influence of tissue conductivity anisotropy on EEG/MEG field and return current computation in a realistic head model: a simulation and visualization study using high-resolution finite element modeling. Neuroimage 30, 813–826. doi: 10.1016/j.neuroimage.2005.10.014
World Medical Association (2013). World Medical Association declaration of Helsinki: ethical principles for medical research involving human subjects. JAMA 310, 2191–2194.
Zavaglia, M., Astolfi, L., Babiloni, F., and Ursino, M. (2006). A neural mass model for the simulation of cortical activity estimated from high resolution EEG during cognitive or motor tasks. J. Neurosci. Methods 157, 317–329. doi: 10.1016/j.jneumeth.2006.04.022
Zhang, X., Guo, J., Li, R., and Lu, Z. (2016). A simulation model and pattern recognition method of electroencephalogram driven by expression. J. Xian Jiaotong Univ. 50, 1–8.
Zhang, X., Lu, Z., Li, H., Zhang, T., Gu, L., and Wang, Y. (2021b). microFE-BCI: an asynchronous practical EEG-based control paradigm assisted by micro-facial-expressions. TechRxiv [Preprint]. doi: 10.36227/techrxiv.14397593.v1
Keywords: BCI, EEG, GAN, modeling, simulation, classification
Citation: Zhang X, Lu Z, Zhang T, Li H, Wang Y and Tao Q (2021) Realizing the Application of EEG Modeling in BCI Classification: Based on a Conditional GAN Converter. Front. Neurosci. 15:727394. doi: 10.3389/fnins.2021.727394
Received: 18 June 2021; Accepted: 04 October 2021;
Published: 11 November 2021.
Edited by:
Mufti Mahmud, Nottingham Trent University, United KingdomReviewed by:
Cosimo Ieracitano, Mediterranea University of Reggio Calabria, ItalyAvinash Kumar Singh, University of Technology Sydney, Australia
Copyright © 2021 Zhang, Lu, Zhang, Li, Wang and Tao. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Zhufeng Lu, zhufeng.luk@qq.com