EEG and epilepsy

Authors: Francois Tadel, Elizabeth Bock, John C. Mosher.

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.

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.

References

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.
Neurology, 62:S2-8.

Standard montages recommended by the American Clinical Neurophysiology Society:

• American Clinical Neurophysiology Society Guidelines

• Guidelines for Standard Electrode Position Nomenclature #5

• A Proposal for Standard Montages to Be Used in Clinical EEG #6

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",

  • "Yes, use one channel file per subject".

Import the anatomy

• Switch to the "anatomy" view of the protocol.

• Right-click on the TutorialEpilepsy folder > New subject > sepi01

  • Leave the default options you set for the protocol

• Right-click on the subject node > Import anatomy folder:

  • Set the file format: "FreeSurfer folder"

  • Select the folder: sample_epilepsy/anatomy

  • Number of vertices of the cortex surface: 15000 (default value)

• In the MRI Viewer, click on "Click here to compute MNI transformation". It compute an affine transformation between the subject space and the MNI ICBM152 space, and sets automatically the anatomical landmarks AC/PC/IH.
  
  

• Set the three fiducial points (indicated here in MRI coordinates):
  • NAS: x=134, y=222, z=74
  • LPA: x=58, y=123, z=69
  • RPA: x=204, y=120, z=75

• Click on [Save] to validate the modifications and continue with the importation of the anatomy. At the end of the process, make sure that the file "cortex_15000V" is selected (downsampled pial surface, that will be used for the source estimation). If it is not, double-click on it 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 following tutorial: Warping default anatomy

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 file "Link to raw file" lets you access directly the contents of the EEG recordings
• The channel file "Deltamed channels" in the (Common files) folder contains the name of the channels, but not their positions. We need to overwrite this file and import manually the positions of the electrodes (either a standard cap or accurate positions digitized with a Polhemus device).

Prepare the channel file

• Import the positions of the electrodes:

  • Right-click the channel file in Common files > Add EEG positions > Import from file:
    

  • Select the file format: "EEG: ANT 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.
• 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
• 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. We need to redefine this manually to get correct groups of sensors, we want to have only real EEG electrodes in the "EEG" category and put everything that we are not going to use a different category.

• Right-click on the channel file > Edit channel file:

  • Note that the EOG, EMG and ECG channels have their type correctly set.
  • 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

  • In many cases you would have to edit manually the type of some channels: select one channel, click on the type and select something in the drop-down list, or select multiple channels > 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

• The channel file we imported contains generic electrodes positions; hence, it cannot be fully aligned with the head surface coming from the MRI. We need to register manually these electrodes positions with the subject anatomy.

• Right-click on the channel file > MRI registration > Edit... You will see several buttons across the toolbar. For instructions, hover your mouse over the buttons to reveal the tooltips.

• 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.
• Click on the button "Project electrodes on scalp surface", to ensure all the electrodes touch the skin surface.

• Click on "OK" and agree to save the modifications.
  
  

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 the "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.

  • In the Record tab, change the default duration that is reviewed to 10s.
    

  • In the figure, click on the [Flip +/-] button to have the negative values pointing up (convention used by most neurologists).
    

• Open a 2D Sensor cap map of the EEG sensor values as well:

  • 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 as well, 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, and disable the automatic positioning of the figures you the figure arrangement doesn't get lost when you open a new figure (Window layout menu at the top-right of the Brainstorm figure > None)

• 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 a balance between the amount of information to display and the ease of use.
  
  

Frequency filters

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
  
  

Time and amplitude resolution

The resolutions of the time and amplitude axes has 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, leave 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.

Recommendations for this dataset are:

• Time axis: 170 pixels/sec (~55 mm/sec)

• Amplitude: 15 μV/pixels (~45 μV/mm)
  
  

User setups

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. This is a good time to use the menu "User setups" to save this window configuration, so that you can 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.

Multiple montages

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). Brainstorm offers a flexible way of doing this.

• 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 as a new "User setup"
• If you don't see the "Longitudinal 3" menu, it is probably because you have using Brainstorm before these predefined montages were made available in the software 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"

  

Mark spikes

Detect heartbeats

When you have a clean ECG signal for your patient, you can detect automatically identify all the heartbeats in the recordings. Because heartbeats can cause sharp waves in some EEG traces, it can help 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 the file).

• 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 with Brainstorm 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 (importing a a text file containing the timing of the markers)

  • Select file sample_epilepsy/data/tutorial_spikes.txt

  • When prompted, enter the event name "SPIKE"

• Note that 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).

Manually marking

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 leave 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 key 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.

Pre-process recordings

Evaluation

Two of the typical pre-processing steps consist in getting rid of the contamination due to the power lines (50 Hz or 60Hz) and of the frequencies we are not interested in (a low-pass filter to remove the high-frequencies and a high-pass filter to remove the very slow components of the signals). Let's start with the 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.

Band-pass filter

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:

• Leave the "Link to raw file" selected in the Process1 list.

• Run the process "Pre-process > Band-pass filter", with the options:

  • Lower cutoff frequency: 0.5 Hz (high-pass filter)

  • Upper cutoff frequency: 80 Hz (low-pass filter)

  • Mirror signal before filtering: On, Sensor types: EEG
    
    

• Note that this new continuous files is saved in your Brainstorm database, while your 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.

Average reference

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 operation is saved in the same way as the SSP and ICA projectors: as 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.
  You can close this figure, and close the recordings.
  
  

Evaluate the modifications

• In Process1, select the filtered file 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 very 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 very low 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, we select a few components and remove all the others. For more information on ICA, please refer to the documentation of the EEGLAB software.

• We will illustrate now how to compute some projectors to remove the contaminations from the eye movements with the Infomax algorithm. To make it run fast on most computers, we will not process the full recordings but only a short segments with a lot of eye movements: the segment [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: A segment where the artifact is present (or all the recordings if you have time)

  • 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).

  • Frequency band: Ignore (0-0Hz), we want to perform the analysis on the broadband signal.

  • Number of ICA components: If this value 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 of the original signals to the ICA algorithm.
    In the current example, the algorithm doesn't converge if we use all the components (25 EEG channels), and we have to enter 24 to make it work.
    [TODO] This issue has to be investigated further.

  • Sensor types: EEG

  • Using existing SSP/ICA: Yes. In this case, using or not the existing re-referencing projector should not make any difference, but it would change the results if you have other projectors already computed.

  • Save averaged artifact: Not applicable (because we are not using events).

  • ICA method: Infomax (calls directly the function runica.m from EEGLAB).

• You can follow the convergence of the Infomax 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). You can notice that the first two components (IC1, IC2) look very correlated with the EOG signals, and their topographies are very typical of eye movements.
  
  

• 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 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.

Epoching

Now we want to extract all the spikes from this filtered file and import them in the database.

• Right-click on the imported and filtered file "Raw (0s,3600s) | high(0.5Hz)" > 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

• Uncheck "Remove DC offset": the DC offset has already been removed with the high-pass filter
  
  

• You should see a new folder "SPIKE" containing 58 epochs in your database.
  
  

• 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 figure open 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: Conditions"

• 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.
  
  

Averaging spikes

• 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).

• You can more views to your reviewing environment (see the tutorial Visual exploration):
  
  

• Explore this average in time with the left and right keyboard arrow.

• 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

Head model

• We are going to use a realistic head model, this requires to calculate some additional surfaces for the subject, to represent the inner skull and outer skull surface. Go to the "Anatomy" view, right-click on the subject > Generate BEM surfaces.

• Use 1922 vertices for each layer (default). Three new surface files are created at the end:
  
  

• Go back to the "Functional data" view, right-click on the channel file > Compute head model. Select the "cortex surface option" and the "OpenMEEG BEM" forward model. Leave all the OpenMEEG options to their defaults except for one: select the option "Use adjoint formulation". Depending on your computer, this may take a long time (10-60 minutes).
  
  

• 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).<<BR>If the OpenMEEG calculation crashes, please refer to the OpenMEEG tutorial.

• An 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 and numerical instabilities).

• If you can compute the model but then later the source maps you get don't make any sense, it is probably 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 after the source estimation 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 chose what is described as option 1b in the section "Epilepsy" in that tutorial: we are going to calculate the noise covariance matrix using 10 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 10s segment between 110s and 160s.

• Right-click on the filtered file "Raw|band" > Noise covariance > Compute from recordings: Baseline = [110,160]ms, leave the other options to the default values.
  
  

Inverse model

• Right-click on the OpenMEEG head model > Compute sources [2016]: dSPM, Unconstrained
  
  

• It creates a shared inversion kernel in the folder (Common files) and a source link for each file in the subject. If you are not familiar with these concepts, please refer to the Source estimation tutorial.

• Double-click on the sources file for the average to open it. Configure its display: surface smoothing, amplitude threshold, colormap. Make sure that all the visualization filters are turned OFF.
  
  

• To display the same information re-interpolated in the volume, right-click on the source file > Cortical activation > Display on MRI
  
  

• You can save movies or contact sheets of the source maps with the Snapshot popup menu
  
  

• Example: Video made with three figures: Download
  
  

• Example: Contact sheet
  
  

• The inverse solutions that you calculated are all shared for all the recordings in this subject, so you can go back to your individual epochs, or even the continuous file, and look at the source maps for these files.

Regions of interest

• The scouts are a very interesting feature for studying the spike spatial propagation. From the displays we generated previously, there is no clear propagation of the spike. Let's try to explore this in more details, at least to illustrate the method.
• Start by placing a small scout at the focus (~7 vertices)
• Place a few other scouts of a similar size at other places around the spike location, in places at which you could suspect a possible propagation.
• Plot the overlayed time series for these scouts. Try with the scout functions MEAN and MAX.
• There is no delay in the peaks for the different regions at -20ms, 0ms and 35ms This temporal exploration seems to confirm that there is no clear spatial propagation of the spike.

• If you need a reminder on how to manipulate scouts: see the Scouts tutorial
  
  

Z-score normalization

A good way to reveal better the source activity at the cortex level is to calculate a Z-score of the source maps with respect with a quiet baseline. We can use the same baseline as for the calculation of the noise covariance matrix.

• Import the block a recordings to use as the baseline. Right-click on the imported and processed file "Raw (0s,3600s) | high(0.5Hz)" > Import in database. Uncheck all the options and select the time window: [120, 130] secondes.
  
  

• Click on the Process2 tab
• In Files A: Drag the sources for the segment you just imported "Raw (120.00s,130.00s)"
• In Files B: Drag the sources for the averaged spike

• Alternatively, select the recordings and manually select the "Process sources" button on both sides
  
  

• Run the process: Standardize > Z-score (static). Double-click on the new file to display it
  
  

• The amplitude in the figure is now displayed without units, it represents a multiple of the standard deviation over the baseline. The interpretation of this new figure is not very different from the previous one, but it helps observing the things with a different angle. This process usually provides smoother maps, with less bias due to the depth of the sources. On the non-standardized minimum norm maps, the values of the deeper sources are usually lower than the ones at the surface. The Z-score helps showing deeper sources that have minimum norm values that are low, but higher than during the baseline.

Source analysis: Full head 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 of being constrained to the surface. More information about this method in the Volume source estimation tutorial.

Head model

• Right-click on the channel file > Compute head model, select "MRI volume" instead of "Cortex surface". This feature is now stable, but still tagged as experimental because not all the Brainstorm features are available for volume-based source models (eg. no regions of interest, no group analysis, etc). Use the same options as before for OpenMEEG.
  
  

Inverse model

• Right-click on the head model > Compute sources. Use all the default options.
  
  

• Right-click on the unconstrained sources > Cortical activation > Display on MRI (3D).

  • To move the slices in the 3D view: right-click + move your mouse in the slice direction.

  • Use all the sliders in the Surface tab to configure the display.
    
    

• Right-click on the same file > Cortical activation > Display on MRI (MRI Viewer)

  • To smooth the display across slices: right-click on the figure > Smooth sources > "3"

  • For a glass brain view: select the option "MIP: Functional" in the MRI Viewer
    
    

Moving dipoles

See tutorial: Computing and displaying dipoles

Note: you can also load dipoles calculated with other programs.

MEM inverse solution

This section is on a different page: ?Link

Additional documentation

Forum posts

• EEG reference: Reference of the gain matrix

• Registration: Match known EEG locations to the MRI image

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
NeuroImage, 30:1161-1170.

Scripting

To reproduce this entire analysis with a script, use the following processes:

• Import recordings > Import anatomy folder

  • Subject name: sepi01

  • Folder to import: sample_epilepsy/anatomy (format = FreeSurfer)

  • Enter the fiducials positions indicated at the beginning of the tutorial
  • Because of the selection of these points, it is usally better to run this first step in interactive mode.

• Import anatomy > Generate BEM surfaces

  • Set all the values to 1922 and the skull thickness to 4mm

• Import recordings > Create link to raw file

  • File to import: sample_epilepsy/data/tutorial_eeg.bin (format = Deltamed Coherence-Neurofile)

  • Uncheck the option "Align sensors using headpoints"

• Import recordings > Set channel file

  • File to import: sample_epilepsy/data/tutorial_electrodes.elc (format = EEG:ANT Xensor)

  • Check the option "Align sensors using headpoints"

• Import recordings > Project electrodes on scalp

• Import recordings > Set channels type

  • Channel types or names: SP1, SP2, RS, PHO, DELR, DELL, QR, QL

• Events > Import from file

  • Event file: sample_epilepsy/data/tutorial_spikes.txt: (format = Array of times)

  • Event name: SPIKE

• Import recordings > Import MEG/EEG: Time

  • Condition name: Empty
  • Time window: default (entire file)
  • Split recordings in time blocks: 0s

• Pre-process > Band-pass filter:

  • Lower cutoff frequency: 0.5 Hz

  • Upper cutoff frequency: 0
  • Sensor types: Empty

  • Check: Overwrite input file

• Sources > Compute noise covariance

  • Baseline: [120s - 130s]

  • Remove DC offset: Block by block
  • Ouptut: Full noise covariance matrix

• Import recordings > Import MEG/EEG: Events

  • Condition name: Empty

  • Event names: SPIKE

  • Time window: default (entire file)

  • Epoch time: [-101.6s, 300.8s]

  • Check all the additional options

• Average > Average files

  • By condition (subject average)
  • Arithmetic average
  • Uncheck: Keep all the events

• File > Save snapshot

  • Snapshot: Sensors/MRI registration
  • Sensor type: EEG
  • Orientation: Left
  • Comment: "Average spike"

• File > Save snapshot

  • Snapshot: Recordings time series
  • Sensor type: EEG

• File > Save snapshot

  • Snapshot: Recordings topography (contact sheet)
  • Sensor type: EEG
  • Contact sheet (start time, stop time): "-40 110"
  • Contact sheet (number of images): 16

• Sources > Compute head model

  • EEG Method: OpenMEEG BEM
  • Edit OpenMEEG options
    • BEM Layers: Select all (scalp, skull, brain)
    • Check: Use adjoint formulation
    • Check: Use adaptive integration

• Sources > Compute sources

  • wMNE
  • Sensor types: EEG
  • Kernel onle: Shared

• File > Save snapshot

  • Snapshot: Sources (one time)
  • Orientation: top
  • Time: 0s

The following script from the Brainstorm distribution reproduces the analysis presented in this tutorial page: brainstorm3/toolbox/script/tutorial_epilepsy.m