正如前文所述
here
,
point.in.poly
已被弃用,取而代之的是
sf::st_intersection
并且在版本2.0-2中被删除。因此,您无法使用该函数,因为它没有定义。对此有两种解决方案:
-
代替
point.in.poly()
通过
sf::st_intersection()
.
-
定义
点in.poly
你自己。代码为
点in.poly
是:
point.in.poly <- function(dt, poly_shape) {
dt_pp <- dt[!is.na(sp::over(dt, sp::geometry(poly_shape))), ]
dt_pp@data <- data.frame(dt_pp@data, sp::over(dt_pp, poly_shape))
dt_pp@proj4string <- dt@proj4string
dt_pp
}