zi = griddata(x, y, z, xi, yi)
fig, ax = plt.subplots()
im = ax.contourf(xi, yi, zi)
ax.scatter(data.Lon, data.Lat, c=data.Z, s=100, vmin=zi.min(), vmax=zi.max())
fig.colorbar(im)
plt.show()
Контуры с наложением карты на неравномерную сетку в python