带图层的spplot未对齐

时间:2015-03-19 13:59:50

标签: r raster

我正在尝试在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)

enter image description here

为了使我遵循本教程的图层居中:http://www.r-bloggers.com/stamen-maps-with-spplot/

我不明白为什么自动居中不起作用

1 个答案:

答案 0 :(得分:0)

这是与ggmap返回的对象的投影(或不存在)有关的问题。我在this threada tutorial中使用dismo::gmap进行了有趣的讨论。我必须更新我的博文以包含此问题。