r - 将栅格添加到ggmap基本 map :在ggplot2中设置alpha(透明度)并将颜色填充为inset_raster()

我想在ggplot2中绘制带有覆盖GoogleMaps基本地图的栅格的地图。因此，我这样使用get_map()和insert_raster()：

library(ggplot2)
library(ggmap)

bm <- ggmap(get_map(location = "Bangkok", maptype = "hybrid"))

bm + inset_raster(as.raster(r), xmin = r@extent[1], xmax = r@extent[2],
                  ymin = r@extent[3], ymax = r@extent[4])

是否可以设置alpha并更改fill颜色？

结果看起来像这样：

最佳答案

没有fortify甚至更快：

阅读下面的原始帖子以获取更多信息

从this blog entry中，我发现我们可以直接在ggplot::geom_polygon()中使用空间多边形

r <- raster(system.file("external/test.grd", package="raster"))
# just to make it reproducible with ggmap we have to transform to wgs84
r <- projectRaster(r, crs = CRS("+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs"))

rtp <- rasterToPolygons(r)

bm <- ggmap(get_map(location = bbox(rtp), maptype = "hybrid", zoom = 13))
bm +
  geom_polygon(data = rtp,
               aes(x = long, y = lat, group = group,
                   fill = rep(rtp$test, each = 5)),
               size = 0,
               alpha = 0.5)  +
  scale_fill_gradientn("RasterValues", colors = topo.colors(255))

如果您只需要可视化某些内容，如何处理绘制速度

如下所述，在绘制大量像素时，这种绘制可能会变得非常缓慢。因此，您可能会考虑在转换为多边形之前减少像素数量（在大多数情况下，这实际上并不会减少地图中的信息量）。因此，raster::aggregate可用于将像素数量减少到合理的数量。

该示例显示了如何将像素数量减少4数量级（即2 * 2，水平*垂直）。有关更多信息，请参见?raster::aggregate。

r <- aggregate(r, fact = 2)
#  afterwards continue with rasterToPolygons(r)...

原始帖子：

一段时间后，我找到了解决此问题的方法。将栅格转换为多边形！然后，该想法基本上在Marc Needham's blog post之后实施。

但是，有一个缺点：大量的多边形会使ggplot变得非常慢，而您将不可避免地面对它。但是，可以通过绘制到png()（或其他）设备中来加快速度。



这是一个代码示例：

library(raster)
library(ggplot2)
library(ggmap)

r <- raster(....) # any raster you want to plot
rtp <- rasterToPolygons(r)
rtp@data$id <- 1:nrow(rtp@data)   # add id column for join

rtpFort <- fortify(rtp, data = rtp@data)
rtpFortMer <- merge(rtpFort, rtp@data, by.x = 'id', by.y = 'id')  # join data

bm <- ggmap(get_map(location = "Shanghai", maptype = "hybrid", zoom = 10))

bm + geom_polygon(data = rtpFortMer,
                  aes(x = long, y = lat, group = group, fill = layer),
                  alpha = 0.5,
                  size = 0) +  ## size = 0 to remove the polygon outlines
     scale_fill_gradientn(colours = topo.colors(255))

结果是这样的：