使用R中calc函数计算RasterStack像素均值时遇报错求助

我使用的RasterStack属性如下：

class      : RasterStack 
dimensions : 116, 62, 7192, 100  (nrow, ncol, ncell, nlayers)
resolution : 200, 200  (x, y)
extent     : 501829.3, 514229.3, 3838920, 3862120  (xmin, xmax, ymin, ymax)
crs        : +proj=utm +zone=36 +datum=WGS84 +units=m +no_defs 
names      : Realization_1, Realization_2, Realization_3, Realization_4, Realization_5, Realization_6, Realization_7, Realization_8, Realization_9, Realization_10, Realization_11, Realization_12, Realization_13, Realization_14, Realization_15, ... 
min values :     0.2835853,     0.2561096,     0.3094338,     0.2750784,     0.3269123,     0.3061016,     0.3261848,     0.2906344,     0.3147314,      0.2518777,      0.3202900,      0.3056327,      0.3103129,      0.3044279,      0.2978966, ... 
max values :     0.8998597,     0.9295183,     0.8803385,     0.9196700,     0.8673874,     0.8853973,     0.8788281,     0.8999460,     0.8963305,      0.9357309,      0.8737690,      0.8818156,      0.8567144,      0.8692995,      0.8821410, ...

我尝试用以下代码计算每个像素在所有图层中的平均值：

r_FI_mean = raster::calc(rstack_FI, fun = mean);

但出现了如下错误：

Error: unable to find an inherited method for function ‘calc’ for signature ‘x = "RasterStack", fun = "numeric"’

用同样的命令计算最小值或最大值时，未出现报错。

问题原因与解决方法 #

错误的核心是R无法识别fun = mean中的mean为函数——大概率是你之前的代码中给mean赋值过一个数值，导致它被覆盖成了数值类型，而非原本的base::mean函数。

解决方法1：明确指定基础包的mean函数

直接调用base::mean，避免变量名冲突：

r_FI_mean = raster::calc(rstack_FI, fun = base::mean)

解决方法2：用匿名函数包裹mean

通过匿名函数传递计算逻辑，确保R识别为函数：

r_FI_mean = raster::calc(rstack_FI, fun = function(x) mean(x, na.rm = TRUE))

注：建议加上na.rm = TRUE参数，避免因缺失值导致结果异常

替代方案：使用terra包（推荐）

raster包已逐步被terra包替代，后者性能更优，语法更简洁：

library(terra)
# 将RasterStack转换为terra的SpatRaster
rstack_FI_terra = rast(rstack_FI)
# 计算像素平均值
r_FI_mean = app(rstack_FI_terra, mean, na.rm = TRUE)

内容的提问来源于stack exchange，提问作者ΚΩΝΣΤΑΝΤΙΝΟΣ ΠΑΝΑΓΙΩΤΟΥ