Beamformer Source Analysis in MEG
|
|
- Christine Boone
- 6 years ago
- Views:
Transcription
1 Beamformer Source Analysis in MEG Douglas Cheyne, PhD Program in Neurosciences and Mental Health Hospital for Sick Children Research Institute & Department of Medical Imaging University of Toronto CCD Macquarie University, March, 2017
2 Source Imaging using MEG Inverse Problem (Helmholtz, 1853) Sensors To determine the underlying current sources in head from the measured external potentials or fields Sources Sources Solutions are non-unique May be highly underdetermined Solution: Construct forward models of the sources and fit models to data
3 Inverse solutions in MEG Main approaches to source modeling in MEG: Dipole fitting (parametric) methods overdetermined solution, well understood need to know number of independent sources difficult to combine results across subjects Linear inverse (minimum-norm estimates) no need to specify number of sources can model distributed sources model based - requires pseudoinverse of large matrices (regularization) requires artifact-free data Spatial filtering (beamforming) adaptive ( data-driven ) stable solutions suppresses artifacts in data automatically performance degrades for highly synchronous sources
4 Spatial filtering methods A spatial filter is the weighted output of the MEG sensor array that reflects activity at a specific brain location over time Signal Space Projection (SSP) W(r) source virtual sensor d(t) = W T M(t) source strength time W(r) = L(r) = forward model (lead field) of source
5 Spatial filtering methods A spatial filter is the weighted output of the MEG sensor array that reflects activity at a specific brain location over time Signal Space Projection (SSP) noise sources source W(r) virtual sensor d(t) = W T M(t) source strength time W(r) = L(r) = forward model (lead field) of source
6 Spatial filtering methods A spatial filter is the weighted output of the MEG sensor array that reflects activity at a specific brain location over time Signal Space Projection (SSP) noise sources source W(r) virtual sensor d(t) = W T M(t) source strength time W(r) = L(r) = forward model (lead field) of source Adaptive Beamformer W(r) noise sources source source strength time W(r) = L (r) C m -1 / L (r) C m -1 L (r)
7 Spatial filtering methods Source images can be produced by generating beamformer weights at each voxel throughout brain volume at arbitrary resolution and plotting source amplitude or power Vector beamformers orthogonal current sources at each voxel* Scalar beamformers compute optimal current direction at each voxel *for EEG # directions = 3, for MEG spherical model, # directions = 2
8 Introduction of Beamforming Methods in EEG/MEG Linearly Constrained Minimum-Variance (LCMV) beamformer (Van Veen et al., 1997) first application of beamformer method to EEG/MEG inverse problem, EEG only Synthetic Aperture Magnetometry (SAM) (Robinson & Vrba, 1999) introduction of scalar beamformer and differential imaging Dynamic Imaging of Coherent Sources (DICS) (Gross et al., 2001) frequency domain beamformer, coherence analysis Eigenspace / spatiotemporal beamformer (Sekihara et al., 2001, 2002; Dalal et al., 2004) eigenspace beamformer (averaged data), NutMEG toolbox Event-related SAM (ersam) (Cheyne et al., 2004, 2006) adapted SAM algorithm to image evoked responses SAMerf (Robinson, 2004) SAM algorithm modified for short time windows and averaging ( not same as ersam ) Event-related beamformer (ERB) (Cheyne et al., 2008) ersam with optimized orientation calculation 5-D beamformer (Dalal et al., 2010) computes beamformer over separate frequency bands and time windows
9 Differential imaging using the SAM beamformer T trials Single trial data Compute beamformer Compute Source Power M channels M x M Active Covariance matrix (C A ) M x M Baseline Covariance matrix (C B ) compute beamformer weights W(r) = C -1 B/B T C -1 B Pseudo-T = W(r) T C A W(r) W(r) T C B W(r) baseline active N samples Forward solution (B) for source r source 3.0 pseudo-t 1.0 Narrow-band Filter (e.g., Hz) and segment into active and control windows repeat for all voxels SAM = Synthetic Aperture Magnetometry (Robinson & Vrba, 1999) Alternatively, compute pseudo-f = F - 1 if F > 1. = 1-1/F if F < 1.0 where F = W(r) T C A W(r) / W(r) T C B W(r)
10 SAM versus Event-related Beamformer
11 Event-related (spatiotemporal) beamformer T trials Single trial data Compute beamformer source activity for voxel j S(j,t) = W(j) x m(t) M channels M x M Covariance matrix (C) compute beamformer weights W(j) = C -1 B/B T C -1 B N samples Forward solution (B) for source j source t = 20 ms Bandpass filter (e.g., 1 50 Hz) Take absolute value and repeat for all voxels Map power at latency of evoked response 3.0 pseudo-z 1.0
12 Beamformer localization of premovement motor field (MF) Time course at peak (virtual sensor) DC-15 Hz button press Self-paced right index finger movement Recording time 10 min ( movements) ERB image in MRIViewer (t = -40 ms, threshold = FWHM)
13 Introduction to Beamformers Beamformer Units and Weight Normalization: Beamformer filter can only suppress noise sources that are correlated across sensors Uncorrelated noise (e.g., random system noise) will be amplified by the weights in a non-uniform manner, with increasing distortion with increasing distance from the sensors This gain error must be removed from image data
14 Spatial filtering methods weight normalization virtual sensor at peak (na-m) s(r,t) = w(r) T m(t) Source location Non-normalized (units = na-m) virtual sensor at peak (pseudo-z) Source activity (peak= 20 nam) + Gaussian noise s(r,t) = wn(r) T m(t) Normalized (units = pseudo-z)
15 Spatial filtering methods weight normalization The non-uniform distortion of beamformer images can be removed by normalizing the weight vector for each voxel w(r), to have unit length ( unit-gain or Borgiatti-Kaplan beamformer) wn M ( ( )) 2 ( r) = w(r) wi r i=0 Units of w(r) = A-m / Tesla (output of beamformer is A-m) Units of w n (r) are arbitrary units (output of beamformer in arbitrary units) The neural activity index (Van Veen, 1997) and pseudo-z (Robinson and Vrba, 1999) include an additional scaling of the unit-gain beamformer to units of specified noise n w wn ( ) 2 ( r) = w(r) wi( r)n w M i=0 Output of beamformer is scaled to units of noise variance. Robinson termed this pseudo-z
16 Introduction to Beamformers Interpreting the spatial extent of activity peaks in beamformer images Do bigger blobs mean more (stronger) activation?
17 Effect of SNR on beamformer resolution Simulated bilateral auditory cortex sources X = 0.0 cm Y = 5.5 and -5.5 cm Z = 6.0 cm + Gaussian noise Q= 60 na-m Q= 40 na-m Q= 20 na-m Stronger source (higher SNR) = more focal source in image
18 Introduction to Beamformers Beamformer solutions depend on: 1) Accuracy of forward solution (L) 2) Ability to obtain stable inverse of covariance matrix (C m ) W(r) = L (r) C m -1 / L (r) C m -1 L (r) Depends on: condition number of covariance matrix (ratio of max / min eigenvalues) number of channels (less is better?) bandwidth of data (smaller bandwidth increases covariance error) number of time samples in data (about 30 s recommended) From: Brookes M. et al., (2008) Optimizing experimental design for MEG beamformer imaging. NeuroImage 39:
19 How much data is required? Decrease in FWHM of beamformer peaks / increase in pseudo-z power with increasing amount of data used to compute covariance matrix (from Brookes et al., 2008)
20 Covariance Matrix Regularization: Beamformers - Regularization In cases where data covariance is rank deficient (e.g., PCA or SSS has been applied to the raw data) or ill-conditioned (e.g., too few samples, averaged data) regularization of covariance matrix may be required [ ] 1 L θ L T θ ( C + µσ) 1 L θ W θ = C + µσ # $ % & 1 Diagonal regularization As µ à As µ à 0 W θ = L θ! " L T θ L θ # $ 1 W θ = C 1 L θ " # L T θ C 1 L θ $ % 1 (minimum noise sensitivity) (maximum spatial resolution) L q = forward solution for target location C = data covariance matrix µ = regularization parameter S = diagonal matrix of sensor noise Σ = σ 11 σ 22 0 σ σ MM
21 Introduction to Beamformers Do I need to clean my data? Advantages: reduces amount of noise sources beamformer has to suppress may be required if noise sources are large, close to brain sources Disadvantages: noise reduction techniques will reduce rank of covariance matrix requiring regularization (reduces spatial resolution)
22 Beamformer suppression of eye-movement artifacts adult performing self-paced movements (KIT-Macquarie MEG 160) Sensor data (all channels) (0 15 Hz) During fixation During visual scanning Motor cortex activation (t = -75 ms) Motor cortex time-course 40 BJ_rightIndex_voxel_0.0_3.0_9.0_(0.944_ 0.264_0.198) with fixation with eye movement 10 Moment (nam) Scanning Fixating Time (sec)
23 Suppression of ferromagnetic artifacts Motor Field (MF) localization in adult subject with metal retainer Average (frontal sensors) 200 ft ERB source image (2 mm) t =- 40 ms (DC 15 Hz) Right index finger movement L R
24 Limitations of beamforming methods The Source Correlation Problem: Beamformers will fail to detect highly synchronous (zero-phase correlated) sources
25 Spatial filtering methods effects of source correlation Simulation (3 sources) Source power Source power Source activity Source power Source power
26 Spatial filtering methods effects of source correlation Simulation (auditory evoked fields) Source 1 (0, 5.5, 6.0) Source 2 (0, -5.5, 6.0) Gaussian noise (10-20 ft / Hz 1/2 ) 150 trials Beamformer source reconstruction left right Right hemisphere source jittered trial-by-trial by 6 ms Average (all sensors) Trial-by-trial latency jitter of 6 ms reduces effect of correlation
27 Spatial filtering methods effects of source correlation From: Cheyne D., Bostan AC., Gaetz W, Pang EW. (2007) Event-related beamforming: A robust method for presurgical functional mapping using MEG Clinical Neurophysiology, Vol 118 pp
28 (a) Whole-Array Beamforming (b) Half-Array Beamforming pseudo Z pseudo Z pseudo Z (c) Left Hemisphere 4 3 N1m P1m -4 P2m Time (ms) Right Hemisphere Whole Array Half Array Time (ms) From Herdman A. and Cheyne D. A Practical Guide to MEG and Beamforming. In: T. Handy (Ed) Brain Signal Analysis: Advances in Neuroelectric and Neuromagnetic Methods. MIT Press, Massachusetts pp
29 Correcting for correlations in beamformer solutions Coherent Source Region Suppression method (Dalal et al., 2006) identify region where correlated source is know to exist can be done iteratively to image (e.g., left and right hemisphere sources) suppression ROI Image 1 (RH suppression) Image 2 (RH suppression) corrected image = max(image1,image2)
30 Limitations of beamforming methods Source Polarity Problem: Scalar beamformer gives optimal current direction but true polarity of the source is unknown, since beamformer weights derived from source power (covariance) This means source polarity can flip randomly form voxel to voxel Solutions: plot absolute power (LCMV, ERB imaging) flip source waveforms to have uniform polarity across subjects constrain source orientation to cortical surface
31 Adaptive (Minimum-variance) Beamforming Source Polarity Problem: Source waveforms (virtual sensors) in brainwave can be rectified or flipped Alternatively vector (LCMV) beamformer can be used waveforms will be rectified!
32 Some recent publications using BrainWave Mersov A., Jobst C., Cheyne D., De Nil, L. (2016) Sensorimotor oscillations prior to speech onset reflect altered motor networks in adults who stutter. Frontiers in Human Neuroscience 10: 443. Pu Y., Cornwell BR., Cheyne D., Johnson BW. (2016) The functional role of human right hippocampal theta rhythm in environmental encoding during spatial navigation. Human Brain Mapping 38: Isabella S., Ferrari P., Jobst, C., Cheyne JA., Cheyne D. (2015) Complementary roles of cortical oscillations in automatic and controlled processing during a rapid serial task. NeuroImage 118:
33 The BrainWave Toolbox What it does: computes event-related (evoked) and SAM beamformer source images computes source activity time courses ( virtual sensors ) average or time-frequency plots volumetric source images in MEG or MNI space (latter requires SPM8 or SPM12 installed) warping of template MRI to head shapes (if MRI not available, requires SPM) surface based images using Freesurfer or CIVET cortical meshes group analysis with contrasts and permutation test group averaging of source activity waveforms or their time-frequency transformations built-in 4D viewers for volumes or surfaces
34 The BrainWave Toolbox Software requirements requires MATLAB w/ standard toolboxes runs on Linux, Mac OS X or Windows platforms optimized using multi-threaded compiled C++ mex functions (more cores = faster ) supports CTF, Yokogawa-KIT, Elekta-Neuromag MEG data import and co-register DICOM MRI fit spherical models to MRI surfaces (requires FSL, Mac and Linux versions)
35 Sample Datasets: CTF_adult: self-paced motor responses (left and right) in adult subject 151-channel CTF MEG system MRI and co-registered Freesurfer and CIVET surfaces KIT_child: motor responses (left and right) in preschooler 112 channel Yokogawa-KIT Child MEG No MRI CTF_group: 8 subjects performing motor task (left and right) 151-channel CTF MEG system pre-calculated group images (volume and surface) Elekta_SEF: right median-nerve stimulation in adult subject 306-channel VectorView MEG system no-sss preprocessing
36
37 BrainWave: Data Organization and Preprocessing Steps 1. Import and epoch MEG data epoch continuous datasets into one dataset per condition, should be same length (e.g., for comparing time courses) each dataset name must begin with <subjectid> code followed by underscore
38 Suggested Reading Beamforming Methods: Barnes, G.R., Hillebrand, A., Statistical flattening of beamformer images. Hum. Brain Mapp. 18, Brookes M., Vrba J., Robinson S., Stevenson C., Peters A., Barnes A., Morris P. (2008) Optimising experimental design for MEG beamformer imaging. NeuroImage 39: Cheyne D., Bakhtazad L. and Gaetz W. (2006) Spatiotemporal mapping of cortical activity accompanying voluntary movements using an event-related beamforming approach Human Brain Mapping 27: Cheyne D., Bostan AC., Gaetz W, Pang EW. (2007) Event-related beamforming: A robust method for presurgical functional mapping using MEG Clinical Neurophysiology, Vol 118 pp Gross, J., Kujala, J., Hamalainen, M., Timmermann, L., Schnitzler, A., Salmelin, R., Dynamic imaging of coherent sources: studying neural interactions in the human brain. Proc. Natl. Acad. Sci. U. S. A. 98 (2), Herdman A., Cheyne D. (2009) A practical guide for MEG and Beamforming. In. Handy T. Brain Signal Analysis MIT Press pp Hillebrand, A., Singh, K.D., Holliday, I.E., Furlong, P.L., Barnes, G.R., A new approach to neuroimaging with magnetoencephalography. Hum. Brain Mapp. 25, Hillebrand, A., Barnes, G.R., The use of anatomical constraints with MEG beamformers. NeuroImage 20, Robinson, S. and Vrba J. (1999). Functional neuroimaging by synthetic aperture magnetometry. Nenonen J, Ilmoniemi RJ, Katila T, editors. Proceedings of Biomag 2000 Conference. p Sekihara, K., Nagarajan, S.S., Poeppel, D., Marantz, A., Miyashita, Y. (2001) Reconstructing spatio-temporal activities of neural sources using an MEG vector beamformer technique. IEEE Trans Biomed Eng 48, Van Veen, B.D., van Drongelen, W., Yuchtman, M., Suzuki, A., (1997) Localization of brain electrical activity via linearly constrained minimum variance spatial filtering. IEEE Trans Biomed Eng 44,
Spatial Filtering Methods in MEG. Part 3: Template Normalization and Group Analysis"
Spatial Filtering Methods in MEG Part 3: Template Normalization and Group Analysis" Douglas Cheyne, PhD" Program in Neurosciences and Mental Health" Hospital for Sick Children Research Institute " &" Department
More informationUser Documentation. BrainWave v User s Documentation -- Toronto, Canada
User Documentation BrainWave v1.5 2012 -- User s Documentation -- Toronto, Canada http://cheynelab.utoronto.ca Table of Contents Getting Started... 5 Introduction... 5 System Requirements... 5 Disclaimer
More informationUser s Documentation 2.3. BrainWave v User s Documentation -- Toronto, Canada
User s Documentation 2.3 BrainWave v2.2 2013 -- User s Documentation -- Toronto, Canada http://cheynelab.utoronto.ca Table of Contents Getting Started... 5 Introduction... 5 System Requirements... 6 Disclaimer
More informationUser s Documentation 2.0. BrainWave v User s Documentation -- Toronto, Canada
User s Documentation 2.0 BrainWave v2.0 2012 -- User s Documentation -- Toronto, Canada http://cheynelab.utoronto.ca Table of Contents Getting Started... 5 Introduction... 5 System Requirements... 5 Disclaimer
More informationUSER S DOCUMENTATION v.3.1
USER S DOCUMENTATION v.3.1 BrainWave 2015 Hospital for Sick Children, Toronto, Ontario http://cheynelab.utoronto.ca Table of Contents GETTING STARTED... 4 INTRODUCTION... 4 SYSTEM REQUIREMENTS... 5 DISCLAIMER
More informationUSER S DOCUMENTATION v.3.0 (beta)
USER S DOCUMENTATION v.3.0 (beta) BrainWave 2015 Toronto, Ontario http://cheynelab.utoronto.ca Table of Contents GETTING STARTED... 4 INTRODUCTION... 4 SYSTEM REQUIREMENTS... 5 DISCLAIMER & LICENSE...
More informationTHE SPATIAL filter is a popular method used to reconstruct
1358 IEEE TRANSACTIONS ON BIOMEDICAL ENGINEERING, VOL. 57, NO. 6, JUNE 2010 Array-Gain Constraint Minimum-Norm Spatial Filter With Recursively Updated Gram Matrix For Biomagnetic Source Imaging Isamu Kumihashi
More informationUSER S DOCUMENTATION v.3.4
USER S DOCUMENTATION v.3.4 BrainWave 2018 Hospital for Sick Children, Toronto, Ontario http://cheynelab.utoronto.ca Table of Contents VERSION UPDATES... 4 GETTING STARTED... 5 INTRODUCTION... 5 SYSTEM
More informationWe are IntechOpen, the world s leading publisher of Open Access books Built by scientists, for scientists. International authors and editors
We are IntechOpen, the world s leading publisher of Open Access books Built by scientists, for scientists 3,900 116,000 120M Open access books available International authors and editors Downloads Our
More information- Graphical editing of user montages for convenient data review - Import of user-defined file formats using generic reader
Data review and processing Source montages and 3D whole-head mapping Onset of epileptic seizure with 3D whole-head maps and hemispheric comparison of density spectral arrays (DSA) Graphical display of
More informationMEG and EEG analysis using
MEG and EEG analysis using François Tadel Biomag 2012 Paris August 31, 2012 Foretaste Graphic interface 2 Foretaste Scripting environment Rapid selection of files and processes to apply Automatic generation
More informationSource Reconstruction in MEG & EEG
Source Reconstruction in MEG & EEG ~ From Brain-Waves to Neural Sources ~ Workshop Karolinska Institutet June 16 th 2017 Program for today Intro Overview of a source reconstruction pipeline Overview of
More informationBESA Research. CE certified software package for comprehensive, fast, and user-friendly analysis of EEG and MEG
BESA Research CE certified software package for comprehensive, fast, and user-friendly analysis of EEG and MEG BESA Research choose the best analysis tool for your EEG and MEG data BESA Research is the
More informationAnalysis of Functional MRI Timeseries Data Using Signal Processing Techniques
Analysis of Functional MRI Timeseries Data Using Signal Processing Techniques Sea Chen Department of Biomedical Engineering Advisors: Dr. Charles A. Bouman and Dr. Mark J. Lowe S. Chen Final Exam October
More informationOHBA M/EEG Analysis Workshop. Mark Woolrich Diego Vidaurre Andrew Quinn Romesh Abeysuriya Robert Becker
OHBA M/EEG Analysis Workshop Mark Woolrich Diego Vidaurre Andrew Quinn Romesh Abeysuriya Robert Becker Workshop Schedule Tuesday Session 1: Preprocessing, manual and automatic pipelines Session 2: Task
More informationImproving the Interpretability of All-to-All Pairwise Source Connectivity Analysis in MEG With Nonhomogeneous Smoothing
r Human Brain Mapping 32:426 437 (2011) r Improving the Interpretability of All-to-All Pairwise Source Connectivity Analysis in MEG With Nonhomogeneous Smoothing Jan-Mathijs Schoffelen* and Joachim Gross
More informationMEG and EEG analysis using
MEG and EEG analysis using François Tadel Dalhousie University, Halifax February 16, 2013 Foretaste Graphic interface 2 Foretaste Scripting environment Rapid selection of files and processes to apply Automatic
More informationMEG and EEG analysis using
MEG and EEG analysis using François Tadel & Sylvain Baillet May 2014 Foretaste Graphic interface 2 Foretaste Scripting environment Rapid selection of files and processes to apply Automatic generation of
More informationBioEST ver. 1.5beta User Manual Chang-Hwan Im (Ph.D.)
BioEST ver. 1.5beta User Manual Chang-Hwan Im (Ph.D.) http://www.bioest.com E-mail: ichism@elecmech.snu.ac.kr imxxx010@umn.edu 1. What is BioEST? BioEST is a special edition of SNUEEG(http://www.bioinverse.com),
More informationThe organization of the human cerebral cortex estimated by intrinsic functional connectivity
1 The organization of the human cerebral cortex estimated by intrinsic functional connectivity Journal: Journal of Neurophysiology Author: B. T. Thomas Yeo, et al Link: https://www.ncbi.nlm.nih.gov/pubmed/21653723
More informationGroup Sta*s*cs in MEG/EEG
Group Sta*s*cs in MEG/EEG Will Woods NIF Fellow Brain and Psychological Sciences Research Centre Swinburne University of Technology A Cau*onary tale. A Cau*onary tale. A Cau*onary tale. Overview Introduc*on
More informationCOHERENT MEG/EEG SOURCE LOCALIZATION IN TRANSFORMED DATA SPACE
Biomedical Engineering: Applications, Basis and Communications, Vol. 22, No. 5 (21) 351 365 DOI: 1.415/S11623721211 COHERENT MEG/EEG SOURCE LOCALIZATION IN TRANSFORMED DATA SPACE Junpeng Zhang,, Sarang
More informationTheory: modeling, localization and imaging
Electromagnetic Brain Mapping with MEG/EEG Theory: modeling, localization and imaging Sylvain Baillet Imaging Group Cognitive Neuroscience & Brain Imaging Lab. Hôpital de la Salpêtrière CNRS UPR640 - LENA
More informationNew Approaches for EEG Source Localization and Dipole Moment Estimation. Shun Chi Wu, Yuchen Yao, A. Lee Swindlehurst University of California Irvine
New Approaches for EEG Source Localization and Dipole Moment Estimation Shun Chi Wu, Yuchen Yao, A. Lee Swindlehurst University of California Irvine Outline Motivation why EEG? Mathematical Model equivalent
More informationM/EEG pre-processing 22/04/2014. GUI Script Batch. Clarification of terms SPM speak. What do we need? Why batch?
22/04/2014 Clarification of terms SPM speak GUI Script Batch M/EEG pre-processing Vladimir Litvak Wellcome Trust Centre for Neuroimaging UCL Institute of Neurology Why batch? What do we need? As opposed
More informationA comparison of minimum norm and MUSIC for a combined MEG/EEG sensor array
Adv. Radio Sci., 10, 99 104, 2012 doi:10.5194/ars-10-99-2012 Author(s) 2012. CC Attribution 3.0 License. Advances in Radio Science A comparison of minimum norm and MUSIC for a combined MEG/EEG sensor array
More informationStatistical Analysis of Neuroimaging Data. Phebe Kemmer BIOS 516 Sept 24, 2015
Statistical Analysis of Neuroimaging Data Phebe Kemmer BIOS 516 Sept 24, 2015 Review from last time Structural Imaging modalities MRI, CAT, DTI (diffusion tensor imaging) Functional Imaging modalities
More informationHST.583 Functional Magnetic Resonance Imaging: Data Acquisition and Analysis Fall 2006
MIT OpenCourseWare http://ocw.mit.edu HST.583 Functional Magnetic Resonance Imaging: Data Acquisition and Analysis Fall 2006 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
More informationA Simple Generative Model for Single-Trial EEG Classification
A Simple Generative Model for Single-Trial EEG Classification Jens Kohlmorgen and Benjamin Blankertz Fraunhofer FIRST.IDA Kekuléstr. 7, 12489 Berlin, Germany {jek, blanker}@first.fraunhofer.de http://ida.first.fraunhofer.de
More informationBiomagnetic inverse problems:
Biomagnetic inverse problems: Magnetic resonance electrical property tomography (MREPT) and magnetoencephalography (MEG) 2018 Aug. 16 The University of Tokyo Takaaki Nara 1 Contents Measurement What is
More informationData Visualisation in SPM: An introduction
Data Visualisation in SPM: An introduction Alexa Morcom Edinburgh SPM course, April 2010 Centre for Cognitive & Neural Systems/ Department of Psychology University of Edinburgh Visualising results remembered
More informationEffect of age and dementia on topology of brain functional networks. Paul McCarthy, Luba Benuskova, Liz Franz University of Otago, New Zealand
Effect of age and dementia on topology of brain functional networks Paul McCarthy, Luba Benuskova, Liz Franz University of Otago, New Zealand 1 Structural changes in aging brain Age-related changes in
More informationModeling versus Accuracy in EEG and MEG Data
Modeling versus Accuracy in EEG and MEG Data John C. Mosher, Richard M. Leahy, Ming-xiong Huang, Michael E. Spencer Los Alamos National Laboratory Group P- MS D Los Alamos, NM 7 mosher@lanl.gov, () -7
More informationAttention modulates spatial priority maps in human occipital, parietal, and frontal cortex
Attention modulates spatial priority maps in human occipital, parietal, and frontal cortex Thomas C. Sprague 1 and John T. Serences 1,2 1 Neuroscience Graduate Program, University of California San Diego
More informationCS 229 Final Project Report Learning to Decode Cognitive States of Rat using Functional Magnetic Resonance Imaging Time Series
CS 229 Final Project Report Learning to Decode Cognitive States of Rat using Functional Magnetic Resonance Imaging Time Series Jingyuan Chen //Department of Electrical Engineering, cjy2010@stanford.edu//
More informationData Visualisation in SPM: An introduction
Data Visualisation in SPM: An introduction Alexa Morcom Edinburgh SPM course, April 2015 SPMmip [-30, 3, -9] 3 Visualising results remembered vs. fixation contrast(s) < < After the results table - what
More informationJournal of Articles in Support of The Null Hypothesis
Data Preprocessing Martin M. Monti, PhD UCLA Psychology NITP 2016 Typical (task-based) fmri analysis sequence Image Pre-processing Single Subject Analysis Group Analysis Journal of Articles in Support
More informationIntroduction to fmri. Pre-processing
Introduction to fmri Pre-processing Tibor Auer Department of Psychology Research Fellow in MRI Data Types Anatomical data: T 1 -weighted, 3D, 1/subject or session - (ME)MPRAGE/FLASH sequence, undistorted
More informationMEG & PLS PIPELINE: SOFTWARE FOR MEG DATA ANALYSIS AND PLS STATISTICS
MEG & PLS PIPELINE: SOFTWARE FOR MEG DATA ANALYSIS AND PLS STATISTICS USER DOCUMENTATION VERSION: 2.00 SEPT. 16, 2014 MEG & PLS Pipeline ([MEG]PLS) Copyright 2013-2014, Michael J. Cheung & Natasa Kovacevic
More informationVisualization of Magnetoencephalographic Data Using Minimum Current Estimates
NeuroImage 10, 173 180 (1999) Article ID nimg.1999.0454, available online at http://www.idealibrary.com on Visualization of Magnetoencephalographic Data Using Minimum Current Estimates K. Uutela,* M. Hämäläinen,*
More informationFMRI data: Independent Component Analysis (GIFT) & Connectivity Analysis (FNC)
FMRI data: Independent Component Analysis (GIFT) & Connectivity Analysis (FNC) Software: Matlab Toolbox: GIFT & FNC Yingying Wang, Ph.D. in Biomedical Engineering 10 16 th, 2014 PI: Dr. Nadine Gaab Outline
More informationTask Level Hierarchical System for BCI-enabled Shared Autonomy Iretiayo Akinola, Boyuan Chen, Jonathan Koss, Aalhad Patankar, Jake Varley & Peter
Task Level Hierarchical System for BCI-enabled Shared Autonomy Iretiayo Akinola, Boyuan Chen, Jonathan Koss, Aalhad Patankar, Jake Varley & Peter Allen Columbia University Shared Autonomy Agent 1 Agent
More informationSection 9. Human Anatomy and Physiology
Section 9. Human Anatomy and Physiology 9.1 MR Neuroimaging 9.2 Electroencephalography Overview As stated throughout, electrophysiology is the key tool in current systems neuroscience. However, single-
More informationGroup Sta/s/cs with BESA and BrainWave
24/04/12 Group Sta/s/cs with BESA and Graciela Tesan, PhD CCD From raw data to p- values 1 MEG and EEG data preprocessing and analysis BESA Sta/s/cs: analyses of EEG ac/vity and MEG distributed sources
More informationNA-MIC National Alliance for Medical Image Computing fmri Data Analysis
NA-MIC fmri Data Analysis Sonia Pujol, Ph.D. Wendy Plesniak, Ph.D. Randy Gollub, M.D., Ph.D. Acknowledgments NIH U54EB005149 Neuroimage Analysis Center NIH P41RR013218 FIRST Biomedical Informatics Research
More informationResearch Article MEG/EEG Source Reconstruction, Statistical Evaluation, and Visualization with NUTMEG
Hindawi Publishing Corporation Computational Intelligence and Neuroscience Volume 211, Article ID 758973, 17 pages doi:1.1155/211/758973 Research Article MEG/EEG Source Reconstruction, Statistical Evaluation,
More informationBasic Introduction to Data Analysis. Block Design Demonstration. Robert Savoy
Basic Introduction to Data Analysis Block Design Demonstration Robert Savoy Sample Block Design Experiment Demonstration Use of Visual and Motor Task Separability of Responses Combined Visual and Motor
More informationSupplementary Figure 1
Supplementary Figure 1 BOLD and CBV functional maps showing EPI versus line-scanning FLASH fmri. A. Colored BOLD and CBV functional maps are shown in the highlighted window (green frame) of the raw EPI
More informationEstimating Noise and Dimensionality in BCI Data Sets: Towards Illiteracy Comprehension
Estimating Noise and Dimensionality in BCI Data Sets: Towards Illiteracy Comprehension Claudia Sannelli, Mikio Braun, Michael Tangermann, Klaus-Robert Müller, Machine Learning Laboratory, Dept. Computer
More informationSelected topics on fmri image processing
fmri Symposium, 24 January, 2006, Ghent Selected topics on fmri image processing Jan Sijbers The fmri folks of the Vision Lab UZA Bio-Imaging Lab Vision Lab Univ. Maastricht T.U. Delft 1 Overview generalized
More informationB degrees of freedom are known as partially adaptive arrays
IEEE TRANSACTIONS ON ANTENNAS AND PROPAGATION, VOL. 36. NO. 3. MARCH 1988 357 Eigenstructure Based Partially Adaptive Array Design Abstract-A procedure is presented for designing partially adaptive arrays
More informationADJUST: An Automatic EEG artifact Detector based on the Joint Use of Spatial and Temporal features
ADJUST: An Automatic EEG artifact Detector based on the Joint Use of Spatial and Temporal features A Tutorial. Marco Buiatti 1 and Andrea Mognon 2 1 INSERM U992 Cognitive Neuroimaging Unit, Gif sur Yvette,
More information7/15/2016 ARE YOUR ANALYSES TOO WHY IS YOUR ANALYSIS PARAMETRIC? PARAMETRIC? That s not Normal!
ARE YOUR ANALYSES TOO PARAMETRIC? That s not Normal! Martin M Monti http://montilab.psych.ucla.edu WHY IS YOUR ANALYSIS PARAMETRIC? i. Optimal power (defined as the probability to detect a real difference)
More informationAUTOMATIC REMOVAL OF SPARSE ARTIFACTS IN ELECTROENCEPHALOGRAM
AUTOMATIC REMOVAL OF SPARSE ARTIFACTS IN ELECTROENCEPHALOGRAM Petr Tichavský, Miroslav Zima Institute of Information Theory and Automation, Pod vodárenskou věží, Prague, Czech Republic Faculty of Nuclear
More informationArtifact detection and repair in fmri
Artifact detection and repair in fmri Paul K. Mazaika, Ph.D. Center for Interdisciplinary Brain Sciences Research (CIBSR) Division of Interdisciplinary Behavioral Sciences Stanford University School of
More informationGroup (Level 2) fmri Data Analysis - Lab 4
Group (Level 2) fmri Data Analysis - Lab 4 Index Goals of this Lab Before Getting Started The Chosen Ten Checking Data Quality Create a Mean Anatomical of the Group Group Analysis: One-Sample T-Test Examine
More informationFMRI Pre-Processing and Model- Based Statistics
FMRI Pre-Processing and Model- Based Statistics Brief intro to FMRI experiments and analysis FMRI pre-stats image processing Simple Single-Subject Statistics Multi-Level FMRI Analysis Advanced FMRI Analysis
More informationMultivariate pattern classification
Multivariate pattern classification Thomas Wolbers Space & Ageing Laboratory (www.sal.mvm.ed.ac.uk) Centre for Cognitive and Neural Systems & Centre for Cognitive Ageing and Cognitive Epidemiology Outline
More informationINDEPENDENT COMPONENT ANALYSIS APPLIED TO fmri DATA: A GENERATIVE MODEL FOR VALIDATING RESULTS
INDEPENDENT COMPONENT ANALYSIS APPLIED TO fmri DATA: A GENERATIVE MODEL FOR VALIDATING RESULTS V. Calhoun 1,2, T. Adali, 2 and G. Pearlson 1 1 Johns Hopkins University Division of Psychiatric Neuro-Imaging,
More informationClassification of Subject Motion for Improved Reconstruction of Dynamic Magnetic Resonance Imaging
1 CS 9 Final Project Classification of Subject Motion for Improved Reconstruction of Dynamic Magnetic Resonance Imaging Feiyu Chen Department of Electrical Engineering ABSTRACT Subject motion is a significant
More informationFunctional MRI in Clinical Research and Practice Preprocessing
Functional MRI in Clinical Research and Practice Preprocessing fmri Preprocessing Slice timing correction Geometric distortion correction Head motion correction Temporal filtering Intensity normalization
More informationFunctional MRI data preprocessing. Cyril Pernet, PhD
Functional MRI data preprocessing Cyril Pernet, PhD Data have been acquired, what s s next? time No matter the design, multiple volumes (made from multiple slices) have been acquired in time. Before getting
More informationSingle Subject Demo Data Instructions 1) click "New" and answer "No" to the "spatially preprocess" question.
(1) conn - Functional connectivity toolbox v1.0 Single Subject Demo Data Instructions 1) click "New" and answer "No" to the "spatially preprocess" question. 2) in "Basic" enter "1" subject, "6" seconds
More informationSPM8 for Basic and Clinical Investigators. Preprocessing. fmri Preprocessing
SPM8 for Basic and Clinical Investigators Preprocessing fmri Preprocessing Slice timing correction Geometric distortion correction Head motion correction Temporal filtering Intensity normalization Spatial
More informationExploratory Analysis of High Dimensional Time Series with Applications to Multichannel Electroencephalograms
Exploratory Analysis of High Dimensional Time Series with Applications to Multichannel Electroencephalograms arxiv:6.7684v [stat.ap] 24 Oct 26 Yuxiao Wang, Chee-Ming Ting 2 and Hernando Ombao,3 Abstract
More informationMapping of Hierarchical Activation in the Visual Cortex Suman Chakravartula, Denise Jones, Guillaume Leseur CS229 Final Project Report. Autumn 2008.
Mapping of Hierarchical Activation in the Visual Cortex Suman Chakravartula, Denise Jones, Guillaume Leseur CS229 Final Project Report. Autumn 2008. Introduction There is much that is unknown regarding
More informationMethods for data preprocessing
Methods for data preprocessing John Ashburner Wellcome Trust Centre for Neuroimaging, 12 Queen Square, London, UK. Overview Voxel-Based Morphometry Morphometry in general Volumetrics VBM preprocessing
More informationOptimum Array Processing
Optimum Array Processing Part IV of Detection, Estimation, and Modulation Theory Harry L. Van Trees WILEY- INTERSCIENCE A JOHN WILEY & SONS, INC., PUBLICATION Preface xix 1 Introduction 1 1.1 Array Processing
More informationNIH Public Access Author Manuscript Neuroimage. Author manuscript; available in PMC 2010 February 1.
NIH Public Access Author Manuscript Published in final edited form as: Neuroimage. 2009 February 1; 44(3): 932 946. doi:10.1016/j.neuroimage.2008.05.063. A Distributed Spatio-Temporal EEG/MEG Inverse Solver
More informationAnalysis of fmri data within Brainvisa Example with the Saccades database
Analysis of fmri data within Brainvisa Example with the Saccades database 18/11/2009 Note : All the sentences in italic correspond to informations relative to the specific dataset under study TP participants
More informationMultivariate Pattern Classification. Thomas Wolbers Space and Aging Laboratory Centre for Cognitive and Neural Systems
Multivariate Pattern Classification Thomas Wolbers Space and Aging Laboratory Centre for Cognitive and Neural Systems Outline WHY PATTERN CLASSIFICATION? PROCESSING STREAM PREPROCESSING / FEATURE REDUCTION
More informationVersion. Getting Started: An fmri-cpca Tutorial
Version 11 Getting Started: An fmri-cpca Tutorial 2 Table of Contents Table of Contents... 2 Introduction... 3 Definition of fmri-cpca Data... 3 Purpose of this tutorial... 3 Prerequisites... 4 Used Terms
More informationMEG Laboratory Reference Manual for research use
MEG Laboratory Reference Manual for research use for Ver. 1R007B Manual version 20040224 Index 1. File... 11 1.1 New... 11 1.2 Open... 11 1.3 Transfer...... 11 1.4 Suspended File List... 12 1.5 Save...
More informationTopographic Mapping with fmri
Topographic Mapping with fmri Retinotopy in visual cortex Tonotopy in auditory cortex signal processing + neuroimaging = beauty! Topographic Mapping with fmri Retinotopy in visual cortex Tonotopy in auditory
More informationWhat is a network? Network Analysis
What is a network? Network Analysis Valerie Cardenas Nicolson Associate Adjunct Professor Department of Radiology and Biomedical Imaging Complex weblike structures Cell is network of chemicals connected
More informationFunction-Structure Integration in FreeSurfer
Function-Structure Integration in FreeSurfer Outline Function-Structure Integration Function-Structure Registration in FreeSurfer fmri Analysis Preprocessing First-Level Analysis Higher-Level (Group) Analysis
More informationINDEPENDENT COMPONENT ANALYSIS WITH FEATURE SELECTIVE FILTERING
INDEPENDENT COMPONENT ANALYSIS WITH FEATURE SELECTIVE FILTERING Yi-Ou Li 1, Tülay Adalı 1, and Vince D. Calhoun 2,3 1 Department of Computer Science and Electrical Engineering University of Maryland Baltimore
More informationGLM for fmri data analysis Lab Exercise 1
GLM for fmri data analysis Lab Exercise 1 March 15, 2013 Medical Image Processing Lab Medical Image Processing Lab GLM for fmri data analysis Outline 1 Getting Started 2 AUDIO 1 st level Preprocessing
More informationPreprocessing II: Between Subjects John Ashburner
Preprocessing II: Between Subjects John Ashburner Pre-processing Overview Statistics or whatever fmri time-series Anatomical MRI Template Smoothed Estimate Spatial Norm Motion Correct Smooth Coregister
More informationAutomated MR Image Analysis Pipelines
Automated MR Image Analysis Pipelines Andy Simmons Centre for Neuroimaging Sciences, Kings College London Institute of Psychiatry. NIHR Biomedical Research Centre for Mental Health at IoP & SLAM. Neuroimaging
More informationEPI Data Are Acquired Serially. EPI Data Are Acquired Serially 10/23/2011. Functional Connectivity Preprocessing. fmri Preprocessing
Functional Connectivity Preprocessing Geometric distortion Head motion Geometric distortion Head motion EPI Data Are Acquired Serially EPI Data Are Acquired Serially descending 1 EPI Data Are Acquired
More informationI.e. Sex differences in child appetitive traits and Eating in the Absence of Hunger:
Supplementary Materials I. Evidence of sex differences on eating behavior in children I.e. Sex differences in child appetitive traits and Eating in the Absence of Hunger: Table 2. Parent Report for Child
More informationERPEEG Tutorial. Version 1.0. This tutorial was written by: Sravya Atluri, Matthew Frehlich and Dr. Faranak Farzan.
ERPEEG Tutorial Version 1.0 This tutorial was written by: Sravya Atluri, Matthew Frehlich and Dr. Faranak Farzan. Contact: faranak.farzan@sfu.ca Temerty Centre for Therapeutic Brain Stimulation Centre
More informationUpdate ASA 4.8. Expand your research potential with ASA 4.8. Highly advanced 3D display of single channel coherence
Update ASA 4.8 Expand your research potential with ASA 4.8. The ASA 4.8 software has everything needed for a complete analysis of EEG / ERP and MEG data. From features like (pre)processing of data, co-registration
More informationUnconstrained Beamforming : A Versatile Approach.
Unconstrained Beamforming : A Versatile Approach. Monika Agrawal Centre for Applied Research in Electronics, IIT Delhi October 11, 2005 Abstract Adaptive beamforming is minimizing the output power in constrained
More informationASA Getting Started. ANT BV, Enschede, Netherlands Advanced Neuro Technology
ASA Getting Started ANT BV, Enschede, Netherlands Advanced Neuro Technology www.ant-neuro.com asa@ant-neuro.com TABLE OF CONTENTS DISCLAIMER... 3 NOTICE... 3 INTRODUCTION... 5 ASA-LAB... 6 ASA-LAB RECORDING
More informationMultiVariate Bayesian (MVB) decoding of brain images
MultiVariate Bayesian (MVB) decoding of brain images Alexa Morcom Edinburgh SPM course 2015 With thanks to J. Daunizeau, K. Brodersen for slides stimulus behaviour encoding of sensorial or cognitive state?
More informationSupplementary Figure 1. Decoding results broken down for different ROIs
Supplementary Figure 1 Decoding results broken down for different ROIs Decoding results for areas V1, V2, V3, and V1 V3 combined. (a) Decoded and presented orientations are strongly correlated in areas
More informationEfficient Visual Coding: From Retina To V2
Efficient Visual Coding: From Retina To V Honghao Shan Garrison Cottrell Computer Science and Engineering UCSD La Jolla, CA 9093-0404 shanhonghao@gmail.com, gary@ucsd.edu Abstract The human visual system
More informationBackground: Bioengineering - Electronics - Signal processing. Biometry - Person Identification from brain wave detection University of «Roma TRE»
Background: Bioengineering - Electronics - Signal processing Neuroscience - Source imaging - Brain functional connectivity University of Rome «Sapienza» BCI - Feature extraction from P300 brain potential
More informationRole of Parallel Imaging in High Field Functional MRI
Role of Parallel Imaging in High Field Functional MRI Douglas C. Noll & Bradley P. Sutton Department of Biomedical Engineering, University of Michigan Supported by NIH Grant DA15410 & The Whitaker Foundation
More informationBasic fmri Design and Analysis. Preprocessing
Basic fmri Design and Analysis Preprocessing fmri Preprocessing Slice timing correction Geometric distortion correction Head motion correction Temporal filtering Intensity normalization Spatial filtering
More informationA P300-speller based on event-related spectral perturbation (ERSP) Ming, D; An, X; Wan, B; Qi, H; Zhang, Z; Hu, Y
Title A P300-speller based on event-related spectral perturbation (ERSP) Author(s) Ming, D; An, X; Wan, B; Qi, H; Zhang, Z; Hu, Y Citation The 01 IEEE International Conference on Signal Processing, Communication
More informationExtracting Coactivated Features from Multiple Data Sets
Extracting Coactivated Features from Multiple Data Sets Michael U. Gutmann and Aapo Hyvärinen Dept. of Computer Science and HIIT Dept. of Mathematics and Statistics P.O. Box 68, FIN-4 University of Helsinki,
More informationImproved rejection of artifacts from EEG data using high-order statistics and independent component analysis
Improved rejection of artifacts from EEG data using high-order statistics and independent component analysis Arnaud Delorme 1,2, Tzyy-Ping Jung 2, Terrence Sejnowski 1, Scott Makeig 2 1 Computational Neurobiology
More informationQIBA PET Amyloid BC March 11, Agenda
QIBA PET Amyloid BC March 11, 2016 - Agenda 1. QIBA Round 6 Funding a. Deadlines b. What projects can be funded, what cannot c. Discussion of projects Mechanical phantom and DRO Paul & John? Any Profile
More informationSPM8 for Basic and Clinical Investigators. Preprocessing
SPM8 for Basic and Clinical Investigators Preprocessing fmri Preprocessing Slice timing correction Geometric distortion correction Head motion correction Temporal filtering Intensity normalization Spatial
More informationDetecting Changes In Non-Isotropic Images
Detecting Changes In Non-Isotropic Images K.J. Worsley 1, M. Andermann 1, T. Koulis 1, D. MacDonald, 2 and A.C. Evans 2 August 4, 1999 1 Department of Mathematics and Statistics, 2 Montreal Neurological
More informationAccelerated MRI Techniques: Basics of Parallel Imaging and Compressed Sensing
Accelerated MRI Techniques: Basics of Parallel Imaging and Compressed Sensing Peng Hu, Ph.D. Associate Professor Department of Radiological Sciences PengHu@mednet.ucla.edu 310-267-6838 MRI... MRI has low
More informationA TEMPORAL FREQUENCY DESCRIPTION OF THE SPATIAL CORRELATION BETWEEN VOXELS IN FMRI DUE TO SPATIAL PROCESSING. Mary C. Kociuba
A TEMPORAL FREQUENCY DESCRIPTION OF THE SPATIAL CORRELATION BETWEEN VOXELS IN FMRI DUE TO SPATIAL PROCESSING by Mary C. Kociuba A Thesis Submitted to the Faculty of the Graduate School, Marquette University,
More information