1 def load_img(filename):
   2     im=Image.open(filename)
   3     return resize(fromstring(im.tostring(),uint8),(im.size[1],im.size[0],3))

   1 from pylab import *
   2 
   3 def load_img(filename):
   4     im=Image.open(filename)
   5     return resize(fromstring(im.tostring(),uint8),(im.size[1],im.size[0],3))
   6 
   7 filename='AERONET_Hornsund.2009072.aqua.250m.jpg'
   8 
   9 a=load_img(filename)
  10 
  11 test=a[1500,500:1500,0]
  12 psd(test,NFFT=512,Fs=250,Fc=0,detrend=mlab.detrend_none,window=mlab.window_hanning, noverlap=0)
  13 
  14 show()