tengo la siguiente base de datos, de latitud, longitud y ph de agua "tomadas" en fuentes.
long <- c(-77.0958901385, -76.7059862868, -77.0490377322, -77.0082420569, -77.029088411, -76.8875419285)
lat <- c(-11.702568503, -12.2926462918, -12.0510554271, -11.8780894556, -12.1207933935, -12.2201247629)
ph <- c(7.64, 7.66, 7.71, 8.10, 8.22, 8.29)
df1<- data.frame(long, lat, ph)
library(sf)
library(dplyr)
library(ggplot2)
library(raster)
df<- df1 %>% st_as_sf(coords = c("long","lat"),crs = 4326)
dflima<-getData("GADM",country="PER", level=3) %>% st_as_sf() %>% filter(NAME_2=="Lima")
df %>% ggplot() +
stat_density2d(aes(x=df1$long, y=df1$lat, fill =..level..),bins=5 ,geom = "polygon")+
scale_fill_distiller(palette = "Spectral")+
geom_sf(data=dflima ,col = "#000000",alpha=0)
en el grafico he intentado hacer uno de densidad, el cual no me grafica la densidad. Lo que quiero hacer es un mapa de calor donde el color rojo tome los valores altos de ph y azul los normales y bajos azules por ejemplo, quiero hacer algo parecido al siguiente grafico, el cual es un mapa de calor el cual esta hecho en otro software.
GRACIAS DE ANTEMANO POR SU AYUDA