访问全层名称栅格堆栈 R

Question

我有一个气候数据栅格堆栈，其中每一层代表每日降水量。每一层的名称都很长，包括每一层代表的月、年和日。

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

例如，第一层的数据是在 2015 年 3 月 10 日收集的。但是，当我使用 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 多个图层。有没有一种方法可以访问带有每个层的时间戳信息的完整层名称？

谢谢！

Answer 1

你一个月只展示几天，但是从你提供的例子来看，你似乎可以做到：

# 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"

访问全层名称栅格堆栈 R

Accessing full layer name raster stack R

r

raster

names