r - 如何在R中的MODIS NDVI数据集（MYD13Q1）上执行时间序列分析？

我有来自MODIS（MYD13Q1）的15年数据，该数据每16天提供NDVI值。所有这些数据均为tiff格式。我想通过时间序列分析，通过查看季节和趋势成分来分析这15年中特定区域植被的变化。

我曾尝试使用bfast和bfast-spatial之类的软件包，但并没有取得太大的成功，因为Internet上用于它们的大多数资料都与landsat数据集有关。我尝试遵循这两个教程。 bfast-spatial，bfast

上面提到的教程仅适用于特定像素，但是我想按照我的目标中提到的那样处理大型数据集。如何通过修改这些教程来实现这一目标