ggmap 上的 SpatStat 结果图答案

【问题标题】：SpatStat result plot on ggmapggmap 上的 SpatStat 结果图
【发布时间】：2016-02-17 10:22:01
【问题描述】：

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

例如，如果我计算一个 quadratcount，我如何将结果添加到 ggmap 的顶部？

感谢您的帮助，

阿诺

【问题讨论】：

请阅读How to Ask 以及如何提供reproducible example。
也尝试提供一个来自 quadratcount 的输出示例

标签： r ggmap spatstat

【解决方案1】：

不幸的是，没有简单的方法可以做到这一点，因为 spatstat 使用基本图形，而 ggplot 则完全不同。也许您可以在这种特定情况下使用此答案：https://stackoverflow.com/a/10338467/3341769

【讨论】：

嗨，Edge，我已经看到了这篇文章。不幸的是，我真的很想使用 spatstat。碰巧，您知道我如何在地图顶部显示由 spatstat 计算的结果（甚至不是来自 ggmap 的地图）？我只需要在 OpenStreetMap 等背景底图的顶部显示结果。谢谢
我同意这个答案不是很好，但这是我目前能做的最好的。最重要的事实是，没有简单的方法可以将 spatstat 基础图形与 ggmap 和 ggplot 结合起来。要在 OpenStreetMap 等之上绘制 spatstat 结果。我猜你需要转换为 sp 中的一个类（例如使用 maptools），然后在那里正确投影。很抱歉没有更具体，但我对此没有太多经验。

【解决方案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 函数得到结果后，我使用 raster 函数将其转换为栅格。然后我将其转换为数据框，然后我可以轻松地将其作为 geom_tile 对象添加到 ggmap 对象中。

【讨论】：