作者：CSDN @ _养乐多_

本文将介绍使用python编程语言，进行遥感数据时间序列最大值所对应的影像索引的代码。代码中使用了numpy和gdal，通过numpy广播机制实现时间序列最大值所对应的影像索引的计算，并以NDVI时间序列数据为例。代码方便易运行，逻辑简单，速度快。只需要输入单波段遥感数据，就可输出最大值所对应的索引的影像。

输入输出如下图所示，

---

---

一、使用说明

• 索引是按文件夹中影像的命名顺序排列的。
• 如果你的索引表示年份的前后顺序，且你想把索引换成年份，可以修改代码也可以在ArcGIS中重分类。
• 如果时间序列数据是年内的数据，且时间顺序按照doy排序，索引对应不同的doy，可以修改代码也可以在ArcGIS中重分类。
• 修改像素值的代码网上比较多，本文不做记录。

二、完整代码