Landsatデータでgreenbrownパッケージを使用

Question

greenbrownパッケージの 'PhenologyRaster'機能を使用して、自分のスタディエリアの成長期をモデル化しようとしています。しかし、関数を実行するたびに空の出力が得られます（たとえば、SOS.2016のレイヤーはNAと表示されます）。私の質問は次のとおりです。1年のデータで機能を実行しているか、Landsatの時系列が多少不規則なため（つまり1年あたり〜30シーンの頻度）、問題が発生していますか？私はPhenologyRatser機能を実行するために、次のコードを使用しています

：

PhenoTest = PhenologyRaster(landsat2016,start=c(2016,1,3),end=c(2016,12,20),freq=24,approach="Deriv",min.mean=-0.5,tsgf='TSGFspline',interpolate=TRUE) 

関数は、次の特性を持つラスタスタックに適用されます。

class  : RasterBrick 
dimensions : 526, 591, 310866, 18 (nrow, ncol, ncell, nlayers) 
resolution : 30, 30 (x, y) 
extent  : 604965, 622695, 4208175, 4223955 (xmin, xmax, ymin, ymax) 
coord. ref. : +proj=utm +zone=10 +datum=WGS84 +units=m +no_defs +ellps=WGS84 +towgs84=0,0,0 
data source : in memory 
names  : X2016.01.03, X2016.01.19, X2016.02.04, X2016.03.07, X2016.03.23, X2016.04.24, X2016.05.10, X2016.05.26, X2016.06.27, X2016.07.13, X2016.07.29, X2016.08.14, X2016.08.30, X2016.09.15, X2016.10.01, ... 
min values :  -0.1964,   NA,  -0.5382,   NA,  -0.4696,  -0.2197,  -0.2803,  -0.4274,  -0.4827,  -0.2631,  -0.5256,  -0.4856,  -0.5631,  -0.3204,  -0.5512, ... 
max values :  0.1714,   NA,  0.2425,   NA,  0.2061,  0.5173,  0.4583,  0.2470,  0.3629,  0.5165,  0.2981,  0.2802,  1.6199,  0.5016,  0.3007, ... 

Answer 1

は私も同じでした問題。私がやったのは、3年間ダミーシリーズを作ってからデータが正常に実行されたということでした。

b.1 <- brick(r.1, r.2, r.3, r.4, r.5, r.6, r.7, r.8, r.9, r.10, r.11, r.12) 
b.2 <- stack(b.2, b.2, b.2) 
pheno.test <- PhenologyRaster(b.2, start=c(2016,1), freq=12, approach="White", 
           tsgf="TSGFspline", interpolate=T)