1

我有一个包含经度、纬度和强度变量（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 对象有关，但我无法弄清楚如何让它工作。任何帮助，将不胜感激。

谢谢，

罗宾

4

1 回答 1

4

如果没有可重现的示例，就不可能测试您的代码，但您可能可以将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')
于 2018-01-02T08:57:54.567 回答