2016-08-20 24 views
0

ラスタブリックのインタラクティブなプロットを作成しようとしています。ピクセルをクリックすると、そのピクセルで時系列のデータが得られます。 (私のラスタ煉瓦はおよそ345の画像である。)プロットピクセル時系列

これは私がやったことです:

library(raster) 

EVI <- "D:\\Modis_EVI\\Original\\EVI_Stack_single5000.tif" 
y.EVI <- brick(EVI) 
plot(as.numeric(click(y.EVI)), type="l", lwd=2) 

しかし、それがすべてでプロットしません。そして、私は4つの画像のような小さなスタックで試してみると、このエラーが発生します:

Error in plot.window(…) : need finite 'xlim' values 

助けてください。

+0

なぜ単一の値をプロットするだろうが、座標系とx値の設定に意味がありますか?最初にtifを表示する必要があるようですか? –

答えて

0

ユーザーが一度クリックできるようにする場合は、n=1clickに指定する必要があります。たとえば:

library(raster) 
b <- brick(replicate(10, raster(matrix(runif(100), ncol=10)))) 

plot_ts <- function(x) { 
    plot(x[[1]]) 
    z <- c(click(x, n=1, show=FALSE)) 
    plot(z, type='l', lwd=2, ylab='y', xlab='time', las=1) 
    z 
} 

z <- plot_ts(b) 

ここで、セルをクリックした後にプロット例だ... enter image description here