Differences between revisions 79 and 125 (spanning 46 versions)

Elekta-Neuromag recordings

This tutorial describes how to process continuous Elekta-Neuromag MEG recordings. It is based on median nerve stimulation acquired at the Massachusetts General Hospital in 2005 with a Neuromag Vectorview 306 system. The sample dataset contains the results for one subject for both arms: 18 minutes, 300 stimulations per arm.

This document shows what to do step by step, but does not really explain what is happening, the meaning of the different options or processes, the issues or bugs you can encounter, and does not provide an exhaustive description of the software features. Those topics are introduced in the basic tutorials based on CTF recordings; so make sure that you followed all those initial tutorials before going through this one.

The script file tutorial_neuromag.m in the brainstorm3/toolbox/script folder performs exactly the same tasks automatically, without any user interaction. Please have a look at this file if you plan to write scripts to process recordings in .fif format.

Contents

Download and installation
Import the anatomy
Access the recordings
Pre-processing
Epoching and averaging
Gradiometers & magnetometers
Source estimation
Scripting
Feedback

Download and installation

Requirements: You have already followed all the basic 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_neuromag.zip
Unzip it in a folder that is not in any of the Brainstorm folders (program folder or database folder). This is really important that you always keep your original data files in a separate folder: the program folder can be deleted when updating the software, and the contents of the database folder is supposed to be manipulated only by the program itself.
Start Brainstorm (Matlab scripts or stand-alone version)
Menu File > Create new protocol > TutorialNeuromag, with the following options:
- No, use individual anatomy,
- Yes, use one channel file per subject.
  (we have only one run, it's possible to share the channels positions and head models between different conditions, never select this option if you have more than one run)

Import the anatomy

Right-click on the TutorialNeuromag folder > New subject > Subject01 (leave all the defaults).
Right-click on the subject node > Import anatomy folder:
- Set the file format: "FreeSurfer folder"
- Select the folder: sample_neuromag/anatomy
- Number of vertices of the cortex surface: 15000 (default value)
Set the 6 required fiducial points (indicated in MRI coordinates):
- NAS: x=131, y=232, z=123
- LPA: x=48, y=136, z=74
- RPA: x=204, y=131, z=67
- AC: x=127, y=143, z=118
- PC: x=127, y=117, z=111
- IH: x=127, y=120, z=167 (anywhere on the midsagittal plane)
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.

Access the recordings

Link the FIF file to the database

The sample_neuromag package contains one run for one subject only.
Switch to the "functional data" view, the middle button in the toolbar above the database explorer.
Right click on the Subject01 > Review raw file:
- Select the file format: MEG/EEG Neuromag FIFF (*.fif)
- Pick the file: sample_neuromag/data/mind004_050924_median01_raw.fif
Select "Event channel" when asked about the events.
Then you have to select from which technical tracks you want to read the events.
In this dataset, the tracks of interest are:
- STI 001: Contains the onsets of the electric stimulation of the right arm.
- STI 002: Contains the onsets of the electric stimulation of the left arm.
- Those are not standard settings, they depend on your acquisition setup. In most cases, you would rather import "events" defined with the Neuromag acquisition software, which are saved in tracks STI 014, STI 101, STI 201...
Select STI 001 and STI 002. The full tracks are read, all the changes of values are detected from those two tracks and a list of events is created.
At the end, you're asked if you want to "Refine the registration with the head points". This operation improves the initial MRI/MEG registration by fitting the head points digitized at the MEG acquisition and the scalp surface. Answer YES.
Always make sure that the final registration looks good. The yellow surface represents the inner surface of the MEG helmet and the green dots represent the head points that were saved with the MEG recordings. In this case the registration looks acceptable. If you consider it is not the case, you can fix it manually: right-click on the channel file > MRI registration > Edit.
The files that appear in the folder (Common files) will be shared for all the conditions of Subject01. This mechanism allows to compute the forward and inverse models only once for all the conditions when there is only one run (= one channel file = one head position).

Review the recordings

Right-click on the "Link to raw file" to see what types of channels are available:
- MEG (all): Represents all the MEG sensors in the file (magnetometers and gradiometers)
- MEG (magnetometers): 101 magnetometers (sensor names ending with a "1")
- MEG (gradiometers): 202 planar gradiometers (sensor names ending with "2" or "3")
- ECG: Electro-cardiogram
- EOG: Electro-oculogram
- Stim: All the stimulation lines
Double-click on the file to display all the MEG sensors.
In the Record tab, look at the list of the available default display configurations for the Neuromag Vectorview MEG systems. If you would like to use different combinations of sensors, you can create your own combinations with the entry "Edit montages..."

Pre-processing

Evaluate the recordings

Drag and drop the "Link to raw file" into the Process1 list.
Select the process "Frequency > Power spectrum density", configure it as follows:
In all the processes, the types of sensors can be specified using the following convention:
- MEG: All the MEG sensors
- MEG MAG: All the magnetometers
- MEG GRAD: All the gradiometers
- MEG GRAD2: Only the first axial gradiometer on each chip (name ends with a "2")
- MEG GRAD3: Only the second axial gradiometer on each chip (name ends with a "3")
Double-click on the PSD file to display it.
This power spectrum shows that there is no MEG sensor that is clearly bad, but some sensors show a strong contamination from the power lines (60Hz and harmonics).

Remove: 60Hz and harmonics

Preliminary note: This process is very time and resource consuming. If you are not planning to analyze anything specifically around 60Hz, you may want to skip this section and and the long computation time associated with it. It is described here as a reference for users who would need it.
The standard procedure illustrated for removing the 60Hz power line noise illustrated in the CTF-based introduction tutorial would be simply to run the process "Pre-process > Notch filter" on the continuous file. However, this case is a bit more complicated.
- If you try to run this process without selecting the option "Process the entire file at once", you would get the following error. The file already contains some SSP projectors, created by one of the Neuromag progams. Using those projectors requires all the sensors, and running frequency filters requires all the time samples at once. As a consequence, the entire file must be read in one block in order to apply frequency filters on it.
- Alternatively, if you try with the option "Process the entire file at once", you would probably crash Matlab or freeze your computer because of a lack of memory. The file is already large and this process requires to create a lot of matrices of the same size in memory. Do not try this unless you have at least 24Gb of memory available on your computer.
- The solution is to re-write an intermediate file with the SSP projectors applied, that you will be able to process easily after.
Select the "Link to raw file" in Process1 and run process "Artifacts > Apply SSP and CTF compensation". This creates a new .fif file which is an exact copy of the original .fif file, but with the SSP projectors applied to the recordings values (instead of being just saved as linear operators on the side).
Clear the Process1 list and select the new continuous file "Raw | clean"
Run the process "Pre-process > Notch filter", configured as follows:
(this may take a while complete...)
To evaluate the results of this process, select the last file generated ("Raw | clean | notch") and run again the process "Frequency > Power spectrum density".
You should observe a significant decrease of the contributions of the removed frequencies (60Hz, 120Hz, 180Hz) compared with the original spectrum.

Heartbeats and eye blinks

Signal Space Projection (SSP) is a method for projecting the recordings away from stereotyped artifacts, such as eye blinks and heartbeats.

Double-click on the last continuous file you got to display all the MEG recordings.
("Raw | clean | notch" with the notch filter, or "Link to raw file" without)
Right-click on the link > EOG > Display time series, to look at the eye movements.
Right-click on the link > ECG > Display time series, to look at the heartbeats.
From the SSP menu in the Record tab, run the automatic detections:
- Detect eye blinks: Select channel EOG 062, event name "blink".
- Detect heartbeats: Select channel ECG 063, event name "cardiac".
- Remove simultaneous: To avoid capturing ocular artifacts in the cardiac SSP.
  Remove events named: "cardiac"
  When too close from events: "blink"
  Minimum delay between events: 250ms.
Review the traces of ECG/EOG channels and make sure the events detected make sense.
From the same menu, run the following processes, separately for GRAD and MAG:
- Compute SSP: Heartbeats: Event "cardiac", sensor types="MEG MAG".
  (select component #1)
- Compute SSP: Heartbeats: Event "cardiac", sensor types="MEG GRAD".
  (select component #1)
- Compute SSP: Eyeblinks: Event "blink", sensor types="MEG MAG", use existing SSP.
  (select component #1)
- Compute SSP: Eyeblinks: Event "blink", sensor types="MEG GRAD".
  (select component #1)
For more information regarding this SSP method, refer to this tutorial:
?Detect and remove artifacts.
At this stage, you should also review the entire file (pressing repeatedly F3) to detect any obvious bad segment in the recordings. If you detect any, select the time window to mark as bad, right click on the figure > Reject time segment.

Epoching and averaging

Import the recordings

Right-click on the raw link > Import in database.
The following figure appears, that asks how to import those recordings in the Brainstorm database.
- Time window: Time range of interest, keep all the time definition.
- Split: Useful to import continuous recordings without events. We do not need this here.
- Events selection: Check the "Use events" option, and select both STI 001 and STI 002. The values in the parenthesis represent the number of occurrences of each event.
- Epoch time: Time instants that will be extracted before an after each event, to create the epochs that will be saved in the database. Set it to [-100, +300] ms
- Use Signal Space Projections: Save and use the SSP vectors created by MaxFilter, MNE and Brainstorm. Keep this option selected.
- Remove DC Offset: Check this option, and select: Time range: [-100, 0] ms.
  For each epoch, this will: compute the average of each channel over the baseline (pre-stimulus interval: -100ms to 0ms), then subtract it from the channel at all the times in [-100,+300]ms.
- Resample recordings: Keep this unchecked
- Create new conditions for epochs: If selected, a new condition is created for each event type (here, it will create two conditions "STI 001" and "STI 002"). If not selected, all the epochs from all the selected events are saved in the same folder.
- Click on [Import] and wait.
Rename the condition "STI_001_5" and "STI_002_5" into respectively "Right" and "Left" (Right-click > Rename, or press F2, or click twice waiting a while before the two clicks).
Rename the groups of trials to Left and Right.
At this stage, you should review the trials one after the other (press F3 to jump to the next file), separately for the magnetometers and the gradiometers, to make sure that you don't have any bad trial that have been imported. If you find a bad trial: right-click on the file or on the figure > Reject trial.

Averaging

Drag and drop all the left and right trials to the Process1 tab.
Run the process Average > Average files > By condition:
Double-click on the Left and Right averages to display all the MEG sensors:

Stimulation artifact

We see clearly our evoked responses in those figures, but the stimulation artifact between 0ms and 4ms causes the automatic scaling of the figures to be too high, especially for the right stimulation. This effect would get a lot worse when displaying the sources for those files. To get a better display of these averages, we can cut these peaks at 0ms by replacing the values with linear interpolations.
Select the averages for the Left and Right condition in the Process1 tab.
Run the process Artifacts > Cut stimulation artifact: 0-4ms, overwrite
Now you can display the Left and Right averages correctly: Right-click on the file > MEG (All) > Display time series.
Display the 2D topographies for each average: Right-click on the file > MEG (All) > 2D Sensor cap (or press CTRL+T from the time series figures).
You can also add a low-pass filter at 100Hz for the visualization.
Everything looks ok, the time series figures show a clean signal with a very high SNR, and the topographies at the first response peak (about 22ms) show an activity on the left side for the Right condition, and on the right side for the Left condition.

Gradiometers & magnetometers

You can also display the MEG MAG (top) and MEG GRAD (bottom) sensors separately:
The topographies look very similar, because all the values are re-interpolated on a surface of virtual magnetometers. To display the real values for the different types of sensors, without this magnetic interpolation, use the menus "No magnetic interpolation".
Right-click on Left > MEG (mag) > No magnetic interpolation > 2D Sensor cap.
Repeat the same operation for "MEG (gradiometers)". You now have three options, because it is not possible to display at the same time the values for the two planar gradiometers:
- Grad 2: First gradiometer, sensor names ending with "2"
- Grad 3: Second gradiometer, sensor names ending with "3"
- Grad norm: Norm of the two gradiometers => sqrt(grad2^{2 + grad3}2)
The three options side-by-side at 22ms (grad 2, grad3, norm):

Source estimation

We need now to calculate a source model for these recordings, using a noise covariance matrix calculated from the pre-stimulation baselines. This process is not detailed a lot here because it is very similar to what is shown in the CTF-based introduction tutorials.

Head model

Right-click on the channel file > Compute head model.
Use the overlapping spheres model and keep all of the options at their default values.
For more information: ?Head model tutorial.

Noise covariance

Select all the Left and Right trials, right-click > Noise covariance > Compute from recordings.
Leave all the default options and click [OK].
For more information: ?Noise covariance tutorial.

Inverse model

Right-click on the headmodel > Compute sources.
Select wMNE and all the default options. This will run a source reconstruction grouping the information from the magnetometers and the gradiometers into the same model. This is the recommended procedure, but you can also estimate sources for the two types of sensors separately if needed.
A shared inversion kernel is created in the (Common files) folder, and a link node is now visible for each recording file, single trials and averages.
For more information: ?Source estimation tutorial.

Scouts

You can easily reproduce exactly the same results as the ones presented in the introduction tutorials. Place scouts to capture the activity in the primary and secondary somatosensory areas to track the processing of the electric stimulations in time, at the surface surface of the brain.

Scripting

The operations described in this tutorial can be reproduced from a Matlab script, available in the Brainstorm distribution: brainstorm3/toolbox/script/tutorial_neuromag.m

-  ⇤ ← Revision 79 as of 2014-06-27 15:21:16 → 
  Size: 17544
  Editor: agrippa
  Comment:
+   ← Revision 125 as of 2014-07-28 19:56:53 → ⇥
  Size: 21133
  Editor: agrippa
  Comment:
-Deletions are marked like this.
+Additions are marked like this.
 Line 1:
-= Import and process Neuromag raw recordings =
This tutorial describes how to process continuous Neuromag MEG recordings. It is based on median nerve stimulation acquired at the Massachusetts General Hospital in 2005 with a Neuromag Vectorview 306 system. The sample dataset contains the results for one subject for both arms.
+= Elekta-Neuromag recordings =
This tutorial describes how to process continuous Elekta-Neuromag MEG recordings. It is based on median nerve stimulation acquired at the Massachusetts General Hospital in 2005 with a Neuromag Vectorview 306 system. The sample dataset contains the results for one subject for both arms: 18 minutes, 300 stimulations per arm.
 Line 15:
- * Select the menu File > Create new protocol. Name it "'''!TutorialNeuromag'''" and select the options:
  * "'''No, use individual anatomy'''",
  * "'''Yes, use one channel file per subject'''".<<BR>> (we have only one run, it's possible to share the channels positions  and head models between different conditions, do not select this option  if you have more than one run per subject) <<BR>><<BR>> {{attachment:createProtocol.gif||height="239",width="222"}}
+ * Menu File > Create new protocol > '''!TutorialNeuromag''', with the following options:
  * '''No, use individual anatomy''',
  * '''Yes, use one channel file per subject'''.<<BR>> (we have only one run, it's possible to share the channels positions  and head models between different conditions, never select this option  if you have more than one run) <<BR>><<BR>> {{attachment:createProtocol.gif||height="372",width="347"}}
 Line 20:
- * Right-click on the !TutorialNeuromag folder > New subject > Subject01 (leave all the defaults)<<BR>><<BR>> {{attachment:createSubjectMenu.gif||height="165",width="343"}}  {{attachment:createSubject.gif||height="309",width="296"}}
+ * Right-click on the !TutorialNeuromag folder > New subject > Subject01 (leave all the defaults).<<BR>><<BR>> {{attachment:createSubjectMenu.gif||height="165",width="343"}}  {{attachment:createSubject.gif||height="309",width="296"}}
 Line 35:
+==== Link the FIF file to the database ====
-Line 41:
+Line 42:
- * You have then to select from which technical track you want to read the events. <<BR>><<BR>> {{attachment:eventsChannels.gif||height="165",width="375"}}
+ * Then you have to select from which technical tracks you want to read the events. <<BR>><<BR>> {{attachment:eventsChannels.gif||height="165",width="375"}}
-Line 46:
+Line 47:
- * Select "STI 001" and "STI 002". The full track is read, all the changes of values are detected from those two tracks and a list of events is created.
 * At the end, you're asked if you want to "Refine the registration with the head  points". This operation improves the initial MRI/MEG registration by fitting the head points digitized at the MEG acquisition and the scalp surface. Answer '''YES'''.<<BR>><<BR>> {{attachment:refine.gif||height="170",width="324"}}
 * The files that appear in the folder "''(Common files)''" will be  shared for all the conditions of Subject01. This mechanism allows to  compute the forward and inverse models only once for all the conditions  when there is only one run (= one channel file = one head position) :<<BR>><<BR>> {{attachment:treeAfterImport1.gif||height="116",width="210"}}
+ * Select '''STI 001''' and '''STI 002'''. The full tracks are read, all the changes of values are detected from those two tracks and a list of events is created.
 * At the end, you're asked if you want to "Refine the registration with the head  points". This operation improves the initial MRI/MEG registration by fitting the head points digitized at the MEG acquisition and the scalp surface. Answer '''YES'''.<<BR>><<BR>> {{attachment:refine.gif||height="170",width="324"}} <<BR>>
 * Always make sure that the final registration looks good. The  yellow surface represents the inner surface of the MEG helmet and the  green dots represent the head points that were saved with the MEG  recordings. In this case the registration looks acceptable. If  you consider it is not the case, you can fix it manually: right-click on the channel file > MRI registration > Edit. <<BR>><<BR>> {{attachment:checkRegistrationFig.gif||height="199",width="218"}}
 * The files that appear in the folder ''(Common files)'' will be  shared for all the conditions of Subject01. This mechanism allows to  compute the forward and inverse models only once for all the conditions  when there is only one run (= one channel file = one head position).<<BR>><<BR>> {{attachment:treeAfterImport1.gif||height="116",width="210"}}

==== Review the recordings ====
 * Right-click on the "Link to raw file" to see what types of channels are available:<<BR>><<BR>> {{attachment:listTypes.gif||height="341",width="332"}}
  * '''MEG (all)''': Represents all the MEG sensors in the file (magnetometers and gradiometers)
  * '''MEG (magnetometers)''': 101 magnetometers (sensor names ending with a "1")
  * '''MEG (gradiometers)''': 202 planar gradiometers (sensor names ending with "2" or "3")
  * '''ECG''': Electro-cardiogram
  * '''EOG''': Electro-oculogram
  * '''Stim''': All the stimulation lines
 * Double-click on the file to display all the MEG sensors.
 * In the Record tab, look at the list of the available default display configurations for the Neuromag Vectorview MEG systems. If you would like to use different combinations of sensors, you can create your own combinations with the entry "Edit montages..."<<BR>><<BR>> {{attachment:montages.gif||height="333",width="457"}}
-Line 51:
+Line 64:
-==== Evalue the recordings ====
+==== Evaluate the recordings ====
-Line 53:
+Line 66:
- * Select the process "Frequency > Power spectrum density", configure it as follows:<<BR>><<BR>> {{attachment:processPsd.gif||height="314",width="426"}}
 * Double-click on the PSD file to display it.
 * This power spectrum shows that there is no MEG sensor that is clearly bad. Some sensors show a strong contamination from the power lines (60Hz and harmonics).
+ * Select the process "'''Frequency > Power spectrum density'''", configure it as follows:<<BR>><<BR>> {{attachment:processPsd.gif||height="358",width="487"}}
 * In all the processes, the types of sensors can be specified using the following convention:
  * '''MEG''': All the MEG sensors
  * '''MEG MAG''': All the magnetometers
  * '''MEG GRAD''': All the gradiometers
  * '''MEG GRAD2''': Only the first axial gradiometer on each chip (name ends with a "2")
  * '''MEG GRAD3''': Only the second axial gradiometer on each chip  (name ends with a "3")
 * Double-click on the PSD file to display it.<<BR>><<BR>> {{attachment:processPsdView.gif||height="179",width="383"}}
 * This power spectrum shows that there is no MEG sensor that is clearly bad, but some sensors show a strong contamination from the power lines (60Hz and harmonics).
-Line 57:
+Line 76:
-== Importing MEG recordings ==
 * The following figure appears, that asks how to import those recordings in the Brainstorm database.<<BR>><<BR>> {{attachment:importOptions.gif}}
  * '''Time window''': Time range of interest. Now we want to keep the whole time definition, we are interested by all the stimulations.
+==== Remove: 60Hz and harmonics ====
 * __Preliminary note__: This process is very time and resource consuming. If you are not planning to analyze anything specifically around 60Hz, you may want to skip this section and and the long computation time associated with it. It is described here as a reference for users who would need it.
 * The standard procedure illustrated for removing the 60Hz power line noise illustrated in the CTF-based introduction tutorial would be simply to run the process "Pre-process > Notch filter"  on the continuous file. However, this case is a bit more complicated.
  * If you try to run this process without selecting the option "Process the entire file at once", you would get the following error. The file already contains some SSP projectors, created by one of the Neuromag progams. Using those projectors requires all the sensors, and running frequency filters requires all the time samples at once. As a consequence, the entire file must be read in one block in order to apply frequency filters on it. <<BR>><<BR>> {{attachment:sinError.gif||height="234",width="658"}}
  * Alternatively, if you try with the  option "Process the entire file at once", you would probably crash Matlab or freeze your computer because of a lack of memory. The file is already large and this process requires to create a lot of matrices of the same size in memory. Do not try this unless you have at least 24Gb of memory available on your computer.
  * The solution is to re-write an intermediate file with the SSP projectors applied, that you will be able to process easily after.
 * Select the "Link to raw file" in Process1 and run process "'''Artifacts > Apply SSP and CTF compensation'''". This creates a new .fif file which is an exact copy of the original .fif file, but with the SSP projectors applied to the recordings values (instead of being just saved as linear operators on the side).<<BR>><<BR>> {{attachment:processApplySsp.gif||height="145",width="270"}}
 * Clear the Process1 list and select the new continuous file "'''Raw | clean'''"
 * Run the process  "'''Pre-process > Notch filter'''", configured as follows:<<BR>>(this may take a while complete...)<<BR>><<BR>> {{attachment:processSin.gif||height="274",width="336"}}
 * To evaluate the results of this process, select the last file generated ('''"Raw | clean | notch"''') and run again the process "'''Frequency > Power spectrum density'''".<<BR>><<BR>> {{attachment:processPsd2.gif||height="351",width="396"}}
 * You should observe a significant decrease of the contributions of the removed frequencies (60Hz, 120Hz, 180Hz) compared with the original spectrum. <<BR>><<BR>> {{attachment:processPsdView2.gif||height="202",width="425"}}

==== Heartbeats and eye blinks ====
Signal Space Projection (SSP) is a method for projecting the recordings away from stereotyped artifacts, such as eye blinks and heartbeats.

 * Double-click on the last continuous file you got to display all the '''MEG '''recordings.<<BR>>("Raw | clean | notch" with the notch filter, or "Link to raw file" without)
 * Right-click on the link > '''EOG '''> Display time series, to look at the eye movements.
 * Right-click on the link > '''ECG '''> Display time series, to look at the heartbeats.
 * From the SSP menu in the Record tab, run the automatic detections:
  * '''Detect eye blinks:''' Select channel '''EOG 062''', event name "blink".<<BR>><<BR>> {{attachment:processDetectEog.gif||height="213",width="309"}}
  * '''Detect  heartbeats:''' Select channel '''ECG 063''', event name "cardiac".<<BR>><<BR>> {{attachment:processDetectEcg.gif||height="214",width="310"}}
  * '''Remove simultaneous:''' To avoid capturing ocular artifacts in the cardiac SSP.<<BR>>Remove events named: "'''cardiac'''"<<BR>>When too close from events: "'''blink'''"<<BR>>Minimum delay between events: '''250ms'''.<<BR>><<BR>> {{attachment:processRemoveSimult.gif||height="234",width="309"}}
 * Review the traces of ECG/EOG channels and make sure the events detected make sense.<<BR>><<BR>> {{attachment:processDetectAfter.gif||height="240",width="229"}}
 * From the same menu, run the following processes, separately for GRAD and MAG:
  * '''Compute SSP: Heartbeats''': Event "cardiac", sensor types="'''MEG MAG'''".<<BR>>(select component #1) <<BR>><<BR>> {{attachment:sspEcgMag.gif||height="218",width="560"}}
  * '''Compute SSP: Heartbeats''': Event "cardiac", sensor types="'''MEG GRAD'''".<<BR>>(select component #1) <<BR>><<BR>> {{attachment:sspEcgGrad.gif||height="219",width="561"}}
  * '''Compute SSP: Eyeblinks''': Event "blink", sensor types="'''MEG MAG'''", use existing SSP.<<BR>>(select component #1) <<BR>><<BR>> {{attachment:sspEogMag.gif||height="218",width="560"}}
  * '''Compute SSP: Eyeblinks''': Event "blink", sensor types="'''MEG GRAD'''".<<BR>>(select component #1)<<BR>><<BR>> {{attachment:sspEogGrad.gif||height="218",width="560"}}
 * For more information regarding this SSP method, refer to this tutorial:<<BR>>[[Tutorials/TutRawSsp|Detect and remove artifacts]].
 * At this stage, you should also review the entire file (pressing repeatedly F3) to detect any obvious bad segment in the recordings. If you detect any, select the time window to mark as bad, right click on the figure > Reject time segment.

== Epoching and averaging ==
==== Import the recordings ====
 * Right-click on the raw link > Import in database.<<BR>><<BR>> {{attachment:importPopup.gif||height="298",width="193"}}
 * The following figure appears, that asks how to import those recordings in the Brainstorm database.<<BR>><<BR>> {{attachment:importOptions.gif||height="364",width="580"}}
  * '''Time window''': Time range of interest, keep all the time definition.
-Line 61:
+Line 113:
-  * '''Events selection''': Check the "Use events" option, and select both STI001 and STI002. The value in the parenthesis represents the number of occurrences of this event.
+  * '''Events selection''': Check the "Use events" option, and select both '''STI 001 '''and '''STI 002'''. The values in the parenthesis represent the number of occurrences of each event.
-Line 63:
+Line 115:
-  * '''Use Signal Space Projections''': Save and use the SSP vectors created by !MaxFilter or during other pre-processing steps. Keep this option selected.
  * '''Remove DC Offset''': Check this option, and select: Time range: [-100, 0] ms. <<BR>>For each epoch, this will: compute the average of each channel over the baseline (pre-stimulus interval: -100ms to 0ms), and subtract it from the channel at all the times in [-100,+300]ms.
+  * '''Use Signal Space Projections''': Save and use the SSP vectors created by !MaxFilter, MNE and Brainstorm. Keep this option selected.
  * '''Remove DC Offset''': Check this option, and select: Time range: [-100, 0] ms. <<BR>>For each epoch, this will: compute the average of each channel over the baseline (pre-stimulus interval: -100ms to 0ms), then subtract it from the channel at all the times in [-100,+300]ms.
-Line 66:
+Line 118:
-  * '''Create new conditions for epochs''': If selected, a new condition is created for each event type (here, it will create a condition "Event_5"). If not selected, all the epochs from all the selected events are saved in the same condition, the one that was selected in the database explorer (if no condition selected, create a condition called "New condition").
  * Click on Import and wait.
 * Rename the condition "STI_001_5" and "STI_002_5" into respectively "Right" and "Left" (Right-click > Rename, or press F2, or click twice waiting a while before the two clicks).<<BR>><<BR>>
+  * '''Create new conditions for epochs''': If selected, a new condition is created for each event type (here, it will create two conditions "STI 001" and "STI 002"). If not selected, all the epochs from all the selected events are saved in the same folder.
  * Click on [Import] and wait.
 * Rename the condition "'''STI_001_5'''" and "'''STI_002''''''_5'''" into respectively "'''Right'''" and "'''Left'''" (Right-click > Rename, or press F2, or click twice waiting a while before the two clicks).
 * Rename the groups of trials to '''Left''' and '''Right'''.<<BR>><<BR>> {{attachment:rename.gif||height="105",width="225"}}
 * At this stage, you should review the trials one after the other (press F3 to jump to the next file), separately for the magnetometers and the gradiometers, to make sure that you don't have any bad trial that have been imported. If you find a bad trial: right-click on the file or on the figure > Reject trial.
-Line 70:
+Line 124:
-== Pre-processing ==
=== Select the files to process ===
We want to apply the same pre-processing operations to all the trials of both conditions. Drag and drop the "Subject01" node from the database explorer to the tab ''Process1''. Select on the "Recordings" button, then click ''Run''.
+==== Averaging ====
 * Drag and drop all the left and right trials to the Process1 tab.
 * Run the process '''Average > Average files > By condition''': <<BR>><<BR>> {{attachment:processAverage.gif||height="452",width="451"}}
 * Double-click on the Left and Right averages to display all the MEG sensors:<<BR>><<BR>> {{attachment:averageArtifact.gif||height="196",width="480"}}
-Line 74:
+Line 129:
- . {{attachment:cutStimList.gif}}

=== Create a processing pipeline ===
 * '''Process > Cut stimulation artifact''': The electric stimulation of the median nerve induces a strong artifact  right after 0ms. This process uses a simple trick to remove this  artifact: re-interpolate the values for the artifact time window (linear  interpolation). It doesn't affect much the data but will make all the  displays much better. Set the time window to [0, 3.9]ms. Select "Overwrite initial file".

 * '''Process > Band-pass filter''': Set the frequency range to [2Hz -90Hz]. Select "Overwrite initial file".

 * '''Detect > Detect bad channels: peak-to-peak''': This process tests the peak-to-peak amplitude of each channel, and defines it is good or bad based on a threshold value that is defined individually for each channel type.Set the following options:
  * '''Time window''': All the epoch (keep the default value)

  * '''Meg gradio''': [0, 3500]

  * '''Meg magneto''': [0, 3500]

  * '''Reject entire trial'''
 * '''Average > Average files''': Select "By condition".<<BR>><<BR>> {{attachment:pipeline.gif}}

=== Run the pipeline ===
Click on run. Wait. With all those process tags added, the comments of the files are getting a little too long. Rename them respectively Right and Left, by renaming directly the list (the node that contains all the epochs in each condition). It takes a while but will improve a lot the readability later.

 . {{attachment:treeAfterPreprocess.gif}} {{attachment:treeAfterPreprocess2.gif}}

Expand the list of trials Left, and notice that the icons of a few trials have a red flag.

 * Those are the trials that were identified as bad by the process "Detect bad channels".
 * Now if you clear the Process1 list and drag and drop again Subject01, you would see the following: there are now only 598 recording files detected, instead of 624 previously. This number is: number total of trials (624) - number of bad trials (28) + the averages (2).
 * On the right side of the panel, a list of checkboxes appeared. They represent the list of tags that were found in the collection of files in Subject01. Tags "bandpass" and "cutstim" are set for all the trials, and the "none" check box represents the files that do not have tags (the two averages).
 * If you uncheck "none", it removes the 2 averaged files from the selection. If you uncheck both "bandpass" and "cutstim", it removes all the pre-processed trials from the selection.

 . {{attachment:treeBad.gif}} {{attachment:treeBadProcess.gif}}

=== Review the epochs manually ===
It is always very important to keep an eye on the quality of the data at the different steps of the analysis. There are always a few epochs that are too artifacted or noisy to be used, or one bad sensor. An automatic detection was already applied, but here is the procedure if you need to do it manually.

 * Double-click on the Right epochs list to expand the list of files contained in it. Then display the time series for the first Right epoch by double-clicking on it.<<BR>><<BR>> {{attachment:reviewEpochs.gif}}
 * Then press F3 to go to the next epoch (or menu "Navigator > Next data file" in the main Brainstorm window). Do that until you reach the last epoch. For each epoch:
  * Check that the amplitude range in consistant with the previous epochs (in this study: between 1000fT and 2500fT).
  * '''Tag channels as bad''': select a sensor by clicking on it and press the Delete key (or right-click > Channel > Mark selected as bad)
  * '''Tag the entire trial as bad''': You can mark a trial as bad, it will be completely ignored from all the following computations. Right-click 1) on the trial in the explorer or 2) in the time series figure, and select menu ''Reject trial''. Or use the keyboard shortcut ''Ctrl+B''.
 * In the case of this study, there are so many trials for each condition and the recorded signal is so strong that you do not really need to spend too much time on the selection of the bad channels and bad epochs. If you just keep everything as it is, it is just fine.

=== Review the average ===
 * Display the time series for each average: Right-click on the file > MEG (All) > Display time series.
 * Display the 2D topographies for each average: Right-click on the file > MEG (All) > 2D Sensor cap (or press CTRL+T from the time series figures)<<BR>><<BR>> {{attachment:reviewAverages.gif}}
+==== Stimulation artifact ====
 * We see clearly our evoked responses in those figures, but the stimulation artifact between 0ms and 4ms causes the automatic scaling of the figures to be too high, especially for the right stimulation. This effect would get a lot worse when displaying the sources for those files. To get a better display of these averages, we can cut these peaks at 0ms by replacing the values with linear interpolations.
 * Select the averages for the Left and Right condition in the Process1 tab.
 * Run the process '''Artifacts > Cut stimulation artifact''': 0-4ms, overwrite<<BR>><<BR>> {{attachment:processCut.gif||height="239",width="471"}}
 * Now you can display the Left and Right averages correctly: Right-click on the file > MEG (All) > Display time series.
 * Display the 2D topographies for each average: Right-click on the file > MEG (All) > 2D Sensor cap (or press CTRL+T from the time series figures).
 * You can also add a low-pass filter at 100Hz for the visualization.<<BR>><<BR>> {{attachment:reviewAverages.gif||height="378",width="791"}}
-Line 120:
+Line 138:
-== Forward model ==
First step of the source estimation process: establishing a model that describes the way the brain electric activities influence the magnetic fields that are recorded by the MEG sensors. This model can be designated in the documentation by the following terms: head model, forward model, lead field matrix.

=== MEG / MRI registration ===
An accurate forward model requires first of all an accurate registration of the anatomy files (MRI+surfaces) and functional recordings (position of the MEG sensors and EEG electrodes). A basic registration is provided by the alignment of the fiducials (Nasion, LPA, RPA), that were both located before the acquisition of the recordings and marked on the MRI in Brainstorm. This registration based on three points only can be very inaccurate, because it is sometimes hard to identify clearly those points, and not everybody identify them the same way. Two methods described in the [[Tutorials/TutImportRecordings|introduction tutorial #3]] may help you to improve this registration.

 * Right-click on the channel file > MRI Registration > Check<<BR>><<BR>> {{attachment:checkRegistration.gif}}  {{attachment:checkRegistrationFig.gif}}
 * The yellow surface represents the inner surface of the MEG helmet, and the green dots represent the head points that where saved with the MEG recordings.
 * In this case the registration looks acceptable. If you consider it is not your case, you can try the two other menus: "Edit..." and "Refine using head points"

=== Compute head model ===
Right-click on any node that contains the channel file (including the channel file itself), and select: "Compute head model". When the computation is done, close the "Check spheres" figure. The lead  field matrix is saved in file "Overlapping spheres" in "Common files",  in a the "Gain" field.

 . {{attachment:headmodelStart.gif}}

== Noise covariance matrix ==
To estimate the sources properly, we need an estimation of the noise level for each sensor. A good way to do this is to compute the covariance matrix of the concatenation of the baselines from all the trials in both conditions.

 * Select at the same time the two groups of trials (Right and Left). To do this: hold the Control key and click successively on the Right and the Left trial lists.
 * Right-click on one of them and select: Noise covariance > Compute from recordings. Leave all the options unchanged, and click on Ok.
 * This operation computes the noise covariance matrix based on the  baseline, [-100, 0]ms, of all the __good__ trials (596 files). The result is stored in a new file in the  "Common files" folder.

 . {{attachment:noisecovMenu.gif}} {{attachment:noisecovOptions.gif}}
+== Gradiometers & magnetometers ==
 * You can also display the MEG MAG (top) and MEG GRAD (bottom) sensors separately: <<BR>><<BR>> {{attachment:averageMagGrad.gif||height="319",width="467"}}
 * The topographies look very similar, because all the values are re-interpolated on a surface of virtual magnetometers. To display the real values for the different types of sensors, without this magnetic interpolation, use the menus "No magnetic interpolation".
 * Right-click on Left > MEG (mag) > No magnetic interpolation > 2D Sensor cap. <<BR>><<BR>> {{attachment:averageMagNointerp.gif||height="196",width="235"}}
 * Repeat the same operation for "MEG (gradiometers)". You now have three options, because it is not possible to display at the same time the values for the two planar gradiometers:
  * Grad 2: First gradiometer, sensor names ending with "2"
  * Grad 3: Second gradiometer, sensor names ending with "3"
  * Grad norm: Norm of the two gradiometers => sqrt(grad2^2 + grad3^2)
 * The three options side-by-side at 22ms (grad 2, grad3, norm): <<BR>><<BR>> {{attachment:averageGradNointerp2.gif||height="165",width="542"}}
-Line 145:
+Line 149:
-Reconstruction of the cortical sources with a weighted minimum norm estimator (wMNE).
+We need now to calculate a source model for these recordings, using a noise covariance matrix calculated from the pre-stimulation baselines. This process is not detailed a lot here because it is very similar to what is shown in the CTF-based introduction tutorials.
-Line 147:
+Line 151:
- * Right-click on "Common files", (or on the head model, or on the subject node, this is all equivalent), and select "Compute sources". <<BR>><<BR>> {{attachment:inverseMenu.gif}}
 * Select "Yes" to ignore the warning about noise covariance used at the same time for averages and single trials (if you do not see such a warning, that's fine too). In this case, this warning does not really apply, the results are the same if we compute separately the minimum norm estimator for the single trials and the averages.
 * Select "Normal mode" (not "Expert mode"), wMNE, and both gradiometers and magnetometers, as in the figure below. Click on Run.<<BR>><<BR>> {{attachment:inverseOptions.gif}} {{attachment:treeAfterInverse.gif}}
 * A shared inversion kernel was created in the ''(Common files)'' folder, and a link node is now visible for each recording file, single trials and averages. For more information about what those links mean and the operations performed to display them, please refer to the [[Tutorials/TutSourceEstimation|tutorial #4 "Source estimation"]].
 * For the average of condition Right: display the time series, the 2D sensor topography, the sources on the cortex, the sources on 3D orthogonal slices, and the sources in the MRI Viewer.<<BR>><<BR>> {{attachment:inverseAllDisplay.gif}}
 * Create a scout called LeftSS that covers the most active region of the cortex at 33ms for condition Right. Display its time series for conditions Left and Right at the same time (select both Left and Right source files, right click > Cortical activations > Scouts time series (check "Overlay: conditions").<<BR>><<BR>> {{attachment:inverseScoutLeft.gif}}
+==== Head model ====
 * Right-click on the channel file > Compute head model.<<BR>><<BR>> {{attachment:headmodelMenu.gif||height="196",width="228"}}
 * Use the overlapping spheres model and keep all of the options at their default values.<<BR>><<BR>> {{attachment:headmodelOptions.gif||height="247",width="222"}}
 * For more information: [[Tutorials/TutHeadModel|Head model tutorial]].

==== Noise covariance ====
 * Select all the Left and Right trials, right-click > Noise covariance > Compute from recordings.<<BR>><<BR>> {{attachment:noisecovPopup.gif||height="188",width="420"}}
 * Leave all the default options and click [OK]. <<BR>><<BR>> {{attachment:noisecovOptions.gif||height="374",width="382"}}
 * For more information: [[Tutorials/TutNoiseCov|Noise covariance tutorial]].

==== Inverse model ====
 * Right-click on the headmodel > Compute sources.<<BR>><<BR>> {{attachment:inversePopup.gif||height="131",width="224"}}
 * Select wMNE and all the default options. This will run a source reconstruction grouping the information from the magnetometers and the gradiometers into the same model. This is the recommended procedure, but you can also estimate sources for the two types of sensors separately if needed. <<BR>><<BR>> {{attachment:inverseOptions.gif||height="233",width="207"}}
 * A shared inversion kernel is created in the ''(Common files)'' folder, and a link node is now visible for each recording file, single trials and averages. <<BR>><<BR>> {{attachment:inverseFiles.gif||height="239",width="229"}}
 * For more information: [[Tutorials/TutSourceEstimation|Source estimation tutorial]].

==== Scouts ====
You can easily reproduce exactly the same results as the ones presented in the introduction tutorials. Place scouts to capture the activity in the primary and secondary somatosensory areas to track the processing of the electric stimulations in time, at the surface surface of the brain.

{{attachment:final.gif||height="419",width="842"}}

== Scripting ==
The operations described in this tutorial can be reproduced from a Matlab script, available in the Brainstorm distribution: '''brainstorm3/toolbox/script/tutorial_neuromag.m '''
-Line 155:
+Line 176:
-<<EmbedContent(http://neuroimage.usc.edu/brainstorm3_register/get_feedback.php?Tutorials/TutMindNeuromag)>>
+<<EmbedContent(http://neuroimage.usc.edu/bst/get_feedback.php?Tutorials/TutMindNeuromag)>>

Software

Users

Development