度々すみません。
http://aidiary.hatenablog.com/entry/20111001/1317441171
に習って、とりあえず動かしてみました。(__.waveは適当な音源)

import numpy
import wave
from pylab import *

if __name__ == "__main__":
wf = wave.open("___.wav", "rb")
data = wf.readframes(wf.getnframes())
data = frombuffer(data, dtype="int16")
length = float(wf.getnframes()) / wf.getframerate() # 波形長さ（秒）

N = 512
hammingWindow = np.hamming(N)

pxx, freqs, bins, im = specgram(data, NFFT=N, Fs=wf.getframerate(), noverlap=0, window=hammingWindow)
axis([0, length, 0, wf.getframerate() / 2])
xlabel("time [second]")
ylabel("frequency [Hz]")
show()

当然、上記ウェブサイトと同様にスペクトログラムは表示されるのですが、
任意の時間、周波数における信号強度を画像（色）ではなく数値として取り出すにはどうしたら良いでしょうか。
もっというと、周波数をある値に固定して、その時間-強度変化を取得したいのです。

pxxが信号強度かと思いprintしてみましたが、二次元配列が出力されるだけで、よく理解できず・・・。

ご教示いただければ幸いです。
よろしくお願いいたします。