本文整理汇总了Python中mne.io.Raw._data[meg_channels,:]方法的典型用法代码示例。如果您正苦于以下问题：Python Raw._data[meg_channels,:]方法的具体用法？Python Raw._data[meg_channels,:]怎么用？Python Raw._data[meg_channels,:]使用的例子？那么恭喜您, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类mne.io.Raw的用法示例。
在下文中一共展示了Raw._data[meg_channels,:]方法的1个代码示例，这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞，您的评价将有助于系统推荐出更棒的Python代码示例。
示例1: clean_data
# 需要导入模块: from mne.io import Raw [as 别名]
# 或者: from mne.io.Raw import _data[meg_channels,:] [as 别名]
def clean_data(self, fn_raw, W_orig, A_orig, idx_keep,
meg_channels=[]):
# import necessary modules
from mne import pick_types
from mne.io import Raw
# load data and get indices of MEG channel
meg_raw = Raw(fn_raw, preload=True)
if not np.any(meg_channels):
meg_channels = pick_types(meg_raw.info, meg=True, eeg=False,
eog=False, stim=False, exclude='bads')
# transform data into FFT- and ICA-space
meg_data = meg_raw._data[meg_channels, :]
fft_data = np.fft.fft(meg_data.astype(np.complex64))
# normalize data
fftsize = fft_data.shape[1]
fft_data = (fft_data - np.dot(np.ones((fftsize, 1)), self.dmean).transpose()) / \
np.dot(np.ones((fftsize, 1)), self.dstd).transpose()
act = np.dot(W_orig, fft_data)
# remove ICs related to signal-of-interest
A_clean = A_orig.copy()
idx_zero = range(A_clean.shape[1])
idx_zero = np.setdiff1d(idx_zero, idx_keep)
A_clean[:, idx_zero] = 0
# transform data back into FFT- and MEG-space
# fft_data = np.dot(A_clean, act)
fft_data = np.dot(A_clean, act) * np.dot(np.ones((fftsize, 1)), self.dstd).transpose() + \
np.dot(np.ones((fftsize, 1)), self.dmean).transpose()
meg_raw._data[meg_channels, :] = np.fft.ifft(fft_data).real
# save cleaned data object
fn_fourierICA = fn_raw[:fn_raw.rfind('-raw.fif')] + ',fourierICA-raw.fif'
meg_raw.save(fn_fourierICA, overwrite=True)