Una vera mappa di calore in r

Question

Vorrei fare una vera mappa di calore in R, proprio come una mappa meteorologica, tranne per il fatto che i miei dati sono molto più semplici.

Considera questi dati 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)

La z sarebbe colore.

Ecco cosa vorrebbe una mappa di calore dall'aspetto discreto per questi dati:

Come posso creare una mappa di calore in modo tale che i colori siano lisci e lo spazio 2D completo è riempito con colori levigati in base ai valori Z.

Includi il codice di esempio, non solo un link che probabilmente mi confonderà ancora di più, e probabilmente ho già visitato quel sito comunque. Grazie :)