Я пытаюсь пересечь два SpatialPolygonsDataFrames и получить SpatialPolygonsDataFrame в результате. К сожалению, используя gIntersection
функцию from rgeos
(которая работает впечатляюще быстро для пересечения полигонов), я не могу извлечь соответствующие кадры данных. Рассмотрим следующий пример:
> fracPoly <- gIntersection( toSingle, fromSingle )
> class(toSingle)
[1] "SpatialPolygonsDataFrame"
> class(fromSingle)
[1] "SpatialPolygonsDataFrame"
> class(fracPoly)
[1] "SpatialPolygons"
Я могу написать функцию-обертку, которая обрабатывает передачу data.frames
, но это будет небольшая боль, чтобы сделать всю проверку правильной, и до того, как я это сделал, я надеялся, что кто-то может подтвердить, что нет лучшего способа или указать мне на другую функцию (или вариант за gIntersection
) который позволил бы мне сохранить связанный data.frames
.
Обновить
При дальнейшем размышлении, это может быть очень преднамеренным поведением gIntersection
. В конце концов, из двух SPDF, чей data.frame вы передаете? Поэтому мне, возможно, придется написать обертку, которая объединит их.