來自mp3的劇情譜圖

Question

我想直接從python 2.7.3（使用ubuntu）的mp3文件中繪製一個譜圖。我可以從wav文件中做到如下。

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

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

什麼是從一個MP3文件，而不是一個WAV做同樣的事情最乾淨的方法？如果我能避免它，我不想將所有的mp3文件轉換爲wav。

Answer 1

我安裝的Debian/Ubuntu的軟件包libav-tools並調用avconv到MP3到臨時wav文件進行解碼：

---

編輯：您的其他問題被關閉，所以在這裏我要擴大我的回答有一個簡單的帶通濾波例子。在你鏈接的文件中，它看起來像大多數鳥鳴集中在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)