スペクトログラム（PythonのSciPyで描く）をあるWAVファイルに施して、何の曲か絵を見ただけで当てる。

Python in Excelでスペクトログラムが描けるようになったのでその続き。

SciPyは音声のWAVファイルもそのまま読める。

そこであるファイルを持ってきてスペクトログラムを掛けてみよう。

まずは広帯域。

これじゃわからん。。。低い周波数だけを取り出し。

もうこれでわかる人もいるだろう。

答えは、学校のチャイムでした。

こちらの素材を使わせてもらいました。

https://g-miya.net/wordpress/se14/

使ったコード：

import numpy as np

from scipy import signal

from scipy.io import wavfile

import matplotlib.pyplot as plt

fs, x = wavfile.read("se_14_chime.wav")

f, t, Sxx = signal.spectrogram(x[:,0], fs, window=signal.windows.hamming(5000), nperseg=5000)

n = int(f.shape[0]/50)

plt.figure(figsize=(12,6))

plt.pcolormesh(t, f[:n], 10*np.log10(Sxx[:n,:]), vmin = -20, shading='gouraud', cmap='jet')

plt.ylabel('Frequency [Hz]')

plt.xlabel('Time [sec]')

plt.colorbar().set_label("[dB]")

plt.show()