给定这两个栅格:
library(raster)
rr1 <- raster(ncol=10, nrow=10)
rr1[] <- runif(ncell(rr1))
rr1[1]=2
rr1[10]=2
rr1[50]=2
rr2 <- raster(ncol=250, nrow=360)
rr2[] <- runif(ncell(rr2))
第一个栅格 rr1 中的像素比第二个栅格 rr2 中的像素大。
我需要的是对 rr2 中与 rr1 =2 的像素重叠(或存在于其中)的小像素的所有值求和
你可以运行
terra::resample(rr2, rr1, sum)
。这将产生与 rr1
(范围、分辨率)完全对齐的输出,并将 rr2
的值与新像素相加(在重叠不完全的情况下进行加权)。除了求和之外还有其他选项,您可以参考帮助页面了解可能的功能。