1 / 10

Improved Artifact correction for Combined EEG/fMRI

Improved Artifact correction for Combined EEG/fMRI K.J. Mullinger 1 , P.S. Morgan 2 , R.W. Bowtell 1 1 Sir Peter Mansfield Magnetic Resonance Centre, School of Physics and Astronomy, University of Nottingham, Nottingham, UK 2 Academic Radiology, University of Nottingham, Nottingham, UK

betty_james
Download Presentation

Improved Artifact correction for Combined EEG/fMRI

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. Improved Artifact correction for Combined EEG/fMRI K.J. Mullinger1, P.S. Morgan2, R.W. Bowtell1 1 Sir Peter Mansfield Magnetic Resonance Centre, School of Physics and Astronomy, University of Nottingham, Nottingham, UK 2Academic Radiology, University of Nottingham, Nottingham, UK

  2. INTRODUCTION Simultaneous EEG and fMRI has been made possible by the development of EEG hardware and correction methods that allow artifacts generated by the scanner gradients and cardiac pulse to be characterised and then subtracted from the EEG recording [1]. Correction methods generally rely on the generation of reproducible gradient artifacts and accurate detection of cardiac R-peaks. Here, we describe the implementation of two methodological developments aimed at improving the reliability of artifact correction: synchronization of EEG sampling to the MR scanner clock [2] and use of the scanner’s physiological logging in identifying cardiac R-peaks. References: [1] Allen et al. Neuroimage 8:229-239,1998. [2] Mandelkow et al. Neuroimage, 32(3):1120-1126,2006 [3] Chia et al. JMRI, 12:678-688,2000.

  3. Figure 1: Philips VCG system METHODS: • All data were acquired with approval from the local research ethics committee. • fMRI and EEG data were acquired simultaneously. • Each study was carried out on 3 healthy volunteers • Experiments were conducted using: MRI • Philips Achieva 3.0 T MR scanner • Standard EPI sequence • 64×64 matrix • 20 slices • 3.25×3.25×3.00 mm3 voxels. • Cardiac and respiratory cycles were simultaneously recorded using the scanner’s physiological monitoring system (vector cardiogram (VCG) [3] (Figure 1) and respiratory belt) whose output is sampled at 500 Hz. EEG • BrainAmp MR EEG amplifier, Brain Vision Recorder software (Brain Products, Munich) • BrainCap MR electrode cap with 32 electrodes (5 kHz sampling rate). • The ECG electrode was placed at the base of the subject’s back. • Triggers marking the beginning of each volume acquisition were recorded on the EEG system.

  4. Study 1: • Data were recorded for 6 minutes (180 volumes) for 3 different situations: • the EEG sampling and imaging gradient waveforms were synchronised by driving the BrainAmp clock using a 5 kHz signal derived from the 10 MHz MR scanner clock (TR = 2 s); • the EEG sampling was not synchronised to the scanner clock (TR = 2 s); • the EEG and MR clocks were synchronised, but a TR of 2.0001 s which is not a multiple of the scanner clock period, was employed. • In each experiment, the time between repetitions of scanner waveforms was TR/20=100 ms, so that gradient artifacts occurred at multiples of 10 Hz. Study 2: • A visual stimulus consisting of a flashing checkerboard was presented • at 7.5 Hz or 10 Hz • for 30 cycles • with 5 s on and 5 s off. • Data were recorded both with and without synchronisation of the MR scanner clock and EEG sampling. • A TR of2.2 s was employed so that the dominant gradient artifacts occurred at multiples of 9.09 Hz, thus avoiding the frequency of the visual stimulus.

  5. ANALYSIS: • Off-line EEG signal correction for both studies was based on averaging and then subtracting gradient and pulse artifacts, as implemented in Brain Vision Analyzer [1]. • Gradient artifact correction employed an artifact template formed from the average over all TR-periods, using the scanner-generated markers. • Pulse artifact correction was based on R-peak markers derived from either the ECG or VCG traces. • The VCG is formed from a four-lead measurement [3] on the chest, which is relatively insensitive to gradient artifact and allows orthogonalisation of the pulse artifact and R-wave. • After artifact correction, data were down-sampled to 500 Hz sampling rate using cardinal splines and an anti-aliasing filter.

  6. 0.14 Synchronised, TR=2s 0.17 Not Synchronised, TR=2s Synchronised, TR=2.0001s 0.18 Table 1: The ratio of standard deviation in the signal corrected:uncorrected, averaged over all channels and subjects. Figure 2: Attenuation of signal power after gradient artifact correction averaged over all channels. Error bars show variation in attenuation across all channels. RESULTS: Study 1: Gradient Artifact Correction • Figure 2 shows the attenuation (-20log10(corrected/uncorrected)) of the EEG voltage at multiples of 10 Hz after gradient correction for one subject, indicating that synchronisation of the EEG sampling to the MR scanner clock, leads to significantly improved correction of gradient artifacts. • This reduction in residual artifact at high frequencies, which was manifested in the data from all three subjects, will be particularly advantageous for measuring gamma band activity during MR scanning. • Figure 2 also shows that even with synchronisation it is imperative to employ a TR which is a multiple of the scanner clock period in order to achieve good gradient artifact correction. • Table 1 shows the ratio of the standard deviation of the EEG signal before and after gradient artifact correction averaged in the time-domain over all channels and subjects. It shows that the improved artifact correction that can be achieved with synchronisation significantly reduces the EEG signal variance.

  7. Figure 4: VCG trace obtained straight from the scanner without any correction. Figure 3: Typical ECG trace before (top) and after (bottom) gradient correction. Pulse Artifact Correction: • ECG traces produced before and after gradient artifact correction, from the electrode placed on the back, are shown in Figure 3. • This shows that the gradient artifact must be removed from the ECG trace before it can be used for detection of R peaks. • Figure 4 shows a corresponding VCG trace which can be used immediately in pulse artifact correction and is never saturated by the gradient artifacts.

  8. Figure 5 indicates that using the VCG rather than ECG trace to define R-peak markers provides a slightly improved level of pulse artifact correction. Use of the VCG consequently offers an alternative approach to pulse artifact correction where difficulties in obtaining an adequate quality ECG trace occur. This may be a particular problem at high field due to the increased magnitude of the pulse artifact. Figure 5: EEG artifacts on TP10 due to the cardio-ballistic effect before (black) and after correction using markers derived from VCG (red) and ECG (green) traces.

  9. Study 2 • Figure 6 shows the results obtained from one representative subject during visual stimulation (similar signal behavior was found in all three subjects). • It shows the Fourier transform of the EEG signal from channel Pz acquired with and without synchronisation, using VCG-derived markers for pulse artifact correction. • Peaks corresponding to electrical activity at multiples of the stimulus frequency (arrowed) are evident in traces obtained with synchronisation (A and C), but are obscured by gradient artifact peaks when synchronisation is not employed (B and D). A B Figure 6: A and B FFT of the EEG trace, from channel Pz, after correction. Data from stimulus “on” (blue) and “off” periods (red) are shown for synchronised (A) and unsynchronised (B) data acquisition. C and D show the difference in the spectra recorded in “on” and “off” periods for synchronised (C) and unsynchronised (B) acquisition. Arrows identify frequencies of significant neuronal activity. The frequency of visual stimulation was 10Hz for this subject. C D

  10. DISCUSSION AND CONCLUSIONS • Synchronisation of the scanner and EEG clocks significantly improves gradient artifact correction, by ensuring that the gradient artifact is identically sampled across TR periods. • It is imperative to choose a TR that is a multiple of the scanner clock period to achieve the best gradient artifact correction [2]. • The VCG trace from the scanner can be used for correction of the pulse artifact and provides a useful alternative when difficulties in obtaining an adequate quality ECG trace occur. • Implementation of these methodological developments has allowed the reproducible detection of 60 Hz electrical activity which was evoked by visual stimulation during combined EEG/fMRI experiments. Acknowledgements:EPSRC and Philips Medical Systems for supporting a studentship for KM.

More Related