` `sp::over()`等效于` `terra`‘

Question

包terra中是否有相当于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

或者你也可以这样做(用迷幻药来改变解决方案)

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