2016-02-01 90 views
0

我有三個柵格:如何從柵格繪圖?

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

哪個包定義了'光柵()'函數?你能否讓這個例子可重現? (也就是說,試試'plot(st1 [6,5 ,, drop = FALSE])或類似的。) – krlmlr

+0

謝謝Krlmlr,我的意思是時間序列不是圖像 –

+0

Coerce'st1 [...]'使用矢量'as.vector()'。不知道你是如何從這個時間序列中獲得的。 – krlmlr

回答

1

目前尚不清楚你想實現什麼。但是,你可以畫出這樣的:

time <- as.Date(Sys.time())+1:3 
plot(time, st1[6,5]) 

或者

s <- ts(as.vector(st1[6,5])) 
plot(s)