本文记录了使用python语言和Pettitt突变检测算法检测遥感数据30年的NDVI时间序列中的突变时间(年份)的代码。
Pettitt可以检测时间序列中最显著的突变点,其检测结果只有一个,相对MK突变检测更准确。p为显著性阈值,大于最大阈值的被视为最显著突变点。
下图是本文实现的Pettitt突变检测方法检测得到的突变点以及满足的置信水平,第一栏为原始NDVI时间序列数据及突变点,从图中可以看出明显的突变位置,第2栏为统计量Ut及显著性水平p,可以看出该突变满足显著性水平p=0.001。
Pettitt突变检测示意图:
一、读取遥感时间序列数据
读取时间序列数据的代码请查看本专栏《python:处理遥感时间序列(代码框架),并保存结果》一文。