У меня есть данные пространственных координат, и я хотел бы построить карту плотности этих координат с наложением карты. Пакет hexbin R выглядит многообещающе, но я изо всех сил пытаюсь понять, как нанести данные hexbin на карту. Я также предпочел бы придерживаться базовой графики, а не ggplot.
Вот простой пример того, что я пытаюсь сделать (это не работает):
library(hexbin)
library(sf)
library(maps)
x <- rnorm(10000, mean = 0, sd = 40)
y <- rnorm(10000, mean = 0, sd = 20)
bin <- hexbin(x, y)
pts <- st_as_sf(cbind.data.frame(x,y), coords = 1:2, crs = "+init=epsg:4326")
map(col = gray(0.95), fill=T, lwd=0.2)
plot(pts, add=TRUE, pch = 3, cex=0.5)
Здесь данные гексбина не отображаются в координатном пространстве карты.
map(col = gray(0.95), fill=T, lwd=0.2)
plot(bin, style = "colorscale", newpage=FALSE)
Есть предложения? Спасибо!