时间序列可视化新思路：5分钟用Python的pyts库生成格拉姆角场(GAF/GADF)热图-编程阁

时间序列可视化新思路：5分钟用Python的pyts库生成格拉姆角场(GAF/GADF)热图

当面对一维时间序列数据时，折线图往往是我们首选的展示方式。但你是否想过，这些看似简单的数据点背后，可能隐藏着用传统方法难以察觉的模式？格拉姆角场（Gramian Angular Field）技术为我们打开了一扇新窗——它能将时间序列转化为二维热图，让数据的自相关性、周期性和突变特征以全新的视觉形式呈现。

对于数据分析师和科研人员来说，这种转换意义重大。想象一下，当你的股票价格数据或传感器读数被转化为彩色热图后，原本需要复杂统计检验才能发现的特征，现在可能只需一眼就能识别。更棒的是，借助Python的pyts库，整个过程只需5分钟就能完成。

1. 格拉姆角场核心原理速览

格拉姆角场（GAF）包含两种主要形式：求和型（GASF）和差分型（GADF）。它们通过巧妙的数学变换，保留了时间序列的全部信息，同时将其重新编码为图像矩阵。这种转换的核心优势在于：

信息无损：转换过程完全可逆，不会丢失原始数据特征
视觉友好：热图形式更符合人类对图像模式的识别习惯
兼容CNN：为后续使用卷积神经网络处理时间序列铺平道路

转换过程主要包含三个关键步骤：

分段聚合：通过PAA(Piecewise Aggregation Approximation)降低时间序列维度
极坐标转换：将缩放后的数据映射到极坐标系
格拉姆矩阵构建：计算角度关系的余弦值形成最终图像

提示：虽然数学原理看似复杂，但pyts库已经封装了所有底层计算，用户只需关注如何解读结果图像。

2. 快速上手：5分钟实战演示

让我们通过一个真实案例快速体验GAF的威力。假设我们有一组来自工业设备的振动传感器数据，采样频率为100Hz，持续10秒（共1000个数据点）。

首先安装必要的库：

pip install pyts numpy matplotlib

然后运行以下代码生成GAF图像：

from pyts.image import GramianAngularField import numpy as np import matplotlib.pyplot as plt # 生成模拟振动信号（含周期性冲击） t = np.linspace(0, 10, 1000) signal = np.sin(2 * np.pi * 5 * t) # 5Hz基础振动 signal[300:305] += 2 # 添加瞬时冲击 signal[700:710] += 1.5 # 添加持续扰动 # 初始化GAF转换器 gaf = GramianAngularField(image_size=100, method='summation') # 转换并绘制图像 X_gaf = gaf.fit_transform(signal.reshape(1, -1)) plt.figure(figsize=(10, 8)) plt.imshow(X_gaf[0], cmap='rainbow', origin='lower') plt.colorbar() plt.title("GASF Visualization of Vibration Signal") plt.show()

这段代码会输出一个100×100像素的彩色热图，其中：

对角线：代表时间序列的原始值（颜色越暖表示值越大）
非对角线：反映不同时间点之间的相互关系
局部模式：周期性振动会形成规则的条纹图案
异常点：冲击和扰动会表现为明显的颜色突变区域

3. 图像解读：从热图中挖掘洞见

学会生成图像只是第一步，更重要的是理解如何从中提取有价值的信息。以下是几种典型模式的解读指南：

热图特征	对应时间序列特征	实际意义
对角线的渐变	趋势成分	设备老化、温度漂移等缓慢变化
垂直于对角线的条纹	周期性波动	旋转机械的固有振动频率
局部色块突变	瞬时异常	设备碰撞、传感器故障等突发事件
对称/反对称区域	自相关特性	系统记忆效应或反馈机制强度

以之前的振动信号为例，我们可以清晰看到：

图像中均匀间隔的暖色带对应5Hz的基础振动
300ms处的红色斑点对应瞬时冲击
700ms处的黄色区域反映持续扰动

这种可视化方式比原始波形图更能凸显异常事件的时空关联性。

4. 高级技巧：参数调优与组合应用

为了获得最佳可视化效果，需要根据数据特性调整几个关键参数：

image_size：控制输出图像分辨率。较大的值保留更多细节但增加计算量，通常设置为原始数据点数的1/5到1/10。

method：选择'summation'(GASF)或'difference'(GADF)。两者的区别在于：

GASF更适合突出幅度变化和整体趋势
GADF对相位变化和局部波动更敏感

# 比较GASF和GADF的不同表现 gaf_gadf = GramianAngularField(image_size=100, method='difference') X_gadf = gaf_gadf.fit_transform(signal.reshape(1, -1)) fig, (ax1, ax2) = plt.subplots(1, 2, figsize=(16, 6)) ax1.imshow(X_gaf[0], cmap='rainbow', origin='lower') ax1.set_title('GASF') ax2.imshow(X_gadf[0], cmap='rainbow', origin='lower') ax2.set_title('GADF') plt.show()

对于金融时间序列分析，可以结合移动窗口技术生成动态GAF图像序列，观察市场模式演变：

# 滚动窗口GAF分析 window_size = 100 # 每个窗口100个数据点 n_windows = len(signal) // window_size plt.figure(figsize=(12, 4 * n_windows)) for i in range(n_windows): segment = signal[i*window_size : (i+1)*window_size] seg_gaf = gaf.fit_transform(segment.reshape(1, -1)) plt.subplot(n_windows, 1, i+1) plt.imshow(seg_gaf[0], cmap='rainbow', origin='lower', extent=[0, window_size, 0, window_size]) plt.title(f'Window {i+1}') plt.tight_layout() plt.show()