Differences between revisions 16 and 17
Revision 16 as of 2010-01-11 14:38:59
Size: 1054
Comment:
Revision 17 as of 2010-01-18 14:37:15
Size: 1573
Comment:
Deletions are marked like this. Additions are marked like this.
Line 43: Line 43:








from scipy import *
import Image
def profmean(img,y1=0,ylen=5,x1=0,xlen=512):
    """
        
        
    """
    y2=y1+ylen
    xlist=range(x1,x1+xlen)
    print xlist
    p=[]
    for x in xlist:
        p.append(a[y1:y2,x,0].mean())

    return p

def load_img(filename):
    im=Image.open(filename)
    return resize(fromstring(im.tostring(),uint8),(im.size[1],im.size[0],3))
    
filename='AERONET_Hornsund.2009072.terra.250m.jpg'

a=load_img(filename)

profmean(a)

Modis-Daten:

MODIS

http://rapidfire.sci.gsfc.nasa.gov/subsets/?subset=AERONET_Hornsund.2009072.terra.250m

http://rapidfire.sci.gsfc.nasa.gov/subsets/?subset=AERONET_Hornsund.2009072&altdates

Satellitenbilder (Aqua Modis) Fram Straße : Modis Images

Bilder als Dateianhänge

  • 30.03.2003 10:40
  • 31.03.2009 11:35

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

Spektrum erzeugen mit der Funktion pylab.psd. Spektrum über eine Zeile:

   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()

from scipy import * import Image def profmean(img,y1=0,ylen=5,x1=0,xlen=512):

  • """ """ y2=y1+ylen xlist=range(x1,x1+xlen) print xlist p=[] for x in xlist:
    • p.append(a[y1:y2,x,0].mean())
    return p

def load_img(filename):

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

filename='AERONET_Hornsund.2009072.terra.250m.jpg'

a=load_img(filename)

profmean(a)

LehreWiki: OpenSource2010/Project/Project Idea2010/Konvektionszellen (last edited 2011-01-17 09:45:59 by anonymous)