栏目分类:
子分类:
返回
名师互学网用户登录
快速导航关闭
当前搜索
当前分类
子分类
实用工具
热门搜索
名师互学网 > IT > 软件开发 > 后端开发 > Python

MNE SSP learning recording——Example: Environmental noise reduction from empty-room recordings

Python 更新时间: 发布时间: IT归档 最新发布 模块sitemap 名妆网 法律咨询 聚返吧 英语巴士网 伯小乐 网商动力

MNE SSP learning recording——Example: Environmental noise reduction from empty-room recordings

MNE SSP learning recording——Example: Environmental noise reduction from empty-room recordings(https://mne.tools/stable/auto_tutorials/preprocessing/50_artifact_correction_ssp.html#tut-artifact-ssp) 记录以下在这次学习中出现的问题(未解决),欢迎指出问题! 1. pycharm上编写如下代码:
import os
import numpy as np
import matplotlib.pyplot as plt
import mne
from mne.preprocessing import (create_eog_epochs, create_ecg_epochs,
                               compute_proj_ecg, compute_proj_eog)
sample_data_folder = mne.datasets.sample.data_path()
sample_data_raw_file = os.path.join(sample_data_folder, 'MEG', 'sample',
                                    'sample_audvis_raw.fif')
# sample_data_raw_file = sample_data_raw_file.replace("\", "/")
raw = mne.io.read_raw_fif(sample_data_raw_file)

system_projs = raw.info['projs']
raw.del_proj()
empty_room_file = os.path.join(sample_data_folder, 'MEG', 'sample',
                               'ernoise_raw.fif')
empty_room_raw = mne.io.read_raw_fif(empty_room_file)

empty_room_raw.del_proj()

for average in (False, True):
    empty_room_raw.plot_psd(average=average, dB=False, xscale='log')
    plt.show()

empty_room_projs = mne.compute_proj_raw(empty_room_raw, n_grad=3, n_mag=3)
mne.viz.plot_projs_topomap(empty_room_projs, colorbar=True, vlim='joint',
                           info=empty_room_raw.info)
plt.show()

fig, axs = plt.subplots(2, 3)
for idx, _projs in enumerate([system_projs, empty_room_projs[3:]]):
    mne.viz.plot_projs_topomap(_projs, axes=axs[idx], colorbar=True,
                               vlim='joint', info=empty_room_raw.info)
	plt.show()

mags = mne.pick_types(raw.info, meg='mag')
for title, projs in [('system', system_projs),
                     ('subject-specific', empty_room_projs[3:])]:
    raw.add_proj(projs, remove_existing=True)
    fig = raw.plot(proj=True, order=mags, duration=1, n_channels=2)
    fig.subplots_adjust(top=0.9)  # make room for title
    fig.suptitle('{} projectors'.format(title), size='xx-large', weight='bold')
    plt.show()
运行结果如下——(运行结果与官网一致,这里引用官网原图)


问题出现在下面这里:


我的代码在pycharm里面运行后就只有第一排的三个图,第二排是空白的,之后还以一个空白的图片,如下所示,


尝试切换pycharm sci-view设置选项也无法解决,后注意到官网这里的mne版本

果断使用pip3 install -U mne==0.24.0安装新版本的mne后运行提示没用import
pooch package,又在pip3 install pooch且imoprt pooch后运行还是没有改变,因该不是mne版本问题!做到这,我怀疑是编译器pycharm的问题,于是使用anaconda安装好mne0.24.0版本和pooch package后重新编写上述代码,还是一样的效果——
与pycharm不同在于(顺序和官网一样了)

但是这里怎么变这么多了

老问题

想着之前好像安装pydot模块也遇到过这样的伪版本问题,怕不是清华镜像的问题,上https://anaconda.org/查看搜索mne看看,

下载一大堆库…


再试试 还是如上一样的运行结果…

但是不管怎样,最终结果还是和官网一样的——

做到这里后又仔细看了几遍官网关于的该example的讲解,还是没发现端倪,这里希望是换源的影响,这里我用的是清华源,不知道有没有遇到同样问题的并成功解决该问题或者成功做出该example的大佬指点一下,感激不尽!!!
转载请注明:文章转载自 www.mshxw.com
本文地址:https://www.mshxw.com/it/503667.html
我们一直用心在做
关于我们 文章归档 网站地图 联系我们

版权所有 (c)2021-2022 MSHXW.COM

ICP备案号:晋ICP备2021003244-6号