Electrophysiological brain imaging based on simulation-driven deep learning in the context of epilepsy.

Identifying the location, the spatial extent and the electrical activity of distributed brain sources in the context of epilepsy through ElectroEncephaloGraphy (EEG) recordings is a challenging task because of the highly ill-posed nature of the underlying Electrophysiological Source Imaging (ESI) problem. To guarantee a unique solution, most existing ESI methods pay more attention to solve this inverse problem by imposing physiological constraints. This paper proposes an efficient ESI approach based on simulation-driven deep learning. Epileptic High-resolution 256-channels scalp EEG (Hr-EEG) signals are simulated in a realistic manner to train the proposed patient-specific model. More particularly, a computational neural mass model developed in our team is used to generate the temporal dynamics of the activity of each dipole while the forward problem is solved using a patient-specific three-shell realistic head model and the boundary element method. A Temporal Convolutional Network (TCN) is considered in the proposed model to capture local spatial patterns. To enable the model to observe the EEG signals from different scale levels, the multi-scale strategy is leveraged to capture the overall features and fine-grain features by adjusting the convolutional kernel size. Then, the Long Short-Term Memory (LSTM) is used to extract temporal dependencies among the computed spatial features. The performance of the proposed method is evaluated through three different scenarios of realistic synthetic interictal Hr-EEG data as well as on real interictal Hr-EEG data acquired in three patients with drug-resistant partial epilepsy, during their presurgical evaluation. A performance comparison study is also conducted with two other deep learning-based methods and four classical ESI techniques. The proposed model achieved a Dipole Localization Error (DLE) of 1.39 and Normalized Hamming Distance (NHD) of 0.28 in the case of one patch with SNR of 10 dB. In the case of two uncorrelated patches with an SNR of 10 dB, obtained DLE and NHD were respectively 1.50 and 0.28. Even in the more challenging scenario of two correlated patches with an SNR of 10 dB, the proposed approach still achieved a DLE of 3.74 and an NHD of 0.43. The results obtained on simulated data demonstrate that the proposed method outperforms the existing methods for different signal-to-noise and source configurations. The good behavior of the proposed method is also confirmed on real interictal EEG data. The robustness with respect to noise makes it a promising and alternative tool to localize epileptic brain areas and to reconstruct their electrical activities from EEG signals.

Selective enhancement of low-gamma activity by tACS improves phonemic processing and reading accuracy in dyslexia.

The phonological deficit in dyslexia is associated with altered low-gamma oscillatory function in left auditory cortex, but a causal relationship between oscillatory function and phonemic processing has never been established. After confirming a deficit at 30 Hz with electroencephalography (EEG), we applied 20 minutes of transcranial alternating current stimulation (tACS) to transiently restore this activity in adults with dyslexia. The intervention significantly improved phonological processing and reading accuracy as measured immediately after tACS. The effect occurred selectively for a 30-Hz stimulation in the dyslexia group. Importantly, we observed that the focal intervention over the left auditory cortex also decreased 30-Hz activity in the right superior temporal cortex, resulting in reinstating a left dominance for the oscillatory response. These findings establish a causal role of neural oscillations in phonological processing and offer solid neurophysiological grounds for a potential correction of low-gamma anomalies and for alleviating the phonological deficit in dyslexia.

Modelling acute and lasting effects of tDCS on epileptic activity.

Transcranial Direct brain stimulation (tDCS) is commonly used in order to modulate cortical networks activity during physiological processes through the application of weak electrical fields with scalp electrodes. Cathodal stimulation has been shown to decrease brain excitability in the context of epilepsy, with variable success. However, the cellular mechanisms responsible for the acute and the long-lasting effect of tDCS remain elusive. Using a novel approach of computational modeling that combines detailed but functionally integrated neurons we built a physiologically-based thalamocortical column. This model comprises 10,000 individual neurons made of pyramidal cells, and 3 types of gamma-aminobutyric acid (GABA) -ergic cells (VIP, PV, and SST) respecting the anatomy, layers, projection, connectivity and neurites orientation. Simulating realistic electric fields in term of intensity, main results showed that 1) tDCS effects are best explained by modulation of the presynaptic probability of release 2) tDCS affects the dynamic of cortical network only if a sufficient number of neurons are modulated 3)VIP GABAergic interneurons of the superficial layer of the cortex are especially affected by tDCS 4) Long lasting effect depends on glutamatergic synaptic plasticity.

Identification of Interictal Epileptic Networks from Dense-EEG.

Epilepsy is a network disease. The epileptic network usually involves spatially distributed brain regions. In this context, noninvasive M/EEG source connectivity is an emerging technique to identify functional brain networks at cortical level from noninvasive recordings. In this paper, we analyze the effect of the two key factors involved in EEG source connectivity processing: (i) the algorithm used in the solution of the EEG inverse problem and (ii) the method used in the estimation of the functional connectivity. We evaluate four inverse solutions algorithms (dSPM, wMNE, sLORETA and cMEM) and four connectivity measures (r 2, h 2, PLV, and MI) on data simulated from a combined biophysical/physiological model to generate realistic interictal epileptic spikes reflected in scalp EEG. We use a new network-based similarity index to compare between the network identified by each of the inverse/connectivity combination and the original network generated in the model. The method will be also applied on real data recorded from one epileptic patient who underwent a full presurgical evaluation for drug-resistant focal epilepsy. In simulated data, results revealed that the selection of the inverse/connectivity combination has a significant impact on the identified networks. Results suggested that nonlinear methods (nonlinear correlation coefficient, phase synchronization and mutual information) for measuring the connectivity are more efficient than the linear one (the cross correlation coefficient). The wMNE inverse solution showed higher performance than dSPM, cMEM and sLORETA. In real data, the combination (wMNE/PLV) led to a very good matching between the interictal epileptic network identified from noninvasive EEG recordings and the network obtained from connectivity analysis of intracerebral EEG recordings. These results suggest that source connectivity method, when appropriately configured, is able to extract highly relevant diagnostic information about networks involved in interictal epileptic spikes from non-invasive dense-EEG data.

Whole-brain simulation of interictal epileptic discharges for patient-specific interpretation of interictal SEEG data.

In patients with refractory epilepsy, the clinical interpretation of stereoelectroencephalographic (SEEG) signals is crucial to delineate the epileptogenic network that should be targeted by surgery. We propose a pipeline of patient-specific computational modeling of interictal epileptic activity to improve the definition of regions of interest. Comparison between the computationally defined regions of interest and the resected region confirmed the efficiency of the pipeline. This result suggests that computational modeling can be used to reconstruct signals and aid clinical interpretation.

Multiscale neuro-inspired models for interpretation of EEG signals in patients with epilepsy.

OBJECTIVE: The aim is to gain insight into the pathophysiological mechanisms underlying interictal epileptiform discharges observed in electroencephalographic (EEG) and stereo-EEG (SEEG, depth electrodes) recordings performed during pre-surgical evaluation of patients with drug-resistant epilepsy. METHODS: We developed novel neuro-inspired computational models of the human cerebral cortex at three different levels of description: i) microscale (detailed neuron models), ii) mesoscale (neuronal mass models) and iii) macroscale (whole brain models). Although conceptually different, micro- and mesoscale models share some similar features, such as the typology of neurons (pyramidal cells and three types of interneurons), their spatial arrangement in cortical layers, and their synaptic connectivity (excitatory and inhibitory). The whole brain model consists of a large-scale network of interconnected neuronal masses, with connectivity based on the human connectome. RESULTS: For these three levels of description, the fine-tuning of free parameters and the quantitative comparison with real data allowed us to reproduce interictal epileptiform discharges with a high degree of fidelity and to formulate hypotheses about the cell- and network-related mechanisms underlying the generation of fast ripples and SEEG-recorded epileptic spikes and spike-waves. CONCLUSIONS: The proposed models provide valuable insights into the pathophysiological mechanisms underlying the generation of epileptic events. The knowledge gained from these models effectively complements the clinical analysis of SEEG data collected during the evaluation of patients with epilepsy. SIGNIFICANCE: These models are likely to play a key role in the mechanistic interpretation of epileptiform activity.

Signal processing and computational modeling for interpretation of SEEG-recorded interictal epileptiform discharges in epileptogenic and non-epileptogenic zones.

Objective.In partial epilepsies, interictal epileptiform discharges (IEDs) are paroxysmal events observed in epileptogenic zone (EZ) and non-epileptogenic zone (NEZ). IEDs' generation and recurrence are subject to different hypotheses: they appear through glutamatergic and gamma-aminobutyric acidergic (GABAergic) processes; they may trigger seizures or prevent seizure propagation. This paper focuses on a specific class of IEDs, spike-waves (SWs), characterized by a short-duration spike followed by a longer duration wave, both of the same polarity. Signal analysis and neurophysiological mathematical models are used to interpret puzzling IED generation.Approach.Interictal activity was recorded by intracranial stereo-electroencephalography (SEEG) electrodes in five different patients. SEEG experts identified the epileptic and non-epileptic zones in which IEDs were detected. After quantifying spatial and temporal features of the detected IEDs, the most significant features for classifying epileptic and non-epileptic zones were determined. A neurophysiologically-plausible mathematical model was then introduced to simulate the IEDs and understand the underlying differences observed in epileptic and non-epileptic zone IEDs.Main results.Two classes of SWs were identified according to subtle differences in morphology and timing of the spike and wave component. Results showed that type-1 SWs were generated in epileptogenic regions also involved at seizure onset, while type-2 SWs were produced in the propagation or non-involved areas. The modeling study indicated that synaptic kinetics, cortical organization, and network interactions determined the morphology of the simulated SEEG signals. Modeling results suggested that the IED morphologies were linked to the degree of preserved inhibition.Significance.This work contributes to the understanding of different mechanisms generating IEDs in epileptic networks. The combination of signal analysis and computational models provides an efficient framework for exploring IEDs in partial epilepsies and classifying EZ and NEZ.

Localization of extended brain sources from EEG/MEG: the ExSo-MUSIC approach.

We propose a new MUSIC-like method, called 2q-ExSo-MUSIC (q ≥ 1). This method is an extension of the 2q-MUSIC (q ≥ 1) approach for solving the EEG/MEG inverse problem, when spatially-extended neocortical sources ("ExSo") are considered. It introduces a novel ExSo-MUSIC principle. The novelty is two-fold: i) the parameterization of the spatial source distribution that leads to an appropriate metric in the context of distributed brain sources and ii) the introduction of an original, efficient and low-cost way of optimizing this metric. In 2q-ExSo-MUSIC, the possible use of higher order statistics (q ≥ 2) offers a better robustness with respect to Gaussian noise of unknown spatial coherence and modeling errors. As a result we reduced the penalizing effects of both the background cerebral activity that can be seen as a Gaussian and spatially correlated noise, and the modeling errors induced by the non-exact resolution of the forward problem. Computer results on simulated EEG signals obtained with physiologically-relevant models of both the sources and the volume conductor show a highly increased performance of our 2q-ExSo-MUSIC method as compared to the classical 2q-MUSIC algorithms.

PET imaging of brain 5-HT1A receptors in the preoperative evaluation of temporal lobe epilepsy.

[(18)F]MPPF PET has previously been used to identify the epileptic lobe in temporal lobe epilepsy (TLE) patients at the group level. This study aims to validate the visual analysis of [(18)F]MPPF PET in the assessment of individual TLE patients for their suitability to undergo temporal lobe resection. Forty-two patients suffering from TLE and 18 control subjects matched for age and gender were prospectively enrolled for [(18)F]MPPF PET. Four subtypes were defined according to the presurgical evaluation: mesio-TLE (MTLE, 32 patients), temporal neocortical epilepsy (NC, five patients), temporo-perisylvian epilepsy (T+, three patients) and temporal epilepsy without further information (t, two patients). Parametric binding potential (BP(ND)) images were obtained using a simplified reference tissue model. Three examiners, who were blinded to other data, visually interpreted each scan and delineated areas of decreased [(18)F]MPPF BP(ND). Statistical parametric mapping (SPM) analysis of MPPF BP(ND) images was also performed. Visual analysis showed a low rate of disagreement between the three examiners (7%). PET scans were considered normal in four patients (9.5%). In the remaining 38 patients (90.5%), areas of focal BP(ND) decrease were identified. A specific pattern was encountered in the MTLE subgroup, consisting of a BP(ND) decrease involving hippocampus, amygdala and temporal pole altogether. Combining the results from the presurgical investigations and the surgical outcome, we estimated that the area of BP(ND) decrease coincided with the epileptogenic zone in 40% of patients in the MTLE subgroup and 33% in the other TLE subtypes. This relatively low precision was due to 47% of patients who showed BP(ND) decreases in the insula ipsilateral to the epileptogenic lobe. The SPM analysis had much lower sensitivity (67%) to detect BP(ND) decreases in the epileptogenic temporal lobe, but revealed areas of increased BP(ND) outside the epileptogenic zone and bitemporal BP(ND) decreases of undetermined clinical significance, which were undetectable by visual analysis, in 29% of patients. In conclusion, visual analysis of [(18)F]MPPF BP(ND) images helps in the correct identification of the epileptogenic temporal lobe in all patients showing BP(ND) decreases, with a false negative rate inferior to 10% and no false positives in control subjects. All TLE patients with [(18)F]MPPF BP(ND) decreases involving hippocampus, amygdala and temporal pole together, with or without extension to the ipsilateral insula, were good candidates for anterior temporal lobectomy. All these patients became seizure free after surgery, even when the clinical presentation was not that of a typical MTLE, or when MRI failed to detect hippocampal atrophy.

COALIA: A Computational Model of Human EEG for Consciousness Research.

Understanding the origin of the main physiological processes involved in consciousness is a major challenge of contemporary neuroscience, with crucial implications for the study of Disorders of Consciousness (DOC). The difficulties in achieving this task include the considerable quantity of experimental data in this field, along with the non-intuitive, nonlinear nature of neuronal dynamics. One possibility of integrating the main results from the experimental literature into a cohesive framework, while accounting for nonlinear brain dynamics, is the use of physiologically-inspired computational models. In this study, we present a physiologically-grounded computational model, attempting to account for the main micro-circuits identified in the human cortex, while including the specificities of each neuronal type. More specifically, the model accounts for thalamo-cortical (vertical) regulation of cortico-cortical (horizontal) connectivity, which is a central mechanism for brain information integration and processing. The distinct neuronal assemblies communicate through feedforward and feedback excitatory and inhibitory synaptic connections implemented in a template brain accounting for long-range connectome. The EEG generated by this physiologically-based simulated brain is validated through comparison with brain rhythms recorded in humans in two states of consciousness (wakefulness, sleep). Using the model, it is possible to reproduce the local disynaptic disinhibition of basket cells (fast GABAergic inhibition) and glutamatergic pyramidal neurons through long-range activation of vasoactive intestinal-peptide (VIP) interneurons that induced inhibition of somatostatin positive (SST) interneurons. The model (COALIA) predicts that the strength and dynamics of the thalamic output on the cortex control the local and long-range cortical processing of information. Furthermore, the model reproduces and explains clinical results regarding the complexity of transcranial magnetic stimulation TMS-evoked EEG responses in DOC patients and healthy volunteers, through a modulation of thalamo-cortical connectivity that governs the level of cortico-cortical communication. This new model provides a quantitative framework to accelerate the study of the physiological mechanisms involved in the emergence, maintenance and disruption (sleep, anesthesia, DOC) of consciousness.

On the origin of epileptic High Frequency Oscillations observed on clinical electrodes.

OBJECTIVE: In this study we aim to identify the key (patho)physiological mechanisms and biophysical factors which impact the observability and spectral features of High Frequency Oscillations (HFOs). METHODS: In order to accurately replicate HFOs we developed virtual-brain/virtual-electrode simulation environment combining novel neurophysiological models of neuronal populations with biophysical models for the source/sensor relationship. Both (patho)physiological mechanisms (synaptic transmission, depolarizing GABAA effect, hyperexcitability) and physical factors (geometry of extended cortical sources, size and position of electrodes) were taken into account. Simulated HFOs were compared to real HFOs extracted from intracerebral recordings of 2 patients. RESULTS: Our results revealed that HFO pathological activity is being generated by feed-forward activation of cortical interneurons that produce fast depolarizing GABAergic post-synaptic potentials (PSPs) onto pyramidal cells. Out of phase patterns of depolarizing GABAergic PSPs explained the shape, entropy and spatiotemporal features of real human HFOs. CONCLUSIONS: The terminology "high-frequency oscillation" (HFO) might be misleading as the fast ripple component (200-600 Hz) is more likely a "high-frequency activity" (HFA), the origin of which is independent from any oscillatory process. SIGNIFICANCE: New insights regarding the origins and observability of HFOs along depth-EEG electrodes were gained in terms of spatial extent and 3D geometry of neuronal sources.

[Not Available]. / Démences médicamenteuses : étude cas/non-cas dans la Banque Nationale de Pharmacovigilance.

The increased incidence of dementia on the aging population makes this disease a major public health problem. Among known causes of dementia, drug etiology is under considered. We investigated the relationship between exposure to drug therapy and dementia with a case/non-case study using reports of the French Pharmacovigilance database. Among 263 962 adverse effects recorded between 1985 and 2005, 79 (0.03%) are dementia. Median age is 66 (range 3-91). There was 41 women and 37 men. The therapeutic drug class associated with dementia were anticonvulsants, antiparkinsonians, antidepressants, anxiolytics, hypnotics, antipsychotics and morphinics. An association between reporting of dementia and non neurotropic drugs were also found, i.e. interferon alfa-2B, vancomycin and allopurinol. The term "Dementia" is only mentioned in the summary of the characteristics of valproate, but it may concern other drugs. Drug etiology for dementia is a reality but is not necessarily attributed as a cause in aging population, in particular.

[Drug-induced dementia: a case/non-case study in the French Pharmacovigilance database]. / Démences médicamenteuses: étude cas/non-cas dans la Banque Nationale de Pharmacovigilance.

The increased incidence of dementia on the aging population makes this disease a major public health problem. Among known causes of dementia, drug etiology is under considered. We investigated the relationship between exposure to drug therapy and dementia with a case/non-case study using reports of the French Pharmacovigilance database. Among 263 962 adverse effects recorded between 1985 and 2005, 79 (0.03%) are dementia. Median age is 66 (range 3-91). There was 41 women and 37 men. The therapeutic drug class associated with dementia were anticonvulsants, antiparkinsonians, antidepressants, anxiolytics, hypnotics, antipsychotics and morphinics. An association between reporting of dementia and non neurotropic drugs were also found, i.e. interferon alfa-2B, vancomycin and allopurinol. The term "Dementia" is only mentioned in the summary of the characteristics of valproate, but it may concern other drugs. Drug etiology for dementia is a reality but is not necessarily attributed as a cause in aging population, in particular.

A Penalized Semialgebraic Deflation ICA Algorithm for the Efficient Extraction of Interictal Epileptic Signals.

As a noninvasive technique, electroencephalography (EEG) is commonly used to monitor the brain signals of patients with epilepsy such as the interictal epileptic spikes. However, the recorded data are often corrupted by artifacts originating, for example, from muscle activities, which may have much higher amplitudes than the interictal epileptic signals of interest. To remove these artifacts, a number of independent component analysis (ICA) techniques were successfully applied. In this paper, we propose a new deflation ICA algorithm, called penalized semialgebraic unitary deflation (P-SAUD) algorithm, that improves upon classical ICA methods by leading to a considerably reduced computational complexity at equivalent performance. This is achieved by employing a penalized semialgebraic extraction scheme, which permits us to identify the epileptic components of interest (interictal spikes) first and obviates the need of extracting subsequent components. The proposed method is evaluated on physiologically plausible simulated EEG data and actual measurements of three patients. The results are compared to those of several popular ICA algorithms as well as second-order blind source separation methods, demonstrating that P-SAUD extracts the epileptic spikes with the same accuracy as the best ICA methods, but reduces the computational complexity by a factor of 10 for 32-channel recordings. This superior computational efficiency is of particular interest considering the increasing use of high-resolution EEG recordings, whose analysis requires algorithms with low computational cost.

Automatic Detection and Classification of High-Frequency Oscillations in Depth-EEG Signals.

GOAL: Interictal high-frequency oscillations (HFOs [30-600 Hz]) have proven to be relevant biomarkers in epilepsy. In this paper, four categories of HFOs are considered: Gamma ([30-80 Hz]), high-gamma ([80-120 Hz]), ripples ([120-250 Hz]), and fast-ripples ([250-600 Hz]). A universal detector of the four types of HFOs is proposed. It has the advantages of 1) classifying HFOs, and thus, being robust to inter and intrasubject variability; 2) rejecting artefacts, thus being specific. METHODS: Gabor atoms are tuned to cover the physiological bands. Gabor transform is then used to detect HFOs in intracerebral electroencephalography (iEEG) signals recorded in patients candidate to epilepsy surgery. To extract relevant features, energy ratios, along with event duration, are investigated. Discriminant ratios are optimized so as to maximize among the four types of HFOs and artefacts. A multiclass support vector machine (SVM) is used to classify detected events. Pseudoreal signals are simulated to measure the performance of the method when the ground truth is known. RESULTS: Experiments are conducted on simulated and on human iEEG signals. The proposed method shows high performance in terms of sensitivity and false discovery rate. CONCLUSION: The methods have the advantages of detecting and discriminating all types of HFOs as well as avoiding false detections caused by artefacts. SIGNIFICANCE: Experimental results show the feasibility of a robust and universal detector.

Spatial projection as a preprocessing step for EEG source reconstruction using spatiotemporal Kalman filtering.

The reconstruction of brain sources from non-invasive electroencephalography (EEG) or magnetoencephalography (MEG) via source imaging can be distorted by information redundancy in case of high-resolution recordings. Dimensionality reduction approaches such as spatial projection may be used to alleviate this problem. In this proof-of-principle paper we apply spatial projection to solve the problem of information redundancy in case of source reconstruction via spatiotemporal Kalman filtering (STKF), which is based on state-space modeling. We compare two approaches for incorporating spatial projection into the STKF algorithm and select the best approach based on its performance in source localization with respect to accurate estimation of source location, lack of spurious sources, computational speed and small number of required optimization steps in state-space model parameter estimation. We use state-of-the-art simulated EEG data based on neuronal population models, for which the number and location of sources is known, to validate the source reconstruction results of the STKF. The incorporation of spatial projection into the STKF algorithm solved the problem of information redundancy, resulting in correct source localization with no spurious sources, and decreased the overall computational time in STKF analysis. The results help make STKF analyses of high-density EEG, MEG or simultaneous MEG-EEG data more feasible.

Source reconstruction via the spatiotemporal Kalman filter and LORETA from EEG time series with 32 or fewer electrodes.

The clinical routine of non-invasive electroencephalography (EEG) is usually performed with 8-40 electrodes, especially in long-term monitoring, infants or emergency care. There is a need in clinical and scientific brain imaging to develop inverse solution methods that can reconstruct brain sources from these low-density EEG recordings. In this proof-of-principle paper we investigate the performance of the spatiotemporal Kalman filter (STKF) in EEG source reconstruction with 9-, 19- and 32- electrodes. We used simulated EEG data of epileptic spikes generated from lateral frontal and lateral temporal brain sources using state-of-the-art neuronal population models. For validation of source reconstruction, we compared STKF results to the location of the simulated source and to the results of low-resolution brain electromagnetic tomography (LORETA) standard inverse solution. STKF consistently showed less localization bias compared to LORETA, especially when the number of electrodes was decreased. The results encourage further research into the application of the STKF in source reconstruction of brain activity from low-density EEG recordings.

A 18F-MPPF PET normative database of 5-HT1A receptor binding in men and women over aging.

UNLABELLED: Neurotransmission imaging studies require normative data for the statistical assessment of neurophysiologic dysfunctions. 2'-Methoxyphenyl-(N-2'-pyridinyl)-p-18F-fluoro-benzamidoethylpiperazine (18F-MPPF) is a specific serotonin 5-HT1A antagonist PET tracer recently characterized, modeled, and used for clinical research to explore abnormalities in the serotoninergic system. Our study reports, to our knowledge, the first large normative imaging database of 18F-MPPF binding potential (BP) over aging, for both males and females. METHODS: Fifty-three healthy volunteers (27 females, 26 males; age, 20-70 y) were selected to undergo structural MRI and single-injection 18F-MPPF multiframe dynamic PET. 18F-MPPF BP values were computed using a nonlinear modeling method with tissue reference. The statistical assessment of the effect of age and sex was performed both at the anatomic structure level, using regions of interest drawn manually on individual MR images, and at the voxel level, using normalized BP parametric images in different statistical parametric mapping designs. RESULTS: A negative linear correlation between age and 18F-MPPF binding (3.6% decrease by decade) was found in females but not in males and involved most of the limbic and paralimbic regions; on the other hand, males in their 30s showed decreased binding in most cerebral regions. CONCLUSION: A comparison of males and females revealed higher BP values independent of age in females in the right hemisphere and a different evolution of BP over aging. These results confirm the necessity of a database for further statistical analysis in individuals or groups with pathology.

Denoising of Ictal EEG Data Using Semi-Blind Source Separation Methods Based on Time-Frequency Priors.

Removing muscle activity from ictal ElectroEncephaloGram (EEG) data is an essential preprocessing step in diagnosis and study of epileptic disorders. Indeed, at the very beginning of seizures, ictal EEG has a low amplitude and its morphology in the time domain is quite similar to muscular activity. Contrary to the time domain, ictal signals have specific characteristics in the time-frequency domain. In this paper, we use the time-frequency signature of ictal discharges as a priori information on the sources of interest. To extract the time-frequency signature of ictal sources, we use the Canonical Correlation Analysis (CCA) method. Then, we propose two time-frequency based semi-blind source separation approaches, namely the Time-Frequency-Generalized EigenValue Decomposition (TF-GEVD) and the Time-Frequency-Denoising Source Separation (TF-DSS), for the denoising of ictal signals based on these time-frequency signatures. The performance of the proposed methods is compared with that of CCA and Independent Component Analysis (ICA) approaches for the denoising of simulated ictal EEGs and of real ictal data. The results show the superiority of the proposed methods in comparison with CCA and ICA.