2017-03-01 13 views
0

私はRasterBrickから時系列を抽出することにより、動物園のオブジェクトを作成しました:時系列プロット:プロットからギャップを削除[動物園]

library(zoo) 
s <- RasterBrick 

target_cell <- 23265 

ss <- zooExtract(s, target_cell) 

windows() 
plot.zoo(ss, type = "o") 

私の結果は以下の通りです。私はすべてのポイントを通過する行をしたいが、私はそれを動作させるように見えることはできません、私はこの質問のための以前の答えを見てきました。私はlines()を試して、それは同じ結果を与えます。

+0

はあなたのデータの値NA持つように見えますが、再現性の例 –

+0

せずに多くを言うことができないはいI多くのデータギャップとその不規則な時系列がある、これは問題ですか? – dtanon

+1

これらは 'ss_filled = zoo :: na.locf(ss);で動作するはずです。 ss_complete_cases = ss [complete.cases(ss)]; plot.zoo(ss_complete_cases、type = "o") 'また、行方不明の日付を保存することをお勧めしますデータ – OdeToMyFiddle

答えて

1

@Richardテルフォードは、プロットのギャップを示唆したように欠落したデータによるものです。 zooから基本機能complete.casesna.locfで ギャップ が削除または下記のように先の観察を充填することができます。

#to identity periods with missing data 

missingPeriod = as.Date(index(ss[!complete.cases(ss),])) 


#to retain only periods with no missing data 
ssComplete = ss[complete.cases(ss),] 


#to retain all periods with gaps filled with previous value 
#maxgap parameter controls number of missing data replaced with prev observation 
N = 5 
ssFilled = zoo::na.locf(ss,maxgap=N); 

#plots 
plot.zoo(ssComplete , type = "o") 

plot.zoo(ssFilled, type = "o") 
関連する問題