0

我有一些经度位置数据，我想显示几个不同研究个体中每一个的随时间变化。我也想在它上面做一个边缘瓦片状的密度图，但我需要它来显示这个密度在地图上的位置，所以我需要一个覆盖在它上面的地理地图。

我的数据看起来像这样：

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()

但是，我一直坚持绘制瓷砖密度并将地图叠加到它上面。结果应该是这样的。

期望的结果

如果您能想出如何将地图叠加到密度图上，那将有很大帮助。非常感谢。

4

0 回答 0