`sp::over()` 相当于 `terra`

Question

terra 包中是否有 sp::over() 的等价物？获取显示 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

请使用 terra 库

在下方找到一种可能的解决方案

Reprex

代码

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 package (v2.0.1)}

创建于 2022-01-17

Answer 2

您的示例数据

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

或者像这样

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

`sp::over()` 相当于 `terra`

`sp::over()` equivalent in `terra`

r

spatial

terra