1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
|
import os.path as op
from mne import Annotations
from mne.io import read_raw_fif
from mne.preprocessing.eog import find_eog_events
data_path = op.join(op.dirname(__file__), '..', '..', 'io', 'tests', 'data')
raw_fname = op.join(data_path, 'test_raw.fif')
event_fname = op.join(data_path, 'test-eve.fif')
proj_fname = op.join(data_path, 'test-proj.fif')
def test_find_eog():
"""Test find EOG peaks."""
ch_name = 'EOG 061'
raw = read_raw_fif(raw_fname)
raw.set_annotations(Annotations([14, 21], [1, 1], 'BAD_blink'))
events = find_eog_events(raw, ch_name=ch_name)
assert len(events) == 4
assert not all(events[:, 0] < 29000)
events = find_eog_events(raw, reject_by_annotation=True, ch_name=ch_name)
assert all(events[:, 0] < 29000)
# threshold option
events_thr = find_eog_events(raw, thresh=100e-6, ch_name=ch_name)
assert len(events_thr) == 5
# test different ways to specify the EOG channel(s)
events = find_eog_events(raw, ch_name=None)
assert len(events) == 4
events = find_eog_events(raw, ch_name=['EEG 060', 'EOG 061'])
assert len(events) == 4
|
