ポリゴン内のポイントまたはポイントがあるポリゴンデータを抽出したい場合がある． 特に難しいことはないが，メモとして残しておく． library(sf) set.seed(1) nc <- system.file("shape/nc.shp", package="sf") %>% st_read nc1 <- nc[1:50, ] nc_points <…