N = int(len(z) ** 0.5)
z = z.reshape(N, N)
plt.imshow(z + 10, extent=(np.amin(x), np.amax(x), np.amin(y), np.amax(y)),
    cmap=cm.hot, norm=LogNorm())
plt.colorbar()
plt.show()
Как построить карту плотности в питоне?