使用terra中多边形的图层边界随机采样点

Question

我只是想知道是否有一种方法可以类似于QGIS功能基于图层边界随机采样点？

我有一个非常大的栅格和一个由跨一定范围内的多个断开连接的多边形组成的矢量。我想仅在多边形的图层边界内对栅格上的随机点进行采样。

我曾经通过使用多边形剪切栅格来做到这一点，但由于我现在正在使用更大的栅格，它需要太多的计算能力，不值得再这样做。

谢谢。

Answer 1

示例数据

library(raster)
r <- rast(system.file("ex/elev.tif", package="terra"))
v <- vect(system.file("ex/lux.shp", package="terra"))[c(3,6)]

你可以做

s <- spatSample(v, 20)
e <- extract(r, s)

看

plot(r)
lines(v, lwd=3, col="gray")
points(s, col="red")