Avoir une conceptualisation de temps difficile celui-ci. J'ai actuellement un fichier netcdf qui contient des moyennes climatiques mensuelles de 1979-2016. Je voudrais faire des piles raster annuelles séparées de cette plus grande (c'est-à-dire passer de 1 pile raster à 38). Je ne peux pas pour la vie de moi comprendre cela. Toute suggestion sera appréciée! Les données peuvent être trouvées ici: http://nimbus.cos.uidaho.edu/abatz/DATA/vpd19792016.ncCréer de nombreux sous-ensembles de piles raster à partir d'une plus grande pile raster
Ce que j'ai commencé à ce jour est la suivante:
library(raster)
library(lubridate)
library(rgdal)
library(tidyverse)
library(tools)
file = "../data/raw/climate/aet_19792016.nc"
file_split <- file %>%
basename %>%
strsplit(split = "_") %>%
unlist
var <- file_split[1]
year <- substr(file_split[2], start = 1, stop = 4)
endyear <- substr(file_split[2], start = 5, stop = 8)
start_date <- as.Date(paste(year, "01", "01", sep = "-"))
end_date <- as.Date(paste(ifelse(year == endyear, year, endyear), "12", "31", sep = "-"))
date_seq <- seq(start_date, end_date, by = "1 month")
month_seq <- month(date_seq)
nc <- nc_open(file)
nc_att <- attributes(nc$var)$names
ncvar <- ncvar_get(nc, nc_att)
tvar <- aperm(ncvar, c(3,2,1))
proj <- "+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs +towgs84=0,0,0 "
raster <- brick(tvar, crs= proj)
extent(raster) <- c(-124.793, -67.043, 25.04186, 49.41686)
names(raster) <- paste(var, year(date_seq),
unique(month(date_seq, label = TRUE)),
sep = "_")
unstack(raster)
Merci d'avance pour vos suggestions!
This is it! Super utile - merci beaucoup d'avoir répondu. –