如何用ggplot2剪断或裁剪或填充紧围绕多边形外部的矩形

Question

我只是想用白色填充一个简单多边形之外的区域。出于某种原因，它搞砸了，在中心画了一个奇怪的木桩，好像它认为它是吸血鬼杀手什么的。

我尝试关注this post，但有些东西不见了。我本以为这样会更容易，但事实证明它是一个非常暴躁的小恶魔。

如何在不弄乱多边形内部区域的情况下，将适合投影的多边形外部区域填充为白色？谢谢

# reproducible example
library(rgeos)
library(maptools)

shpct.tf <- tempfile() ; td <- tempdir()

download.file( 
    "ftp://ftp2.census.gov/geo/pvs/tiger2010st/09_Connecticut/09/tl_2010_09_state10.zip" ,
    shpct.tf ,
    mode = 'wb'
)

shpct.uz <- unzip( shpct.tf , exdir = td )

# read in connecticut
ct.shp <- readShapePoly( shpct.uz[ grep( 'shp$' , shpct.uz ) ] )

# box outside of connecticut
ct.shp.env <- gEnvelope( ct.shp )

# difference between connecticut and its box
ct.shp.diff <- gDifference( ct.shp.env , ct.shp )

# prepare both shapes for ggplot2
f.ct.shp <- fortify( ct.shp )
outside <- fortify( ct.shp.diff )


library(ggplot2)

# create all layers + projections
plot <- ggplot(data = f.ct.shp, aes(x = long, y = lat))  #start with the base-plot 
layer1 <- geom_polygon(data=f.ct.shp, aes(x=long,y=lat), fill='black')
layer2 <- geom_polygon(data=outside, aes(x=long,y=lat), fill='white')
co <- coord_map( project = "merc" )

# this works
plot + layer1 

# this does not
plot + layer1 + layer2

# this also does not
plot + layer1 + layer2 + co

Answer 1

ct.shp.diff 由四个多边形组成：

R> length(ct.shp.diff@polygons[[1]]@Polygons)
# 4

或

R> nlevels(outside$group) 
# 4

因此，您需要layer2 中的群体美学（否则 ggplot 会尝试绘制单个多边形，这会导致各部分之间出现奇怪的连接）：

layer2 <- geom_polygon(data=outside, aes(x=long, y=lat, group=group), fill='white')
plot + layer1 + layer2 + co