J'essaie d'intersecter deux SpatialPolygonsDataFrames et d'obtenir un SpatialPolygonsDataFrame comme résultat. Malheureusement, en utilisant la gIntersection
fonction de rgeos
(qui fonctionne de manière impressionnante rapidement pour intersecter les polygones), je n'arrive pas à récupérer les trames de données associées. Prenons l'exemple suivant:
> fracPoly <- gIntersection( toSingle, fromSingle )
> class(toSingle)
[1] "SpatialPolygonsDataFrame"
> class(fromSingle)
[1] "SpatialPolygonsDataFrame"
> class(fracPoly)
[1] "SpatialPolygons"
Je peux écrire une fonction wrapper qui gère le transfert de data.frames
, mais ce sera une douleur mineure de bien vérifier et avant j'espérais que quelqu'un pourrait soit confirmer qu'il n'y a pas de meilleur moyen ou me diriger vers une autre fonction (ou option pour gIntersection
) ce qui me permettrait de conserver les associés data.frames
.
Mise à jour
Après réflexion, cela peut être un comportement très délibéré de la part de gIntersection
. Après tout, des deux SPDF, dont vous transmettez le data.frame? Il se peut donc que je doive écrire un wrapper qui fusionne les deux.