如何用Python处理时间序列数据?resample重采样

使用pandas的resample方法进行时间序列数据处理及聚合的核心步骤如下:1. 确保dataframe或series具有datetimeindex,这是resample操作的前提;2. 使用resample(‘freq’)指定目标频率,如’d’(日)、’w’(周)、’m’(月)等;3. 应用聚合函数如.mean()、.sum()、.ohlc()等对每个时间区间内的数据进行汇总;4. 可通过label和closed参数控制时间区间的标签位置和闭合端点;5. 对缺失值使用fillna()方法进行填充或保留nan;6. 结合apply()、agg()、pipe()实现自定义聚合、多列多函数聚合和后处理流程;7. 区分resample()与asfreq(),前者用于聚合降频,后者用于频率转换而不聚合。

如何用Python处理时间序列数据?resample重采样

用Python处理时间序列数据,特别是涉及频率转换和数据聚合时,Pandas库里的resample方法简直是神器。它能让你轻松地将数据从一个时间粒度转换到另一个,比如把每分钟的数据汇总成每日、每周甚至每月的数据,同时还能对这些新的时间段内的数据进行各种统计计算。

如何用Python处理时间序列数据?resample重采样

解决方案

要使用resample,你首先得确保你的DataFrame或Series有一个DatetimeIndex。这是基础,也是关键。一旦有了这个索引,resample方法就能像变魔术一样工作。

想象一下,你有一组高频的传感器数据,每秒记录一次,但你现在只想看每天的平均值。

立即学习“Python免费学习笔记(深入)”;

如何用Python处理时间序列数据?resample重采样

import pandas as pdimport numpy as np# 创建一个示例时间序列数据# 假设从2023年1月1日0点开始,每分钟一个数据点start_time = pd.to_datetime('2023-01-01 00:00:00')time_index = pd.date_range(start=start_time, periods=1000, freq='T') # 'T' for minute frequencydata = np.random.rand(1000) * 100df = pd.DataFrame({'value': data}, index=time_index)print("原始数据(前5行):")print(df.head())print("n原始数据信息:")df.info()# 将数据重采样为每日平均值daily_avg = df.resample('D').mean()print("n重采样为每日平均值(前5行):")print(daily_avg.head())# 重采样为每周总和weekly_sum = df.resample('W').sum()print("n重采样为每周总和(前5行):")print(weekly_sum.head())# 重采样为每月OHLC(开高低收)monthly_ohlc = df['value'].resample('M').ohlc()print("n重采样为每月OHLC(前5行):")print(monthly_ohlc.head())# 如果想指定标签的位置(比如用区间的结束时间作为标签)daily_avg_end_label = df.resample('D', label='right').mean()print("n重采样为每日平均值(标签为区间结束时间,前5行):")print(daily_avg_end_label.head())# 还可以用loffset来调整时间戳daily_avg_shifted = df.resample('D', loffset='-12H').mean() # 将标签向后移12小时print("n重采样并调整标签(前5行):")print(daily_avg_shifted.head())

resample的字符串参数非常灵活,支持各种频率,比如'D'(日)、'W'(周)、'M'(月)、'Q'(季度)、'A'(年)、'H'(小时)、'T''min'(分钟)、'S'(秒)等等。后面跟着的聚合函数(如.mean(), .sum(), .first(), .last(), .ohlc(), .count())决定了每个新时间段内的数据如何被汇总。

resample背后的数据聚合逻辑是什么?

在我看来,resample的核心逻辑,其实就是一种智能的“分箱”(binning)操作。它根据你指定的频率,在时间轴上划分出一个个连续的、互不重叠的时间区间。然后,它会把落在同一个时间区间内的所有原始数据点“收集”起来,形成一个临时的组。最后,你调用 .mean(), .sum() 或者其他聚合函数时,就是对这些组内的数据进行计算,得到每个新时间区间的最终结果。

如何用Python处理时间序列数据?resample重采样

举个例子,如果你把数据从分钟级别重采样到日级别(resample('D')),那么它会识别出所有属于1月1日的数据,把它们归到一起;再识别出所有属于1月2日的数据,又归到一起,以此类推。每个“日”就形成了一个独立的箱子。聚合函数就是对这些箱子里的数据进行操作。

这个过程很巧妙,因为它不仅仅是简单地选择数据,而是真正地重新组织数据。如果某个时间段内没有数据,resample默认会为这个时间段创建一个条目,并填充NaN,这在处理不连续或稀疏的时间序列时非常有用,因为它能帮你清晰地看到数据的缺失情况。我觉得这种“显式缺失”比隐式缺失要好得多,至少你知道哪里没数据。

如何处理resample操作中的缺失值和边界情况?

缺失值和边界情况在实际数据处理中是常态,resample也提供了参数来应对。

缺失值处理:正如前面提到的,如果某个重采样后的时间区间内没有原始数据点,resample会填充NaN。处理这些NaN,你可以在resample之后链式调用Pandas的缺失值处理方法:

# 假设我们有一个稀疏的时间序列sparse_time_index = pd.to_datetime(['2023-01-01', '2023-01-03', '2023-01-06'])sparse_data = [10, 20, 30]sparse_df = pd.DataFrame({'value': sparse_data}, index=sparse_time_index)# 重采样为每日,会有NaNdaily_sparse = sparse_df.resample('D').mean()print("n稀疏数据重采样为每日(含NaN):")print(daily_sparse)# 填充缺失值daily_filled_ffill = daily_sparse.fillna(method='ffill') # 用前一个有效值填充print("n缺失值前向填充:")print(daily_filled_ffill)daily_filled_bfill = daily_sparse.fillna(method='bfill') # 用后一个有效值填充print("n缺失值后向填充:")print(daily_filled_bfill)daily_filled_zero = daily_sparse.fillna(0) # 填充0print("n缺失值填充0:")print(daily_filled_zero)

选择哪种填充方式取决于你的业务逻辑和数据特性。有时候,NaN本身就是一种信息,表示那个时间点确实没有数据,不需要填充。

边界情况(closedlabel):resample默认情况下,区间的左边界是包含的,右边界是不包含的(closed='left'),并且新生成的时间戳标签是区间的左边界(label='left')。但你可以在特定场景下调整它们:

closed: 决定区间的哪一端是闭合的(包含)。'left' (默认): [start, end)'right': (start, end]label: 决定新时间段的标签是使用区间的开始时间还是结束时间。'left' (默认): 使用区间的开始时间作为标签。'right': 使用区间的结束时间作为标签。

# 原始数据,方便观察边界data_for_boundary = pd.DataFrame({'value': range(5)}, index=pd.to_datetime(['2023-01-01 00:00', '2023-01-01 23:59', '2023-01-02 00:00', '2023-01-02 23:59', '2023-01-03 00:00']))print("n原始数据用于边界观察:")print(data_for_boundary)# 默认行为:closed='left', label='left'# 2023-01-01的数据包括00:00到23:59# 2023-01-02的数据包括00:00到23:59default_resample = data_for_boundary.resample('D').sum()print("n默认重采样(closed='left', label='left'):")print(default_resample)# closed='right', label='right'# 2023-01-01的数据包括(2022-12-31 23:59, 2023-01-01 23:59]# 标签是结束时间right_closed_right_label = data_for_boundary.resample('D', closed='right', label='right').sum()print("n重采样(closed='right', label='right'):")print(right_closed_right_label)

理解这些参数对于确保数据聚合的准确性至关重要,特别是当你的业务对时间边界有严格要求时。

除了简单的聚合,resample还能做哪些高级操作?

resample的强大之处远不止于简单的聚合函数。它能与Pandas的其他功能结合,实现更复杂的分析。

自定义聚合函数 (apply):如果内置的mean, sum等不能满足你的需求,你可以传入一个自定义函数到apply方法中。这个函数会接收每个重采样区间的Series数据。

# 计算每个月数据的范围 (max - min)monthly_range = df['value'].resample('M').apply(lambda x: x.max() - x.min())print("n每月数据范围(自定义聚合):")print(monthly_range.head())# 甚至可以返回多个值,但需要确保返回的是Series或DataFrame,或者每个区间返回一个标量# 比如返回每个月的Q1和Q3def get_quantiles(x):    return pd.Series({'Q1': x.quantile(0.25), 'Q3': x.quantile(0.75)})monthly_quantiles = df['value'].resample('M').apply(get_quantiles)print("n每月Q1和Q3(自定义聚合返回多值):")print(monthly_quantiles.head())

这给了你极大的灵活性,可以根据业务需求实现任何复杂的聚合逻辑。

多列聚合:如果你DataFrame有多列,并且想对不同的列应用不同的聚合函数,可以使用agg方法。

df_multi_col = df.copy()df_multi_col['another_value'] = np.random.rand(len(df)) * 50# 对'value'列求平均,对'another_value'列求和daily_agg_multi = df_multi_col.resample('D').agg({    'value': 'mean',    'another_value': 'sum'})print("n多列聚合(不同函数):")print(daily_agg_multi.head())# 对所有数值列应用多个聚合函数daily_agg_all_cols = df_multi_col.resample('D').agg(['mean', 'sum', 'std'])print("n所有数值列应用多个聚合函数:")print(daily_agg_all_cols.head())

agg方法让多维度分析变得异常便捷。

pipe结合:pipe方法允许你将resample的结果直接传递给一个函数或一系列函数,实现更流畅的工作流。

def process_resampled_data(resampled_df):    # 假设这里有一些后处理逻辑,比如填充缺失值并计算一个新列    processed_df = resampled_df.fillna(0)    processed_df['value_ratio'] = processed_df['value'] / processed_df['another_value']    return processed_df# 将resample的结果通过pipe传递给自定义处理函数processed_output = df_multi_col.resample('D').mean().pipe(process_resampled_data)print("n使用pipe进行后处理:")print(processed_output.head())

这有助于构建清晰、可读的数据处理管道。

asfreq()resample()区别虽然不是resample本身的高级操作,但理解asfreq()resample()的区别对于时间序列处理至关重要。resample()是聚合操作,它会改变数据点的数量,并对每个区间进行聚合。而asfreq()仅仅是改变频率,它不会进行聚合,如果新的频率点上没有数据,它会填充NaN,或者直接取最近的一个点(取决于参数)。

# df是原始分钟级数据# resample到每小时,取平均hourly_resample = df.resample('H').mean()print("nresample到每小时(平均):")print(hourly_resample.head())# asfreq到每小时,不聚合,直接取对应时间点的值,没有则NaN# 注意:asfreq通常用于升采样,或在降采样时仅取特定点hourly_asfreq = df.asfreq('H')print("nasfreq到每小时(不聚合):")print(hourly_asfreq.head())

asfreq在需要保持原始数据点结构,只是改变观察频率时非常有用,比如从日数据变为月数据,但只关心每个月的第一天。

这些高级用法让resample不仅仅是一个简单的聚合工具,更是时间序列分析中不可或缺的利器。通过灵活运用,你可以解决各种复杂的时间序列数据挑战。

以上就是如何用Python处理时间序列数据?resample重采样的详细内容,更多请关注创想鸟其它相关文章!

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。
如发现本站有涉嫌抄袭侵权/违法违规的内容, 请发送邮件至 chuangxiangniao@163.com 举报,一经查实,本站将立刻删除。
发布者:程序猿,转转请注明出处:https://www.chuangxiangniao.com/p/1363704.html

(0)
打赏 微信扫一扫 微信扫一扫 支付宝扫一扫 支付宝扫一扫
Python怎样处理文本数据?字符串操作完整指南
上一篇 2025年12月14日 03:34:28
Python如何分析数据相关性—热力图与相关系数矩阵
下一篇 2025年12月14日 03:34:38

相关推荐

  • composer require-dev和require有什么不同_Composer Require与Require-Dev区别解析

    require用于声明项目运行必需的依赖,如框架、数据库组件和第三方SDK,这些包会随项目部署到生产环境;2. require-dev用于声明仅在开发和测试阶段需要的工具,如PHPUnit、PHPStan、Faker等,不会默认部署到生产环境;3. 安装时composer install根据环境决定…

    2026年5月10日
    1000
  • Matplotlib 地图中多类型图例的创建与优化

    Matplotlib 地图中多类型图例的创建与优化Matplotlib 地图中多类型图例的创建与优化Matplotlib 地图中多类型图例的创建与优化Matplotlib 地图中多类型图例的创建与优化

    本教程旨在解决matplotlib地图可视化中,如何在一个图例中同时展示颜色块(如区域分类)和自定义标记(如特定兴趣点)的问题。文章详细介绍了当传统`patch`对象无法正确显示标记时,如何利用`matplotlib.lines.line2d`创建标记图例句柄,并将其与颜色块图例句柄合并,从而生成一…

    2026年5月10日 用户投稿
    100
  • Golang JSON序列化:控制敏感字段暴露的最佳实践

    本教程探讨golang中如何高效控制结构体字段在json序列化时的可见性。当需要将包含敏感信息的结构体数组转换为json响应时,通过利用`encoding/json`包提供的结构体标签,特别是`json:”-“`,可以轻松实现对特定字段的忽略,从而避免敏感数据泄露,确保api…

    2026年5月10日
    000
  • 利用海象运算符简化条件赋值:Python教程与最佳实践

    本文旨在探讨Python中海象运算符(:=)在条件赋值场景下的应用。通过对比传统if/else语句与海象运算符,以及条件表达式,分析海象运算符在简化代码、提高可读性方面的优势与局限性。并通过具体示例,展示如何在列表推导式等场景下合理使用海象运算符,同时强调其潜在的复杂性及替代方案,帮助开发者更好地掌…

    2026年5月10日
    100
  • Debian syslog性能优化技巧有哪些

    提升Debian系统syslog (通常基于rsyslog)性能,关键在于精简配置和高效处理日志。以下策略能有效优化日志管理,提升系统整体性能: 精简配置,高效加载: 在rsyslog配置文件中,仅加载必要的输入、输出和解析模块。 使用全局指令设置日志级别和格式,避免不必要的处理。 自定义模板: 创…

    2026年5月10日
    000
  • 比特币新手教程 比特币交易平台有哪些

    比特币是一种去中心化的数字货币,基于区块链技术实现点对点交易,具有匿名性、有限发行和不可篡改等特点;新手可通过交易所购买,P2P交易获得比特币,常用平台包括Binance、OKX和Huobi;交易流程包括注册账户、实名认证、绑定支付方式、充值法币并下单购买,可选择市价单或限价单;比特币存储方式有交易…

    2026年5月10日
    000
  • c++中的SFINAE技术是什么_c++模板编程中的SFINAE原理与应用

    SFINAE 是“替换失败不是错误”的原则,指模板实例化时若参数替换导致错误,只要存在其他合法候选,编译器不报错而是继续重载决议。它用于条件启用模板、类型检测等场景,如通过 decltype 或 enable_if 控制函数重载,实现类型特征判断。尽管 C++20 引入 Concepts 简化了部分…

    2026年5月10日
    000
  • Go语言mgo查询构建:深入理解bson.M与日期范围查询的正确实践

    本文旨在解决go语言mgo库中构建复杂查询时,特别是涉及嵌套`bson.m`和日期范围筛选的常见错误。我们将深入剖析`bson.m`的类型特性,解释为何直接索引`interface{}`会导致“invalid operation”错误,并提供一种推荐的、结构清晰的代码重构方案,以确保查询条件能够正确…

    2026年5月10日
    100
  • RichHandler与Rich Progress集成:解决显示冲突的教程

    在使用rich库的`richhandler`进行日志输出并同时使用`progress`组件时,可能会遇到显示错乱或溢出问题。这通常是由于为`richhandler`和`progress`分别创建了独立的`console`实例导致的。解决方案是确保日志处理器和进度条组件共享同一个`console`实例…

    2026年5月10日
    000
  • Golang goroutine与channel调试技巧

    使用go run -race检测数据竞争,结合runtime.NumGoroutine监控协程数量,通过pprof分析阻塞调用栈,利用select超时避免永久阻塞,有效排查goroutine泄漏、死锁和数据竞争问题。 Go语言的goroutine和channel是并发编程的核心,但它们也带来了调试上…

    2026年5月10日
    000
  • 使用 Jupyter Notebook 进行探索性数据分析

    Jupyter Notebook通过单元格实现代码与Markdown结合,支持数据导入(pandas)、清洗(fillna)、探索(matplotlib/seaborn可视化)、统计分析(describe/corr)和特征工程,便于记录与分享分析过程。 Jupyter Notebook 是进行探索性…

    2026年5月10日
    000
  • 《魔兽世界》将于6月11日开启国服回归技术测试

    《魔兽世界》将于6月11日开启国服回归技术测试《魔兽世界》将于6月11日开启国服回归技术测试《魔兽世界》将于6月11日开启国服回归技术测试《魔兽世界》将于6月11日开启国服回归技术测试

    《%ign%ignore_a_1%re_a_1%》官方宣布,将于6月11日开启国服回归技术测试,时间为7天,并称可以在6月内正式开服,玩家们可以访问官网下载战网客户端并预下载“巫妖王之怒”客户端,技术测试详情见下图。 WordAi WordAI是一个AI驱动的内容重写平台 53 查看详情 以上就是《…

    2026年5月10日 用户投稿
    200
  • php常量怎么用_PHP常量(define/const)定义与使用方法

    PHP中可通过define函数和const关键字定义常量,用于存储不可变值。define适用于全局作用域,支持动态名称和条件定义,如define(‘SITE_NAME’, ‘MyWebsite’);const在编译时生效,语法简洁但限制多,只能在类或全…

    2026年5月10日
    000
  • 如何在HTML中插入表单元素_HTML表单控件与输入类型使用指南

    HTML表单通过标签构建,包含action和method属性定义数据提交目标与方式,常用input类型如text、password、email等适配不同输入需求,配合label、required、placeholder提升可用性,结合textarea、select、button等控件实现完整交互,是…

    2026年5月10日
    000
  • 网站标题关键词更新后,搜索引擎为何仍显示旧标题?

    网站标题更新后,搜索引擎为何显示旧标题? 网站SEO优化中,站长常修改网站标题关键词,期望搜索结果显示自定义标题。然而,即使更新标签、meta keywords、meta description和结构化数据中的name属性后,搜索结果仍显示旧标题,这令人费解。本文将对此进行解释。 问题:站长修改了网…

    2026年5月10日
    100
  • 创建指定大小并填充特定数据的Golang文件教程

    本文将介绍如何使用Golang创建一个指定大小的文件,并用特定数据填充它。我们将使用 `os` 包提供的函数来创建和截断文件,从而实现快速生成大文件的目的。示例代码展示了如何创建一个10MB的文件,并将其填充为全零数据。掌握这些方法,可以方便地在例如日志系统或磁盘队列等场景中,预先创建测试文件或初始…

    2026年5月10日
    000
  • Python命令怎样使用profile分析脚本性能 Python命令性能分析的基础教程

    使用Python的cProfile模块分析脚本性能最直接的方式是通过命令行执行python -m cProfile your_script.py,它会输出每个函数的调用次数、总耗时、累积耗时等关键指标,帮助定位性能瓶颈;为进一步分析,可将结果保存为文件python -m cProfile -o ou…

    2026年5月10日
    000
  • 使用 WebCodecs VideoDecoder 实现精确逐帧回退

    本文档旨在解决在使用 WebCodecs VideoDecoder 进行视频解码时,实现精确逐帧回退的问题。通过比较帧的时间戳与目标帧的时间戳,可以避免渲染中间帧,从而提高用户体验。本文将提供详细的解决方案和示例代码,帮助开发者实现精确的视频帧控制。 在使用 WebCodecs VideoDecod…

    2026年5月10日
    000
  • 如何插入查询结果数据_SQL插入Select查询结果方法

    如何插入查询结果数据_SQL插入Select查询结果方法如何插入查询结果数据_SQL插入Select查询结果方法如何插入查询结果数据_SQL插入Select查询结果方法如何插入查询结果数据_SQL插入Select查询结果方法

    使用INSERT INTO…SELECT语句可高效插入数据,通过NOT EXISTS、LEFT JOIN、MERGE语句或唯一约束避免重复;表结构不一致时可通过别名、类型转换、默认值或计算字段处理;结合存储过程可提升可维护性,支持参数化与动态SQL。 将查询结果数据插入到另一个表中,可以…

    2026年5月10日 用户投稿
    000
  • Discord.py 交互按钮超时与持久化解决方案

    本教程旨在解决Discord.py中交互按钮在一段时间后出现“This Interaction Failed”错误的问题。我们将深入探讨视图(View)的超时机制,并提供通过正确设置timeout参数以及利用bot.add_view()方法实现按钮持久化的具体方案,确保您的机器人交互功能稳定可靠,即…

    2026年5月10日
    000

发表回复

登录后才能评论
关注微信