Python's standard library makes it easy to read the wav file, use struct to unpack the binary data and matplotlib conveniently comes equipped with a spectrogram plot to carry out the fourier transform.

import wave, struct, pylab, matplotlib
from numpy import array

f = wave.open('Spectrogram.wav')
frames = f.readframes(f._nframes)
v = array(struct.unpack("%sh"% f._nframes*f._nchannels,
          frames), dtype='float' ).flatten()

pylab.specgram(v,Fs=f._framerate,NFFT=2048,noverlap=1024)
matplotlib.pyplot.bone()
pylab.ylim(0,6000)
pylab.show()

Good coloring. That's bone. And the lettering is something called Silian Rail.

Just for fun

I ran aphex twin's "equation" from the album windowlicker through my script just to compare the results.

There should be clearly visible face there, but I'm putting the loss in quality down to the terrible sample rate more than anything else. I'll update when I can get a better copy ripped.