我有三个栅格:
library (raster)
r <- raster(nrows=10, ncols=10); r <- setValues(r, 1:ncell(r))
r1 <- raster(nrows=10, ncols=10);r1 <- setValues(r1, 1:ncell(r))
r2 <- raster(nrows=10, ncols=10);r2 <- setValues(r1, 1:ncell(r))
st1=stack(r,r1,r2)
所有栅格的getvalues第6行和第5列
> st1[6,5,]
layer.1 layer.2 layer.3
[1,] 55 55 55
我试图将其绘制为时间序列,但情节不正确:
plot(st1[6,5,])
答案 0 :(得分:1)
目前尚不清楚您想要实现的目标。但你可以这样画:
time <- as.Date(Sys.time())+1:3
plot(time, st1[6,5])
或者
s <- ts(as.vector(st1[6,5]))
plot(s)