Доступ к растровому стеку полного имени слоя R

У меня есть растровый стек для климатических данных, где каждый слой представляет собой меру суточных осадков. Имена каждого слоя длинные и включают месяц, год и день, которые каждый слой представляет.

prcpSmall
class       : RasterStack 
dimensions  : 790, 812, 641480, 6  (nrow, ncol, ncell, nlayers)
resolution  : 1000, 1000  (x, y)
extent      : 1590250, 2402250, 149500, 939500  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=lcc +lat_1=25 +lat_2=60 +lat_0=42.5 +lon_0=-100 +x_0=0 +y_0=0 +ellps=WGS84 +units=m +no_defs 
names       : FIA_DAYME//03.2015.10, FIA_DAYME//03.2015.11, FIA_DAYME//03.2015.12, FIA_DAYME//03.2015.13, FIA_DAYME//03.2015.14, FIA_DAYME//03.2015.15

Например, данные в первом слое были собраны 10 марта 2015 года. Но когда я использую names(), отметка времени исчезает.

names(prcpSmall)
[1] "FIA_DAYMET_prcp_2003.2015.10" "FIA_DAYMET_prcp_2003.2015.11"
[3] "FIA_DAYMET_prcp_2003.2015.12" "FIA_DAYMET_prcp_2003.2015.13"
[5] "FIA_DAYMET_prcp_2003.2015.14" "FIA_DAYMET_prcp_2003.2015.15"

Полный набор данных включает более 4000 слоев. Есть ли способ получить доступ к полному имени слоя с информацией о времени для каждого слоя?

Спасибо!


r names raster
person NR503    schedule 29.01.2018    source источник
comment
Добро пожаловать в SO! Не могли бы вы опубликовать образцы данных через dput(prcpSmall).   -  person Ralf Stubner    schedule 29.01.2018

Ответы (1)


arrow_upward
0
arrow_downward

Вы показываете только несколько дней в течение одного месяца, но из приведенного вами примера может показаться, что вы можете:

# n <- names(prcpSmall)
n <- c("FIA_DAYMET_prcp_2003.2015.10", "FIA_DAYMET_prcp_2003.2015.11", "FIA_DAYMET_prcp_2003.2015.12", "FIA_DAYMET_prcp_2003.2015.13", "FIA_DAYMET_prcp_2003.2015.14", "FIA_DAYMET_prcp_2003.2015.15")

d <- gsub("FIA_DAYMET_prcp_20", "", n)

as.Date(d, "%m.%Y.%d")
# [1] "2015-03-10" "2015-03-11" "2015-03-12" "2015-03-13" "2015-03-14" "2015-03-15"
person Robert Hijmans    schedule 29.01.2018
comment
Этот код отлично работал. Спасибо, что поделились функцией gsub, я ее раньше не видел. - person NR503; 30.01.2018