我如何计算sd？ as.double（x）中的错误：无法将类型“ S4”强制转换为类型“ double”的向量]

Question

有人知道我的代码有什么问题吗？我编辑了帖子，因为我没有给您数据。我要计算sd。计算均值。

这里是裁剪数据的链接：

https://drive.google.com/drive/folders/1ljT1fzaDlSmn_3j7zHshS5lrV1wBvVQD

library(raster)
library(rgdal)

#load data
setwd("C:/Users/cathe/Documents/GEOTIFF_offsetcorrected")
dbase = "C:/Users/cathe/Documents/GEOTIFF_offsetcorrected"

#polygon with crop-extend
shape_data <- readOGR("C:/Users/cathe/OneDrive/Documents/ArcGIS/Projects/Cologne/WGS_1984_UTM_Zone_32N/Cologne.shp", stringsAsFactors=FALSE)

#load tif files
ALL_FILES <- list.files(path = dbase, pattern = "*\\.tif$|*.TIF$")

#output
outfiles = file.path("C:/Users/cathe/Documents/Output"),paste0(basename(tools::file_path_sans_ext(ALL_FILES)),".tif"))


#crop
for (i in seq_along(ALL_FILES)) { 
r = crop(stack(ALL_FILES[i]), shape_data 
writeRaster(r, filename=outfiles[i],
bylayer=FALSE,
format="GTiff",
datatype= "INT1U",
options="COMPRESS=ZIP",
overwrite=TRUE)
}

#mean
mean <- mean(r)
plot(mean)
writeRaster(x = mean, filename = "mean1.tif", driver = "GeoTiff", overwrite=TRUE)

#sd
standard_dev <- sd(r)
standard_dev2 <- sd(r, na.rm =TRUE)
standard_deviation <- calc(r, sd)

有人知道我的代码有什么问题吗？我编辑了帖子，因为我没有给您数据。我要计算sd。均值的计算有效！这是裁剪的链接...

Answer 1

功能：sd(x, na.rm = FALSE)

Answer 2

您希望为RasterBrick中的每个像元计算sd。