Я хочу создать контурную карту, похожую на использование stat_density2d, но основанную на третьей переменной интенсивности. Таким образом, контуры будут относиться к интенсивности, а не к плотности точек.
Вот пример работы с использованием плиток:
install.packages("ggmap")
library(ggmap)
# Create map of Dublin
dub_map <- get_map(location = "dublin", zoom = 11)
## Define spatial point of interest (Dublin)
onepnt <- data.frame(lon=-6.270089, lat=53.341819)
#### Generate longitude and latitude points based on random distributions
with mean of 'onepnt' and given st.dev.
rlon <- rnorm(500, mean=onepnt$lon, sd=0.08)
rlat <- rnorm(500, mean=onepnt$lat, sd=0.08)
df <- cbind.data.frame(lon=rlon, lat=rlat)
# Generate intensity variable "E"
df$E <- runif(500, 0.5, 3.5)
# Generate the map
dev.off()
ggmap(dub_map, extent = "device") +
stat_summary_2d(data = df, aes(x = lon, y = lat,
z = E), fun = mean, alpha = 0.6, bins = 30) +
scale_fill_gradient(name = "E (units)", low = "green", high =
"red")
Но я ищу что-то вроде этого https://sromalewski.files.wordpress.com/2011/09/bike-share-station-suggestions-091911-913x1024.png