怎么在Python中利用NumPy滑动窗口实现高效时序特征提取

云枫吖_6893

云枫吖_6893

2026-09-23

848人浏览

原创

应选用numpy.lib.stride_tricks.sliding_window_view,它通过调整strides和shape构造内存友好的视图,不复制数据,速度比python循环快10–100倍,输出形状为(n-w+1, w),需显式导入。

怎么在python中利用numpy滑动窗口实现高效时序特征提取

滑动窗口函数选 numpy.lib.stride_tricks.sliding_window_view 而不是手写循环

NumPy 1.20+ 内置的 sliding_window_view 是目前最直接、内存友好的选择。它不复制数据,而是通过调整数组的 stridesshape 构造视图,速度比 Python 循环快 10–100 倍,也比 np.array([arr[i:i+w] for i in range(len(arr)-w+1)]) 少大量内存开销。

常见错误是忽略其返回结果的维度变化:输入一维数组 arr(长度 N),窗口大小 w,输出是二维数组 shape (N-w+1, w),不是三维或自动降维。

  • 必须显式导入:from numpy.lib.stride_tricks import sliding_window_view
  • 若 NumPy 版本 numba.jit 加速的手写循环,或升级 NumPy
  • 对 float32 数据,sliding_window_view 返回视图仍保持原 dtype,无需额外转换

提取统计特征时优先用 axis=1 向量化计算

拿到窗口视图后,别再用 for 遍历每个窗口算均值/标准差。直接在第二轴(即每个窗口内部)调用 NumPy 的聚合函数,利用底层 C 实现批量处理。

例如计算每个窗口的均值、标准差、最大值差(max-min):

import numpy as np
from numpy.lib.stride_tricks import sliding_window_view
<p>x = np.random.randn(1000)
windows = sliding_window_view(x, window_shape=50)  # shape (951, 50)</p><p>means = windows.mean(axis=1)      # shape (951,)
stds  = windows.std(axis=1)       # shape (951,)
ranges = windows.max(axis=1) - windows.min(axis=1)  # shape (951,)</p>

注意:std 默认是样本标准差(ddof=0),若需无偏估计,显式写 windows.std(axis=1, ddof=1)axis=1 不能省略,否则会压缩整个二维数组成单个标量。

A Python CLI skill for Cutout.Pro visual APIs — background removal, face cutout, and photo enhancement. Supports file upload & image URL input.
A Python CLI skill for Cutout.Pro visual APIs — background removal, face cutout, and photo enhancement. Supports file upload & image URL input.

调用 Cutout.Pro 视觉处理 API 进行背景移除、人像抠图和照片增强,支持文件上传与图片 URL 输入。

下载

处理边界与非整除窗口时,明确用 pad 还是截断

sliding_window_view 默认只生成完整窗口,末尾不足 window_shape 的元素直接丢弃(即“valid”模式)。这在时序建模中常导致特征长度变短,影响后续对齐。

  • 需要等长输出?先用 np.pad 补零或边缘值:x_padded = np.pad(x, (0, 49), mode='edge'),再套窗口
  • 补全后注意:补的是原始序列末端,不是每个窗口单独补;mode='reflect''wrap' 在某些周期信号中更合理
  • 避免用 scipy.signal.convolve 替代——它适合线性滤波,但不直接支持任意统计量(如分位数、偏度)

高阶特征(如偏度、分位数)慎用 axis=1 + scipy.stats

scipy.stats.skewscipy.stats.mstats.mquantiles 等函数不支持 axis 参数,无法直接作用于二维窗口视图。强行 np.apply_along_axis 会退化为 Python 循环,丧失向量化优势。

可行方案:

  • 偏度:用公式手动实现(基于三阶中心矩),复用 windows.mean(axis=1)windows.var(axis=1, ddof=0)
  • 分位数:用 np.quantile(windows, q, axis=1)(NumPy ≥ 1.15),q 可为标量或数组,返回对应维度结果
  • 自定义函数且必须用 scipy?用 numba.jit(nopython=True) 包一层,输入 windows,内部用 for 但编译后接近 C 速度

真正容易被忽略的是:当窗口数超 10⁵ 且窗口宽超 100 时,即使视图不复制数据,windowsstrides 结构也可能触发某些 NumPy 函数的隐式拷贝(如部分 np.sort 变体),建议用 np.may_share_memory(x, windows) 检查是否仍是视图。

Python免费学习笔记(深入):立即使用
在学习笔记中,你将探索 Python 的核心概念和高级技巧!

相关文章

PHP速学视频免费教程(入门到精通)
PHP速学视频免费教程(入门到精通)

PHP怎么学习?PHP怎么入门?PHP在哪学?PHP怎么学才快?不用担心,这里为大家提供了PHP速学教程(入门到精通),有需要的小伙伴保存下载就能学习啦!

下载

本站声明:本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系admin@php.cn

相关专题

更多
python打包成可执行文件
python打包成可执行文件

本专题为大家带来python打包成可执行文件相关的文章,大家可以免费的下载体验。

2023.07.20

1551

4

python能做什么
python能做什么

python能做的有:可用于开发基于控制台的应用程序、多媒体部分开发、用于开发基于Web的应用程序、使用python处理数据、系统编程等等。本专题为大家提供python相关的各种文章、以及下载和课程。

2023.07.25

3704

7

format在python中的用法
format在python中的用法

Python中的format是一种字符串格式化方法,用于将变量或值插入到字符串中的占位符位置。通过format方法,我们可以动态地构建字符串,使其包含不同值。php中文网给大家带来了相关的教程以及文章,欢迎大家前来阅读学习。

2023.07.31

1569

3

python教程
python教程

Python已成为一门网红语言,即使是在非编程开发者当中,也掀起了一股学习的热潮。本专题为大家带来python教程的相关文章,大家可以免费体验学习。

2023.08.03

21157

23

python环境变量的配置
python环境变量的配置

Python是一种流行的编程语言,被广泛用于软件开发、数据分析和科学计算等领域。在安装Python之后,我们需要配置环境变量,以便在任何位置都能够访问Python的可执行文件。php中文网给大家带来了相关的教程以及文章,欢迎大家前来学习阅读。

2023.08.04

2607

5

python eval
python eval

eval函数是Python中一个非常强大的函数,它可以将字符串作为Python代码进行执行,实现动态编程的效果。然而,由于其潜在的安全风险和性能问题,需要谨慎使用。php中文网给大家带来了相关的教程以及文章,欢迎大家前来学习阅读。

2023.08.04

2667

5

scratch和python区别
scratch和python区别

scratch和python的区别:1、scratch是一种专为初学者设计的图形化编程语言,python是一种文本编程语言;2、scratch使用的是基于积木的编程语法,python采用更加传统的文本编程语法等等。本专题为大家提供scratch和python相关的文章、下载、课程内容,供大家免费下载体验。

2023.08.11

1083

5

python合并两个列表
python合并两个列表

Python是一种强大的编程语言,具有许多方便的功能和工具。在Python中,有多种方法可以合并两个列表。php中文网给大家带来了相关的教程以及文章,欢迎大家前来学习阅读。

2023.08.10

576

4

python是前端还是后端
python是前端还是后端

Python属于前端也属于后端,其灵活性和丰富的生态系统使得开发人员能够在不同的领域中灵活运用。本专题为大家提供python相关的文章、下载、课程内容,供大家免费下载体验。

2023.08.11

2063

5

热门下载

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

精品课程

更多
相关推荐
/
热门推荐
/
最新课程