# NOT RUN {
pts = boost(sf::st_sfc(sf::st_point(c(.5,.5)),
sf::st_point(c(1.5, 1.5)),
sf::st_point(c(2.5, 2.5))))
pol = boost(sf::st_sfc(sf::st_polygon(
list(rbind(c(0,0), c(2,0), c(2,2), c(0,2), c(0,0))))))
bg_intersects(pol, pts)
# }
Run the code above in your browser using DataLab