我想使用terra包,通过仅选择栅格值大于特定阈值的日期来子集SpatRaster。我通过生成具有3个图层的SpatRaster重新编写了代码(我的原始栅格是时态栅格)。不幸的是,我遇到了这个错误:[subset] no(valid)layer selected。然而,我也尝试使用subset()函数,同样的错误。
library(terra)
d <- rast(ncol=10, nrow=10, nlyr = 3)
values(d[[1]]) <- 50:ncell(d)
values(d[[2]]) <- 60:ncell(d)
values(d[[3]]) <- 70:ncell(d)
set.seed(0)
d <- d[[d>60]]## error [subset] no (valid) layer selected
1条答案
按热度按时间vltsax251#
示例数据
基于时间戳值的子集。
或者
要根据某些统计信息(如最大值)选择图层,可以先计算统计信息,然后再进行子集划分。