|
Size: 16925
Comment:
|
← Revision 128 as of 2022-08-22 14:15:40 ⇥
Size: 0
Comment:
|
| Deletions are marked like this. | Additions are marked like this. |
| Line 1: | Line 1: |
| = Review continuous recordings and edit markers = This tutorial describes how to review a continuous file and add time markers. It is based on a median nerve stimulation experiment recorded at the Montreal Neurological Institute in 2011 with a CTF MEG 275 system. The sample dataset contains 6 minutes of recordings at 1200Hz for one subject and includes 100 stimulations of each arm. The segmentation of the T1 MRI of the subject was performed using [[Tutorials/LabelFreeSurfer|FreeSurfer]]. <<TableOfContents(2,2)>> == 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 [[http://neuroimage.usc.edu/brainstorm3_register/download.php|Download]] page of this website, and dowload the file: '''sample_raw.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) * Select the menu File > Create new protocol. Name it "!TutorialRaw" and select the options: "No, use individual anatomy", "Yes, use one channel file per subject". == Import the anatomy == * Create a new subject Subject01 * In the anatomy view, right-click on the subject > Import !FreeSurfer folder, and select the folder "sample_raw/Anatomy". This imports automatically all the available files in a segmentation folder generated with !FreeSurfer, as illustrated in [[Tutorials/LabelFreeSurfer|this page]]. * Leave the default number of vertices of the cortex surface (15000) * When asked to place the fiducials, you can use the file "sample_raw/Anatomy/fiducials.jpg". <<BR>> || {{attachment:fiducials.jpg}} ||The picture shows the three head position coils used in the CTF systems. The positions of the fiducials were digitized with a Polhemus device at the middle of each of those coils. In the end, you should have the following MRI coordinates (+/- a few millimeters): <<BR>><<BR>>NAS: x=126, y=215, z=138<<BR>>LPA: x=58, y=136, z=121<<BR>>RPA: x=197, y=138, z=118<<BR>>AC: x=128, y=137, z=157<<BR>>PC: x=128, y=113, z=157<<BR>>IH: x=128, y=125, z=216 (anywhere on the midsagittal plane) || * Click on Save when you're done, and wait until the process is over. In the end, a figure shows the surfaces that were generated. Close it.<<BR>><<BR>> {{attachment:anat.gif|fiducials.jpg}} * Four files are now available as the anatomy of Subject01: a T1 MRI, a head surface, and two cortex surfaces (high and low resolution). The anatomy of the subject is ready. == Access the raw file == The basic tutorials you read before explain how to import recordings in the database: this operation creates a copy of all the data in Matlab .mat files in the Brainstorm database folders. You could process continuous recordings in the same way, but the .mat format has this limitation that the entire file has to be read even when you want to access just a portion of it. Long recordings usually cannot fit in memory and have to be split in small blocks of a few seconds, which makes it very difficult to review and process. Brainstorm offers the possibility to visualize continuous MEG/EEG recordings in any of the supported file formats without having to fully "import" them. A link to the native file is created in the database, which can be then manipulated almost like the "imported" recording blocks. Only the description of the file is saved in the database, and when displaying it the values are read directly from the native file. In addition, an interface allows to edit the time markers that are saved in the file. Those markers can then be used to import the recordings in the database (ie. to do the segmentation of the continous recordings in epochs/trials). Then the imported epochs/trials (hard copies in .mat format) can be pre-processed and averaged. * Select the "Functional data (sorted by subject)" exploration mode (second button in the toolbar above the database explorer). * Right-click on the subject node, and select: "Review raw file". Select the "MEG: CTF" file type, and pick the ds folder in "/sample_raw/Data".<<BR>><<BR>> {{attachment:menuReview.gif}} * Then 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 before the MEG acquisition on the scalp surface with an ICP algorithm. Answer yes. Even if the result is not perfect, it usually improves the positioning of the head in the MEG helmet.<<BR>><<BR>> {{attachment:refine.gif}} <<BR>> {{attachment:refineBefore.gif}} {{attachment:refineAfter.gif}} * Two new files appeared in the database explorer:<<BR>><<BR>> {{attachment:linkInTree.gif}} * The channel file contains the definition of the sensors, exactly as when importing the files in the database with the "Import MEG/EEG" menu. It is saved in the folder ''(Common files)'', because the subject was created using the option "Yes, use one channel file per subject". Therefore, the same channel file will be used for all the folders of Subject01. * The node named "Link to raw file" contains all the information that was read from the continous file (file format, time vector, sampling frequency, events, bad channels, path to the original file, etc.), but no recordings. The MEG and EEG values recorded will be read directly from the native file. == Review the recordings == === Open the file === Right-click on the data file > MEG (all) > Display time series. {{attachment:displayTsMenu.gif}} You can see a new tab "Event" and a figure showing the recordings. {{attachment:rawPanel.gif}} === Navigate in time === As described in the basic tutorials, you can set the current time by using either the time panel (buttons and text field), or the figure (click on the white or grey areas of the figure). But you can notice that only a few seconds are visible in the figure, while the time panel (top left of the previous figure), indicates that we have 360s of recordings. Only a small block of the continuous file has been loaded in memory. This small time window can be configured with the tab '''Event/Page settings''', with the text boxes '''Start''' and '''Duration'''. The time series figure is similar to the ones that were presented in the previous tutorials, with a few new elements. The navigation bar at the bottom represents the time of the entire raw file, where the events are also represented by dots. The ''''<<<'''' and ''''>>>'''' buttons are the same as the ones in the time panel, and jump to the previous/next segment in the file. Clicking on the bar or dragging the red cursor change the current time window as well. === Sensor selection === Let's switch to a nicer representation of the recordings time series: click on the "Display mode" button in the toolbar of the main Brainstorm window. {{attachment:tsColumn.gif}} Now the traces are displayed in columns, but All the channels are displayed in the same figure, which makes it unreadable. Select a subset of channels by right-clicking on the figure > Display setup, or by using one of the keyboard shortcuts (Shift+A, B, C...). Default groups of sensors re available for some MEG systems, but you can also create your own groups of sensors with the menu "Edit selections". {{attachment:tsChannelSelection.gif}} === Amplitude scale === In this display mode, the amplitude scale is represented on the right of the figure. You can adjust this vertical scale: * Use the buttons "'''^'''" and "'''v'''" on the right side of the figure. The shortcuts for those buttons are indicated in the tooltips (leave the mouse for a short while over a button) * Hold the '''Shift key''' and move the mouse wheel, or use the keys "'''+'''" and "'''-'''".'''''' * Use the button "'''...'''" on the right side of the figure ("Set scale manually") to set the scale to a precise level. When scrolling in time to a different page, the amplitude scale is by default kept. You can change this behavior to re-evaluate automatically an optimal scale each time you change the current time window. This option is called "Auto-scale amplitude" and is disabled by default. To activate it: click on the "'''AS'''" button on the right of the figure, or check the menu "'''Display > Auto-scale amplitude'''" in the Event tab. == Adding/removing events == The panel "Recordings" in the Event tab lists the events that were detected on the "STI 001" (right arm) and "STI 002" (left arm) input channel: * on the left, the groups of events and the number of occurrence for each group (323 stimulations on the right, 301 on the left); * on the right, all the occurrences of the selected event group, described by the time instant at which they occur. the events occurring in the current time window are shown at the top of the figure: if you click on a marker, it gets selected in the Event tab as well; First create a new group of events, define its name (test) and its color (red), with the "Events" menu in the Event tab: {{attachment:addGroup.gif}} {{attachment:addGroup_done.gif}} Then set the current time were you want to add a new Test even, by clicking on the figure (current time = where the vertical red line is). Select the group event "Test" and add a few occurrences with any of the three methods: * In the Event tab: select the menu Events > Add event * In the time series figure: right-click > Add event * In the time series figure: press Ctrl + A {{attachment:addEvent_done.gif}} Now remove all the events occurrences, but not the group "Test": * In the Event tab: select one or more event occurrences (mouse click + Shift or Control), and press the "Delete" key. * In the time series figure: click on the event to delete (on the blue dot), and right-click > Delete event * In the time series figure: click on the event to delete, and press Control + A You can also use this interface to create events that have a temporal extension, ie. they last for more than one time sample. This is usually used to define bad/artifacted segments in the data. * In the time series window, select a time range (click + move) instead of just setting the current time. * Note that you can click outside of the white area to select the time (on top of the figure), in order not to select channels instead of moving. * Add an event (menus or Control+A): note the way it is represented in the figure and in the Event panel. * The first occurrence you add in an event group defines the event type: single time point, or time range. You cannot mix different types of events in a group. {{attachment:extEventSel.gif}} {{attachment:extEvent.gif}} === Saving modifications === When you close the raw file viewer, or the last figure that shows a part of the raw file, the dataset is unloaded, the file is released and the memory is freed. If you edited the events for this file, you would be asked whether to save the modifications or not. If you answer "Yes", the modifications to the events are only saved in the Brainstorm definition of the file, not in the raw file itself. So you would see your changes the next time you open the file with Brainstorm, but not if you open the raw file again with an external program. {{attachment:saveModif.gif}} Note that events you edit are not saved automically until that moment. As you would do with any other type of work, save your work regularly, to limit the damages of a program or computer crash. In the Event tab, menu File > Save modifications. {{attachment:saveModifMenu.gif}} == Shortcuts == === Keyboard shortcuts === * '''Left / right arrows''': * No other key: Change current time, sample by sample * With '''Control '''key: Jump to previous/next time segment (same as the "<<<" and ">>>" buttons) * With '''Shift '''key: Jump to next event of the selected group * On MacOS, these shortcuts are different: please read the tooltip from the the ">", ">>", and ">>>" buttons in the time panel to get the appropriate shortucuts. * '''Page-up / page-down''': * Same as left/right arrows, but faster. * If epochs are defined in the file: '''Control + page-up/page-down''' jumps to the next/previous epoch. * '''F3/Shift+F3''': Jump to the next/previous epoch or page * '''Plus / minus''': Adjust the vertical scale of the time series * '''Control + A''': Add / delete event occurrence * '''Control + T''': Open a 2D topography window at the current time * '''Shift + Letter''': Changes the set of electrodes currently displayed in the figure (list available by right-clicking on the figure > Display setup > ...) * '''Enter''': Display the selected channels in a separate figure (selected channels = lines on which you clicked, that are shown in red) * '''Escape''': Un-select all the selected channels * '''Delete''': Mark the selected channels as bad === Mouse shortcuts === * '''Mouse click''': Change current time * '''Mouse click on a channel''': Select the channel * '''Mouse click + move''': Select time range * '''Mouse wheel''': Zoom around current time * '''Control + mouse wheel''': Zoom vertically * '''Shift + mouse wheel''': Adjust the vertical scale of the time series * '''Right-click''': Display popup menu * '''Right-click + move''': Move in a zoomed figure * '''Double click''': Restore initial zoom settings (but do not restore the vertical scale of the time series) == Display options == When a block of recordings is loaded in memory, a few basic operations can be applied on it on the fly. === Event tab / Display menu === . {{attachment:displayOptions.gif}} * '''Remove baseline''': When selected, remove from each channel its average value over the entire current time window. This means that if you change the length of the time window, the value that is removed from each channel may change. * '''Update amplitude scale''': When selected, each time the time window changes, the vertical scale of the time series is adapted to the maximum value over the time window. When not selected: the vertical scales keeps its last value when you jump to another part of the file. * '''Use CTF 3rd gradient compensation''': Enable/disable the CTF noise-correction based on the reference sensors. * '''Use Signal Space Projections''': Apply/ignore the SSP vectors to the recordings, if some are available in the file. === Online filter === With the Filter tab in the main Brainstorm window, you can apply a band-pass filter to the recordings, or remove a set of specific frequencies (example: the high frequencies that are used to track the position of the head, or the power lines contamination). {{attachment:onlineFilter.gif}} === EEG Average reference === When reviewing EEG recordings, you can switch between the original values and the average reference by clicking on the "AVG REF" button in the toolbar of the main Brainstorm window. Average reference computation: at each time point, the average of all the electrodes values is removed to each electrode. {{attachment:avgRef.gif}} == Other views == === Topography views === Exactly as introduced in the [[Tutorials/TutExploreRecodings|tutorial #4 "Exploring the recordings"]], you can display a variety of 2D/3D mappings of those recordings. Right-click on the node "Link to raw file" > MEG (all) > ... {{attachment:topo.gif}} === Cortical sources === As presented in tutorial #5 to #7, you can compute successively for this raw file: a head model, a noise covariance matrix, and an inverse model. * '''Head model''': The goal of this tutorial is only to illustrate what you can do with the raw file viewer in Brainstorm, not to localize anything for real. Keep all the default options''' '''for the head model computation. * '''Noise covariance matrix''': Set as the baseline section the entire window before the first electric stimuation (the first occurrence of Event #5): [27.799, 41.170] seconds. Leave the other options to the default values. * '''Sources''': Keep the default options for the minimum norm method (wMNE). At the end you can see three new files in the database. Right-click on the sources and try the different visualization menus. '''__Important__''': Remember that you probably still have the online filter on; the band-pass options you specify in the Filter panel are kept until you close Brainstorm. {{attachment:sourcesTree.gif}} {{attachment:sourcesView.gif}} |
Differences between revisions 72 and 128 (spanning 56 versions)