我有一堆数据,我在县一级绘制,without borders . 我想加入州界 . 我有一个状态shapefile(多边形),但 spplot
似乎没有任何方法可以添加到上一个 Map 之上 . 有没有办法做到这一点,而不重写面板功能,以采取两个SPDF(这似乎非常专业的可能是其他人的问题)?
这是一个可重复的例子:
library(sp)
Srs1 = Polygons(list(Polygon(cbind(c(2,4,4,1,2),c(2,3,5,4,2)))), "s1")
Srs2 = Polygons(list(Polygon(cbind(c(5,4,2,5),c(2,3,2,2)))), "s2")
county <- SpatialPolygonsDataFrame( SpatialPolygons(list(Srs1,Srs2)),
data.frame( z=1:2, row.names=c("s1","s2") ) )
SrsA <- Polygons(list(Polygon(cbind(c(3,5,5,1,3),c(3,4,6,5,3)))),"sA")
state <- SpatialPolygonsDataFrame( SpatialPolygons(list(SrsA)),
data.frame( z=1,row.names="sA" ))
spplot( county, zcol="z",col=NA )
spplot( state, add=TRUE ) # Note the add=TRUE does nothing here, but that's the spirit of what I want to accomplish
1 回答
要使用
spplot
函数进行绘图,可以使用sp.layout
参数 . 例如,创建相应布局项的列表,例如然后绘制,将上面的列表项作为列表传递给
sp.layout
参数:如果两个空间data.frames不完全重叠,则x和y限制可能不正确 . 如有必要,您可以通过从
bbox(obj)
中提取适当的限制来更正此问题 .例如,