[英]Accessing full layer name raster stack 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
例如,第一層中的數據是在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多個圖層。 有沒有辦法訪問每個圖層的完整名稱和時間戳信息?
謝謝!
您只顯示幾天一個月的時間,但是從您提供的示例看來,您可以這樣做:
# 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"
聲明:本站的技術帖子網頁,遵循CC BY-SA 4.0協議,如果您需要轉載,請注明本站網址或者原文地址。任何問題請咨詢:yoyou2525@163.com.