小波去噪python代码_一维数据小波阈值去噪-python实现

简介

处理的一维数据类型：.asc微震文件

小波层数：5 小波基：sym8

阈值公式：

,cD1为第一层分解的细节系数，N为数据长度

阈值函数：软硬阈值折中的方法

2.代码

#模块调用

import numpy as np

import pandas as pd

import matplotlib.pyplot as plt

import math

import pywt

#封装成函数

def sgn(num):

if(num > 0.0):

return 1.0

elif(num == 0.0):

return 0.0

else:

return -1.0

def wavelet_noising(new_df):

data = new_df

data = data.values.T.tolist() # 将np.ndarray()转为列表

w = pywt.Wavelet('sym8')

# [ca3, cd3, cd2, cd1] = pywt.wavedec(data, w, level=3) # 分解波

[ca5, cd5, cd4, cd3, cd2, cd1] = pywt.wavedec(data, w, level=5) # 分解波

length1 = len(cd1)

length0 = len(data)