本文介绍如何利用 Pandas 的 resample() 方法，将原始时间分辨率为 0.25 秒的频谱 CSV 数据高效下采样至 1 秒分辨率——即每 4 列时间序列合并为 1 列，按行（频率通道）保持结构不变，通过时间索引驱动聚合，避免手动切片与循环。
本文介绍如何利用 pandas 的 `resample()` 方法，将原始时间分辨率为 0.25 秒的频谱 csv 数据高效下采样至 1 秒分辨率——即每 4 列时间序列合并为 1 列，按行（频率通道）保持结构不变，通过时间索引驱动聚合，避免手动切片与循环。
在射频信号处理、光谱分析等场景中，原始采集数据常具有高时间分辨率（如每 0.25 秒一帧），但后续建模或可视化往往无需如此精细的粒度。此时需对时间轴进行有语义的下采样（downsampling）：不是简单丢弃列，而是将相邻时间点按物理意义聚合（如取均值以保留能量特征）。Pandas 提供了专为此类时序重采样设计的 .resample() 方法，其核心优势在于：自动对齐时间边界、支持多种聚合函数、可扩展性强，且完全向量化，性能远超 groupby 或循环切片。
✅ 正确做法：转置 + 时间重采样 + 转回
由于 .resample() 默认作用于 DataFrame 的行方向（即时间索引在 index 上），而你的原始数据中时间信息位于列名（columns），因此需先将 DataFrame 转置（.T），使时间成为索引，再执行重采样，最后转回原始结构：
# 假设已读入数据，且列名为时间戳字符串（如 '0.0', '0.25', '0.5', ...）或可解析为 datetime
# 第一步：将列名转换为 datetime 索引（起始时间 + 250ms 频率）
import pandas as pd
import numpy as np
# 示例：构造带时间索引的原始数据（3600 列 × 200 行）
count = 3600
df = pd.DataFrame({
f"freq_{i}MHz": np.random.randn(count) * 100
for i in range(1, 201)
}, index=pd.date_range("2025-04-01", periods=count, freq="250ms")).T
# 第二步：转置 → 重采样（1s 窗口均值）→ 再转置还原结构
df_downsampled = df.T.resample("1S").mean().T
print(f"原始形状: {df.shape}") # (200, 3600)
print(f"下采样后形状: {df_downsampled.shape}") # (200, 900)? 关键说明：
"1S"表示 1 秒（注意大小写，"1s"亦可，但"S"是 Pandas 官方推荐的秒单位缩写）；.resample(...).mean()会严格按日历秒对齐（如00:00:00,00:00:01, …），自动将每连续 4 个 250ms 时间点（如00:00:00.000,00:00:00.250,00:00:00.500,00:00:00.750）归入同一窗口并计算均值。
⚠️ 注意事项与常见陷阱
-
列名必须是合法时间索引：若原始 CSV 的列名为纯数字（如
0,1,2, …）或字符串（如't0','t1'），需先映射为pd.date_range：# 若原列为 [0, 0.25, 0.5, ..., 899.75]（单位：秒） original_seconds = np.arange(0, 3600 * 0.25, 0.25) # 长度 3600 df.columns = pd.to_datetime("2025-04-01") + pd.to_timedelta(original_seconds, unit='s') 避免使用
axis=1参数：旧版 Pandas 支持df.resample("1S", axis=1)，但自 v2.0.0 起已被弃用，官方明确要求使用.T.resample().T模式。-
缺失值处理：若某 1 秒窗口内存在 NaN，
.mean()默认跳过 NaN 计算（等效skipna=True）；若需严格要求 4 个有效值才聚合，可改用：df.T.resample("1S").apply(lambda x: x.mean() if len(x.dropna()) == 4 else np.nan).T 其他聚合方式：除
.mean()外，还可灵活替换为.max(),.std(),.first(), 或自定义函数（如加权平均、中位数等）。
✅ 验证结果正确性（小规模测试）
建议抽取前几行/列手工验证逻辑：
# 取前2频率 × 前8时间点（即2秒内8个0.25s样本）
df_mini = df.iloc[:2, :8].T # shape: (8, 2)
print("原始8点（转置后）：")
print(df_mini)
resampled_mini = df_mini.resample("1S").mean()
print("\n重采样后（每秒1行）：")
print(resampled_mini)
# 输出将清晰显示：第1行 = 前4点均值，第2行 = 后4点均值综上，df.T.resample("1S").mean().T 是处理此类“列即时间”型频谱数据下采样的最简洁、最鲁棒、最符合 Pandas 设计哲学的标准解法。它不仅代码简短，更具备时间对齐精度高、可读性强、易于扩展（如改为 2 秒或加窗滤波）等优势，应作为首选方案纳入你的信号预处理流水线。
