Rの逆クリッピング（消去）？

他のオブジェクトの内側にある 部分を保存する通常のクリップとは対照的に、リバースクリップは、別のオブジェクトの境界の外側にある空間オブジェクトの部分のみを保存します。

ArcMapで逆クリップを実行していますか？ArcMapでそれを行う方法を示します。

Rでこれを行うにはどうすればよいですか？

再現可能な例（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")

ここでやりたいことは、ロンドンを除くすべての英国を救うことです。視覚的には、結果の画像の黒い形状が穴になるようにします。

シンプルな機能に対する回答：

sfパッケージはGeometry Engine Open Sourceを利用するため、st_withinなどのコマンドのリストにアクセスできます。

そのようなコマンドの1つであるst_differenceは、仕事をします：

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))

また、この記事の下部に向かって参照してください

また、st_as_sfを使用してSpをsfに強制することでも実行できます。属性は管理が難しい場合があるため、警告に注意してください！

パッケージgDifferenceからの単純なアプリケーションのようですrgeos：

> 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)

投影警告は、LondonBoroughsシェープファイルにファイルがないため.prjです。

アウトラインや別のソリッドポリゴンではなく、穴であることを確認するためだけに：

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

パーティーに少し遅れましたが、「逆」引数を使用してマスクでこれを行う簡単な方法があります。

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