パッケージのgIntersect
関数を使用してrgeos
、2つのSpatialPolygonsDataFrames(spdf)で交差するポリゴンを特定しました。結果は、sfdf#1のポリゴンIDを含む行名と、spdf#2のポリゴンIDを列名する論理行列です。この行列を使用して、関数を使用する前に、2つの元のspdfをサブセット化し、重複のないポリゴンを省略したいと思いますgIntersection
。論理マトリックスを、少なくともいくつかのオーバーラップがあるポリゴンを表すIDの2つのベクトル(spdfごとに1つ)に簡単に変換できます。
私の質問は-ポリゴンIDのベクトルを使用してspdfをサブセットできますか?もしそうなら、それは@data
スロットの値を正しくサブセットしますか?私は一般的にこの構文を使用してspdfをサブセット化できることを知っています:
spdf[spdf@data$myVariable[someVector], ]
これはIDには機能しません。IDがspdfリスト構造の奥深くに埋まっているためです。
spdf[r, ]
)