import numpy as np
import matplotlib
matplotlib.use('Agg')
from scipy.interpolate import griddata
from mpl_toolkits.basemap import Basemap, maskoceans
import matplotlib.pyplot as plt
from numpy.random import seed
# Define map bounds
xMin, xMax = 2.5, 8.0
yMin, yMax = 50.6, 53.8
# Create map
m = Basemap(projection='merc',llcrnrlon=xMin,llcrnrlat=yMin,urcrnrlon=xMax,urcrnrlat=yMax,resolution='i')
m.drawcountries(linewidth=1)
m.drawcoastlines(linewidth=1)
# Load data
x = np.array([54.325666666667,52.36,53.269444444444,55.399166666667,54.116666666667,53.614444444444,53.491666666667,53.824130555556,52.918055555556,54.03694,52.139722,52.462242867998,52.926865008825,54.853888888889,52.317222,53.240026656696,52.632430667762,52.642696895243,53.391265948394,52.505333893732,52.648187308904,52.098821802977,52.896643913235,52.457270486008,53.223000488316,52.701902388132,52.0548617826,53.411581103636,52.434561756559,52.749056395511,53.123676213651,52.067534268959,53.194409573306,52.27314817052,51.441334059998,51.224757511326,51.990941918858,51.447744494043,51.891830906739,51.960667359998,51.969031121385,51.564889021961,51.857593837453,51.449772459909,51.658528382201,51.196699902606,50.905256257898,51.497306260089])
y = np.array([2.93575,3.3416666666667,3.6277777777778,3.8102777777778,4.0122222222222,4.9602777777778,5.9416666666667,2.9452777777778,4.1502777777778,6.04167,4.436389,4.5549006792363,4.7811453228565,4.6961111111111,4.789722,4.9207907082729,5.1734739738872,4.9787572406902,5.3458010937365,4.6029300588208,5.4003881262577,5.1797058644882,5.383478899702,5.5196324030324,5.7515738887123,5.8874461671401,5.8723225499118,6.1990994508938,6.2589770334531,6.5729701105864,6.5848470019087,6.6567253619722,7.1493220605216,6.8908745111116,3.5958241584686,3.8609657214986,4.121849767852,4.342014,4.3126638323991,4.4469005114756,4.9259216999194,4.9352386335384,5.1453989235756,5.3770039280214,5.7065946674719,5.7625447234516,5.7617834850481,6.1961067840608])
z = np.array([42.7,41.84,41.8,48.35,43.4,44,40.5,45.7,50.9,42.46,-1.15,4.4,1.22,50.6,-3.35,10.79,0.8,-2.4,0.73,8.5,7.25,1.9,-1.3,-3.66,1.22,-3.35,48.16,2.9,3.6,15.82,5.18,29.07,-0.2,34.75,8.03,1.68,11.86,19.2,3.5,-4.27,-0.71,14.94,0.66,22.56,21.95,30,114.3,19.5])
x,y = m(x,y)
# target grid to interpolate to
xi = np.arange(x.min(),x.max(),50000)
yi = np.arange(y.min(),y.max(),50000)
xi,yi = np.meshgrid(xi,yi)
# interpolate
zi = griddata((x,y),z,(xi,yi),method='cubic')
# pyplot
cs = m.pcolormesh(xi,yi,zi,latlon=False,cmap=plt.cm.jet)
cbar = m.colorbar(cs)
plt.savefig("test2.png");