在ggmap上的SpatStat结果图

Question

我想知道是否可以将从spatstat（https://cran.r-project.org/web/packages/spatstat/index.html）计算的结果添加到从ggmap创建的绘图的顶部。

例如，如果我计算一个quadratcount，我怎样才能在ggmap的顶部添加结果。？

感谢您的帮助，

阿诺

Answer 1

不幸的是，没有简单的方法可以做到这一点，因为spatstat使用基本图形和ggplot根本不同。也许你可以在这个具体案例中使用这个答案：https://stackoverflow.com/a/10338467/3341769

Answer 2

最后，我已经能够做我想要的了。 我是这样做的：

library(raster)
library(ggmap)
library(spatstat)

px  = spatstat::pixellate(myDataset, eps=0.004, what=c("number"))
traj.raster = raster::raster(px)
raster_df = as.data.frame(rasterToPoints(traj.raster))

从pixellate函数获取结果后，我使用栅格函数将其转换为栅格。然后我将其转换为数据框，然后我可以轻松地将其作为geom_tile对象添加到ggmap对象。