ECoG/sEEG tutorial

Authors: Francois Tadel, Marcel Heers.

This tutorial introduces some concepts that are specific to the management of ECoG/sEEG recordings in the Brainstorm environment. It is based on a clinical case from the Epilepsy Center at the University Hospital of Freiburg, Germany.

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.

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.

Download and installation

• Requirements: You have already followed all the introduction tutorials and you have a working copy of Brainstorm installed on your computer.

• SPM: The coregistration between pre- and post-implantation MRI volumes requires the SPM12 toolbox. If you are running Brainstorm from a Matlab installation, SPM can be installed or configured as a plugin: menu Plugins > SPM12 > Install. If you are using the stand-alone compiled version of Brainstorm: all the needed SPM scripts have been compiled and included in this executable.

• Download the dataset:

  • Go to the Download page of this website, and get the file: sample_ecog.zip

  • Unzip it outside of any of the Brainstorm folders (program folder or database folder)

• Brainstorm:

  • Start Brainstorm (Matlab scripts or stand-alone version).

  • Select the menu File > Create new protocol. Name it "TutorialEcog" and select:
    "No, use individual anatomy",
    "No, use one channel file per acquisition run".

Import the anatomy

The entire dataset, anatomy and recordings, is organized following the BIDS specification. It could be loaded automatically with the menu Import BIDS dataset, as illustrated in the Resting/OMEGA tutorial. However, in this tutorial we will detail the step-by-step instructions to load the anatomy and ECoG recordings manually.

Pre-implantation MRI

The pre-implantation T1 MRI is available as a .nii file (sample_ecog/sub-ecog01/ses-postimp/anat). This volume was processed with FreeSurfer 7.1.1, and the output segmentation is available in the derivatives folder (sample_ecog/derivatives/freesurfer/sub-ecog01_ses-preimp). We don't need to import both to the Brainstorm database: in order to minimize the number of operations, we will import only the FreeSurfer folder.

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

• Right-click on the TutorialEcog folder > New subject > ecog01

  • Keep the default options you set for the protocol

• Right-click on the subject node > Import anatomy folder (auto): More information

  • Set the file format: FreeSurfer + Volume atlases

  • Select the folder: sample_ecog/derivatives/freesurfer/sub-ecog01_ses-preimp

• At the end, make sure the file "cortex_15000V" is selected (downsampled pial surface), otherwise double-click to select it as the default cortex surface.
  
  

• Rename the file MRI into T1pre, for better readability in the rest of the tutorial.
  
  

• When processing your own datasets: Instead of using FreeSurfer, you can decide to run the segmentation of the pre-implantation MRI directly from Brainstorm, using CAT12.

• For help with using the volume atlases, refer to the introduction tutorial Explore the anatomy.

Post-implantation MRI

• The pre-implantation MRI will be used as the anatomical reference for this subject. We will now import a second scan done after the sEEG/ECoG implantation, on which we can see the recording contacts. In this dataset, the post-implantation volume is another T1 MRI scan (contacts hyposignal appears in black), but it is more commonly a CT scan (contacts hypersignal appears in white).

• Right-click on the subject node > Import MRI:
  Select: sample_ecog/sub-ecog01/ses-postimp/anat/sub-ecog01_ses-postimp_T1w.nii.gz

• Do you want to apply the transformation to the MRI file? YES
  This will reorient the MRI in Brainstorm's standard orientation, so you can see the coronal/sagittal/axial views correctly oriented.

• How to register the new volume? SPM
  For more information about these options, see the section Volume coregistration. Note that you need to have the SPM12 toolbox installed on your computer and added to your Matlab path for using this registration method.

• Reslice the volume? NO
  This operation rewrites the volume with the orientation and resolution of the pre-implantation MRI, so that the two volumes can be overlaid in the MRI viewer. In this tutorial, we are going to place the ECoG/sEEG contacts manually on this post-implantation volume. For a better image quality, it is better to do it using the non-resliced volume. We do want the resliced volume, at least to check that pre- and post-implantation MRI volumes are correctly registered, but we can compute the resliced volume separately.

• Rename the new volume sub-ecog01_ses-postimp_T1w_spm into T1post_spm.
  
  

• Right-click on the T1post volume > Register with default MRI > Reslice / world coordinates.
  
  

• Right-click on the resliced T1post volume > Display > Overlay on default MRI
  
  

• The MRI viewer opens, showing the post-implantation volume as a colored layer on top of the previous volume. Adjust the transparency and amplitude threshold of this layer in the section "Data options" of the Surface tab, adjust its colormap with the popup menu of the figure. Use this display to validate that the coregistration of the two volume is correct, all the parts of the head must align well.
  
  

Generate skull surfaces

• Skull surfaces are useful for defining the position of ECOG contacts, and displaying them in 3D views. The contacts are projected on the inner skull surfaces, and then they are represented as cylinders oriented tangentially to the skull.

• Right-click on the subject folder > Generate BEM surfaces > Brainstorm. Keep the default options.

• When the computation is done, double-click on the the original head surface ("head mask") to select it as the default instead of the new BEM version, which is not estimated well because the volume was defaced.
  
  
  
  

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: "ECOG: Deltamed/Micromed/..."

  • Select the recordings: sample_ecog/sub-ecog01/ses-postimp/ieeg/..._ieeg.eeg

• The new "Link to raw file" entry lets you access directly the contents of the original file. The menu "Review raw file" does not actually copy any data to the database. More information.
  
  

Edit the channel types

At the previous step, we explicitly defined that we were importing ECoG recordings. However, this file contains both sEEG and ECoG recordings, therefore we need to modify manually the type of some channels.

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

• Select all the TA and TB channels, right-click > Set channel type: enter SEEG.
  
  

Edit the contacts positions

For 3D mapping and to generate epileptogenicity maps, we need accurate 3D positions for the contacts of the grids and sEEG electrodes. Placing the contacts requires a good understanding of the implantation scheme reported by the neurosurgeon, and some skills in reading MRI scans. This tutorial will now explain how to do this manual marking with the MRI Viewer from the post-implantation MRI, which shows clear hyposignal artifacts around the ECoG and sEEG contacts. We prefer using the non-resliced volume because it offers a clearer view of these artifacts.

If you have no interest in this or if you want to skip faster to the analysis of the signals, you can load the 3D positions from the BIDS dataset: read the advanced section Read BIDS electrodes.tsv.

ECoG grid: G

• Right-click on the channel file > MRI registration > Edit (MRI Viewer: T1post_spm).

• Adjust the brightness and contrast of the colormap in order to see the contacts artifacts correctly:
  • right-click on the figure and move your mouse up/down, or
  • left-click on the colormap and move your mouse up/down (brightness) or left/right (contrast), or

  • right-click on the figure > Colormap:Anatomy (figures below: contrast=100, brightness=100)

• For help with all the options iEEG tab, refer to the Epileptogenicity tutorial.
  Related tutorial: Edit implantation without recordings.

• In the iEEG tab, select "G", and select the information we know about this grid:

  • Type: ECoG

  • Number of contacts: 8 8 (8x8 grid, for a total of 64 contacts)

  • Contact spacing: 10 mm

  • Contact height: 1 mm
    (the grid is thinner, but giving it some thickness makes the figures easier to read)

  • Contact diameter: 4 mm
    (Ad-Tech contacts have 4mm diameters, with a 2.3mm exposure)

  • Wire width: Keep the default value, or adjust it later to your taste
• ECoG grids can be defined by pointing successively at their 4 corners. In this example, the contacts of the grid G are labelled sequentially: A1 A2 .. A8 B1 B2 .. B8 .. .. H1 H2 .. H8. You should point at the 4 corners starting by the first two ones, then finish rotating in the same direction: A1, A8, H8, H1.

• Navigate in the volume and position the cursor at the center of the artifact corresponding to contact G_A1 (look at the implantation schemes above, it is the more frontal and mesial in the grid), then click on button [#1] in the iEEG tab.
  If you can't find this point, you could enter manually the coordinates of the point: click on the button "Set current coordinates" at the bottom of the MRI Viewer window, and copy the MRI coordinates from the screen capture below.
  
  

• Zoom in/out with the buttons or the associated shortcuts (Ctrl+Scroll or +/-) and explore the volume with (Shift+)x/y/z. Additional display options are available in the popup menu for this figure. All the shortcuts are listed in this tutorial.

• Place the cursor at contact G_A8 and click on button [#2].
  
  

• Place the cursor at contact G_H8 and click on button [#3].
  
  

• Place the cursor at contact G_H1 and click on button [#4].
  
  

• After placing the last corner, it should start computing for a while. The positions of the 64 contacts are distributed in the polygon you defined, and then adjusted to the shape of the brain. To get the figure below, with all the contacts projected on all slices, select the checkbox "MIP: Functional" in the MRI Viewer.
  
  

• In this window, the contacts can be moved by simple drag-and-drop of the dots representing them. If you move them by mistake, there is unfortunately no undo button. But you can cancel your modifications in two ways: close the figure and discard the modifications, or select "G" in the iEEG tab and then use the menu Contacts > Use default positions, to recompute the initial positions from the 4 four corners of the ECoG grid.

• In the same menu Contacts, you can find additional projection options: Project on inner skull and Project on cortex. All these menus work well only when skull surfaces are available for the patient (see the section above that explains how to use the menu Generate BEM surfaces).

• From time to time, you can click on the menu Contacts > Save modifications in the iEEG tab. This would prevent you from losing your work if anything unexpected happens.

ECoG strips: IHA, IHB, IHC

IHA

• In the iEEG tab, select strip IHA and edit the options:

  • Type: ECOG-mid (for strips or grids that are not in aligned with the skull surface)

  • Number of contacts: 4

  • Contact spacing: 10 mm

  • Contact height: 1 mm

  • Contact diameter: 4 mm

• ECoG strips are defined only with two points, the first and the last contact. The intermediate contacts are uniformly distributed between the first and last point. In the convention used here, the first contact is always the deepest one.

• If you want to hide the contacts of previously edited devices: select the item in the list at the top of the iEEG tab, then click on the button Show/hide selected electrodes just above the list.

• Place the cursor on IHA1 (the deepest contact on the most posterior strip) and click [#1]. To jump directly to the coordinates of this point: click on the button Set current coordinates at the bottom of the MRI viewer, and report the MRI coordinates displayed in the screen capture below.
  
  

• Place the cursor on IHA4 and click [#2].
  
  

• Because you have selected the type ECOG-mid, the contacts are not be aligned on the skull or the cortex automatically. The contacts between the first and the last one are by default distributed uniformly on a straight line, which probably doesn't match their actual positions. Move slightly contacts IHA2 and IHA3 so that they appear in the middle of the artifacts.

IHB

• Repeat the same steps with IHB: Set it as ECOG-mid, contact spacing=10mm.

• Place the cursor on IHB1 (the deepest contact on the middle strip) and click [#1].
  
  

• Place the cursor on IHB4 and click [#2].
  
  

• Move the IHB2 and IHB3 contacts on the slices to center them on the artifacts.

IHC

• Repeat the same steps with IHC: Set it as ECOG-mid, contact spacing=10mm.

• Place the cursor on IHC1 (the deepest contact on the anterior strip) and click [#1].
  
  

• Place the cursor on IHC4 and click [#2].
  
  

• Move the IHC2 and IHC3 contacts on the slices to center them on the artifacts.

sEEG electrodes: TA, TB

TA

• In the iEEG tab, select electrode TA and edit the options:

  • Type: SEEG

  • Number of contacts: 10

  • Contact spacing: 10 mm

  • Contact length: 2.5 mm

  • Contact diameter: 1 mm

  • Electrode diameter: 0.9 mm

  • Electrode length: 120 mm

• The contacts actually have the same diameter as the electrode itself, but we use a slightly higher diameter for display, so that they render well in most 3D and 2D views.
• sEEG electrodes are defined with two points, the tip of the electrode (deepest contact) and a second point anywhere along the axis the electrode. One typical solution is to use the entry point in the skull as this second point, because the screws are always easy to identify. In this example, the electrodes are curved, and the entry point is not aligned with most of the contact, therefore we suggest to pick the last contact on the electrode instead. The intermediate contacts are placed along electrode, separated according to the "contact spacing" option, starting from the tip.

• Place the cursor on TA1 (the deepest contact on the more mesial probe) and click [Set tip]. To jump directly to the coordinates of this point: click on the button Set current coordinates at the bottom of the MRI viewer, and report the MRI coordinates displayed in the screen capture below.
  
  

• Place the cursor on TA10 (the last contact on this electrode) and click [Set skull entry].
  This is obviously not the location of the screw in the skull, but as the electrode is bent, this is a better point to approximate the contact locations on a straight line.
  
  

• The electrode is curved, therefore you need to move most of the intermediate contacts TA2-TA9 to more appropriate positions, at the center of the matching artifact. Unfortunately, Brainstorm does not handle well curved electrodes yet: if you move the contacts, the electrode remains represented as a straight line. This is a display issue only, which has no impact on any computation, but we hope to improve it at some point.

TB

• Repeat the same steps with TB. with the same configuration:

  • Type: SEEG

  • Number of contacts: 10

  • Contact spacing: 10 mm

  • Contact length: 2.5 mm

  • Contact diameter: 1 mm

  • Electrode diameter: 0.9 mm

  • Electrode length: 120 mm

• Place the cursor on TB1 (the deepest contact on the more mesial probe) and click [Set tip].
  
  

• Place the cursor on TB10 (the last contact on this electrode) and click [Set skull entry].
  
  

• Move the intermediate contacts TB2-TB9 to the centers of the artifacts.

• When you are done with marking all the electrodes, click on [Save] in the MRI viewer.
  
  

3D display

• Right-click on the channel file > Display sensors > Explore all the options available.
  For help with the display options in the iEEG tab: Epileptogenicity tutorial.
  You can add more anatomy elements to the figure with the button "Add a surface" at the top-right of the Surface tab. For more help: Display the anatomy.
  
  

Project on brain surface

You probably noticed that many the contacts of the ECoG grid appear under the surface of the cortex in these 3D figures. This is mostly due to the brain shift between the volumes we used to mark the contacts (post-implantation MRI) and the volume from which the cortex surface was reconstructed (pre-implantation MRI). One solution to compensate for this brain shift and produce more readable figures is to project the contacts on the pre-implantation surfaces.

• In the iEEG tab, select the grid G, then use menu Contacts > Project on inner skull.
  
  

• You could also test what you obtain by projecting the contacts on the inner skull. It produces similar results, but with a different method, projecting the contacts possibily a bit further away from the brain. Depending on the quality of the surfaces, the shape of the brain, and the position of the contacts, one solution may produce weird results, try the other one in that case.

Alternative: Read BIDS electrodes.tsv

A BIDS dataset may contain information about the localization of the iEEG contacts in different systems of coordinates. These files are documented in the in BIDS specification. In this dataset, we can find the following files in the folder sample_ecog/sub-ecog01/ses-postimp/ieeg. Files

• sub-ecog01_ses-postimp_space-IXI549Space_coordsystem.json: MNI space metadata

• sub-ecog01_ses-postimp_space-IXI549Space_electrodes.tsv: XYZ coordinates, MNI space

• sub-ecog01_ses-postimp_space-ScanRAS_coordsystem.json: Subject space metadata

• sub-ecog01_ses-postimp_space-ScanRAS_electrodes.tsv: XYZ coordinates, subject space

The tag space- indicates how the coordinates should be interpreted (list of possible values):

• IXI549Space: Refers to a MNI normalization computed by SPM12 (see Coordinate systems). You should be able to load these locations on any volume for which MNI coordinates are available, but with the approximations due to the volume-based normalization and possible errors due to the use of different types of normalization. In this particular dataset, the MNI coordinates were computed from Brainstorm with the SPM12 Segment non-linear normalization. When possible, always prefer importing coordinates in subject space.

• ScanRAS: Refers to the scanner coordinates (named world coordinates in Brainstorm) of a specific anatomical image (MRI or CT). The interpretation of the coordinates in the file electrodes.tsv requires the information available in the adjacent file coordsystem.json. Typically, this pair of files defines coordinates associated to a specific volume in the BIDS dataset.

For the subject space (space-ScanRAS), the contents of the json file is the following:

 "iEEGCoordinateSystem": "ScanRAS",
 "iEEGCoordinateUnits": "mm",
 "iEEGCoordinateSystemDescription": "Scanner-based RAS coordinates of the post-implantation T1 MRI",
 "iEEGCoordinateProcessingDescription": "SEEG contacts manually placed with Brainstorm on the post-implanation T1 MRI",
 "iEEGCoordinateProcessingReference": "https://neuroimage.usc.edu/brainstorm/Tutorials/ECoG",
 "IntendedFor": "sub-ecog01/ses-postimp/anat/sub-ecog01_ses-postimp_T1w.nii.gz"

The field IntendedFor points at the file sub-ecog01_ses-postimp_T1w.nii.gz. This means that the coordinates in the files sub-ecog01_ses-postimp_space-ScanRAS_electrodes.tsv refer to the post-implantation T1 MRI.

To add the coordinates in subject space (be careful: if you already marked the contacts manually, you would lose all your work - if you want to try this, duplicate the subject first):

• Right-click on the channel file > Add EEG positions > Import from file

• Select format "EEG: BIDS electrodes.tsv, subject space mm"
  (or adjust depending on the tsv file you select and the information defined in the json file)

• Select file sub-ecog01_ses-postimp_space-ScanRAS_electrodes.tsv
  
  

Note that the electrodes.tsv file does not include all the information we defined manually in the previous section (ECOG vs. ECOG-mid, contact diameter...), therefore for optimal display you might need to go through the configuration of each implanted device in the iEEG tab and edit the values as previously.

Display the recordings

ECoG/sEEG time series

• Right-click on any "Link to raw file" > ECOG+SEEG > Display time series.
  
  

• For instructions on how to use this viewer, please refer to these tutorials:

  • Review continuous recordings

  • EEG and epilepsy: Review EEG recordings

• Typical display configuration includes:

  • Set the page duration to 20s or 30s: Record tab > Page settings.

  • Change the display mode to "columns": Record tab > Toolbar > First button.

  • Change the montage (bipolar/monopolar): Record tab > Drop-down menu.

  • Set visualization filters: Filter tab.
  • Disable the automatic scaling when changing page: Button [AS] in the time series figure.
  • Adjust the amplitude scale: Right-click+move up/down, or buttons [^] and [v].

  • Set the resolution precisely: Right-click on the figure > Figure > Set axes resolution.

• Montage options include, for each group separately or all the channels together:

  • Bipolar 1: Difference of pairs of adjacent sensors (sparse), eg. A1-A2, A3-A4, A5-A6

  • Bipolar 2: Difference of pairs of adjacent sensors (exhaustive), eg. A1-A2, A2-A3, A3-A4

  • Local average reference: Average reference computed independently for each group

2D/3D topography

With the menus 2DLayout, 2D Electrodes and 3D Electrodes you can visualize the ECoG/sEEG values at the current time point. Move the time cursor and it will update the values represented on the contacts.

Interpolate on the anatomy

From the same ECOG+SEEG display menu, you can also project the recordings on the MRI or the cortex. You can edit the colormap for the displayed values like any other colormap.

MRI volume: All the voxels in the neighborhood of a contact are attributed the value associated with this contact. In the figure below, both "MIP: functional" and "MIP: anatomy" options are selected (maximum intensity projection along each axis = "glass brain" view).

Cortex surface: The values of the vertices are interpolated from the nearest SEEG contacts (the magnitude decreases with the distance to the contact).

Power spectrum

We recommend you always start your data analysis with a spectral evalution of the recordings, it may help you identify bad channels. This is described in tutorials Power spectrum and EEG and epilepsy.

• In Process1, select the continuous file.

• Run process Frequency > Power spectrum density (Welch): All file, Length=10s, ECOG,SEEG.

• Double-click on the PSD files to display them, change to a log-scale for the X axis, zoom in. None of the channels seem to have noise levels that are higher or lower than the others, all the channels are possibly good.
  
  

• If the spectrum plots are too noisy and difficult to read, you can recompute them with lower frequency resolution. Decrease the duration of the estimator window, possibly down to 1s depending on the sampling frequency. This will estimate the PSD by averaging the power of more time windows with less data bins, the result will be smoother.

Import epochs of interest

At this point of the analysis, we are still looking at the original files, no data was copied to the database. The montages are saved in the Brainstorm preferences, the bad channels and new events are saved in the links of the database but not reported to the original files. If you delete your protocol at this point, you would only lose the event marking and bad channel selection.

We are now going to import a segment a recordings around each of the three seizures marked in this file (10s before and 20s after the onset). This will make real copies of the data in the database, so we can run additional processes on them.

• Double-click on the on the "Link to raw file" to open the time series figure (ECOG+SEEG).

• In the Record tab: double-click on the event category seizure_begin and rename it seizure.

• Close the file viewer to save the modification to the events.

• Right-click on Link to raw file > Import in database:

  • Time window: All
  • Split in time blocks: Disabled

  • Use events: Enabled, select seizure

  • Epoch time: [-10000, +20000] ms (imports -10s to 20s around the event "seizure")

  • Remove DC offset: Disabled
  • Resample: Disabled

  • Create a separate folder for each event type: Disabled
    
    

• At the end, you should have three new folders SZ1/SZ2/SZ3, the same name as the original .TRC files, but without the tag "raw" on top. These new folders contain copies of the SEEG recordings, if you delete these folders from the database explorer, you lose the recordings they contain.

• The imported epochs are saved with a new timing: for the seizure onsets, the reference time t=0s is now the event "Onset", which has been removed from the list. You can still see the other marker Seizure, which should be slightly before 0s. For the baselines, the reference t=0 is the beginning of the segment.
  
  

Anatomical labelling

As illustrated in the introduction tutorial Explore the anatomy, all the anatomical atlases available in the subject (e.g. DKT, Desikan-Killiany, Destrieux) are available in the MRI Viewer. Select the current atlas from the figure popup menu, and the anatomical label under the cursor is displayed at the top-right corner. You can overlay all the labels on the MRI with the menu Show atlas.

We can use these anatomical parcellations, and also the surface parcellations of the cortex available as scouts, in order to label automatically all the ECOG/SEEG contacts.

If accurate MNI coordinates are expected, it is recommended to use a non-linear MNI normalization instead of the initial linear registration computed at the beginning of this tutorial, as illustrated in the SEEG epileptogenicity tutorial. In the anatomy view, right-click on the pre-implantation MRI (first file in the subject anatomy) > MNI normalization > Segment. This requires SPM12.

Switch to the functional view. Right-click on the channel file > iEEG atlas labels > Select all the available options: XYZ coordinates in various coordinate systems, volume parcellations, surface parcellations.

When using surface atlases, select preferencially the high-resolution surfaces (~300000V): a higher density of labelled points on the surface will increase the chances to label correctly the contacts. Selecting the type of surface ("white", "mid" or "cortex") depends on multiple factors, and we don't have enough experience yet to give clear recommendations. To help you decide, here is what these labels mean:

• White: Represents the interface between the white matter and the gray matter

• Cortex: Inflated from the white surface, represents the interface between the gray matter and the CSF.

• Mid: In the FreeSurfer case, this surface is computed by Brainstorm as the layer half-way between the white and the cortex surfaces. If this is available, this might be your best guess. Try it yourself and please share your results and comments on the user forum.

More options and details in the tutorial SEEG epileptogenicity maps.

Additional documentation

Related tutorials

• EEG/Epilepsy tutorial

• sEEG epileptogenicity maps

• EEG-Video synchronization

Forum

• Automated screen captures of electrodes: https://neuroimage.usc.edu/forums/t/34628/2

• Set colors for each contact: https://neuroimage.usc.edu/forums/t/35057