如何在光栅砖上运行计算？答案

【问题标题】：How can i run calculations on a rasterbrick?如何在光栅砖上运行计算？
【发布时间】：2020-07-08 22:12:55
【问题描述】：

我有一个包含每日时间序列和温度数据的光栅砖（总结如下）。如何从此栅格砖创建单个栅格图层，显示温度为

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

【问题讨论】：

标签： r conditional-statements netcdf r-raster

【解决方案1】：

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

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

要获取每个单元格中值低于零的层的比例，您可以这样做

r <- mean(x < 0)

【讨论】：