基于Python GDAL为长时间序列遥感图像绘制时相变化曲线图

本文介绍基于Python中gdal模块，对大量多时相栅格图像，批量绘制像元时间序列折线图的方法。

  首先，明确一下本文需要实现的需求：现有三个文件夹，其中第一个文件夹存放了某一研究区域原始的多时相栅格遥感影像数据（每一景遥感影像对应一个时相，文件夹中有多景遥感影像），每一景遥感影像都是.tif格式；第二个文件夹与第三个文件夹则分别存放了前述第一个文件夹中原始遥感影像基于2种不同滤波方法处理后的遥感影像（同样是每一景遥感影像对应一个时相，文件夹中有多景遥感影像），每一景遥感影像同样也都是.tif格式。我们希望分别针对这三个文件夹中的多张遥感影像数据，随机绘制部分像元对应的时间序列曲线图（每一个像元对应一张曲线图，一张曲线图中有三条曲线）；每一张曲线图的最终结果都是如下所示的类似的样式，X轴表示时间节点，Y轴就是具体的像素值。

  知道了需求，我们便开始代码的书写。具体代码如下：

# -*- coding: ut