ax.add_patch(sea_patch)
x = [3.395833, 31.233333, 28.045556, 45.35]
y = [6.453056, 30.05, -26.204444, 2.033333]
z = [0, 90, 180, 270]
map = Basemap(projection='ortho', lat_0=0, lon_0=20, resolution='l')
map.plot(x, y, '.', latlon=1)
interpol_func = interpolate.interp2d(x, y, z, kind='linear')
newx = np.linspace(min(x), max(x))
newy = np.linspace(min(y), max(y))
X, Y = np.meshgrid(newx, newy)
Z = interpol_func(newx, newy)
map.pcolormesh(X, Y, Z, latlon=1, zorder=3)
filloceans(map, color='blue')
color matplotlib map using bicubic interpolation