Recorte reverso (apagamento) em R?

Um clipe reverso salva apenas a parte do seu objeto espacial que está fora dos limites de outro objeto, em oposição a um clipe comum que salva as partes que estão dentro do outro objeto.

Executando clipe reverso no ArcMap? mostra como fazer isso no ArcMap.

Como faço isso em R?

Exemplo reproduzível (em máquinas Linux):

system("wget 'https://github.com/Robinlovelace/Creating-maps-in-R/archive/master.zip' -P /tmp/")
unzip("/tmp/master.zip", exdir = "/tmp/master")
uk <- readOGR("/tmp/master/Creating-maps-in-R-master/data/", "ukbord")
lnd <- readOGR("/tmp/master/Creating-maps-in-R-master/data/", "LondonBoroughs")
plot(uk)
plot(lnd, add = T, col = "black")

O que eu quero fazer aqui é salvar todo o Reino Unido, exceto Londres. Visualmente, quero que a forma preta na imagem resultante seja um buraco.

Resposta para recursos simples:

O pacote sf baseia-se no código-fonte aberto do Geometry Engine e, como tal, pode acessar a lista de comandos como st_within etc.

Um desses comandos, st_difference, fará o trabalho:

require(sf)

# make a square simple feature
s <- rbind(c(1,1), c(1,5), c(5,5), c(5,1), c(1,1))
s.sf <-st_sfc(st_polygon(list(s)))
s.pol = st_sf(ID = "sq", s.sf)

# make a smaller triangle simple feature
s2 <- rbind(c(2,2), c(3,4), c(4,2), c(2,2))
s2.sf <-st_sfc(st_polygon(list(s2)))
s2.pol = st_sf(ID = "tr", s2.sf)

# find the 'difference', i.e. reverse of st_intersection
t <- st_difference(s.pol,s2.pol)

plot(t)

# have a look at the new geometry, a well known text format with exterior followed by hole
st_geometry(t)[[1]]
POLYGON((1 1, 1 5, 5 5, 5 1, 1 1), (2 2, 4 2, 3 4, 2 2))

veja também no final deste artigo

também pode ser feito coagindo Sp a sf com st_as_sf. Preste atenção aos avisos, pois os atributos podem ser difíceis de gerenciar!

Parece ser uma aplicação simples gDifferencedo rgeospacote:

> require(rgeos)
> ukhole = gDifference(uk, lnd)
Warning message:
In RGEOSBinTopoFunc(spgeom1, spgeom2, byid, id, "rgeos_difference") :
  spgeom1 and spgeom2 have different proj4 strings
> plot(ukhole)

O aviso de projeção ocorre porque o LondonBoroughsshapefile não possui um .prjarquivo.

Apenas para ter certeza de que é um buraco e não um contorno ou outro polígono sólido:

> gArea(lnd) + gArea(ukhole) - gArea(uk)
[1] 0

Um pouco atrasado para a festa, mas existe uma maneira simples de fazer isso com a máscara usando o argumento 'inverso';

ukhole <- mask(uk, lnd, inverse = TRUE)