Python 读取WAV文件

import wave

import struct

from scipy import *

from pylab import *

#读取wav文件,我这儿读了个自己用python写的音阶的wav

filename = '1.wav'

wavefile = wave.open(filename, 'r') # open for writing

#读取wav文件的四种信息的函数。期中numframes表示一共读取了几个frames,在后面要用到滴。

nchannels = wavefile.getnchannels()

sample_width = wavefile.getsampwidth()

framerate = wavefile.getframerate()

numframes = wavefile.getnframes()

print("channel",nchannels)

print("sample_width",sample_width)

print("framerate",framerate)

print("numframes",numframes)

#建一个y的数列,用来保存后面读的每个frame的amplitude。

y = zeros(numframes)

#for循环,readframe(1)每次读一个frame,取其前两位,是左声道的信息。右声道就是后两位啦。

#unpack是struct里的一个函数,用法详见http://docs.python.org/library/struct.html。简单说来就是把#packed的string转换成原来的数据,无论是什么样的数据都返回一个tuple。这里返回的是长度为一的一个

#tuple,所以我们取它的第零位。

for i in range(numframes):

val = wavefile.readframes(1)

left = val[0:2]

#right = val[2:4]

v = struct.unpack('h', left )[0]

y[i] = v

#framerate就是44100,文件初读取的值。然后本程序最关键的一步!specgram!实在太简单了。。。

Fs = framerate

specgram(y, NFFT=1024, Fs=Fs, noverlap=900)

show()

补充知识:matlab生成wav文件并用python验证

在进行频谱分析时,发现MATLAB和python读取wav文件的波形不一致,导致不能得出正确结果,为了验证MATLAB和python哪部分有问题,于是有了这篇博客。

1、需求分析

用MATLAB生成一个正弦波并保存为wav文件,然后用python读取这个wav文件画出波形,查看python读取出来的波形和matlab生成的波形是否一致。

2、代码实现

2.1 MATLAB生成wav文件

Df=5; %频率间隔

fs = 8000; %采样频率

T = 1/fs; %采样周期

N=fs/Df; %序列点数

time = (N-1)./fs; %第一个声音片段的总时长

t=0:T:(N-1)./fs;

y = sin(2*pi*200*t); %生成第一个声音片段,注意需要用db2mag()函数把dB转换成magnitude。

sound(y,fs) %可以播放声音的函数 sound()

filename = ('sinwave.wav'); %给文件取名

plot(t,y)

audiowrite(filename,y,fs) %存储.wav音频文件,在这里文件名为sinwave.wav

其中生成的正弦波波形如下图所示:

2.2 python读取wav文件

import wave

import numpy as np

import pylab as plt

#打开wav文件 ,open返回一个的是一个Wave_read类的实例,通过调用它的方法读取WAV文件的格式和数据。

f = wave.open(r"E:\练习\音频信号处理\spectrum.m\sinwave.wav","rb")

#读取格式信息

#一次性返回所有的WAV文件的格式信息,它返回的是一个组元(tuple):声道数, 量化位数(byte单位), 采

#样频率, 采样点数, 压缩类型, 压缩类型的描述。wave模块只支持非压缩的数据,因此可以忽略最后两个信息

params = f.getparams()

nchannels, sampwidth, framerate, nframes = params[:4]

#读取波形数据

#读取声音数据,传递一个参数指定需要读取的长度(以取样点为单位)

str_data = f.readframes(nframes)

f.close()

#将波形数据转换成数组

#需要根据声道数和量化单位,将读取的二进制数据转换为一个可以计算的数组

wave_data = np.fromstring(str_data,dtype=np.short)

#通过取样点数和取样频率计算出每个取样的时间。

time=np.arange(0,nframes)/framerate

print(params)

print(wave_data.max())#找出正弦波的峰值

#time 也是一个数组,与wave_data[0]或wave_data[1]配对形成系列点坐标

plt.plot(time,wave_data/32768)

plt.show()

python读取wav文件显示出来的正弦波如下图所示:

python读取wav文件的信息及求出波形峰值的结果如下图所示:

我们可以看出峰值为32767(该wav文件有16bit,且带一个符号位。215=32768215=32768 ),我们进行归一化处理后的到波形如下图所示:

由上图可以看出MATLAB生成的正弦波保存为wav文件后,python读取该wav文件。两种语言下正弦波相同。

python音频频谱_Python读取WAV音频文件 画频谱的方法相关推荐

  1. python音频频谱_Python 读取WAV音频文件 画频谱的实例

    Python 读取WAV文件 import wave import struct from scipy import * from pylab import * #读取wav文件,我这儿读了个自己用p ...

  2. python绘制音频频谱_Python 读取wav画频谱

    Python 读取wav音频文件,并用specgram画频谱 import wave import struct from scipy import * from pylab import * #读取 ...

  3. python dicom图像分割_python读取DICOM头文件的实例

    这篇文章主要介绍了关于python 读取DICOM头文件的实例,有着一定的参考价值,现在分享给大家,有需要的朋友可以参考一下 用dicompyler软件打开dicom图像,头文件如图所示: 当然也可以 ...

  4. python读取音频文件_python 读取wav 音频文件的两种方式

    JavaScript算法(冒泡排序.选择排序与插入排序) 冒泡排序.选择排序与插入排序复杂度都是二次方级别的,放在一起说吧. 介绍一些学习这三个排序方法的比较好的资料.冒泡排序看介绍的冒泡排序,选择排 ...

  5. python 读取wav 音频文件的两种方式

    python 中,常用的有两种可以读取wav音频格式的方法,如下所示: 1 import scipy 2 from scipy.io import wavfile 3 4 import soundfi ...

  6. 如何用C语言编写wav读取函数,C++读取WAV音频文件的头部数据的实现方法

    C++读取WAV音频文件的头部数据的实现方法 前言: 在这里分享一下自己的心得,希望和大家一起分享技术,如果有什么不足,还请大家指正.写出这篇目的,就是希望大家一起成长,我也相信技术之间没有高低,只有 ...

  7. python用os.system打开wav文件_使用python读取wav格式文件

    ** 使用python读取wav格式文件 ** - 基本概念 [采样频率] 即取样频率, 指每秒钟取得声音样本的次数.采样频率越高,声音的质量也就越好,声音的还原也就越真实,但同时它占的资源比较多.由 ...

  8. 音频制作与转换—WAV音频文件

    音频制作与转换-WAV音频文件 在音视频制作或转换中,常需要面临如下选择: 是否选择WAV这种音频文件格式:选择何种编码格式:如何设置采样频率(Hz.kHz).采样精度(bit).码率(kBps.比特 ...

  9. python播放本地视频_python opencv 读取本地视频文件 修改ffmpeg的方法

    Python + opencv 读取视频的三种情况: 情况一:通过摄像头采集视频 情况二:通过本地视频文件获取视频 情况三:通过摄像头录制视频,再读取录制的视频 摄像头采集.本地视频文件的读取.写视频 ...

最新文章

  1. 是什么管脚分配_FPGA是什么?里面有什么?能做什么?
  2. win8 html文件怎么打开,技术员研习win8系统html文件图标变成空白的技巧
  3. WebRTC 速成课程
  4. Oracle数据库的导入和导出的两种方式
  5. 【渝粤教育】 国家开放大学2020年春季 1107传感器与测试技术 参考试题
  6. IF-ERRORLEVEL使用方法
  7. python easygui_Python里的easygui库
  8. Python all函数 - Python零基础入门教程
  9. 从经理的角度看领域设计
  10. C#查找指定窗口的子窗口的句柄
  11. 轻松使用EasyRecovery恢复丢失照片
  12. IAR FOR STM8 简单使用教程
  13. JavaCV的摄像头实战之七:推流(带声音)
  14. excel 自定义宏所有excel可用
  15. Tamami教你孕前如何选购防辐射服
  16. The forked VM terminated without properly saying goodbye. VM crash or System.exit called?
  17. 10分钟教会你,如何用模拟量调试电机
  18. 笔记本win10宽带共享wifi热点教程
  19. L84.linux命令每日一练 -- 第11章 Linux系统管理命令 -- rpm和yum
  20. VTK移动立方体法创建多个等值面的透视3D模型

热门文章

  1. 计算机内存不足怎样解决,怎么解决电脑内存不足的问题
  2. 已知两个经纬度求距离C++代码
  3. 房子并发什么意思_梦到买房子是什么意思 梦见房子寓意什么
  4. 线上项目开源 -- WeAre(AR 相册)
  5. android传输大图片,Android中 Activity传递图片
  6. debian8 jessie 更换为国内源
  7. Ubuntu16.04中VIM的安装
  8. WIndows安装vim
  9. 微信小程序app.json文件常用配置说明
  10. CSDN博客写文章如何选择封面,看完这篇博客就够了(数千字手把手教学)