simplefeatures_bits.R

# subset out polygons from x that intersect polygons in y
# should work for other geom types too
get_intersecting <- function(x, y) {
  df     <- as.data.frame(st_intersects(x, y, sparse = FALSE))
  df$k   <- apply(df, MARGIN = 1, FUN = function(x) ifelse(any(x) == TRUE, T, F))
  x$k    <- df$k
  xout   <- filter(x, k == TRUE)
  x$k    <- NULL
  xout$k <- NULL
  return(xout)
}