提问人：Jean-Luc Dupouey 提问时间：11/4/2023 更新时间：11/16/2023 访问量：22
Terra：:p roject 中的加权和：什么是权重？
weighted sum in terra::project: what are weights?
问：
我使用该函数将一个栅格投影到另一个栅格中，该栅格位于不同的坐标系中。文档指出，对于 ，它计算所有非 NA 贡献网格像元的加权和。我不太明白这里使用的权重是什么：是在投影的模板栅格中计算的输入栅格像素的面积吗？还是别的什么？terraterraprojectprojectmethod="sum"
答：
0赞
Jason
11/16/2023
#1
文档不是很清楚，我也有同样的问题。 在这些情况下，在引擎盖下使用 gdal_warp，尽管 GDAL 文档没有提供比 terra 文档更多的信息，但我发现这个带有注释的 Github 拉取请求确实有助于澄清权重是贡献单元的区域。method = "average"terra
以下 reprex 还确认它是面积加权总和：
library(terra)
#> terra 1.7.55
r <- rast(nrows = 1, ncols = 3)
values(r) <- c(3,5,7)
r2 <- rast(nrows = 1, ncols = 2)
plot(r)
text(r, digits = 1)
lines(as.polygons(r2), col = "red")
r_resample <- resample(r, r2, method = "sum")
plot(r_resample)
text(r_resample, digits = 1)
#this agrees with manual area weighted sum calculation:
#left hand pixel
(3*1 + 5*0.5)
#> [1] 5.5
#right hand pixel
(5*0.5 + 7*1)
#> [1] 9.5
创建于 2023-11-16 with reprex v2.0.2
评论