从mp3绘制频谱图 [英] Plot spectogram from mp3

问题描述

我正在尝试直接从python 2.7.3中的mp3文件中绘制频谱图(使用ubuntu).我可以按照下面的方式从wav文件中完成此操作.

I am trying to plot a spectogram straight from an mp3 file in python 2.7.3 (using ubuntu). I can do it from a wav file as follows.

#!/usr/bin/python
from scikits.audiolab import wavread
from pylab import *

signal, fs, enc = wavread('XC124158.wav')
specgram(signal)
show()

从mp3文件而不是wav中做同一件事的最干净的方法是什么?如果可以避免的话，我不想将所有的mp3文件都转换为wav.

What's the cleanest way to do the same thing from an mp3 file instead of a wav? I don't want to convert all the mp3 files to wav if I can avoid it.

推荐答案

我将安装Debian/Ubuntu软件包

I'd install the Debian/Ubuntu package libav-tools and call avconv to decode the mp3 to a temporary wav file:

您的另一个问题已经关闭，因此我将在这里用一个简单的带通滤波示例来扩展我的答案.在您链接的文件中，似乎大多数鸟鸣声都集中在4 kHz-5.5 kHz中.

Your other question was closed, so I'll expand my answer here a bit with a simple bandpass filtering example. In the file you linked it looks like most of the birdsong is concentrated in 4 kHz - 5.5 kHz.

import os
from subprocess import check_call
from tempfile import mktemp
from scikits.audiolab import wavread, play
from scipy.signal import remez, lfilter
from pylab import *

# convert mp3, read wav
mp3filename = 'XC124158.mp3'
wname = mktemp('.wav')
check_call(['avconv', '-i', mp3filename, wname])
sig, fs, enc = wavread(wname)
os.unlink(wname)

# bandpass filter
bands = array([0,3500,4000,5500,6000,fs/2.0]) / fs
desired = [0, 1, 0]
b = remez(513, bands, desired)
sig_filt = lfilter(b, 1, sig)
sig_filt /=  1.05 * max(abs(sig_filt)) # normalize

subplot(211)
specgram(sig, Fs=fs, NFFT=1024, noverlap=0)
axis('tight'); axis(ymax=8000)
title('Original')
subplot(212)
specgram(sig_filt, Fs=fs, NFFT=1024, noverlap=0)
axis('tight'); axis(ymax=8000)
title('Filtered')
show()

play(sig_filt, fs)

这篇关于从mp3绘制频谱图的文章就介绍到这了，希望我们推荐的答案对大家有所帮助，也希望大家多多支持IT屋！