как я могу выполнять вычисления на растровом кирпиче - PullRequest
0 голосов
/ 09 июля 2020

У меня есть растровый блок, содержащий ежедневные временные ряды и данные о температуре (обобщенные ниже). Как я могу создать один растровый слой из этого растрового блока, показывающий среднее количество дней (в году), когда температура <0? </p>

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 

спасибо за любую помощь

1 Ответ

0 голосов
/ 09 июля 2020

Задавая вопрос R, всегда включайте минимальный, самодостаточный, воспроизводимый пример (его создание - также лучший способ ответить на ваши собственные вопросы!). Например,

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

Чтобы получить для каждой ячейки долю слоев со значением ниже нуля, вы можете сделать

r <- mean(x < 0)
...