Terra R 中的遮罩和裁剪功能

Question

当栅格包仍然处于活动状态时，我通常必须执行以下操作将栅格图层裁剪为不规则多边形（针对 terra 更新了一些函数）。我想知道有没有更快/更有效的方法来做同样的事情？:

kml <- vect("pathtokml.kml")
r <- rast("pathtoraster.tif")

crp <- mask(r, project(kml, r)) #fix for different CRS's
crp2 <- crop(crp, ext(kml))
plot(crp2)

Answer 1

先裁剪然后遮罩的效率更高，使用 terra 您可以一步完成此操作。

x <- crop(r, project(kml, r), mask=TRUE)