MEG: Basic Data Processing Analysis

Similar documents
Lauri Parkkonen. Jyväskylä Summer School 2013

Modern Tools for Noninvasive Analysis of Brainwaves. Advances in Biomaterials and Medical Devices Missouri Life Sciences Summit Kansas City, March 8-9

NOISE REDUCTION TECHNIQUES IN ECG USING DIFFERENT METHODS Prof. Kunal Patil 1, Prof. Rajendra Desale 2, Prof. Yogesh Ravandle 3

Non-Invasive EEG Based Wireless Brain Computer Interface for Safety Applications Using Embedded Systems

ELECTROMYOGRAPHY UNIT-4

X Space X Time X Condition

BOLD fmri: signal source, data acquisition, and interpretation

EC209 - Improving Signal-To-Noise Ratio (SNR) for Optimizing Repeatable Auditory Brainstem Responses

FYS3240 PC-based instrumentation and microcontrollers. Signal sampling. Spring 2017 Lecture #5

Student s Learning Activities (Lichtenberg and the electrophorus)

Multihop Routing in Ad Hoc Networks

Ensemble Empirical Mode Decomposition: An adaptive method for noise reduction

Metrics for Assistive Robotics Brain-Computer Interface Evaluation

Neural Coding of Multiple Stimulus Features in Auditory Cortex

Pulse Sequences: Rapid Gradient Echo

Tektronix. Identifying Objects. A Case Study and Class Exercise. Ward Cunningham Computer Research Laboratory Tektronix, Inc. Object-Oriented Design

1Ph_FW_AC-Controller_R-L_Load -- Overview

Sensors, Signals and Noise

21/01/2014. Fundamentals of the analysis of neuronal oscillations. Separating sources

Modeling, Architectures and Signal Processing for Brain Computer Interfaces

Classification of Four Class Motor Imagery and Hand Movements for Brain Computer Interface

Magnitude & Intensity

A Gentle Introduction to Image Processing and Reconstruction in FMRI

Physiological signal(bio-signals) Method, Application, Proposal

Biomechanical Instrumentation Considerations in Data Acquisition ÉCOLE DES SCIENCES DE L ACTIVITÉ PHYSIQUE SCHOOL OF HUMAN KINETICS

EE 529 Remote Sensing Techniques. Radar

Dynamic Data-Driven Adaptive Sampling and Monitoring of Big Spatial-Temporal Data Streams for Real-Time Solar Flare Detection

COMPARISON OF VARIOUS FILTERING TECHNIQUES USED FOR REMOVING HIGH FREQUENCY NOISE IN ECG SIGNAL

Single channel noise reduction

Optical Correlator for Image Motion Compensation in the Focal Plane of a Satellite Camera

H 2 O and fat imaging

Practical Image and Video Processing Using MATLAB

Biomedical Sensor Systems Laboratory. Institute for Neural Engineering Graz University of Technology

Superposition Coding in the Downlink of CDMA Cellular Systems

EDL Group #3 Final Report - Surface Electromyograph System

Digital Halftoning. Sasan Gooran. PhD Course May 2013

Introduction to Computational Neuroscience

PET: New Technologies & Applications, Including Oncology

Dynamic Spectrum Access in Cognitive Radio Networks. Xiaoying Gan 09/17/2009

TOROIDAL ALFVÉN EIGENMODES

Large-Signal Measurements Going beyond S-parameters

Design of an Optimal High Pass Filter in Frequency Wave Number (F-K) Space for Suppressing Dispersive Ground Roll Noise from Onshore Seismic Data

40 Hz Event Related Auditory Potential

Outline. Noise and Distortion. Noise basics Component and system noise Distortion INF4420. Jørgen Andreas Michaelsen Spring / 45 2 / 45

Amplitude balancing for AVO analysis

Optical lever for KAGRA

Page 1/10 Digilent Analog Discovery (DAD) Tutorial 6-Aug-15. Figure 2: DAD pin configuration

Experimental study of broadband trailing edge noise of a linear cascade and its reduction with passive devices

A Study on the Effect of Performing VLF Withstand Tests on Field Aged Degraded Joints

1 Trigonometric Identities

THE OFFICINE GALILEO DIGITAL SUN SENSOR

Product data sheet: OT FIT 15/ /500 LT2 S

EMG. The study of muscle function through the investigation of the electrical signal the muscles produce

Improving ECG Signal using Nuttall Window-Based FIR Filter

MATLAB: Plots. The plot(x,y) command

Off-axis response of Compton and photoelectric polarimeters with a large field of view

Collaborative transmission in wireless sensor networks

Arkansas Tech University MATH 1203: Trigonometry Dr. Marcel B. Finan. Review Problems for Test #3

EKG De-noising using 2-D Wavelet Techniques

Exploring the Effects of Image Persistence in Low Frame Rate Virtual Environments

Open Access Research of Dielectric Loss Measurement with Sparse Representation

STRETCHED WIRE OFFSET MEASUREMENTS: 40 YEARS OF PRACTICE OF THIS TECHNIQUE AT CERN

ece BRADLEY DEPARTMENT of ELECTRICAL & COMPUTER ENGINEERING

Wavelet Based Classification of Finger Movements Using EEG Signals

Efficient Characterization and Testing of MRAM Devices. Parametric Testing for In-Line Monitoring

PhysioTel Digital Signal Simulator Manual

Fundamentals of neuronal oscillations and synchrony

Advanced Digital Motion Control Using SERCOS-based Torque Drives

Time-of-arrival estimation for blind beamforming

2.5D Finite Element Simulation Eddy Current Heat Exchanger Tube Inspection using FEMM

Lecture 4 Biopotential Amplifiers

REPORT DOCUMENTATION PAGE

COE/EE152: Basic Electronics. Lecture 5. Andrew Selasi Agbemenu. Outline

Magnetoencephalography and Auditory Neural Representations

Channel Modeling ETIN10. Wireless Positioning

Human Vision, Color and Basic Image Processing

Application of the new algorithm ISAR- GMSA to a linear phased array-antenna

Design of Instrumentation Systems for Monitoring Geo-Hazards in Transportation. By Barry R. Christopher Christopher Consultants Roswell, Ga.

MAT 115: Finite Math for Computer Science Problem Set 5

ECE 359 Spring 2003 Handout # 16 April 15, SNR for ANGLE MODULATION SYSTEMS. v(t) = A c cos(2πf c t + φ(t)) for FM. for PM.

Automation for a Changing World. Delta Power Meter DPM Series VARITEL

Localization of microscale devices in vivo using addressable transmitters operated as magnetic spins

Math 1111 Math Exam Study Guide

INTERPRETATION of IGEC RESULTS

Notes on Optical Amplifiers

Today s menu. Last lecture. Series mode interference. Noise and interferences R/2 V SM Z L. E Th R/2. Voltage transmission system

ni.com Sensor Measurement Fundamentals Series

A Finite Impulse Response (FIR) Filtering Technique for Enhancement of Electroencephalographic (EEG) Signal

INSTANTANEOUS FREQUENCY ESTIMATION FOR A SINUSOIDAL SIGNAL COMBINING DESA-2 AND NOTCH FILTER. Yosuke SUGIURA, Keisuke USUKURA, Naoyuki AIKAWA

ELT Receiver Architectures and Signal Processing Fall Mandatory homework exercises

EE482: Digital Signal Processing Applications

INTEGRATED APPROACH TO ECG SIGNAL PROCESSING

Validation of the Happify Breather Biofeedback Exercise to Track Heart Rate Variability Using an Optical Sensor

Thermochromic Solid-State Emission of Dipyridyl Sulfoxide Cu(I) Complexes

Design and Performance Simulation of a Ku-Band Rotating Fan-Beam Scatterometer

Chapter 6. The Production Function. Production Jargon. Production

from signals to sources asa-lab turnkey solution for ERP research

BME 701 Lecture 1. Measurement and Instrumentation

Image Processing. Image Processing. What is an Image? Image Resolution. Overview. Sources of Error. Filtering Blur Detect edges

RECOMMENDATION ITU-R P Prediction of sky-wave field strength at frequencies between about 150 and khz

Transcription:

MEG: Basic Data Processing and Time-frequency Analysis Stephan Grimault, PhD November 22, 2006

General outline 1) Basic Pre-processing and processing of MEG data basic ERF (ERP) analysis and activation map Pre-processing and processing: definitions and purpose Pre-processing steps Processing steps 2) Another way to analyze the data: Time-frequency Analysis

Basic ERF (ERP) analysis Protocol example Stimuli

Activation map (1)

Activation map (2)

Pre-processing & Processing: definition and purpose Pre-processing : - Clean the data to keep only the functional signals of interest (involves detection and correction of noise and acquisition artifacts) - optimize data processing Processing: - analysis of functional signals of interest - to evaluate to cognitive hypotheses

General outline 1) Basic Pre-processing and processing of MEG data basic ERF (ERP) analysis and activation map Pre-processing and processing: definitions and purpose Pre-processing steps Processing steps 2) Another way to analyze the data: Time-frequency Analysis

Pre-processing is not trivial: MEG measures magnetic fields from femtotesla (10-15 T) to picotesla (10-12 T) Earth s magnetic field: 4,710-5 T. small magnetic field measurements lead to artifacts

Pre-processing steps (1): Environmental Noise Reduction

Pre-processing steps (2): DC offset removal

Pre-processing (3): Frequency filtering: High-pass filter Low-pass filter

6 f 4 f f

6 f 4 f f

6 f 4 f f 0.1 f

6 f 4 f f

6 f 4 f f 0.05 f

6 f 4 f f 0.1 f 0.05 f

γ β α θ 24-40-100Hz 12-45Hz 8.5-12Hz 4-8.5Hz δ 0.05-4Hz

Sampling rate = 600Hz t = 0.0016 s Signal Theory frequency MAX = Sampling/2 Low Pass Filter = Sampling rate /4 : 600Hz 150Hz

Pre-processing (3): Frequency filtering example High Pass Filter Low Pass Filter

Pre-processing (4): Eye artifact (blinks, motion) 20pT

Eye blink correction ex. 1

Eye blink correction ex. 2

Pre-processing (5): Current Noise at 60 Hz (50 Hz Europe)

Pre-processing (6): ECG or watch artifact

Pre-processing (7): Dental artifacts

Pre-processing (8): Subject motion artifact

General outline 1) Basic Pre-processing and processing of MEG data basic ERF (ERP) analysis and activation map Pre-processing and processing: definitions and purpose Pre-processing steps Processing steps 2) Another way to analyze the data: Time-frequency Analysis

Processing (1): single trial recording

Processing (1): Why record N trial SNR N 1 10 100 200

Processing (2) : select condition

Processing (3) : averaging

Processing ends work of interpretation begins! Protocol example Bonne chance! Stimuli

Basic Pre-processing and Processing of MEG data: Summary Pre processing Environmental Noise Reduction Frequency filtering : High Pass, Low Pass Eye blink correction 60 Hz or 50 Hz (Europe) ECG or watch artifact Processing Recording N trials by condition Select condition (trigger) Averaging

General outline 1) Basic Pre-processing and processing of MEG data basic ERF (ERP) analysis and activation map Pre-processing and processing: definitions and purpose Pre-processing steps Processing steps 2) Another way to analyze the data: Time-frequency Analysis

Time Frequency map representation

Keep Information from thetrials

Average after keep information from Trials

TF map calculation : the wavelets

TF map calculation : the wavelets

Time frequency analysis: one step more

Time frequency analysis: two step more Localisation? Localisation Discussion?

Thanks to: Franco Lepore (CERNEC) Pierre Jolicoeur (CERNEC) Équipe MEG (Anne Sophie, Christophe,, Jean- Marc, Kevin, Manon, Mihaela)