zoukankan      html  css  js  c++  java
  • One EEG preprocessing pipeline

    The preprocessing pipeline of EEG data from EEG-fMRI paradigm differs from that of regular EEG data, because they are mainly influenced by MRI artefacts. However, after removing the MRI artefacts at the first step, following ones would be the same as the regular pipeline.


    Figure above: an example of raw EEG data with MRI artefacts - the very dense part. Only after removing the dominating artefacts could we go on with other preprocessing and analyses.

    Environment requirement

    • Matlab (R2015b).
    • EEGLab toolbox (v13.6.5b). Plug-ins needed: BERGEN, FMRIB, AAR.

    Testing data

    Acquired using an Brain Products system with 64 scalp channels. Standard procedure to setup the EEG system: fs = 5000 Hz, low-pass hard ware filter at 250 Hz. Using SyncBox to sync the EEG and fMRI clocks.
    Gathering continuous data for about 10-min.

    Pipeline

    1. Removing gradient artefacts: using BERGEN plug-in. Sliding window length for artefact template = 31.
    2. Filtering: first high-pass at 0.5 Hz, then low-pass at 45 Hz - do NOT using band-pass filter, do them separately.
    3. Resample the data to fs = 500 Hz. Resample before QRS detection is beneficial (as far as I can see).
    4. Detect R-peaks, visually check incorrect ones and remove cardiac artefacts: using FMRIB plug-in (optimal basis set: first 3 PCs used). Using Edit -> Event values to delete the incorrect QRS events.
    5. Import channel locations: using 'dipfit2.3standard_BESAstandard-10-5-cap385.elp'.
    6. Select the data points to be analysed - to save time.
    7. Remove ECG and bad channels (Must done before re-reference).
    8. Re-reference: average reference.
    9. Removing EOG (blinks) and EMG artefacts using AAR plug-in directly.

    Note

    • The data is a continuous data. For task data (i.e. you may epoch them first), baseline correction is needed.
    • The AAR plugin is good to detect EOG and EMG artefacts, but sometimes they do not pick up ICs relating to residuals of R-peak artefacts.
    • Considering what mentioned above, an alternative is using 'ICA + ADJUST plug-in' to replace AAR in step 9. This is a better strategy with enough experience on recognising artefact ICs.
    • Removing power interference may be needed sometimes: 'Tools -> Filter the data -> Basic FIR filter -> Lower & Higher: both 50Hz, using Notch filter'.
    • The steps and their order may vary among different groups, while the main parts remain unchanged.
    • Inspired by 52Brain. Thanks a lot to the contributors.
  • 相关阅读:
    用ssh整合时,用sessionfactory的getCurrentSession()获取不到session
    layim+signalr2.0+mongodb在线轻聊版解决方案(可提供演示)
    (绿色)修正版gooflow流程解决方案(源码分享+在线演示+UI地址下载)
    一个开源的可视化的jQuery工作流插件
    自定义流程gooflow.08 demo在线演示
    asp.net mvc 系统操作日志设计
    第三方系统平台如何对接gooflow2.0
    一个供新手把玩的jQueryUI在线文档
    Ninject 自动注册
    JQuery.imgAreaSelect 参数说明
  • 原文地址:https://www.cnblogs.com/minks/p/5990196.html
Copyright © 2011-2022 走看看