mne-python;时频分析;论据可能使用不当;如何使tfr\u morlet图正常工作?

svmlkihl  于 2021-07-13  发布在  Java
关注(0)|答案(0)|浏览(321)

我为这个问题含糊不清而道歉。这反映了我对tfr\u morlet图形输出的行为缺乏理解。
指定。我正在研究64通道脑电图数据。我有一个刺激事件,这是我的时代开始标记。我把我的epoch分为4个较小的epoch,每个1秒,其中“epoch1”包含刺激事件和非刺激eeg数据,包括在它的1秒持续时间内。
我的工作流程如下:
加载.vhdr文件
应用筛选器
应用eeg cap蒙太奇
创建事件从刺激标记开始的纪元(epoch0)
制作从epoch0裁剪的另一个历元(epoch1)(tmin=0,tmax=1,include\u tmax=true)。这给出了我感兴趣的用于tfr\u morlet函数的epoch。
epoch1上的ica
时代一号上的莫雷特(后ica)
我已经阅读了tfrèu morlet函数一段时间,并修改了参数以查看它如何影响结果图。我的工作流程中的每一个步骤似乎都正常工作,最重要的是,通过ica成功地移除了刺激事件。
下面是调用tfr\u morlet函数时的一段代码:(您将注意到我在代码段中对基线参数所做的一些修补工作)


# TF-POWER ANALYSIS

# settings for tf graph

mode = 'percent'
title = file[:-5] + " TF-Power Analysis"
name_save = file[:-5]
freqs = np.linspace(start=1, stop=100, num=100)
cycles = 4
n_cycles = freqs / cycles
vmin_max=0.075
decim = 1
base_1, base_2 = None, None
base_1, base_2 = epoch_1_ica.tmax-0.01, epoch_1_ica.tmax
base_1, base_2 = 0, None

# plot tf graph

power = tfr_morlet(inst=epoch_1_ica, freqs=freqs, n_cycles=n_cycles, use_fft=True, decim=decim, n_jobs=1, average=False, return_itc=False, zero_mean=False, output='complex')

# if tfr_morlet(average=False)

avgpower=power.average()
tf_plot=avgpower.plot(picks=None, baseline=(base_1, base_2), mode=mode, title=title, vmin=-vmin_max, vmax=vmin_max, cmap='RdBu_r')

# else:

# tf_plot = power.plot(picks=None, baseline=(base_1, base_2), mode=mode, title=title, vmin=-vmin_max, vmax=vmin_max, cmap='RdBu_r')

# save tf graph

tf_plot.savefig(fname="EEG Data\\Graphs\\tfPower\\simple_version\\" + name_save + 
    "_decim" + str(decim) + 
    "_mode" + str(mode) + 
    "_vmin" + str(vmin_max) +
    "_base" + str(base_1) +
    "-" + str(base_2) +
    "_" + "tf_power.png",
    facecolor="white")

下面是我修改参数后的图形输出之一:

从实际数据来看,各时代的开始和结束都没有大的信号。这只出现在tfr\u morlet图中。我已经调整了tfr\u morlet基线参数值,它似乎会影响这些奇怪的信号,但它也完全修改了所有其他信号。此外,应该有更多的信号出现在较低的频率范围,但他们基本上是不存在的。这有可能是因为数据,但开头和结尾的大信号肯定不是数据的反映。
总之,我有点困了,在网上找不到什么帮助。如果有人有类似的经验,这些问题出现在你的tf图,请让我知道你的想法,我可能会如何处理这一点。
如果不共享我的其余代码就可以解决这个问题,那就更好了。我很有信心,问题出在我上面提供的代码片段上。
感谢各位在这方面的时间和专业知识。我非常期待你的评论。当然,如果我同时解决这个问题,我也会发布我的解决方案。

暂无答案!

目前还没有任何答案,快来回答吧!

相关问题