Rgeosは、ポリゴンを交差させるときに関連する値をドロップします

2つのSpatialPolygonsDataFramesを交差させて、結果としてSpatialPolygonsDataFrameを取得しようとしています。残念ながら、gIntersection関数rgeos（ポリゴンを交差させるために非常に高速に動作します）を使用すると、関連するデータフレームを取得できないようです。次の例を考えてみましょう。

> fracPoly <- gIntersection( toSingle, fromSingle )
> class(toSingle)
[1] "SpatialPolygonsDataFrame"
> class(fromSingle)
[1] "SpatialPolygonsDataFrame"
> class(fracPoly)
[1] "SpatialPolygons"

の転送を処理するラッパー関数を書くことはできますが、data.framesすべてのチェックを正しく行うのはささいな痛みになります。以下のためにgIntersection）これは私が関連を保持できるようになりますdata.frames。

更新

さらに考察すると、これはによる非常に意図的な動作である可能性がありますgIntersection。結局のところ、2つのSPDFのうち、誰のdata.frameを渡すのですか？そのため、2つをマージするラッパーを作成する必要があります。

一部のプロジェクトでは、同じニーズがありました。を維持するだけでなく、data.frame出力ジオメトリタイプを管理するためのコードを追加し、完全な交差ジオプロセスを得るために、いくつかのクリーニング（ジオメトリコレクションのクリーンなど）に進む必要がありました。Rでこのようなタスクを実行する必要がある場合は、RFigisGeoパッケージを試すことができます。

#install RFigisGeo
require(devtools) 
install_github("RFigisGeo", "openfigis")
require(RFigisGeo)

#compute intersection
result <- getIntersection(features1, features2)

gIntersectionの動作は、交差するデータを設計により渡すことではありません。

交差する空間オブジェクト間には一般的な一致がないため、属性に対する任意の操作には、未知のユーザーの意図に関する仮定が必要です。これが、データスロットを通過させない理由です...

... gIntesection（）の設計は、ユーザーのみがジオメトリが変更されたエンティティの属性をどうするかを知ることができるため、意図的です。ユーザーが異なると仮定することもありますが、row.names（）メカニズムで行われるように、交差するジオメトリのIDを渡す以外の一般的な解決策はありません。

驚いたことに、ラスターパッケージには交差関数があり、これも単純に交差してデータを渡します。

ラスターパッケージには、属性データも処理しようとすることでrgeosを拡張する関数がいくつかあります。この場合は、raster :: intersectおよび関数のリストを参照してください：？ "raster-package"（セクションXIV）

私がこれについて得た完全な情報：http : //r-sig-geo.2731867.n2.nabble.com/Intended-usage-of-gIntersection-td7587120.html

上記の回答が機能しなかった（私のような）人のために、ここのリンクは、ラスターの交差でこの正確なことを実行できることを説明しています。

Rで2つのポリゴンをクリップするときにすべての属性データを保持するにはどうすればよいですか？

これを使用して、SpatialPolygonsシェープファイルでSpatialPointsDataFrameをトリミングしました。元のデータを維持しながら、SpatialPointsDataFrameのトリミング/クリップバージョンを作成します。