如何在R中将SpatialPolygons添加到气泡图中

Question

我想使用sp包中的函数bubble来绘制数据，并添加SpatialPolygons图层。我可以使用spplot轻松地绘制数据，但出于某种原因它与bubble不兼容。例如：

library(sp)
# Create SpatialPolygons
Sr1 <- Polygon(cbind(c(2,4,4,1,2),c(2,3,5,4,2)))
Srs1 <- Polygons(list(Sr1), "s1")
SpP <- SpatialPolygons(list(Srs1))
# Create SpatialPointsDataFrame
pp <- data.frame(x1=2:4,x2=2:4,att=2:4)
coordinates(pp) <- ~x1+x2
# Plot using spplot
spplot(pp, sp.layout=list("sp.polygons", SpP, fill="blue"))
# Plot using bubble
bubble(pp, sp.layout=list("sp.polygons", SpP, fill="blue")) 

一些人已经通过将SpatialPolygons转换为SpatialLines来实现此操作，而且看起来我可以通过ggplot2（请参见帖子）来完成它，但我不明白为什么气泡图不能与SpatialPolygons一起使用。有什么技巧可以让它起作用吗？

- Marie Auger-Methe

Answer 1

一种解决方法是使用自定义面板函数：

bubble(pp, "att",
       panel=function(...) {
         sp.polygons(SpP, fill="blue")
         sp:::panel.bubble(...)
       }) 

- rcs

Answer 2

我不确定为什么你的 spplot 没有起作用，但这正是 ggplot 的强项：

library(sp)
library(ggplot2)

Sr1 <- Polygon(cbind(c(2,4,4,1,2),c(2,3,5,4,2)))
Srs1 <- Polygons(list(Sr1), "s1")
SpP <- SpatialPolygons(list(Srs1))

pp <- data.frame(x1=2:4,x2=2:4,att=2:4)

spdf <- SpatialPolygonsDataFrame(SpP, data=data.frame(row.names="s1", val=1))
sp_map <- fortify(spdf)

gg <- ggplot()
gg <- gg + geom_map(data=sp_map, map=sp_map,
                    aes(x=long, y=lat, map_id=id), 
                    color="black", fill="blue")
gg <- gg + geom_point(data=pp, aes(x=x1, y=x2, size=att), color="orange")
gg <- gg + scale_size_continuous(range=c(5,10))
gg <- gg + coord_equal(ylim=c(1.5,5), xlim=c(1,4.5))
gg <- gg + labs(x=NULL, y=NULL)
gg <- gg + theme_bw()
gg <- gg + theme(panel.grid=element_blank())
gg <- gg + theme(panel.border=element_blank())
gg

（我可能会混淆泡泡的x和y坐标）

这样做可以更好地控制美学效果。

- hrbrmstr