我如何从R中的堆叠栅格图层计算标准偏差？

Question

我不知道我的代码有什么问题。我使用R版本4.0.0（2020-04-24）平台：x86_64-w64-mingw32 / x64（64位）运行于：Windows 10 x64（内部版本18362）。我堆叠并裁剪了86 GeoTiff，以计算每个像素的平均值。一切都很好，但是我无法计算标准差（请参见下面的错误）。

我使用RasterStack r做

library(raster)
mean <- mean(r)

standard_dev <- sd(r)
#Fehler in as.double(x) : 
# cannot coerce type 'S4' to vector of type 'double'

standard_dev2 <- sd(r, na.rm =TRUE)
#Fehler in as.double(x) : 
#  cannot coerce type 'S4' to vector of type 'double'

Answer 1

始终包括一个独立的，可重复的示例。您可以从所使用的软件包手册中的示例开始，例如

library(raster)
b <- brick(system.file("external/rlogo.grd", package="raster"))

解决方案

x <- calc(b, sd)

x
class      : RasterLayer 
dimensions : 77, 101, 7777  (nrow, ncol, ncell)
resolution : 1, 1  (x, y)
extent     : 0, 101, 0, 77  (xmin, xmax, ymin, ymax)
crs        : +proj=merc +lon_0=0 +k=1 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs 
source     : memory
names      : layer 
values     : 0, 38.5746  (min, max)