0

0

pandas 如何在 rolling().apply() 中使用 numba jit 加速

舞姬之光

舞姬之光

发布时间:2026-01-29 14:01:37

|

157人浏览过

|

来源于php中文网

原创

@njit函数不能直接用于rolling().apply(),因其不支持pandas.Series输入;正确做法是设raw=True使窗口传入ndarray,再用@njit处理一维数组并返回标量。

pandas 如何在 rolling().apply() 中使用 numba jit 加速

为什么 rolling().apply() 不能直接用 @njit

因为 pandasrolling().apply() 要求传入函数能接收一个 Seriesndarray(取决于 raw 参数),而 numba@njit 编译后函数默认不支持 pandas.Series —— 它只认纯 numpy.ndarray,且要求类型在编译时可推断。直接把 @njit 函数塞进 apply() 会报类似 TypingError: Failed in nopython mode pipeline... cannot determine Numba type of <class></class> 的错误。

正确做法:设 raw=True + 用 @njit 处理 ndarray

关键点是让 rolling 把窗口数据以 numpy.ndarray 形式传给你的函数,而不是 Series。这靠设置 raw=True 实现。此时传入函数的参数就是一维 ndarraynumba 可以顺利编译和运行。

  • raw=True 是必须的,否则传进来的是 Series@njit 直接拒收
  • 函数签名必须只接受一个 ndarray,返回标量(如 float64);不能有可选参数、不能调用 pandas 方法
  • 如果窗口含 NaNraw=True 下仍会传入 nan 值,需在 @njit 函数内显式处理(比如用 np.isnan() 判断)
  • 示例:计算滚动窗口中位数(np.mediannumba 中不可用,得手写)
<pre class="brush:php;toolbar:false;">@njit
def rolling_median(arr):
    n = len(arr)
    if n == 0:
        return np.nan
    # 简单排序取中位(仅作示意;生产环境建议用 partial sort)
    temp = arr.copy()
    for i in range(n):
        for j in range(i + 1, n):
            if temp[i] > temp[j]:
                temp[i], temp[j] = temp[j], temp[i]
    if n % 2 == 1:
        return temp[n // 2]
    else:
        return (temp[n // 2 - 1] + temp[n // 2]) / 2.0
<h1>使用方式</h1><p>s = pd.Series([1.0, 3.0, 2.0, 5.0, 4.0])
result = s.rolling(3).apply(rolling_median, raw=True)</p>

性能对比与注意事项

加速效果取决于函数复杂度。对简单操作(如 np.sum),用 <code>@njit 反而可能更慢(启动开销+小数组无优势);但对自定义逻辑(如条件聚合、分位数、迭代计算),提速常达 2–10 倍。不过要注意:

  • raw=True 会禁用 rolling 内置的 min_periods 对齐逻辑:若窗口内有效值不足,apply 仍会传入含 NaN 的数组,你得自己检查有效长度
  • @njit 不支持 axis 参数或高维数组;rolling 在 DataFrame 上使用时,需逐列处理或改用 apply(axis=1) + raw=True,但后者传入的是行向量(1D),不是二维切片
  • 调试困难:@njit 报错信息晦涩,建议先用纯 Python 版本验证逻辑,再加装饰器

替代方案:用 numba 手写滚动循环(更可控)

rolling().apply(raw=True) 无法满足需求(比如要跳过 NaN、动态窗口、或需要索引信息),更稳的方式是绕过 pandas 的 rolling,直接用 numba 写一个带循环的函数,输入整个数组,输出结果数组。这样完全掌控内存布局、缺失值策略和边界行为。

Spirit Me
Spirit Me

SpiritMe允许用户使用数字化身制作视频,这些化身可以模拟用户的声音和情感

下载

例如实现等价于 s.rolling(3).mean() 的 jit 版本:

<code>@njit
def rolling_mean_jit(arr, window):
    n = len(arr)
    out = np.full(n, np.nan)
    for i in range(window - 1, n):
        window_arr = arr[i - window + 1:i + 1]
        # 过滤 NaN(numba 支持 np.nansum / np.nanmean 从 0.58+,但旧版需手动)
        total, count = 0.0, 0
        for j in range(window):
            if not np.isnan(window_arr[j]):
                total += window_arr[j]
                count += 1
        if count > 0:
            out[i] = total / count
    return out
<h1>调用</h1><p>result = rolling_mean_jit(s.to_numpy(), 3)</p>

这种写法失去 pandas 的自动对齐和日期索引保留能力,但换来最大灵活性和可预测性能。真正高频、长序列、定制逻辑强的场景,这是更可靠的选择。

热门AI工具

更多
DeepSeek
DeepSeek

幻方量化公司旗下的开源大模型平台

豆包大模型
豆包大模型

字节跳动自主研发的一系列大型语言模型

WorkBuddy
WorkBuddy

腾讯云推出的AI原生桌面智能体工作台

腾讯元宝
腾讯元宝

腾讯混元平台推出的AI助手

文心一言
文心一言

文心一言是百度开发的AI聊天机器人,通过对话可以生成各种形式的内容。

讯飞写作
讯飞写作

基于讯飞星火大模型的AI写作工具,可以快速生成新闻稿件、品宣文案、工作总结、心得体会等各种文文稿

即梦AI
即梦AI

一站式AI创作平台,免费AI图片和视频生成。

ChatGPT
ChatGPT

最最强大的AI聊天机器人程序,ChatGPT不单是聊天机器人,还能进行撰写邮件、视频脚本、文案、翻译、代码等任务。

相关专题

更多
Python 时间序列分析与预测
Python 时间序列分析与预测

本专题专注讲解 Python 在时间序列数据处理与预测建模中的实战技巧,涵盖时间索引处理、周期性与趋势分解、平稳性检测、ARIMA/SARIMA 模型构建、预测误差评估,以及基于实际业务场景的时间序列项目实操,帮助学习者掌握从数据预处理到模型预测的完整时序分析能力。

82

2025.12.04

Python 数据清洗与预处理实战
Python 数据清洗与预处理实战

本专题系统讲解 Python 在数据清洗与预处理中的核心技术,包括使用 Pandas 进行缺失值处理、异常值检测、数据格式化、特征工程与数据转换,结合 NumPy 高效处理大规模数据。通过实战案例,帮助学习者掌握 如何处理混乱、不完整数据,为后续数据分析与机器学习模型训练打下坚实基础。

34

2026.01.31

go语言 数组和切片
go语言 数组和切片

本专题整合了go语言数组和切片的区别与含义,阅读专题下面的文章了解更多详细内容。

56

2025.09.03

TypeScript类型系统进阶与大型前端项目实践
TypeScript类型系统进阶与大型前端项目实践

本专题围绕 TypeScript 在大型前端项目中的应用展开,深入讲解类型系统设计与工程化开发方法。内容包括泛型与高级类型、类型推断机制、声明文件编写、模块化结构设计以及代码规范管理。通过真实项目案例分析,帮助开发者构建类型安全、结构清晰、易维护的前端工程体系,提高团队协作效率与代码质量。

76

2026.03.13

Python异步编程与Asyncio高并发应用实践
Python异步编程与Asyncio高并发应用实践

本专题围绕 Python 异步编程模型展开,深入讲解 Asyncio 框架的核心原理与应用实践。内容包括事件循环机制、协程任务调度、异步 IO 处理以及并发任务管理策略。通过构建高并发网络请求与异步数据处理案例,帮助开发者掌握 Python 在高并发场景中的高效开发方法,并提升系统资源利用率与整体运行性能。

116

2026.03.12

C# ASP.NET Core微服务架构与API网关实践
C# ASP.NET Core微服务架构与API网关实践

本专题围绕 C# 在现代后端架构中的微服务实践展开,系统讲解基于 ASP.NET Core 构建可扩展服务体系的核心方法。内容涵盖服务拆分策略、RESTful API 设计、服务间通信、API 网关统一入口管理以及服务治理机制。通过真实项目案例,帮助开发者掌握构建高可用微服务系统的关键技术,提高系统的可扩展性与维护效率。

345

2026.03.11

Go高并发任务调度与Goroutine池化实践
Go高并发任务调度与Goroutine池化实践

本专题围绕 Go 语言在高并发任务处理场景中的实践展开,系统讲解 Goroutine 调度模型、Channel 通信机制以及并发控制策略。内容包括任务队列设计、Goroutine 池化管理、资源限制控制以及并发任务的性能优化方法。通过实际案例演示,帮助开发者构建稳定高效的 Go 并发任务处理系统,提高系统在高负载环境下的处理能力与稳定性。

62

2026.03.10

Kotlin Android模块化架构与组件化开发实践
Kotlin Android模块化架构与组件化开发实践

本专题围绕 Kotlin 在 Android 应用开发中的架构实践展开,重点讲解模块化设计与组件化开发的实现思路。内容包括项目模块拆分策略、公共组件封装、依赖管理优化、路由通信机制以及大型项目的工程化管理方法。通过真实项目案例分析,帮助开发者构建结构清晰、易扩展且维护成本低的 Android 应用架构体系,提升团队协作效率与项目迭代速度。

109

2026.03.09

JavaScript浏览器渲染机制与前端性能优化实践
JavaScript浏览器渲染机制与前端性能优化实践

本专题围绕 JavaScript 在浏览器中的执行与渲染机制展开,系统讲解 DOM 构建、CSSOM 解析、重排与重绘原理,以及关键渲染路径优化方法。内容涵盖事件循环机制、异步任务调度、资源加载优化、代码拆分与懒加载等性能优化策略。通过真实前端项目案例,帮助开发者理解浏览器底层工作原理,并掌握提升网页加载速度与交互体验的实用技巧。

108

2026.03.06

热门下载

更多
网站特效
/
网站源码
/
网站素材
/
前端模板

精品课程

更多
相关推荐
/
热门推荐
/
最新课程
最新Python教程 从入门到精通
最新Python教程 从入门到精通

共4课时 | 22.5万人学习

Django 教程
Django 教程

共28课时 | 5万人学习

SciPy 教程
SciPy 教程

共10课时 | 2万人学习

关于我们 免责申明 举报中心 意见反馈 讲师合作 广告合作 最新更新
php中文网:公益在线php培训,帮助PHP学习者快速成长!
关注服务号 技术交流群
PHP中文网订阅号
每天精选资源文章推送

Copyright 2014-2026 https://www.php.cn/ All Rights Reserved | php.cn | 湘ICP备2023035733号