r - Rで隣接する領域をマージします(空間データを集約します)?

Question

ひどく言葉遣いされた以前の質問を言い換える必要があったと思います（削除しました）。ここで別の試みです。隣接する領域に結合して、共通の境界線が消え、外側の線だけが見えるようにしたいと考えています。

再現可能な例を次に示します。

require(shapefiles)
require(sp)

xx <- readShapeSpatial(system.file("shapes/sids.shp", package="maptools")[1],
                   IDvar="FIPSNO", proj4string=CRS("+proj=longlat +ellps=clrk66"))

# show all the subregions
plot(xx)

ここで、領域 3 と 5 のみを考えてみましょう

plot(xx[c(3,5),])

これらの領域を集約するにはどうすればよいですか。実際には、私がやりたいことは、すべての国を示す大陸全体の地図を作成し、北アメリカと南アメリカを示す地図を作成するようなものです。

私には、これはかなり一般的なタスクのように見えますが、これを行うための適切な機能が今のところ見つかりません。機能を見落としているだけですか、それとも単に手動で機能させることができますか?

Answer 1

このパッケージは、この場合に使用できるデータrgeosを処理するための優れたツールを多数提供します。Spatial*

例えば：

library(rgeos)
regionOfInterest <- gUnion(xx[3,], xx[5,])

これも同じ結果になり、複数のポリゴンの場合により便利な場合があります。

regionOfInterest <- gUnionCascaded(xx[c(3,5), ])

からの結果plot(regionOfInterest):