def main():
mne.utils.set_log_level('ERROR')
layout = load_layout(MEG)
print "Reading and processing data from files.."
filenames = sys.argv[1:]
new_components = []
for idx, fname in enumerate(filenames):
print "Handling " + str(idx+1) + ". subject"
raw = mne.io.read_raw_fif(fname, preload=True)
picks = mne.pick_types(raw.info, eeg=True, meg=True)
raw.drop_channels([ch_name for ix, ch_name in
enumerate(raw.info['ch_names']) if ix not in picks])
raw.drop_channels(raw.info['bads'])
fica = FourierICA(wsize=WSIZE, n_components=COMPONENTS, maxiter=7000,
sfreq=raw.info['sfreq'], hpass=BAND[0], lpass=BAND[1])
fica.fit(raw._data)
components = get_components(fica, raw.info, len(raw.times), layout)
handle = plot_components(components, layout)
input_ = int(raw_input("Which component to use: "))
if input_ == -1:
continue
new_components.append(components[input_ - 1])
import pdb; pdb.set_trace()
print "kissa"
def main():
mne.utils.set_log_level('ERROR')
if MEG:
layout = None
else:
layout = load_layout()
filenames = sys.argv[1:]
components = []
for fname in filenames:
print "Opening " + fname
part = pickle.load(open(fname, "rb"))
components.extend(part)
plot_components(components, layout)
import pdb; pdb.set_trace()
print "kissa"
