作者：CSDN @ _养乐多_

本文记录了使用gdal、ras、numpy 库实现遥感时间序列数据逐像素处理的代码。并以求NADVI时间序列最大值为例。代码可扩展到其他多种对时间序列的处理，比如MK趋势分析，求时间序列中值、众数、标准差、和，时间序列拟合、异常检测、机器学习预测等多种应用中。

代码可以输出时间序列最大值对应的时间，可以是年份（1996，1997，1998, …），也可以是月份(1, 3, 5, 6, …)，也可以是时间(1996.2, 1996.3, 1997.4, 1998.3, …)，也可以是DOY（12， 24， 45，200，…），可以根据需求输出结果。

一、代码框架

请参考《python：处理遥感时间序列（代码框架），并保存结果》。

二、完整代码

第一版本，没有整理的shi山代码，但是结果没问题。