Точки пересечения и многоугольники в R

Если вы выполните overlay(pts, polys), где pts - это объект SpatialPointsDataFrame, а polys - это объект SpatialPolygonsDataFrame, то вы получите вектор той же длины, что и точки, дающие строку фрейма данных многоугольников. Итак, все, что вам нужно сделать, чтобы объединить данные многоугольника во фрейм данных точек, это:

 o = overlay(pts, polys)
 pts@data = cbind(pts@data, polys[o,])

ТЕМ НЕ МЕНИЕ! Если какая-либо из ваших точек выходит за пределы всех ваших полигонов, то overlay возвращает NA, что приведет к сбою polys [o,], поэтому либо убедитесь, что все ваши точки находятся внутри полигонов, либо вам придется подумать о другом способе присвоения значения для точек вне многоугольника ...

С новым пакетом sf это теперь легко и быстро:

library(sf)
out <- st_intersection(points, poly)

Дополнительные опции

Если вы не хотите, чтобы все поля из многоугольника добавлялись к точечному объекту, просто вызовите dplyr::select() для многоугольного объекта перед:

library(magrittr)
library(dplyr)
library(sf)

poly %>% 
  select(column-name1, column-name2, etc.) -> poly

out <- st_intersection(points, poly)

Если у вас возникнут проблемы, убедитесь, что ваш многоугольник действителен:

st_is_valid(poly)

Если вы видите здесь несколько FALSE результатов, попробуйте сделать его действительным:

poly <- st_make_valid(poly) 

Обратите внимание, что эти «допустимые» функции зависят от sf установки, скомпилированной с liblwgeom.

Вы делаете это в одной строке с пакетом point.in.poly из spatialEco.

library(spatialEco)

new_shape <- point.in.poly(pts, polys)

из документации: point.in.poly «пересекает точечные и полигональные классы пространственных объектов и добавляет к точкам атрибуты полигонов».