有没有办法在R中对时间序列NDVI图像运行savitzky golay过滤器。我已经尝试过使用包中给出的以下代码＆＃39; signal＆＃39;
sg <- sgolayfilt(timeseries,3,5)。
但它返回以下错误;
Error in if (all(is.na(x))) return(x) :
argument is not interpretable as logical
文件＆＃34; timeseries＆＃34;这是堆叠的栅格NDVI图像。在这方面有人可以帮助我。 谢谢你的帮助。
答案 0 :(得分：1)
我有12个栅格图层，然后将它们堆叠成栅格堆栈
Date您可能需要根据数据调整p和n。顺便说一句，n必须是奇数维度，根据我的经验，p小于n
我希望它能提供帮助：）