Skip to content

Instantly share code, notes, and snippets.

@agramfort
Created September 11, 2016 13:17
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save agramfort/638b9a1d38196f6bdd03cb8bace94b8e to your computer and use it in GitHub Desktop.
Save agramfort/638b9a1d38196f6bdd03cb8bace94b8e to your computer and use it in GitHub Desktop.
import os.path as op
import mne
data_path = op.join(mne.datasets.sample.data_path(), 'MEG', 'sample')
raw = mne.io.read_raw_fif(op.join(data_path, 'sample_audvis_raw.fif'))
raw.crop(0., 30.)
raw.load_data()
events = mne.read_events(op.join(data_path, 'sample_audvis_raw-eve.fif'))
picks = mne.pick_types(raw.info, meg='grad')
epochs = mne.Epochs(raw, events, [1, 2], picks=picks)
ica = mne.preprocessing.ICA()
ica.fit(epochs)
ica_fname = 'debug-ica.fif'
ica.save(ica_fname)
ica = mne.preprocessing.read_ica(ica_fname)
# ica.labels_ = dict() #Adding this solves the problem
ecg_epochs = mne.preprocessing.create_ecg_epochs(raw)
ecg_inds, ecg_scores = ica.find_bads_ecg(ecg_epochs, verbose=True)
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment