EEG and epilepsy
Authors: Francois Tadel, Elizabeth Bock, John C. Mosher, Marcel Heers.
This tutorial introduces some concepts that are specific to the management of EEG recordings in the Brainstorm environment. It also describes a standard pipeline for analyzing epilepsy recordings. It is based on a clinical case from the Epilepsy Center at the University Hospital of Freiburg, Germany. The anonymized dataset can be downloaded directly from the Brainstorm download page.
Note that the operations used here are not detailed, the goal of this tutorial is not to introduce Brainstorm to new users. For in-depth explanations of the interface and theoretical foundations, please refer to the introduction tutorials. This page is the new version where shared channel files are not used anymore - if you really want to use shared channel files, see the ?old version.
NOT FOR CLINICAL USE:
The performance characteristics of the methods and software implentation presented in this tutorial have not been certified as medical devices and should be used for research purposes only.
- Dataset description
- Download and installation
- Import the anatomy
- Access the recordings
- Review EEG recordings
- Mark spikes
- Pre-process recordings
- Artifact cleaning with ICA
- Import the recordings
- Source analysis: Surface
- Source analysis: Volume
- Source analysis: MEM inverse solution
- Additional documentation
This tutorial dataset (EEG and MRI data) remains proprietary of the Epilepsy Centre, University Hospital Freiburg, Germany. Its use and transfer outside the Brainstorm tutorial, e.g. for research purposes, is prohibited without written consent from the Epilepsy Centre in Freiburg. For questions please contact A. Schulze-Bonhage, MD, PhD: email@example.com
This tutorial dataset was acquired in a patient who suffered from focal epilepsy with focal sensory, dyscognitive and secondarily generalized seizures since the age of eight years. He does not have any typical risk factors for epilepsy. The high resolution 3T epilepsy MRI including postprocessing was found to be normal. FDG-PET of the brain did not show any pathological changes in the glucose metabolism. Non-invasive telemetry revealed left fronto-central sharp waves, polyspikes and bursts of beta band activity (max. amplitude FC1, Cz) especially during sleep. The tutorial dataset was acquired during one night of the non-invasive telemetry recording at the Epilepsy Center Freiburg, Germany.
Afterwards the patient underwent invasive EEG to identify the epileptogenic area and to map functionally important cortex. Details about invasive EEG and source localization from invasive EEG in this patient are reported in Dümpelmann, et al. (2011). Subsequently a left frontal tailored resection was performed. The histological analysis revealed a focal cortical dysplasia type IIB according to the classification of Palmini, et al. (2004). The postsurgical outcome is Engel 1A with a follow-up of 5 years.
The EEG data distributed here was recorded at 256Hz, using a Neurofile NT digital video-EEG system with 128 channels and a 16-bit A/D converter. The signal was filtered in the recording system with a high-pass filter with a time constant of 1 second (cutoff frequency ~ 0.16Hz) and a low-pass filter with a cutoff frequency of 344 Hz. The spikes were marked with Brainstorm by the epileptologists at the Epilepsy Center in Freiburg.
Overview of the data processing
The proper identification of epileptiform discharges ("spikes") is a complicated topic beyond the scope of this tutorial. Residents in Neurology train intensively on how to identify true interictal and ictal discharges, particularly as distinct from so-called "normal variants" that can be abundant. Sleep and drowsy states of the brain can generate "vertex waves," "K-complexes," "positive occipital sharp transients of sleep" (POSTS), and "wickets," to name a few variants, none of which are epileptic. A good overview of terminology and application can be found, for example, on the Medscape website, Epileptiform Discharges.
With this caveat, we nonetheless give an overview of the processing approach. Although automation has been proposed for decades, "spike hunting" is often done by manual inspection of the recorded EEG waveforms (whether scalp, subdural, or depths). The classic scalp sensor arrangement is the International 10-20 System (Wikipedia Site, Bioelectromagnetism Book, Chapter 13.3) arranged at 10 percent and 20 percent distances about the circumference of the head. Each of these 21 electrodes is acquired with respect to some reference (as all potentials must be). In reviewing the recordings, however, several "montages" are traditionally recommended that digitally "re-reference" the original recordings into new linear combinations of electrodes.
A classic montage is the "double banana" (Google Search) which emphasizes local changes in the scalp EEG by forming sequential bipolar pairs, such as "Fp1-F3", "F3-C3." Because F3 is found twice in this montage and of opposite sign, then an epileptic spike centered under F3 will appear as a reversed polarity in these two channels, a visual cue the trained epileptologist seeks when rapidly scanning through hours of recordings. This montage is more formally known as "LB-18.3" or "Longitudinal Bipolar 3" in the nomenclature of the American Clinical Neurophysiology Society (ACNS) Guidelines (see References Below).
ACNS guidelines suggest using both "longitudinal" and "transversal" bipolar montages to survey your data. For temporal epilepsy cases, you may also add a "temporal ring" montage. European neurologists often prefer to review the recordings using an average reference montage. Brainstorm provides several variations of all these montages, allows the users full flexibility in creating their own, and can run different montages simultaneously in multiple windows.
The user has to step along through the data, look for abnormal brain activity, then use Brainstorm's event markers to tag suspect intervals. With the suspected events marked and saved, the user can return later to perform source analyses on these intervals. With this brief overview, we detail below an exercise with the sample epilepsy data.
Note that the aim of this tutorial is not to train you on how to do spike recognition itself, but rather to illustrate how to manipulate the Brainstorm interface to build an environment adapted to this task. The formal data reviewing should be done by clinical neurophysiologists or other experienced personnel, with specialized training in separating true epileptic activity from other more normal variants of brain activity.
Dümpelmann M, Ball T, Schulze-Bonhage A (2011)
sLORETA allows reliable distributed source reconstruction based on subdural strip and grid recordings. Human Brain Mapping.
Palmini A, Najm I, Avanzini G, Babb T, Guerrini R, Foldvary-Schaefer N, Jackson G, Luders HO, Prayson R, Spreafico R, Vinters HV (2004) Terminology and classification of the cortical dysplasias.
Standard montages recommended by the American Clinical Neurophysiology Society:
Download and installation
Requirements: You have already followed all the introduction tutorials and you have a working copy of Brainstorm installed on your computer.
Go to the Download page of this website, and download the file: sample_epilepsy.zip
Unzip it in a folder that is not in any of the Brainstorm folders (program folder or database folder)
- Start Brainstorm (Matlab scripts or stand-alone version)
Select the menu File > Create new protocol. Name it "TutorialEpilepsy" and select the options:
"No, use individual anatomy",
"No, use one channel file per acquisition run".
Import the anatomy
- Switch to the "anatomy" view of the protocol.
Right-click on the TutorialEpilepsy folder > New subject > sepi01
- Keep the default options you set for the protocol
Right-click on the subject node > Import anatomy folder (auto):
Set the file format: "FreeSurfer"
Select the folder: sample_epilepsy/anatomy
The "auto" import menu computed automatically a linear MNI normalization, default 3D positions for the anatomical fiducials NAS/LPA/RPA, and downsampled the cortex surface to 15000 vertices.
Click on [Save] to validate the modifications and continue with the importation of the anatomy. At the end, make sure the file "cortex_15000V" is selected (downsampled pial surface, that will be used for the source estimation). If it is not, double-click to select it as the default cortex surface.
Without the individual MRI
If you do not have access to an individual MR scan of the subject (or if its quality is too low to be processed with FreeSurfer), but if you have digitized the head shape of the subject using a tracking system, you have an alternative: deform one of the Brainstorm templates (Colin27 or ICBM152) to match the shape of the subject's head. For more information, read the tutorials: Warping default anatomy, Using the anatomy templates
Access the recordings
Link the recordings
- Switch to the "functional data" view (2nd button, on top of the database explorer).
Right-click on the subject folder > Review raw file:
- Select the file format: "EEG: Deltamed Coherence-Neurofile(*.bin)"
Select the file: sample_epilepsy/data/tutorial_eeg.bin
The new files "Link to raw file" let you access directly the contents of the original EEG files. The menu "Review raw file" does not actually copy any data to the database. More information.
- The channel file "Deltamed channels" contains the name of the channels, but not their positions. We need to import separately the positions of the electrodes (either a standard cap or accurate positions digitized with a Polhemus device or any other tracking system).
Prepare the channel file
- Import the positions of the electrodes:
Right-click the channel file > Add EEG positions > Import from file:
- Select the file format: "EEG: ANT ASA/Xensor (*.elc)"
Select the file: sample_epilepsy/data/tutorial_electrodes.elc
- This file contains the default electrodes positions from the ASA software (ANT).
- If you have digitized the position of the electrodes (eg. with a Polhemus device) you have to make sure it contains the position of the anatomical references NAS/LPA/RPA. This is the best way to make sure the electrode cap will be registered automatically with the MRI in Brainstorm.
- The channels are matched by name: the position file you import must include the labels of the electrodes and they must be named exactly in the the same way as in your recordings.
In this example, you will get a question "Apply MRI transformation?". If you answer yes, it would apply the voxel=>world transformation initially available in the MRI volume. Whether you need to answer Yes or No depends on which coordinates system is closer to the one use in your position file. Try both and you'll see which one brings the electrodes closer to the head of the subject. Here, answering NO would give a better immediate result, but as we will run the automatic alignment electrodes/head surface later, the answer to this question doesn't matter much.
- The recordings contain signals coming from different types of electrodes:
- 29 EEG electrodes
- EOG1, EOG2: Electrooculograms
- EMG, ECG: Electromyogram and electrocardiogram
- SP1, SP2: Sphenoidal electrodes
- RS: Electrode on the right shoulder
- PHO: Photo stimulation channel
- DELR, DELL, QR, QL: Not used
Right-click on the channel file > Edit channel file:
- Note that the EOG, EMG and ECG channels have their type correctly detected.
- All the other non-EEG channels were set to "EEG_NO_LOC" when we imported the channel locations: SP1, SP2, RS, PHO, DELR, DELL, QR, QL
It happens commonly that the file format for the electrodes positions does not describe the type of the channels correctly, typically all the signals saved in the files are classified as EEG. You can edit the names and types of the channels with this in interface: select one channel, click on the type and select something in the drop-down list. To edit multiple channels at once: right-click in the window > Set channel type.
For this particular study, we can use the channel file as it is configured now, just close the figure. Discard any modification you may have done.
Register electrodes with MRI
Right-click on the channel file > MRI registration > Edit.
Hover your mouse over the toolbar buttons for instructions. Click on the [Label] button in the toolbar to show the names of the electrodes. The default positions are already quite good, and the head shape is correct, only limited manual registration will be required in the next few steps.
Click on the button "Refine registration using head points" to find a better registration between the head shape defined by the electrodes and the head surface coming from the MRI. Note that this works only if have accurate head shapes and electrodes positions. If you don't, skip this step and edit the registration manually.
Click on the button "Project electrodes on scalp surface", to ensure all the electrodes touch the skin surface. Then click on "OK" and agree to save the modifications.
To see again the electrodes, right-click on the channel file > Display sensors > EEG (head).
To see or edit the positions of the electrodes in the MRI Viewer: right-click on the channel file > Display sensors > EEG (MRI Viewer). Select the menu "MIP: Functional" to see all the electrodes. To edit the channel file: right-click > Electrodes > Set electrode position.
Review EEG recordings
Display the recordings
Open a time series figure with the "Average reference" montage
Right-click on the "Link to raw file" > EEG > Display time series.
- In the Record tab, select the first button in the toolbar (Display mode for time series) to view the signals in columns.
- Select "Average reference" in the drop-down list. This menu allows you to set a predefined montages to linearly re-arranging the waveforms. Note the keyboard shortcuts in the menus.
Change the default duration that is reviewed to 10s.
In the figure, select the display option "Flip Y Axis" to have the negative values pointing up (convention used by many neurologists).
Open a 2D Sensor cap map of the EEG sensor values:
Right-click on "Link to raw file" again > EEG > 2D Sensor cap
In the Record tab, set the Montage to this view to "Average reference"
Open the ECG and EOG traces, to avoid confusing spikes with cardiac or ocular artifacts:
Right-click on the Link to raw file > ECG > Display time series
Right-click on the Link to raw file > EOG > Display time series
- The ECG is almost mandatory. The EOG is optional: it can be helpful for beginners but an experienced reviewer will easily recognize the eye movements directly in the EEG data.
Re-arrange the figures in a convenient way, for example as illustrated below. Then disable the automatic positioning of the figures (layout menu at the top-right of the Brainstorm figure > None), so that your figure arrangement doesn't get lost when you open a new figure.
Having a lot of windows open may slow down significantly the display because each time you change the current time, all the figures have to be updated. A lot of space is also wasted on the screen due to window frames. The number of windows to open has to be balanced between the amount of information to display and the ease of use.
Go to the Filter tab to enable some display frequency filters. General recommendations are:
High-pass filter: 0.5 Hz
Low-pass filter: 80 Hz
Note that if you have filters selected in this panel, the display of the EEG signals will be slower. Each time you will go to the next page of recordings, the filters will be applied on the fly to the recordings. The computation time is not very long at each page, but can become annoying when reviewing a lot of data. For a faster display of filtered signals, you may consider apply the filters to the file (with the process Pre-process > Band-pass filter) and then review the recordings without visualization filters.
Time and amplitude resolution
The resolutions of the time and amplitude axes have a lot of importance for the visual detection of epileptic spikes. The shapes we are looking for are altered by the horizontal and vertical scaling. The distance unit on a screen is the pixel, we can set how much time is represented by one pixel horizontally and how much amplitude is represented by one pixel vertically.
In the Brainstorm interface, this resolution is usually set implicitly: you can set the size of the window, the duration or recordings reviewed at once (text box "duration" in tab Record) and the maximum amplitude to show in the figure (buttons [...] and [AS] on the right of the time series figure). From there, you can also zoom in time ([<], [>], mouse wheel) or amplitude ([^], [v], Shift+mouse wheel). These parameters are convenient to explore the recordings interactively but don't allow us to have reproducible displays with constant time and amplitude resolutions.
To set the figure resolution explicitly: right-click on the figure > Figure > Set axes resolution. Note that this interface does not store the input values, it just modifies the other parameters (figure size, time window, max amplitude) to fit the resolution objectives. If you modify these parameters (resize the figure, keep the button [AS] selected and scroll in time, etc) the resolution is lost, you have to set it again manually. In particular, make sure you disable the auto-scaling ([AS] button in the time series figure) if you want to preserve the aspect ratio while you scroll through the data.
This operation typically has to be repeated everytime you open a new file. For a faster access to this option, use the keyboard shortcut Ctrl+O. The option window offers by default the last values that you entered, just press Enter to apply them again.
Recommendations for this dataset are:
Time axis: 170 pixels/sec (~55 mm/sec)
Amplitude: 15 μV/pixels (~45 μV/mm)
This preparation of the reviewing environment requires a large number of operations, and would become quickly annoying if you have to repeat it every time you open a file. You can use the menu "User setups" to save a window configuration and reload it in one click later. In the menu "Window layout", at the top-right of the Brainstorm window, select User setup > New setup. Enter a name of your choice for this particular window arrangement.
This operation will also disable the automatic window arrangement (Window layout > None). To reload it later, open one figure on the dataset you want to review and then select your new entry in the User setup menu.
It may be interesting for some cases to display different groups of sensors in multiple windows (eg. with an MEG system with 300 sensors), or some complicated epilepsy cases where you would like to review at the same time multiple montages (eg. longitudinal and transversal bipolar montages).
Open your full reviewing environment as described before, where the EEG signals are displayed with the "Average reference" montage.
Open another view on the same data with the "Longitudinal 3" montage ("double-banana" LB-18.3)
Right-click on the "Link to raw file" again > EEG > Display time series
Alternatively, you can right-click on the existing figure > Figure > Clone figure.
- Then set the montage for this new figure to "Longitudinal 3".
- Resize all the figures to make room for the new window.
- Save this window configuration as a new "User setup".
- If you don't see the "Longitudinal 3" menu, it is probably because you have been using Brainstorm before these predefined montages were made available in the distribution. To add them manually:
- In the Record tab, select "Edit montages" in the drop-down menu
- Click on the "Load montage" button.
- Go to the folder "brainstorm3/toolbox/sensors/private/", and select the first file.
- Note that a new entry (probably "Longitudinal 1") is added to the list of available montages.
- Repeat the operation with all the files in the folder "brainstorm3/toolbox/sensors/private/".
- Click on "Save" to close the montage editor and now select "Longitudinal 3".
More information available in the tutorial Montage editor.
Scalp current density
In the example below, see how the montage Scalp current density can enhance the visual detection of spikes. More information.
When you have a clean ECG signal for your patient, you can automatically identify all the heartbeats in the recordings. Because heartbeats can cause sharp waves in some EEG traces, it helps the reviewing process to have all the cardiac events marked in the recordings.
Right-click on the "Link to raw file" > EEG > Display time series (or simply double-click on it).
In the tab Record, menu Artifacts > Detect heartbeats: Channel=ECG, All file.
Import the spike markers
Some spikes were marked by the epileptologists at the Epilepsy Center in Freiburg and saved in an external text file. We are going to import this file manually.
In the tab Record, menu File > Add events from files...
Select format Array of times (text file containing the timing of the markers)
Select file sample_epilepsy/data/tutorial_spikes.txt
When prompted, enter the event name "SPIKE"
A new category SPIKE is visible in the events list, containing 58 markers. Click on a few of them and try to identify the shape of the spike (mostly visible on the channel FC1). Then close the viewer and save the modifications.
- The two other types of events that were present initially in the file (REM/REM_Ende) indicated the beginning and the end of periods or REM sleep (the patient is sleeping during the entire session).
If you are marking the events by yourself, you could follow this procedure:
- Close all the current figures ("Close all" button at the top-right corner of the Brainstorm window).
- Double-click on the "Link to raw file" to open a continuous file viewer, and load your reviewing environment (menu User setups).
Start by creating a group of events (Events > Add group), and select it in the list of events.
Make sure that the time and amplitude resolutions are what you are used to
(right-click on the figure > Figure > Set axes resolution)
Scroll through the recordings using the [<<<] and [>>>] buttons or shortcuts such as F3 or F4 (complete list and descriptions available when you hover your mouse over these buttons).
You can adjust the gain of the electrodes to observe better an event with the buttons [^] and [v], with the keyboard (+/-) or the mouse ([Shift+mouse wheel] or [Right-click+move up/down]).
- When you identify a spike, click in a white area of the figure in order to place the time cursor at the peak of the spike. If you click on the signal itself, it selects the corresponding channel, but you can use the shortcut Shift+Click to prevent this behavior and force the time cursor to be moved instead.
- Press Ctrl+E to add a marker where the time cursor is.
If you are marking multiple types of events, it is convenient to set up some additional keyboard shortcuts. Using the menu Events > Edit keyboard shortcuts, you can associate custom events to the keys 1 to 9 of the keyboard. Define the name of the event type to create for each key, and then simply press the corresponding key to add/delete a marker at the current time position.
- To jump to the next/previous event in the current category: use the keyboard shortcuts [Shift+arrow right] and [Shift+arrow left]
More information on the data viewer, see tutorial: Review continuous recordings.
Two of the typical pre-processing steps consist in removing the power lines artifacts (50 Hz or 60Hz) and the frequencies we are not interested in (a low-pass filter to remove the high-frequency noise and a high-pass filter to remove the very slow components of the signals). Let's start with a spectral evaluation of this file.
- In the Process1 box: Drag and drop the "Link to raw file".
Run process Frequency > Power spectrum density (Welch): All file, Length=10s, Overlap=50%.
Double-click on the new PSD file to display it.
This frequency spectrum does not show any particular peak at 50/60Hz, there is no notch filter to apply on these recordings. If we had to, we would run the process "Pre-processing > Notch filter" as explained in the tutorial Detect and remove artifacts.
The filters we selected for reviewing the recordings were for visualization only, they were not applied to the file. In order to apply these filters permanently to the recordings, we need to do the following:
- Keep the "Link to raw file" selected in the Process1 list.
Run process Pre-process > Band-pass filter: [0.5,80]Hz, 60dB, no mirror, sensors=EEG
- Note that this new continuous file is saved in your Brainstorm database, while the original file is saved in a separate folder (sample_epilepsy). If you delete the link to the original file with the database explorer, it would not delete the actual file. If you delete the link to the filtered file, it would delete the file itself.
The average reference montage we used to review the recordings was for visualization only, it did not modify the file. To apply this transformation permanently, we need to run the re-referencing process.
- Double-click on the filtered file "Raw | band(0.5-80Hz)" to open it.
In the Record tab, select the menu Artifacts > Re-reference EEG: Enter AVERAGE.
The re-referencing is saved in the same way as the SSP and ICA projectors: a linear operator applied dynamically to the recordings. At the end of the process, the window "Select active projectors" is shown, so that you can control the list of projectors currently applied to the file. Close this figure and close the recordings.
Average reference and bad channels
You should be careful with the average reference when adding or removing bad channels. Changing the number of channels modifies the way the average reference is computed.
Montage: When using the montage "Average reference" to change only the visulization, the average reference operator is computed dynamically, it always takes into account the bad channels at the moment you review the recordings. Therefore the average reference for display is always correct, you do not have to modify anything if you change the list of bad channels.
Permanent re-referencing: When using the process "Re-reference EEG" to alter permanently the recordings, the average reference operator is computed using the channels that are marked as bad at the moment of the computation, and saved in the file. If later you modify the list of bad channels, this operator will become incorrect, you will have to delete it and recompute it. To avoid any confusion, always compute the average reference after marking the bad channels.
Evaluate the modifications
- In Process1, select the filtered file "Raw | band(0.5-80Hz)".
Run again the process Frequency > Power spectrum density (Welch), with same options.
Double-click on the new PSD file to display the spectrum of the filtered file. You can observe the effects of the filter (drops on the two sides of the spectrum, zoom in to observe it better) and the average reference (all the signals have now a similar level of power).
Artifact cleaning with ICA
The standard cleaning pipeline presented in the introduction tutorials recommends to use the SSP approach to remove the cardiac and blink artifacts. This is valid in most MEG experiments.
In EEG, there is usually not such a clear artifact caused by the heartbeats: cleaning for cardiac artifacts is not part of the standard processing pipeline. Eye movements typically cause huge artifacts on the frontal channels and need to be marked as bad or corrected. In the present case, the patient is sleeping, therefore there are no blinks, but we can observe a lot of fast eye movements during the periods of REM sleep, and other slower eye movements in the rest of the recordings.
The SSP technique is not adapted for cleaning events for which we cannot set markers, and in general it is not suitable for low-density EEG. The topographies are a lot smoother in EEG than in MEG, it is more difficult to identify a spatial filter that targets one specific artifact with a simple spatial PCA decomposition, especially with a low number of sensors.
The ICA approach (Independent Component Analysis) is much more popular in the EEG community. It relies on the the same logic as the SSP: identifying spatial topographies that are specific from an artifact and then removing them from the recordings. But instead of simply identifying the most powerful spatial components with a PCA decomposition, it tries to identify components that are independent in time. Using the information from the time dimension makes it possible to work with a lower number of sensors, and to identify some components of the signal that are completely uncorrelated with the others.
The ICA analysis can also be used to identify the signals from a few brain sources of interest: in this case, select a few components and remove all the others.
Various algorithms are available from the ICA process. You can explore the documentation of all these methods for reference documentation on ICA:
We will now illustrate how to compute some projectors to remove the contaminations from the eye movements with the Picard algorithm.
To make it run fast on most computers, we will not process the full recordings but only a short segment which contains a lot of eye movements: [500,700]s.
- Double-click on the filtered file "Raw | band(0.5-80Hz)" to open it.
In the Record tab, select the menu Artifacts > ICA components.
Time window: Time segment where the artifact is present.
Note that all the BAD segments are removed from the selected data before executing the ICA algorithm, causing some discontinuities in the signals.
Event name: Empty, we want to use all the recordings from the selected time window.
Event window: Not applicable (because the previous option is empty).
Sensor types: EEG
Band-pass filter: Ignored (0-0Hz), we want to perform the analysis on the broadband signal.
Resample: Disabled (0 Hz)
ICA algorithm: Picard
Number of ICA components: If this value N is higher than 0 and lower than the number of channels, it will first perform a spatial PCA decomposition, keep only the N principal components, and pass this subspace to the ICA algorithm. Enter 0 here.
Sort components based on correlation with: Sort ICA components based on their correlation with a selection of 'reference' channels (default: ECG & EOG). The sorting is based on the maximum correlation of a given ICA time-series with any of the reference channels. If empty, the components are sorted in decreasing order of the EEG variance accounted for by each component.
Save averaged artifact: Not applicable (because we are not using events).
You can follow the convergence of the selected algorithm in the Matlab command window. When it is done, the window "Select active projectors" is shown. Select all the components, display their spatial topography and time series (read the tooltips to see which button does what). The first two components (IC1, IC2) look very correlated with the EOG signals, and their topographies are typical of eye movements.
If the algorithm doesn't converge: try to run the ICA decomposition before the EEG referencing, or limit the number of ICA components. In this example, the rank of the data (28) is lower than the number of signals (29), this is due to re-referencing EEG signals, and causes Infomax and FastICA to fail converging. If you set the parameter "number of ICA components" to 28, it would solve the problem.
Select the first two components to remove them: observe the effect on the EEG signals. It removes very well the artifacts due to the eye movements without affecting the rest of the recordings.
Go to around a SPIKE event, make sure that removing these two components doesn't affect the shape of the spike. We do not want them to alter our events of interest.
- Select the first two components, click on [Save] to keep this selection.
To save the ICA time series in a new file in the database, follow the instructions in the SSP tutorial: Extract the time series.
For practicing the visual identification of ICA components for various types of artifacts: https://labeling.ucsd.edu/tutorial/practice
Examples of results obtained with the four methods available in Brainstorm in this example. From left to right: method name and execution time, topography of the first two components, EOG signals and time series of the first two components.
Import the recordings
Now we want to extract all the spikes from this filtered file and import them in the database.
Right-click on the filtered file (Raw|band) > Import in database.
Select the option "Use events" and select the category SPIKE, with 58 events.
Epoch time: [-300, +500]ms around the spike events.
- Remove DC offset: Disabled, the DC offset has already been removed with the high-pass filter.
- Create a separate folder for each event type: Yes, it will import the epochs into a folder "SPIKE" instead of a folder named after the EEG original filename.
You should see a new folder "SPIKE" containing 58 epochs in your database. The channel file gets duplicated for each new folder.
- We recommend you always review all the trials after importing them. To review the imported trials: double-click on the first trial, switch back to butterfly view (first button in Record tab), then navigate between trials with F3 / Shift+F3. Channels and trials can be marked as bad independently.
You can overlay one or more electrodes from multiple files with the tab "Cluster":
- Display the time series for one trial (you must have at least one file loaded to create a cluster).
In the Brainstorm window: Click on the [+] to add the Cluster tab.
- In the Cluster tab: Click on [New Ind] and enter "FC1" (cluster with one electrode only).
- In the Cluster tab: Double-click to rename the cluster into something more meaningful ("FC1").
In the Cluster tab: Select the option "Overlay: Files"
In the database explorer: Right-click on the group of trials > Clusters time series.
Alternatively, you can display all the trials as an image, for one sensor at a time:
In the database explorer: Right-click on the on the group of trials > Display as image > EEG.
- Drag and drop all the SPIKE files or the SPIKE folder into the Process1 tab.
Run process Average > Average files, By folder (subject average), Arithmetic average.
Double click on the new file "Avg: SPIKE" to review it. Add a "2D sensor cap" view (popup or Ctrl+T).
Go to the Filter tab to disable all the visualization filters: these filters are not working well on short signals, the high-pass filter at 0.5Hz can introduce errors in the display of the average.
- To have all the figures automatically re-arrange themselves again, select "Weighted" or "Tiled" in the "Window layout" menu (top-right corner).
Explore this average in time with the left and right keyboard arrow. You can add more views to your reviewing environment (see the tutorial Visual exploration):
Typically, the analysis could be limited to the ascending part of the peak (from -11ms to 0ms), as it gives a rather clear information on the primary epilepsy focus. For the rest of the tutorial, we will use a larger time window (-40ms to 100ms) in order to illustrate the visualization tools. Just keep in mind that it can be hazardous to interpret waves that are far in time or in space from the primary focus.
To get a synthetic view at the sensor level of the evolution of the 2D map in time:
Right-click on the figure > Snapshot > Time contact sheet: Figure: [-40ms, 110ms], 16 images.
Source analysis: Surface
We are going to use a realistic head model. This requires additional surfaces to model the skull properly. Go to the "Anatomy" view, right-click on the subject > MRI segmentation > Generate BEM surfaces, method = Brainstorm.
Use 1922 vertices for each layer (default). Three new surface files are created:
Go back to the "Functional data" view, right-click on the channel file in the SPIKE folder > Compute head model. Select "cortex surface" and "OpenMEEG BEM". Leave all the OpenMEEG options to their defaults. Depending on your computer, this may take between 10 and 60 minutes.
At the first execution, OpenMEEG will be downloaded and installed as a Brainstorm plugin.
- Open your resource monitor to check the memory usage: if it goes over 98%, your system doesn't have enough RAM to compute this model. You should kill the process "om_gain", kill Matlab and try again with a different configuration.
If the automatic download doesn't work, download and install OpenMEEG manually (menu Help). If the OpenMEEG calculation crashes, please refer to the BEM tutorial.
If you don't have enough memory, try with the option Use adjoint formulation. Another easy way to make this work a lot faster and with a limited amount of memory is to use less vertices for each layer, for example 1082 vertices (scalp) and 642 vertices (outer skull and inner skull). However, with less refined models you may run into other issues (intersecting meshes or numerical instabilities).
- If you can compute the model but the source maps you get don't make any sense, it could be because the BEM surfaces were poorly defined. When some cortex vertices are too close to the inner skull surface, the model can be unstable and you would typically see source maps containing zeros everywhere except for a few vertices. If this happens, delete the existing BEM surfaces and try computing new ones with a higher number of vertices for the inner and outer skull surfaces.
If you cannot get OpenMEEG to work, or if the results definitely do not make sense, try using a different forward model: "3-shell sphere". It's a spherical model, so it would perform better in the regions of the head that are close to the sphere. See the Head model tutorial.
Noise covariance matrix
For minimum norm inverse models, we need to estimate the level of noise at the level of the sensors. Defining what can be considered as "noise" in the middle of continuous brain recordings is a difficult question, discussed in the Noise covariance tutorial.
Here we choose what is described as option 1b in the section "Epilepsy" in that tutorial: we are going to calculate the noise covariance matrix using 50 seconds from the continuous file that are very "quiet", ie. that contain no apparent epileptic activity, no REM and no other special artifacts. We can use for this purpose the segment between 110s and 160s.
Right-click on the filtered file "Raw|band" > Noise covariance > Compute from recordings: Baseline = [110,160]ms, keep the other options to the default values.
The noise covariance gets saved in the same folder as the recordings, you need to copy it manually to the other folders. Right-click > Copy to other folders.
Right-click on the OpenMEEG head model > Compute sources : sLORETA, Unconstrained
It creates a shared inversion kernel and a source link for each file in the folder. If you are not familiar with these concepts, please refer to the Source estimation tutorial.
To display the same information re-interpolated in the volume, right-click on the source file > Cortical activation > Display on MRI. Right-click on the figure to set the smoothing level and other options.
You can save movies or contact sheets of the source maps with the Snapshot popup menu
Regions of interest
- From the displays we generated previously, there is no clear propagation of the spike. The scouts can be an interesting tool to explore further this spatial propagation.
Start by placing a small scout at the focus, then place a few other scouts of a similar size at other places around the spike location, where you could suspect a possible propagation. Plot the overlaid time series for these scouts. If you need help with this: see the Scouts tutorial.
Source analysis: Volume
If the results you obtain are not satisfying with the surface-based source estimation, you can run again the same analysis with a different source space, sampling the entire head volume instead. More information about this method in the Volume source estimation tutorial. We will use this new example to estimate the dipole scanning tools.
Right-click on the channel file > Compute head model, select MRI volume. To go faster this time, you can select a spherical head model: 3-shell sphere. Use the default sphere.
Then you have to define what source grid you want to use. Keep the default options.
Right-click on the new head model > Compute sources : Dipole modeling.
Note that this is not a global inverse solution like the minimum norm. It represents the performance of each individual dipole, ie. its capacity to explain all the recordings at one time sample on its own.
- In Process1, select the Dipoles link associated with the average spike.
Run process Sources > Dipole scanning: [-40,100]ms
Double-click on the new dipoles file. In the Dipoles tab, set the threshold Goodness to 90% to show only the most meaningful dipoles. The color of the dipoles below represent the time at which each dipole was estimated. This representation seems to indicate there are two distinct brain regions: the activity at 40ms is deeper than at the peak of the spike (0ms).
More information about dipole scanning: Dipoles: Scanning and displaying.
The dipole scanning approach presented in the previous section estimates the activity of a fixed grid of dipoles and selects the most significant one at each time. Another option to represent the source activity as a sequence of moving dipoles is to use the FieldTrip function ft_dipolefitting.m to run a non-linear dipole fit at each time point. This process requires that you install FieldTrip on your computer first. Note that it works only with a spherical head model.
- In Process1, select the EEG recordings for the average spike.
Run process Sources > FieldTrip: ft_dipolefitting: Time window=[-40,100]ms, Sensors=EEG.
Displaying the dipoles with a goodness of fit above 95% produces very similar figures. It is a good way to cross-validate the dipole results you obtain: if these two approaches lead to completely different interpretations, you should not trust either one.
More information about dipole fitting: FieldTrip dipole fitting
We will illustrate two ways of computing time-frequency maps. Exploring the continuous recordings in time-frequency space can help with the detection of the epileptic activity, while computing the average of the time-frequency power of a few epochs can help understanding the dynamic of the event.
Right-click on the filtered continuous file (Raw|band) > Import in database. Select a block of 10s between 1890s and 1900s (this segment contains a lot of spikes). Do not use events.
- In Process1, select the new imported file: "Raw (1890.00s,1900.00s)"
Run process Frequency > Time-Frequency (Morlet wavelets):
Sensors=EEG, 1/f compensation, Frequency=Log 2:40:80
Display the new time-frequency file. In the Display tab, select channel FC1. You can notice that the vertical stripes in the time-frequency map match the spikes in this segment of file. The epileptic brain activity can be easier to detect in time-frequency space than in the original signals.
- In Process1, select all the imported spike epochs.
Run process Frequency > Time-Frequency (Morlet wavelets):
Sensors=EEG, Spectral flattening=None, Frequency=Log 2:40:80, Save average.
- In Process1, select the new time-frequency average.
Run process Standardize > Baseline normalization: Baseline=[-200,-100]ms, ESR/ERD.
Right-click on the average > All channels, then adjust the colormap.
Source analysis: MEM inverse solution
The maximum of entropy on the mean (MEM) framework solution gathers several source localization techniques dedicated to estimate accurately the spatial extent of the EEG/MEG generators, a very important property in the context of localizing the generators of epileptic discharges, for which it is known that a relatively extended brain region should be involved (several cm2) to detect these discharges from the scalp. You can read more information on MEM in the Brain entropy in space and time (BEst) in the BEst tutorial.
Here are the steps to perform MEM localization on epileptic spikes.
Head model: Double-click on the head model computed on the cortical surface to select it.
(MEM was not implemented for the volume head models)
Baseline: MEM solution relies on a data-driven parcelization of the brain. To do so, it requires some baseline data that can be taken either from a segment of signal of interest or from another baseline section with no epileptic discharges extracted from the data (typically 2s).
Inverse model: Right-click on the average signal "Avg: SPIKE (58 files)" > Compute sources .
The MEM option panel appears. You have the choice between 3 methods. cMEM is the standard MEM in the time domain. wMEM focuses in the oscillatory patterns by computing MEM in the wavelet domain. rMEM is dedicated in finding synchronous sources using continuous wavelets.
Select cMEM for this tutorial. Select a time window from -0.1s to 0.3s under the option Data definition.
To select the baseline, select within data and select a time window from -0.3s to -0.1s. This will select -0.3s to -0.1s from "Avg: SPIKE (58 files)". There are 2 more options available to select baseline:
o Find- if you want to select an external baseline saved in the protocol, type the name of that file next to option Find and select Find. (e.g., if your baseline is named ‘xyz’, type ‘xyz’ and select find)
o Import-to import a file from a specific brainstorm database.
To learn how to use the option parameters of cMEM, we encourage you to read the MEM tutorial. Here, we will use the default values. Press ok to start the source estimation; MEM is an iterative localization approach, so computation time will depend on the time window selected.
Double click on the source file that was created to display the results
EEG reference: Reference of the gain matrix
EEG reference: Adding missing channel to the recordings
Registration: Match known EEG locations to the MRI image
Registration: Flipping EEG electrodes left-right
Pre-processing: Bad channel replacement
Minimum number of electrodes for EEG source modeling
Song J, Davey C, Poulsen C, Luu P, Turovets S, Anderson E, Li K, Tucker D
EEG source localization: Sensor density and head surface coverage
Journal of Neuroscience Methods, Dec 2015
Michel CM, Brunet D
EEG source imaging: A practical review of the analysis steps
Frontiers in Neurology, Apr 2019
Inverse models in epilepsy
Distributed source models vs. ECD (equivalent current dipole)
Kobayashi K, Yoshinaga H, Ohtsuka Y, Gotman J (2005)
Dipole modeling of epileptic spikes can be accurate or misleading
Epilepsia, 2005 Mar;46(3):397-408.
Distributed source models vs. BOLD
Heers M, Hedrich T, An D, Dubeau F, Gotman J, Grova C, Kobayashi E (2014)
Spatial correlation of hemodynamic changes related to interictal epileptic discharges with electric and magnetic source imaging. Human Brain Mapping, published online 24 Feb 2014.
ECD vs. BOLD
Benar CG, Grova C, Kobayashi E, Bagshaw AP, Aghakhani Y, Dubeau F, Gotman J (2006)
EEG–fMRI of epileptic spikes: Concordance with EEG source localization and intracranial EEG
MEM in clinical epilepsy
Related methodological developments:
Amblard C, Lapalme E, Lina JM (2004)
Biomagnetic source detection by maximum entropy and graphical models, IEEE Trans Biomed Eng.
Lapalme E, Lina JM, Mattout J (2006)
Data-driven parceling and entropic inference in MEG, NeuroImage.
Ability of MEM to recover the spatial extent of epileptic spikes:
Grova C, Daunizeau J, Lina JM, Benar CG, Benali H, Gotman J (2006)
Evaluation of EEG localization methods using realistic simulations of interictal spikes, NeuroImage.
Chowdhury RA, Lina JM, Kobayashi E, Grova C (2013)
MEG Source Localization of spatially extended generators of epileptic activity: Comparing Entropic and Hierarchical Bayesian Approaches, PloS ONE.
Multimodal concordance on clinical cases:
Grova C, Daunizeau J, Kobayashi E, Bagshaw AP, Lina JM, Dubeau F, Gotman J (2008)
Concordance between distributed EEG source localization and simultaneous EEG-fMRI studies of epileptic spikes, NeuroImage.
Heers M, Hedrich T, An D, Dubeau F, Gotman J, Grova C, Kobayashi E (2014)
Spatial correlation of hemodynamic changes related to interictal epileptic discharges with electric and magnetic source imaging, Hum Brain Mapp.
Background readings: Epilepsy and EEG/MEG
Ebersole JS, Husain AM, Nordli DR
Current Practice of Clinical Electroencephalography
Atlas of Pediatric EEG
Anderson CT, Carlson CE, Li Z, Raghavan M
Magnetoencephalography in the preoperative evaluation for epilepsy surgery
Curr Neurol Neurosci, Rep. 2014 May;14(5):446.
Kharkar S, Knowlton R
Magnetoencephalography in the presurgical evaluation of epilepsy
Epilepsy Behav. 2015 May;46:19–26.
The following script from the Brainstorm distribution reproduces the analysis presented in this tutorial page: brainstorm3/toolbox/script/tutorial_epilepsy.m