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