前面分享了一篇论文打算复现研究一下，查了查用到的方法在这里分享Jan Verbeselt写的教程。教程原文请大家点击文末 阅读原文 跳转

经过测试代码能够运行通过。由于需要下載MODIS数据运行时间比较长，我是代码跑起来之后出去吃了个晚饭才运行完的

基于R的MODIS数據下载和分析

安装和加载MODIS数据分析所需的R包

在加载包之前，需要先安装R包使用下面的代码可以把所需的R包都安好。

 

 定义了一个timeser函数用于創建时间序列
 
 

 

 
 

 使用像元可靠性信息清洗MODIS NDVI数据
 
 

 

 下面将数据从栅格中取出，作为一个向量使用timeser函数生成一个时间序列。
 
 

 

 接下来就是栅格断點计算：
 
 

 

 使用click函数选择像元并对像元进行BFAST分析：