使用ggmap中的geom_raster绘制热图

Question

我有一个包含经度，纬度和强度变量（var1.pred）的数据框。

我想在ggmap地图上绘制一个平滑的填充等高线图。我使用geom_tile绘图完成了它，但它看起来并不平滑。我减少了瓷砖尺寸，但是图表占用了大量存储空间。我想要做的是使用geom_raster在地图对象上绘制一个漂亮干净的等高线图。 使用下面的代码我使用ggplot：

完成

raster <- ggplot(idw.output, aes(x = lon, y = lat, z = var1.pred))+
geom_raster(aes(fill = var1.pred), alpha = 0.5) +
scale_fill_gradient(low = "white", high = "blue")+
geom_contour(colour = "white", binwidth = 1) +
labs(fill = "Frequency", title = "Frequency per area", x = 'Longitude', y = 'Latitude') 
raster  

这将返回图：

但是我无法弄清楚如何将它与ggmap对象结合起来。我尝试了各种各样的事情，例如：

ggmap(fr) + #fr is a map from get_map
ggplot(idw.output, aes(x = lon, y = lat, z = var1.pred))+
geom_raster(aes(fill = var1.pred), alpha = 0.5) +
scale_fill_gradient(low = "white", high = "blue")+
geom_contour(colour = "white", binwidth = 1) +
labs(fill = "Frequency", title = "Frequency per area", x = 'Longitude', y = 'Latitude')

但我收到错误：

错误：不知道如何在剧情中添加o

我知道这个问题与组合ggplot和ggmap对象有关，但我无法弄清楚如何让它工作。任何帮助将不胜感激。

谢谢，

罗宾

Answer 1

如果没有可重复的示例，则无法测试您的代码，但您可以将ggplot调用移至base_layer参数，以便继续添加geoms。

ggmap(fr, base_layer = ggplot(idw.output, aes(x = lon, y = lat, z = var1.pred))) + 
    geom_raster(aes(fill = var1.pred), alpha = 0.5) +
    scale_fill_gradient(low = "white", high = "blue")+
    geom_contour(colour = "white", binwidth = 1) +
    labs(fill = "Frequency", title = "Frequency per area", x = 'Longitude', y = 'Latitude')