6
SpatialPointsDataFrameにポリゴンをオーバーレイし、SPDFデータを保存する方法は?
私が持っているSpatialPointsDataFrameいくつかの追加データで。ポリゴン内のこれらのポイントを抽出すると同時に、SPDFオブジェクトとそれに対応するデータを保持したいと思います。 これまでのところ、私はほとんど運がなく、共通のIDを照合してマージすることに頼りましたが、これは、個々のIDSでデータをグリッド化したためにのみ機能します。 ここに簡単な例があります、私は赤い四角の中のポイントを探しています。 library(sp) set.seed(357) pts <- data.frame(x = rnorm(100), y = rnorm(100), var1 = runif(100), var2 = sample(letters, 100, replace = TRUE)) coordinates(pts) <- ~ x + y class(pts) plot(pts) axis(1); axis(2) ply <- matrix(c(-1,-1, 1,-1, 1,1, -1,1, -1,-1), ncol = 2, byrow = TRUE) ply <- SpatialPolygons(list(Polygons(list(Polygon(ply)), ID …