非官方中文翻译 本页译自 MNE-Python 官方文档： 查看英文原文 ↗。 文档版权归 MNE Developers 所有（BSD-3-Clause）。

高频 SEF 数据集

本示例查看高频体感诱发场（SEF）响应。

# Author: Jussi Nurminen (jnu@iki.fi)
#
# License: BSD-3-Clause
# Copyright the MNE-Python contributors.

import os

import mne
from mne.datasets import hf_sef

fname_evoked = os.path.join(hf_sef.data_path(), "MEG/subject_b/hf_sef_15min-ave.fif")

print(__doc__)

读取诱发数据

evoked = mne.Evoked(fname_evoked)

    Found the data of interest:
        t =     -50.00 ...     250.00 ms (SEF)
        0 CTF compensation matrices available
        nave = 2790 - aspect type = 100
No projector specified for this dataset. Please consider the method self.add_proj.

创建高通滤波版本

evoked_hp = evoked.copy()
evoked_hp.filter(l_freq=300, h_freq=None)

Setting up high-pass filter at 3e+02 Hz

FIR filter parameters
---------------------
Designing a one-pass, zero-phase, non-causal highpass filter:
- Windowed time-domain design (firwin) method
- Hamming window with 0.0194 passband ripple and 53 dB stopband attenuation
- Lower passband edge: 300.00
- Lower transition bandwidth: 75.00 Hz (-6 dB cutoff frequency: 262.50 Hz)
- Filter length: 133 samples (0.044 s)

	General
	Filename(s)	hf_sef_15min-ave.fif
	MNE object type	Evoked
	Measurement date	1970-01-01 at 00:00:00 UTC
	Participant	Unknown
	Experimenter	Jussi Nurminen (jussi)
	Acquisition
	Aggregation	average of 2790 epochs
	条件	SEF
	Time range	-0.050 – 0.250 s
	Baseline	off
	Sampling frequency	3000.00 Hz
	Time points	901
	Channels
	Magnetometers	102
	Gradiometers	204
	Head & sensor digitization	93 points
	Filters
	Highpass	300.00 Hz
	Lowpass	990.00 Hz

在单个通道上比较高通滤波与未滤波数据

ch = "MEG0443"
pick = evoked.ch_names.index(ch)
edi = {"HF": evoked_hp, "Regular": evoked}
mne.viz.plot_compare_evokeds(edi, picks=pick)