## page was renamed from OpenSource2010/Project/Project Idea/MODIS_Bild_Darstellen
{{attachment:modis_test.png}}


{{{#!python
import Image
from pylab import *
from mpl_toolkits.basemap import Basemap
import numpy as np

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

import scipy.ndimage as nd

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

a=load_img(filename)

a=nd.zoom(a[:,:,0],0.1)

north=80.2363
south=73.7612
east=40.4156
west=-9.3026

filename='AERONET_Hornsund.2009072.terra.250m.jgw'
cds=array(open(filename).read().split('\n')[:-1]).astype(float)

#north=cds[5]
#west=cds[4]


[yn,xn]=shape(a)

# make grid
x=linspace(west,east,xn)
y=linspace(north,south,yn)

[lons,lats] = meshgrid(x,y)



figure(2)
m = Basemap(width=2400000,height=1600000,projection='stere',lat_ts=77.0,lon_0=15.0,lat_0=77.0,resolution='l')


xm,ym=m(lons,lats)

xi=linspace(m.llcrnrx,m.urcrnrx,xn) # define the new grid
yi=linspace(m.llcrnry,m.urcrnry,yn)
modis_img_nmpg=griddata(xm.flatten(),ym.flatten(),a.flatten(),xi,yi)


#m.imshow(xm,ym,a[:,:,0]) # geht nicht!
m.imshow(modis_img_nmpg,cm.bone, interpolation='bilinear',aspect='auto')
m.drawcoastlines()
m.fillcontinents(color='gray',lake_color='aqua')
m.drawmapboundary(fill_color='aqua')

m.drawmeridians(np.arange(0,360,5),labels=[1,0,0,1])
m.drawparallels(np.arange(-90,90,5),labels=[0,1,0,1])

show()
savefig('modis_test.png')
}}}