Une vraie carte de chaleur dans R

Question

Je voudrais faire une véritable carte de chaleur dans R, un peu comme une carte météo, sauf que mes données sont beaucoup plus simples.

Considérez ces données 3D:

x <- c(1,1,1,1,1,2,2,2,2,2,3,3,3,3,3,4,4,4,4,4)
y <- c(1,2,3,4,5,1,2,3,4,5,1,2,3,4,5,1,2,3,4,5)
z <- rnorm(20)

Le Z serait de couleur.

Voici ce qu'une carte thermique discrète aimerait pour ces données:

Comment puis-je faire une carte thermique de telle sorte que les couleurs sont lisses et que l'espace 2D complet est rempli de couleurs lissées en fonction des valeurs z.

Veuillez inclure un exemple de code, pas seulement un lien qui me confondra probablement encore plus, et j'ai probablement déjà visité ce site de toute façon. Merci :)