如何按时间索引SpatRaster对象的单个图层?
基于时间属性查询SpatRaster图层失败问题
之前偶然发现可基于时间属性查询SpatRaster对象的图层,例如按年份r["2017"]、年月r["2017-10"]查询。但现在尝试复现该操作时,始终报错:[subset] no (valid) layer selected。
我查阅了?time以及rspatial.org的文档,未找到与该索引功能相关的说明,请问我忽略了什么?
library(terra) #> terra 1.6.3 r <- rast(ncols = 10, nrows = 10, nlyr = 365) time(r) <- seq(from = as.POSIXlt("2001-01-01", tz = "UTC"), to = as.POSIXlt("2001-12-31", tz = "UTC"), by = "day") r #> class : SpatRaster #> dimensions : 10, 10, 365 (nrow, ncol, nlyr) #> resolution : 36, 18 (x, y) #> extent : -180, 180, -90, 90 (xmin, xmax, ymin, ymax) #> coord. ref. : lon/lat WGS 84 #> time : 2001-01-01 to 2001-12-31 UTC time(r) |> class() #> [1] "POSIXct" "POSIXt" r["2001"] #> Error: [subset] no (valid) layer selected time(r) <- as.Date("2002-01-01") + 0:364 r #> class : SpatRaster #> dimensions : 10, 10, 365 (nrow, ncol, nlyr) #> resolution : 36, 18 (x, y) #> extent : -180, 180, -90, 90 (xmin, xmax, ymin, ymax) #> coord. ref. : lon/lat WGS 84 #> time (days) : 2002-01-01 to 2002-12-31 time(r) |> class() #> [1] "Date" r["2002"] #> Error: [subset] no (valid) layer selected
内容的提问来源于stack exchange,提问作者dimfalk
相关产品推荐
相关产品推荐

