0

我一直在尝试使用greenbrown包的“PhenologyRaster”功能来模拟我研究区域的生长季节。但是，每次我运行该函数时，我都会得到空输出（例如，SOS.2016 层将显示为 NA）。我的问题如下：我遇到问题是因为我在一年的数据上运行该函数，还是因为 Landsat 时间序列有些不规则（即每年约 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, ...
4

1 回答 1

0

我也有同样的问题。我所做的是，我创建了一个 dummy series 三年，然后成功运行数据。

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)
于 2017-05-20T17:51:11.407 回答