我正在尝试在ggmap上叠加意大利边界地图。我正在使用spplot但层没有对齐。
这是我的代码:
library(raster)
library(ggmap)
library(RColorBrewer)
library(latticeExtra)
setwd("c:\\temp")
gadm<-getData('GADM', country='Italy', level=2)
bbPoints <- bbox(gadm)
gmap <- get_map(c(bbPoints), maptype='watercolor', source='stamen', crop=FALSE)
bbMap <- attr(gmap, 'bb')
height <- with(bbMap, ur.lat - ll.lat)
width <- with(bbMap, ur.lon - ll.lon)
latCenter <- with(bbMap, ll.lat + ur.lat)/2
lonCenter <- with(bbMap, ll.lon + ur.lon)/2
spplot(gadm, "PID",border="green",col=NA)+
layer(grid.raster(gmap,x=lonCenter, y=latCenter,
width=width, height=height,
default.units='native'), under = TRUE)
为了使我遵循本教程的图层居中:http://www.r-bloggers.com/stamen-maps-with-spplot/
我不明白为什么自动居中不起作用
答案 0 :(得分:0)
这是与ggmap
返回的对象的投影(或不存在)有关的问题。我在this thread和a tutorial中使用dismo::gmap
进行了有趣的讨论。我必须更新我的博文以包含此问题。