from pylab import *
#略
pxx, freqs, bins, im = specgram(data, NFFT=N, Fs=wf.getframerate(), noverlap=0, window=hammingWindow)
axis([0, length, 0, wf.getframerate() / 2])
show()
みたいなコードでスペクトログラムが画像で出せるのですが、これを数値データでファイルに書き出して欲しい場合はどうすればいいのでしょうか?
スタック・オーバーフローはプログラマーとプログラミングに熱心な人のためのQ&Aサイトです。すぐ登録できます。
登録してこのコミュニティに参加するpxx
にスペクトログラム、freqs
に周波数の刻み、bins
に時刻の刻みが入っているので、それらを使いやすい形式でファイルに保存すればいいと思います。
これらをプロットするには以下のようにしますnp=numpy
, plt=matplotlib.pyplot
W, T = np.meshgrid(bins, freq)
plt.pcolormesh(W, T, pxx)
このspecgram
はmatplotlib.pyplot.specgramだと思いますが、そもそもこれは図に表示するためのものなので計算するだけならscipy.signal.spectrogramの方適しているかもしれません。こちらはサンプルにあるとおり、
f, t, Sxx = signal.spectrogram(x, fs)
plt.pcolormesh(t, f, Sxx)
で表示できるようです。