解释 import pywt import matplotlib.pyplot as plt import numpy as np import preprocess path = r'F:\biyesheji\4\hebingmat' x_train, y_train, x_valid, y_valid, x_test, y_test = preprocess.prepro( d_path=path, length=784, number=30, normal=True, rate=[0.6, 0.2, 0.2], enc=False, enc_step=28) for i in range(0, len(x_valid)): N = 784 fs = 12000 t = np.linspace(0, 784 / fs, N, endpoint=False) wavename = 'cmor3-3' totalscal = 256 fc = pywt.central_frequency(wavename) cparam = 2 * fc * totalscal scales = cparam / np.arange(totalscal, 1, -1) [cwtmatr, frequencies] = pywt.cwt(x_valid[i], scales, wavename, 1.0 / fs) plt.contourf(t, frequencies, abs(cwtmatr)) plt.axis('off') plt.gcf().set_size_inches(784 / 100, 784 / 100) plt.gca().xaxis.set_major_locator(plt.NullLocator()) plt.gca().yaxis.set_major_locator(plt.NullLocator()) plt.subplots_adjust(top=1, bottom=0, right=1, left=0, hspace=0, wspace=0) plt.margins(0, 0) x = r'./cwt_picture/valid/' + str(i) + '-' + str(y_valid[i]) + '.jpg' plt.savefig(x)

时间: 2023-11-12 15:06:36 浏览: 149
ZIP

pywt-master.zip

这段代码主要是进行连续小波变换(CWT)并可视化结果。首先,使用 `pywt` 库导入所需的函数和模块。然后,使用 `preprocess` 库中的 `prepro` 函数加载数据集并进行预处理。接下来,使用 `for` 循环遍历验证集中的每个样本,对每个样本进行 CWT 变换。在变换过程中,设置了信号长度、采样频率、小波名称、总比例尺数、中心频率和尺度等参数。最后,使用 `matplotlib.pyplot` 库将结果可视化并保存为图片。其中 `plt.contourf` 函数用于绘制等高线图,`plt.axis` 函数用于设置坐标轴,`plt.gcf` 函数用于获取当前的 `figure` 对象,`plt.gca` 函数用于获取当前的 `axes` 对象,`plt.subplots_adjust` 函数用于调整子图的布局,`plt.margins` 函数用于设置边缘大小。最后,将图片保存到指定的文件夹中。
阅读全文

相关推荐

import pandas as pd import matplotlib.pyplot as plt import numpy as np import pywt file_name = 'E:/liuyuan/ceshi/zhongyao/Subject_1_0cmH20_norm_breaths.csv' data = pd.read_csv(file_name, skiprows=1, usecols=[0, 2], names=['Time', 'Flow']) x = list() y = list() for i in range(len(data)): x.append(float(data.values[i][0])) y.append(float(data.values[i][1])) start_index = 0 end_index = 5372 time = np.arange(start_index, end_index) flow = np.arange(start_index, end_index) time = data['Time'][start_index:end_index] flow = data['Flow'] def wavelet_filter(data): wavelet = 'db4' # 选择小波基函数 level = 5 # 小波变换的层数 # 小波变换 coeffs = pywt.wavedec(data, wavelet, level=level) threshold = np.std(coeffs[-level]) * np.sqrt(2 * np.log(len(data))) coeffs[1:] = (pywt.threshold(c, threshold, mode='soft') for c in coeffs[1:]) filtered_data = pywt.waverec(coeffs, wavelet) return filtered_data 对Flow进行小波变换滤波 filtered_flow = wavelet_filter(flow) fig, ax = plt.subplots(figsize=(10, 5)) plt.xlim(0, 60) ax.set_ylim(-0.7, 0.7) ax.set_xlabel('Time(s)', fontsize=10) ax.set_ylabel('Flow(L/s)', fontsize=10) ax.plot(time, filtered_flow, label='Filtered Flow') ax.legend() ax.grid(True, linewidth=0.3, alpha=0.5, color='gray') plt.tight_layout() # 自动调整子图的布局 plt.show()import pandas as pd import matplotlib.pyplot as plt import numpy as np import pywt file_name = 'E:/liuyuan/ceshi/zhongyao/Subject_1_0cmH20_norm_breaths.csv' data = pd.read_csv(file_name, skiprows=1, usecols=[0, 2], names=['Time', 'Flow']) x = list() y = list() for i in range(len(data)): x.append(float(data.values[i][0])) y.append(float(data.values[i][1])) start_index = 0 end_index = 5372 time = np.arange(start_index, end_index) flow = np.arange(start_index, end_index) time = data['Time'][start_index:end_index] flow = data['Flow'] def wavelet_filter(data): wavelet = 'db4' # 选择小波基函数 level = 5 # 小波变换的层数 coeffs = pywt.wavedec(data, wavelet, level=level) threshold = np.std(coeffs[-level]) * np.sqrt(2 * np.log(len(data))) coeffs[1:] = (pywt.threshold(c, threshold, mode='soft') for c in coeffs[1:]) # 逆小波变换 filtered_data = pywt.waverec(coeffs, wavelet) return filtered_data filtered_flow = wavelet_filter(flow) fig, ax = plt.subplots(figsize=(10, 5)) plt.xlim(0, 60) ax.set_ylim(-0.7, 0.7) ax.set_xlabel('Time(s)', fontsize=10) ax.set_ylabel('Flow(L/s)', fontsize=10) ax.plot(time, filtered_flow, label='Filtered Flow') ax.legend() ax.grid(True, linewidth=0.3, alpha=0.5, color='gray') plt.tight_layout() plt.show()在此代码上添加切分代码,并按照合适窗口大小切分完准确显示

优化:import numpy as np import scipy.signal as signal import scipy.io.wavfile as wavfile import pywt import matplotlib.pyplot as plt def wiener_filter(x, fs, cutoff): # 维纳滤波函数 N = len(x) freqs, Pxx = signal.periodogram(x, fs=fs) H = np.zeros(N) H[freqs <= cutoff] = 1 Pxx_smooth = np.maximum(Pxx, np.max(Pxx) * 1e-6) H_smooth = np.maximum(H, np.max(H) * 1e-6) G = H_smooth / (H_smooth + 1 / Pxx_smooth) y = np.real(np.fft.ifft(np.fft.fft(x) * G)) return y def kalman_filter(x): # 卡尔曼滤波函数 Q = np.diag([0.01, 1]) R = np.diag([1, 0.1]) A = np.array([[1, 1], [0, 1]]) H = np.array([[1, 0], [0, 1]]) x_hat = np.zeros((2, len(x))) P = np.zeros((2, 2, len(x))) x_hat[:, 0] = np.array([x[0], 0]) P[:, :, 0] = np.eye(2) for k in range(1, len(x)): x_hat[:, k] = np.dot(A, x_hat[:, k-1]) P[:, :, k] = np.dot(np.dot(A, P[:, :, k-1]), A.T) + Q K = np.dot(np.dot(P[:, :, k], H.T), np.linalg.inv(np.dot(np.dot(H, P[:, :, k]), H.T) + R)) x_hat[:, k] += np.dot(K, x[k] - np.dot(H, x_hat[:, k])) P[:, :, k] = np.dot(np.eye(2) - np.dot(K, H), P[:, :, k]) y = x_hat[0, :] return y # 读取含有噪声的语音信号 rate, data = wavfile.read("shengyin.wav") data = data.astype(float) / 32767.0 # 维纳滤波 y_wiener = wiener_filter(data, fs=rate, cutoff=1000) # 卡尔曼滤波 y_kalman = kalman_filter(data) # 保存滤波后的信号到文件中 wavfile.write("wiener_filtered.wav", rate, np.int32(y_wiener * 32767.0)) wavfile.write("kalman_filtered.wav", rate, np.int32(y_kalman * 32767.0))

import numpy as np import matplotlib.pyplot as plt import pywt from skimage import io, color # 读取灰度图像并转换为RGB图像 img_gray = io.imread('lena.png', as_gray=True) img = color.gray2rgb(img_gray) # 对图像的三个通道进行DWT变换 coeffs_r = pywt.dwt2(img[:, :, 0], 'haar') coeffs_g = pywt.dwt2(img[:, :, 1], 'haar') coeffs_b = pywt.dwt2(img[:, :, 2], 'haar') # 设置压缩比率 compress_ratio = 0.5 # 计算阈值 threshold_r = np.sort(np.abs(coeffs_r[1].ravel()))[::-1][int(compress_ratio * len(coeffs_r[1].ravel()))] threshold_g = np.sort(np.abs(coeffs_g[1].ravel()))[::-1][int(compress_ratio * len(coeffs_g[1].ravel()))] threshold_b = np.sort(np.abs(coeffs_b[1].ravel()))[::-1][int(compress_ratio * len(coeffs_b[1].ravel()))] # 对小于阈值的系数进行置零 coeffs_r = list(coeffs_r) coeffs_r[0] = np.round(coeffs_r[0]) coeffs_r[1] = np.where(np.abs(coeffs_r[1]) < threshold_r, 0, coeffs_r[1]) coeffs_r[2] = np.where(np.abs(coeffs_r[2]) < threshold_r, 0, coeffs_r[2]) coeffs_g = list(coeffs_g) coeffs_g[0] = np.round(coeffs_g[0]) coeffs_g[1] = np.where(np.abs(coeffs_g[1]) < threshold_g, 0, coeffs_g[1]) coeffs_g[2] = np.where(np.abs(coeffs_g[2]) < threshold_g, 0, coeffs_g[2]) coeffs_b = list(coeffs_b) coeffs_b[0] = np.round(coeffs_b[0]) coeffs_b[1] = np.where(np.abs(coeffs_b[1]) < threshold_b, 0, coeffs_b[1]) coeffs_b[2] = np.where(np.abs(coeffs_b[2]) < threshold_b, 0, coeffs_b[2]) # 合并三个通道的系数 coeffs = [np.stack([coeffs_r[i], coeffs_g[i], coeffs_b[i]], axis=-1) for i in range(len(coeffs_r))] # 对图像进行IDWT反变换 img_dwt = pywt.idwt2(coeffs, 'haar') # 显示原始图像和压缩后的图像 fig, axes = plt.subplots(nrows=1, ncols=2, figsize=(8, 4)) ax = axes.ravel() ax[0].imshow(img) ax[0].set_title("Original image") ax[1].imshow(img_dwt) ax[1].set_title("Compressed image (DWT)") ax[1].set_xlabel("Compression ratio: {:.2f}".format(compress_ratio)) plt.tight_layout() plt.show()

最新推荐

recommend-type

WorkerError(解决方案).md

项目中常见的问题,记录一下解决方案
recommend-type

2024-2025第一学期一上U1~3.pdf

2024-2025第一学期一上U1~3.pdf
recommend-type

Redis详解与常见问题解决方案中文最新版本

redis是一个key-value存储系统。和Memcached类似,它支持存储的value类型相对更多,包括string(字符串)、list(链表)、set(集合)、zset(sortedset --有序集合)和hash(哈希类型)。这些数据类型都支持push/pop、add/remove及取交集并集和差集及更丰富的操作,而且这些操作都是原子性的。在此基础上,redis支持各种不同方式的排序。与memcached一样,为了保证效率,数据都是缓存在内存中。区别的是redis会周期性的把更新的数据写入磁盘或者把修改操作写入追加的记录文件,并且在此基础上实现了master-slave(主从)同步。感兴趣的朋友可以过来看看
recommend-type

ASP+ACCESS航班在线定票系统设计(源代码+论文)(源代码+论文+说明文档).zip

【项目资源】:包含前端、后端、移动开发、操作系统、人工智能、物联网、信息化管理、数据库、硬件开发、大数据、课程资源、音视频、网站开发等各种技术项目的源码。包括STM32、ESP8266、PHP、QT、Linux、iOS、C++、Java、python、web、C#、EDA、proteus、RTOS等项目的源码。【项目质量】:所有源码都经过严格测试,可以直接运行。功能在确认正常工作后才上传。【适用人群】:适用于希望学习不同技术领域的小白或进阶学习者。可作为毕设项目、课程设计、大作业、工程实训或初期项目立项。【附加价值】:项目具有较高的学习借鉴价值,也可直接拿来修改复刻。对于有一定基础或热衷于研究的人来说,可以在这些基础代码上进行修改和扩展,实现其他功能。【沟通交流】:有任何使用上的问题,欢迎随时与博主沟通,博主会及时解答。鼓励下载和使用,并欢迎大家互相学习,共同进步。
recommend-type

MATLAB实现小波阈值去噪:Visushrink硬软算法对比

资源摘要信息:"本资源提供了一套基于MATLAB实现的小波阈值去噪算法代码。用户可以通过运行主文件"project.m"来执行该去噪算法,并观察到对一张256x256像素的黑白“莱娜”图片进行去噪的全过程。此算法包括了添加AWGN(加性高斯白噪声)的过程,并展示了通过Visushrink硬阈值和软阈值方法对图像去噪的对比结果。此外,该实现还包括了对图像信噪比(SNR)的计算以及将噪声图像和去噪后的图像的打印输出。Visushrink算法的参考代码由M.Kiran Kumar提供,可以在Mathworks网站上找到。去噪过程中涉及到的Lipschitz指数计算,是基于Venkatakrishnan等人的研究,使用小波变换模量极大值(WTMM)的方法来测量。" 知识点详细说明: 1. MATLAB环境使用:本代码要求用户在MATLAB环境下运行。MATLAB是一种高性能的数值计算和可视化环境,广泛应用于工程计算、算法开发和数据分析等领域。 2. 小波阈值去噪:小波去噪是信号处理中的一个技术,用于从信号中去除噪声。该技术利用小波变换将信号分解到不同尺度的子带,然后根据信号与噪声在小波域中的特性差异,通过设置阈值来消除或减少噪声成分。 3. Visushrink算法:Visushrink算法是一种小波阈值去噪方法,由Donoho和Johnstone提出。该算法的硬阈值和软阈值是两种不同的阈值处理策略,硬阈值会将小波系数小于阈值的部分置零,而软阈值则会将这部分系数缩减到零。硬阈值去噪后的信号可能有更多震荡,而软阈值去噪后的信号更为平滑。 4. AWGN(加性高斯白噪声)添加:在模拟真实信号处理场景时,通常需要对原始信号添加噪声。AWGN是一种常见且广泛使用的噪声模型,它假设噪声是均值为零、方差为N0/2的高斯分布,并且与信号不相关。 5. 图像处理:该实现包含了图像处理的相关知识,包括图像的读取、显示和噪声添加。此外,还涉及了图像去噪前后视觉效果的对比展示。 6. 信噪比(SNR)计算:信噪比是衡量信号质量的一个重要指标,反映了信号中有效信息与噪声的比例。在图像去噪的过程中,通常会计算并比较去噪前后图像的SNR值,以评估去噪效果。 7. Lipschitz指数计算:Lipschitz指数是衡量信号局部变化复杂性的一个量度,通常用于描述信号在某个尺度下的变化规律。在小波去噪过程中,Lipschitz指数可用于确定是否保留某个小波系数,因为它与信号的奇异性相关联。 8. WTMM(小波变换模量极大值):小波变换模量极大值方法是一种小波分析技术,用于检测信号中的奇异点或边缘。该技术通过寻找小波系数模量极大值的变化来推断信号的局部特征。 9. 系统开源:该资源被标记为“系统开源”,意味着该MATLAB代码及其相关文件是可以公开访问和自由使用的。开源资源为研究人员和开发者提供了学习和实验的机会,有助于知识共享和技术发展。 资源的文件结构包括"Wavelet-Based-Denoising-MATLAB-Code-master",表明用户获取的是一套完整的项目文件夹,其中包含了执行小波去噪算法所需的所有相关文件和脚本。
recommend-type

管理建模和仿真的文件

管理Boualem Benatallah引用此版本:布阿利姆·贝纳塔拉。管理建模和仿真。约瑟夫-傅立叶大学-格勒诺布尔第一大学,1996年。法语。NNT:电话:00345357HAL ID:电话:00345357https://theses.hal.science/tel-003453572008年12月9日提交HAL是一个多学科的开放存取档案馆,用于存放和传播科学研究论文,无论它们是否被公开。论文可以来自法国或国外的教学和研究机构,也可以来自公共或私人研究中心。L’archive ouverte pluridisciplinaire
recommend-type

【交互特征的影响】:分类问题中的深入探讨,如何正确应用交互特征

![【交互特征的影响】:分类问题中的深入探讨,如何正确应用交互特征](https://img-blog.csdnimg.cn/img_convert/21b6bb90fa40d2020de35150fc359908.png) # 1. 交互特征在分类问题中的重要性 在当今的机器学习领域,分类问题一直占据着核心地位。理解并有效利用数据中的交互特征对于提高分类模型的性能至关重要。本章将介绍交互特征在分类问题中的基础重要性,以及为什么它们在现代数据科学中变得越来越不可或缺。 ## 1.1 交互特征在模型性能中的作用 交互特征能够捕捉到数据中的非线性关系,这对于模型理解和预测复杂模式至关重要。例如
recommend-type

c语言从链式队列 中获取头部元素并返回其状态的函数怎么写

在 C 语言中,如果你有一个链式队列的数据结构,通常它会包含两个指针,一个指向队首(front),一个指向队尾(rear)。以下是一个简单的链式队列头部元素获取函数的示例,假设 `Queue` 是你的链式队列结构体,并且已经包含了必要的成员变量: ```c typedef struct Queue { void* data; // 存储数据的指针 struct Queue* front; // 队首指针 struct Queue* rear; // 队尾指针 } Queue; // 获取头部元素并检查是否为空(如果队列为空,返回 NULL 或适当错误值) void*
recommend-type

易语言实现画板图像缩放功能教程

资源摘要信息:"易语言是一种基于中文的编程语言,主要面向中文用户,其特点是使用中文关键词和语法结构,使得中文使用者更容易理解和编写程序。易语言画板图像缩放源码是易语言编写的程序代码,用于实现图形用户界面中的画板组件上图像的缩放功能。通过这个源码,用户可以调整画板上图像的大小,从而满足不同的显示需求。它可能涉及到的图形处理技术包括图像的获取、缩放算法的实现以及图像的重新绘制等。缩放算法通常可以分为两大类:高质量算法和快速算法。高质量算法如双线性插值和双三次插值,这些算法在图像缩放时能够保持图像的清晰度和细节。快速算法如最近邻插值和快速放大技术,这些方法在处理速度上更快,但可能会牺牲一些图像质量。根据描述和标签,可以推测该源码主要面向图形图像处理爱好者或专业人员,目的是提供一种方便易用的方法来实现图像缩放功能。由于源码文件名称为'画板图像缩放.e',可以推断该文件是一个易语言项目文件,其中包含画板组件和图像处理的相关编程代码。" 易语言作为一种编程语言,其核心特点包括: 1. 中文编程:使用中文作为编程关键字,降低了学习编程的门槛,使得不熟悉英文的用户也能够编写程序。 2. 面向对象:易语言支持面向对象编程(OOP),这是一种编程范式,它使用对象及其接口来设计程序,以提高软件的重用性和模块化。 3. 组件丰富:易语言提供了丰富的组件库,用户可以通过拖放的方式快速搭建图形用户界面。 4. 简单易学:由于语法简单直观,易语言非常适合初学者学习,同时也能够满足专业人士对快速开发的需求。 5. 开发环境:易语言提供了集成开发环境(IDE),其中包含了代码编辑器、调试器以及一系列辅助开发工具。 6. 跨平台:易语言支持在多个操作系统平台编译和运行程序,如Windows、Linux等。 7. 社区支持:易语言有着庞大的用户和开发社区,社区中有很多共享的资源和代码库,便于用户学习和解决编程中遇到的问题。 在处理图形图像方面,易语言能够: 1. 图像文件读写:支持常见的图像文件格式如JPEG、PNG、BMP等的读取和保存。 2. 图像处理功能:包括图像缩放、旋转、裁剪、颜色调整、滤镜效果等基本图像处理操作。 3. 图形绘制:易语言提供了丰富的绘图功能,包括直线、矩形、圆形、多边形等基本图形的绘制,以及文字的输出。 4. 图像缩放算法:易语言实现的画板图像缩放功能中可能使用了特定的缩放算法来优化图像的显示效果和性能。 易语言画板图像缩放源码的实现可能涉及到以下几个方面: 1. 获取画板上的图像:首先需要从画板组件中获取到用户当前绘制或已经存在的图像数据。 2. 图像缩放算法的应用:根据用户的需求,应用适当的图像缩放算法对获取的图像数据进行处理。 3. 图像重新绘制:处理后的图像数据需要重新绘制到画板上,以实现缩放后的效果。 4. 用户交互:提供用户界面,让用户能够通过按钮、滑块等控件选择缩放比例和模式,以及触发缩放操作。 5. 性能优化:为了确保图像缩放操作流畅,需要考虑代码的执行效率和资源的合理利用。 在易语言社区中,用户可以根据自己的需求修改和扩展画板图像缩放源码,或者根据提供的API进一步开发更多高级图像处理功能,从而丰富软件的功能和用户体验。
recommend-type

"互动学习:行动中的多样性与论文攻读经历"

多样性她- 事实上SCI NCES你的时间表ECOLEDO C Tora SC和NCESPOUR l’Ingén学习互动,互动学习以行动为中心的强化学习学会互动,互动学习,以行动为中心的强化学习计算机科学博士论文于2021年9月28日在Villeneuve d'Asq公开支持马修·瑟林评审团主席法布里斯·勒菲弗尔阿维尼翁大学教授论文指导奥利维尔·皮耶昆谷歌研究教授:智囊团论文联合主任菲利普·普雷教授,大学。里尔/CRISTAL/因里亚报告员奥利维耶·西格德索邦大学报告员卢多维奇·德诺耶教授,Facebook /索邦大学审查员越南圣迈IMT Atlantic高级讲师邀请弗洛里安·斯特鲁布博士,Deepmind对于那些及时看到自己错误的人...3谢谢你首先,我要感谢我的两位博士生导师Olivier和Philippe。奥利维尔,"站在巨人的肩膀上"这句话对你来说完全有意义了。从科学上讲,你知道在这篇论文的(许多)错误中,你是我可以依