0

scipyとmatplotlibを使って音声ファイルのスペクトログラムを表示したいです。しかし下記の関数では0除算エラーが出てしまいました。

def show_spec(fn):
    sample_rate, X = wavfile.read(fn)
    print sample_rate, X[:,0].shape
    specgram(X[:,0], Fs=sample_rate)

実行結果

44100 (442991,)
python2.7/site-packages/matplotlib/axes/_axes.py:7042: RuntimeWarning: divide by zero encountered in log10
  Z = 10. * np.log10(spec)

なぜでしょうか。

2

1 件の回答 1

1

matplotlib.pyplot.specgram はデフォルトでは信号の強さを対数スケールに変換して描画します。もしある周波数の強さが 0 だとすると log10(0) となり、計算結果が -∞ 、加えて0除算例外が発生することになります。

プログラムが間違っているわけではありません。

無音部分が原因ならそれを取り除くか完全な無音ではなくするという方法もあると思いますが、 with np.errstate(divide='ignore') の中で specgram を呼び出せば抑制できるかもしれません。
試してないですし、その対応が適切かどうか私にはわかりませんが・・・。

参考: RuntimeWarning: Divide by Zero error: How to avoid? PYTHON, NUMPY - Stack Overflow

1
  • 回答ありがとうございます。 パラメータを変更したところ例外は発生しなくなりました。 2016年2月19日 0:11

この質問に回答するには、ログインする必要があります。

求めていた回答ではありませんか? のタグが付いた他の質問を参照する。