Python实现小波变换去噪，利用Pywavelets处理脑电信号，代码参考

小波变换及pyWavelets的运用

小波变换的数学原理
PyWavelets中几个函数的API

==Wavelet==
==Waverec==
==threshold==

代码

注意事项以及我也没弄懂的地方

小波变换的数学原理

先略了，日后补上。

PyWavelets中几个函数的API

Wavelet

Waverec

threshold

代码

读取数据

// An highlighted block
import numpy as np
import matplotlib.pyplot as plt
import pywt
import pywt.data
from mne.io import concatenate_raws, read_raw_edf, read_raw_fif
from mne import read_events
import pyedflib
import math as m

file = pyedflib.EdfReader("S001R01.edf")
#number of channels: 64 信号信息
totalSignals = file.signals_in_file #totalSignals = 64
signal_labels = file.getSignalLabels() #label names of electrode in 10-10 system
#把data转为64x19920
''' -80 removes the ending 80 zero data points'''
data = np.zeros((totalSignals, file.getNSamples()[0]-80))
for i in np.arange(totalSignals):
    data[i, :] = file.readSignal(i)[:-80]

jupyternotebook绘图设置问题

plt.rcParams['figure.figsize'] = (20.0, 8.0)

小波变换，选用小波函数sym5

	mode =pywt.Modes.smooth
    w=pywt.Wavelet('sym5')
    a = data[23]
    ca = []
    cd = []
    for i in range(4):
        (a,d) = pywt.dwt(a,w,mode)
        ca.append(a)
        cd.append(d)

    sigmod=(np.median(abs(ca[3]))/0.6475)
    T=sigmod*m.sqrt(2*m.log(len(data[23]),10))
    lamda=T/(m.log(m.exp(1)+2*(len(ca[3])-1)))
    print (lamda)

通过阈值法来对A1部分进行阈值选择，其中阈值的选择如↑

ca[3]=pywt.threshold(data=ca[3], value=lamda, mode='soft', substitute=8*lamda)
coeff_list2 = [np.array(ca[3]),np.array(cd[3]),np.array(cd[2]),np.array(cd[1]),np.array(cd[0])]

newdata=pywt.waverec(coeff_list2,'sym5')
fig = plt.figure()
ax_main = fig.add_subplot(5,1,1)
ax_main.set_title('after')
ax_main.plot(newdata)
ax_main.set_xlim(0,len(newdata)-1)

ax = fig.add_subplot(5,1,2)
ax.set_title('orgin')
ax.plot(data[23])
ax.set_xlim(0,len(data[23])-1)
plt.show()

小波逆变换，并绘图

如图所示

注意事项以及我也没弄懂的地方

这里编译器是jupyter

本文链接：https://blog.csdn.net/h_lonely/article/details/103908737

智能推荐

Spark单机版安装及配置（win10）

写在前面：本篇博文记录在win10上安装spark单机版的过程，需要安装的软件有：Java JDK、Scala、Hadoop、Spark，本次安装的版本以及版本之间的对应要求如下：在安装之前，新建一个spark文件夹，除了java之外其他软件都放在这个文件夹下，方便管理。 1 安装java 1.1 查看java是否安装及版本信息如果已经安装了java,需要查看版本，Java版本最好是1.8或...

ADB命令指令控制Android系统的WIFI连接开启/关闭状态

Android系统的usb调试已开启，连接电脑（有ADB环境） 1.在cmd中，进入shell指令模式 2.进入ROOT指令模式 3. 打开WIFI指令 4. 关闭WIFI指令...

在阿里云部署flask接口，外网无法访问问题解决

1：修改阿里云防火墙打开官网，控制台，自己的服务器，点防火墙，然后添加规则，开放端口（5000） 2：ssh连接服务器修改防火墙，我这里使用finalShell连接到阿里云服务器输入：参考文章：https://blog.csdn.net/qq_39481696/article/details/85098957...

java.io理解（一）

流包括输入流和输出流，何为输入流？何为输出流呢？二者应该怎么使用呢？本人对输入流和输出流的理解是：针对程序而言，将资源记录到程序内，就是输入流；将资源导入到程序外部介质中就是输出流。输入流用read方法读取资源到程序，输出流用write方法将资源从程序内部保存到相应的对象中。举一个比较形象的例子，把大脑比作程序，书本比作资源，人类通过read书本，将知识存入大脑中，这就是输入流的过程；相反，将大脑...

Vue 实践过程中的几个问题

前言本篇是我在使用vue过程中，遇到的几个小问题和之前不了解的东西，记录下来，希望能够帮助各位踩坑。如果喜欢的话可以点波赞，或者关注一下，希望本文可以帮到大家。本文首发于我的个人blog：obkoro1.com 本篇记录个人遇到的问题如下: 路由变化页面数据不刷新问题 setTimeout/setInterval this指向改变，无法用this访问VUe实例 setInterval路由跳转继...