加权平均值返回空栅格

问题描述 投票:0回答:1

我不明白为什么这会返回错误,但这些是我需要使用的值。

library(terra)
#> terra 1.7.78
r1 <- rast(ncols=10, nrows=10, nlyr=6, vals=0.00:0.56)
r2 <- rast(ncols=10, nrows=10, nlyr=6, vals=0.00000000:0.03338733)
r3 <- weighted.mean(r1, r2, na.rm=TRUE)
r terra
1个回答
0
投票

由于所有权重均为零,因此无法计算加权平均值。对于更合理的权重,这似乎有效:

library(terra)
set.seed(1)
r1 <- rast(ncols=10, nrows=10, nlyr=6, vals=sample(10, 600, replace=TRUE))
r2 <- rast(ncols=10, nrows=10, nlyr=6, vals=sample(10, 600, replace=TRUE))
r3 <- weighted.mean(r1, r2, na.rm=TRUE)
r3
#class       : SpatRaster 
#dimensions  : 10, 10, 1  (nrow, ncol, nlyr)
#resolution  : 36, 18  (x, y)
#extent      : -180, 180, -90, 90  (xmin, xmax, ymin, ymax)
#coord. ref. : lon/lat WGS 84 (CRS84) (OGC:CRS84) 
#source(s)   : memory
#name        :      sum 
#min value   : 2.214286 
#max value   : 8.481481 
© www.soinside.com 2019 - 2024. All rights reserved.