4

Rで（GISではとても簡単に）操作を実行する方法を理解しようとしています.

spdepパッケージからポリゴン データ セットの例を見てみましょう

library("spdep")
c <- readShapePoly(system.file("etc/shapes/columbus.shp", package="spdep")[1])
plot(c)

を使用して、論理ステートメントでポリゴンを選択できることがわかりましたsubset。例えば：

cc <- サブセット(c, c@data$POLYID<5) プロット(cc)

ここで、空間データに結合したい別のデータ フレームがあるとします。

POLYID=1:9
TO.LINK =101:109
link.data <- data.frame(POLYID=POLYID, TO.LINK=TO.LINK)

これら 2 つのデータセットを使用して、2 つの空間データ フレームを取得するにはどうすればよいですか。

  1. まず、2 番目のデータ フレームに ID を持つポリゴンで構成されます
  2. 2 つ目は、反対のセット (2 つ目のデータ フレームには存在しないポリゴン) で構成されます。

どうすればこの時点に到達できますか？

4

1 に答える 1