Как отрезать или обрезать или наполнить белым прямоугольником, плотно окружающим полигон с помощью ggplot2

Я просто пытаюсь белым заполнить область за пределами простого полигона. По какой-то причине он завораживает, рисуя странный колл через центр, как будто он думает, что это убийца-вампир или что-то в этом роде.

Я пробовал следовать этому сообщению, но что-то ушло из бананов. Я бы подумал, что это будет проще, но это окажется довольно раздражительным маленьким демоном.

Как я могу заполнить область за пределами удобного для проекции полигона, не завинчивая область внутри полигона? спасибо

# 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

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

licensed under cc by-sa 3.0 with attribution.