如何在Python中用NumPy实现平滑滤波算法

风强大大_3883

风强大大_3883

2026-10-09

661人浏览

原创

numpy 本身不提供开箱即用的平滑滤波器,但可通过 np.convolve 或 scipy.ndimage.convolve1d 高效实现;移动平均需手动归一化核、注意边界处理与 nan 影响。

如何在python中用numpy实现平滑滤波算法

NumPy 本身不提供开箱即用的「平滑滤波器」(如高斯、Savitzky-Golay),但能高效支撑其底层实现——关键在于用 np.convolve 或 scipy.ndimage.convolve1d(后者更稳),而非徒手写循环。

用 np.convolve 做简单移动平均平滑

这是最轻量、纯 NumPy 的方式,适合一维信号(如时间序列),但要注意边界处理和归一化。

  • 移动平均核必须手动归一化:比如长度为 5 的均值滤波,要用 np.ones(5)/5,不能只用 np.ones(5),否则结果整体放大
  • mode 参数决定输出长度:"valid" 最安全(丢弃边界,无填充),"same" 保持原长(但 NumPy 不自动补零,需手动 pad 或改用 scipy.signal.convolve)
  • 示例:
    import numpy as np<br>x = np.array([1, 2, 3, 4, 5, 4, 3, 2, 1])<br>kernel = np.ones(3) / 3<br>y = np.convolve(x, kernel, mode="same")  # 注意:mode="same" 在 NumPy 中会补零,但首尾两元素是近似值

为什么不用 np.mean 滑动窗口?

有人想用 np.mean 配合切片循环实现滑动平均,这在小数组上可行,但实际中极易出错且慢。

testing-python
testing-python

使用pytest编写和评估有效的Python测试。适用于编写测试、审查测试代码、调试测试失败或提高测试覆盖率。

下载
  • 手动管理索引易越界,尤其当窗口大小 > 数组长度时,IndexError 直接中断
  • Python 循环 + 切片触发大量内存拷贝,np.convolve 底层调用优化过的 C 实现,速度通常快 10–100 倍
  • 无法利用 CPU 向量化指令;而 convolve 内部已做 stride 和 cache 友好处理

遇到 NaN 数据怎么办?

原始信号含 np.nan 时,np.convolve 会污染整个结果(只要核覆盖任一 NaN,对应输出即为 NaN)。这不是 bug,是设计行为。

  • 先用 pd.Series(x).rolling(w).mean()(如果允许引入 pandas)——它默认跳过 NaN
  • 纯 NumPy 方案:用 np.nanmean + 手动滑动(仅当数据稀疏 NaN 且窗口小才可接受)
  • 更鲁棒的做法:用 scipy.ndimage.uniform_filter1d(x, size=5, mode="nearest"),它天然支持 mode 控制边界,并跳过 NaN(需设 origin=0 和 mode="constant" 配合 cval=np.nan,但注意版本差异)

真正难的不是写出一个能跑的平滑函数,而是选对核、控住边界、处理缺失值、并意识到 np.convolve 的 "same" 模式在首尾其实做了隐式零填充——这点常被忽略,导致信号起止处失真。如果数据敏感,务必可视化前 5 和后 5 个点的输入/输出对比。

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

相关文章

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

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

下载

相关标签:

python

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

相关专题

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

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

2023.07.20

1671

4

python能做什么
python能做什么

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

2023.07.25

4224

7

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

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

2023.07.31

1669

3

python教程
python教程

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

2023.08.03

24517

23

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

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

2023.08.04

2987

5

python eval
python eval

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

2023.08.04

3027

5

scratch和python区别
scratch和python区别

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

2023.08.11

1163

5

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

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

2023.08.10

596

4

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

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

2023.08.11

2343

5

热门下载

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

精品课程

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