绘制覆盖密度/计数图的边缘地图

Question

我有一些经度位置数据，我希望显示其在几个不同的研究对象中随时间的变化。我还想在它上面做一个边缘的类似瓷砖的密度图，但是我需要它来显示这个密度在地图上的位置，所以我需要一个覆盖在它上面的地理地图。

我的数据看起来像这样：

SO <- data.frame(date = rep(seq(as.Date("2000/1/1"), by = "day", length.out = 365), 3),
             julian = rep(seq(1,365,1),3),          
             ind = c(rep(1,365), rep(2,365), rep(3,365)),
             longitude = c(rnorm(365, mean = 90, sd =5), rnorm(365, mean = 85, sd =2), rnorm(365, mean = 92, sd =3)))    

到目前为止，我已经设法绘制了经过平滑线和翻转坐标的经度，以便在x轴上留下经度。我的代码目前看起来像这样：

ggplot(SO, aes(x = julian, y = longitude)) +
          geom_point(aes(color=factor(ind)), size = 0.1) +
          stat_smooth(aes(group= factor(ind)), se = FALSE)+
          stat_smooth(aes(color=factor(ind))) +
          coord_flip() + scale_x_reverse()    

然而，我一直困扰着绘制瓷砖密度并将地图覆盖到它上面。结果应如下所示。

如果你能想出如何将地图叠加到已经有很大帮助的密度图上。非常感谢你。