我如何运行在 rasterbrick 上进行计算？

Question

我有一个包含每日时间序列和温度数据的 rasterbrick（总结如下）。我如何从这个 rasterbrick 中创建一个单一的栅格层显示温度 <0 的平均天数（每年）？

class      : RasterBrick 
dimensions : 201, 464, 93264, 5479  (nrow, ncol, ncell, nlayers)
resolution : 0.25, 0.25  (x, y)
extent     : -40.5, 75.5, 25.25, 75.5  (xmin, xmax, ymin, ymax)
crs        : +proj=longlat +datum=WGS84 +no_defs 
source     : F:/new maxent/tn_ens_spread_0.25deg_reg_1980-1994_v21.0e.nc 
names      : X1980.01.01, X1980.01.02, X1980.01.03, X1980.01.04, X1980.01.05, X1980.01.06, X1980.01.07, X1980.01.08, X1980.01.09, X1980.01.10, X1980.01.11, X1980.01.12, X1980.01.13, X1980.01.14, X1980.01.15, ... 
Date       : 1980-01-01, 1994-12-31 (min, max)
varname    : tn

Answer 1

这是一个最小的、独立的、可重现的例子

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

要为每个单元格获取值小于零的层的比例，您可以执行

r <- mean(x < 0)

我如何运行在 rasterbrick 上进行计算？

How can i run calculations on a rasterbrick?

r

netcdf

conditional-statements

r-raster

我如何 运行 在 rasterbrick 上进行计算？

How can i run calculations on a rasterbrick?

r

netcdf

conditional-statements

r-raster

我如何运行在 rasterbrick 上进行计算？