无法在 ggmap 上绘制 sf 图层

Question

之前一些问题的主题 - 例如。我无法让 sf 层在 ggmap 层上渲染。

require(ggmap)
require(raster)
require(sf)

target_bbox <- c(left = -0.65, bottom = 51.2, right = 0.45, top = 51.8)
map = get_stamenmap(target_bbox, zoom = 9, maptype = 'toner-2010')

# can test this with: ggmap(map)

e = extent(list(x = c(-0.65, 0.45), y = c(51.2, 51.8)))
r = raster(volcano)
crs(r) = CRS("+proj=longlat +datum=WGS84")
p = raster::rasterToPolygons(r)
sf = st_as_sf(p) %>% st_transform(3857)

# test with: ggplot() + geom_sf(data = sf, aes(fill = layer), col = NA)

ggmap(map) + 
   geom_sf(data = sf, aes(fill = layer), col = NA, alpha = 0.4, inherit.aes = FALSE)
Coordinate system already present. Adding new coordinate system, which will replace the existing one.

坐标系冲突但没有错误，但没有图层可见。非常感谢任何建议。

This github solution - 使用自定义函数 ggmap_bbox() 从 ggmap 对象中提取边界框 - 仅成功地在上面添加了另一个警告：

Warning message:
In st_is_longlat(x) :
  bounding box has potentially an invalid value range for longlat data

Answer 1

看起来您正试图在伦敦上空绘制 sf 火山对象，但您已将其定义为 lon/lat 0,0。除此之外，您对 ggmap 然后 geom_sf 的调用工作正常。问题只是 sf 对象不在 ggmap 的边界框中。

我不得不将其更改为 'Spatial' 对象以使用 maptools::elide 移动它，然后又变回 'sf' 对象，因为我找不到简单的移动方法SF 多边形。偏移是眼球近似值。

sf_moved <- sf %>% 
  st_transform(3857) %>% # <- may be unnecessary
  as('Spatial') %>% 
  maptools::elide(.,shift = c(-80000, 6650000)) %>% 
  as('sf') %>% 
  st_set_crs(3857) %>% 
  st_transform(4326) 

ggmap(map) + 
  geom_sf(data = sf_moved, inherit.aes = FALSE, 
          col = NA, aes(fill = layer), alpha = .4) + 
  scale_fill_viridis_c()

无法在 ggmap 上绘制 sf 图层

Cannot plot sf layer on top of ggmap

r

ggmap

sf