r - `terra` 中的 `sp::over()` 等价物

Question

是否有等效的sp::over()in package terra？要获取显示 SpatVector 的哪些几何图形覆盖另一个 SpatVector 的哪些几何图形的数据框 - 像这样，但仅使用terra：

# get a polygons map:
library(terra)
lux <- vect(system.file("ex/lux.shp", package="terra"))
plot(lux)
text(lux, lux$NAME_2)

# get points that overlay some of those polygons:
pts <- vect(cbind(c(5.8, 6, 6.2), c(49.85, 49.5, 49.6)), crs = crs(lux))
plot(pts, col = "blue", add = TRUE)

# find which points overlay which polygons:
library(sp); library(raster)
over(as(pts, "Spatial"), as(lux, "Spatial"))

#   ID_1     NAME_1 ID_2           NAME_2 AREA    POP
# 1    1   Diekirch    3          Redange  259  18664
# 2    3 Luxembourg    9 Esch-sur-Alzette  251 176820
# 3    3 Luxembourg   10       Luxembourg  237 182607

Answer 1

您的示例数据

library(terra)
lux <- vect(system.file("ex/lux.shp", package="terra"))
pts <- vect(cbind(c(5.8, 6, 6.2), c(49.85, 49.5, 49.6)), crs = crs(lux))

您可以使用extract（也在raster）

extract(lux, pts)[,-1]
#  ID_1     NAME_1 ID_2           NAME_2 AREA    POP
#1    1   Diekirch    3          Redange  259  18664
#2    3 Luxembourg    9 Esch-sur-Alzette  251 176820
#3    3 Luxembourg   10       Luxembourg  237 182607

或者你可以这样做（由lovalery对解决方案的一种变体）

i <- relate(pts, lux, "within") |> apply(2, any)
lux[i,] |> data.frame()
#  ID_1     NAME_1 ID_2           NAME_2 AREA    POP
#1    1   Diekirch    3          Redange  259  18664
#2    3 Luxembourg    9 Esch-sur-Alzette  251 176820
#3    3 Luxembourg   10       Luxembourg  237 182607

在现在更简单的版本1.5-13（当前可以像这样安装的开发版本）中install.packages('terra', repos='https://rspatial.r-universe.dev')

i <- is.related(lux, pts, "intersects")
lux[i,] |> data.frame()
#  ID_1     NAME_1 ID_2           NAME_2 AREA    POP
#1    1   Diekirch    3          Redange  259  18664
#2    3 Luxembourg    9 Esch-sur-Alzette  251 176820
#3    3 Luxembourg   10       Luxembourg  237 182607

Answer 2

请在下面使用该terra库找到一种可能的解决方案

代表

• 代码

library(terra)

values(lux)[apply(terra::relate(lux, pts, "covers"), 2, function(x) which(x == TRUE)),]

• 输出

#>    ID_1     NAME_1 ID_2           NAME_2 AREA    POP
#> 3     1   Diekirch    3          Redange  259  18664
#> 10    3 Luxembourg    9 Esch-sur-Alzette  251 176820
#> 11    3 Luxembourg   10       Luxembourg  237 182607

^{由reprex 包（v2.0.1）于 2022-01-17 创建}