我必须将包含 30 个图层(每个像素 30 个值)(466 x 435 x 30)的“星星”类光栅文件转换为具有单个图层的相同类型(星星)的新光栅文件(466 x 435 x 1)。对于每个像素,其值应等于相应像素处 30 层的平均值。下面我展示了从 R studio 中截取的名为 sim 的光栅文件的屏幕截图
enter image description here
如果 stars 对象命名为 r,以下表达式计算每个像素的平均值(不包括 NA)。 st_apply 用于在维度 mean(即像素)上应用 1:2 函数。
stars
r
NA
st_apply
mean
1:2
st_apply(r,1:2,mean,na.rm = TRUE)
这是另一个示例和更多信息:
https://geobgu.xyz/r/raster-processing.html#pixel-means