基于Numpy.convolve使用Python实现滑动平均滤波的思路详解


Posted in Python onMay 16, 2019

基于Numpy.convolve使用Python实现滑动平均滤波的思路详解基于Numpy.convolve使用Python实现滑动平均滤波的思路详解

1.滑动平均概念

滑动平均滤波法(又称递推平均滤波法),时把连续取N个采样值看成一个队列 ,队列的长度固定为N ,每次采样到一个新数据放入队尾,并扔掉原来队首的一次数据.(先进先出原则)  把队列中的N个数据进行算术平均运算,就可获得新的滤波结果。N值的选取:流量,N=12;压力:N=4;液面,N=4~12;温度,N=1~4

优点:  对周期性干扰有良好的抑制作用,平滑度高  适用于高频振荡的系统 

缺点:  灵敏度低  对偶然出现的脉冲性干扰的抑制作用较差  不易消除由于脉冲干扰所引起的采样值偏差  不适用于脉冲干扰比较严重的场合  比较浪费RAM 

2.解决思路

可以发现滑动平均滤波法计算很类似与一维卷积的工作原理,滑动平均的N就对应一维卷积核大小(长度)。

步长会有些区别,滑动平均滤波法滑动步长为1,而一维卷积步长可以自定义。还有区别就是一维卷积的核参数是需要更新迭代的,而滑动平均滤波法核参数都是一。

我们应该怎么利用这个相似性呢?其实也很简单,只需要把一维卷积核大小(长度)和N相等,步长设置为1,核参数都初始为1就可以了。由于一维卷积具备速度快,然后我们就可以使用一维卷积来实现这个功能了,快速高效。

使用深度学习框架实现这个功能是否有些大材小用了?是有些大材小用了,因为这里使用卷积的核参数不用更新,其实没必要使用复杂的深度学习框架,如果Numpy中可以实现这些功能就更简单方便了。

说干就干,经过查找发现Numpy.convolve可以实现我们想要的功能。

3.Numpy.convolve介绍

numpy.convolve(a, v, mode=‘full')

参数:

a:(N,)输入的一维数组

v:(M,)输入的第二个一维数组

mode:{‘full', ‘valid', ‘same'}参数可选

‘full' 默认值,返回每一个卷积值,长度是N+M-1,在卷积的边缘处,信号不重叠,存在边际效应。

‘same' 返回的数组长度为max(M, N),边际效应依旧存在。

‘valid'  返回的数组长度为max(M,N)-min(M,N)+1,此时返回的是完全重叠的点。边缘的点无效。

和一维卷积参数类似,a就是被卷积数据,v是卷积核大小。

4.算法实现

def np_move_avg(a,n,mode="same"):
  return(np.convolve(a, np.ones((n,))/n, mode=mode))

原理说明

运行平均值是卷积数学运算的一个例子。对于运行平均值,沿着输入滑动窗口并计算窗口内容的平均值。对于离散的1D信号,卷积是相同的,除了代替计算任意线性组合的平均值,即将每个元素乘以相应的系数并将结果相加。那些系数,一个用于窗口中的每个位置,有时称为卷积核。现在,N值的算术平均值是(x_1 + x_2 + ... + x_N) / N,所以相应的内核是(1/N, 1/N, ..., 1/N),这正是我们通过使用得到的np.ones((N,))/N。

边缘处理

该mode的参数np.convolve指定如何处理边缘。在这里选择了same模式,这样可以保证输出长度一种,但你可能还有其他优先事项。这是一个说明模式之间差异的图:

import numpy as np
import matplotlib.pyplot as plt
def np_move_avg(a,n,mode="same"):
  return(np.convolve(a, np.ones((n,))/n, mode=mode))
modes = ['full', 'same', 'valid']
for m in modes:
  plt.plot(np_move_avg(np.ones((200,)), 50, mode=m));
plt.axis([-10, 251, -.1, 1.1]);
plt.legend(modes, loc='lower center');
plt.show() ​

基于Numpy.convolve使用Python实现滑动平均滤波的思路详解基于Numpy.convolve使用Python实现滑动平均滤波的思路详解

5.参考

1. https://stackoverflow.com/questions/13728392/moving-average-or-running-mean

总结

以上所述是小编给大家介绍的Python实现滑动平均滤波的思路详解(基于Numpy.convolve),希望对大家有所帮助,如果大家有任何疑问欢迎给我留言,小编会及时回复大家的!

Python 相关文章推荐
Python对列表去重的多种方法(四种方法)
Dec 05 Python
numpy向空的二维数组中添加元素的方法
Nov 01 Python
Python3 导入上级目录中的模块实例
Feb 16 Python
django配置连接数据库及原生sql语句的使用方法
Mar 03 Python
Django中信号signals的简单使用方法
Jul 04 Python
python爬虫-模拟微博登录功能
Sep 12 Python
python scrapy重复执行实现代码详解
Dec 28 Python
使用python turtle画高达
Jan 19 Python
Python logging日志库空间不足问题解决
Sep 14 Python
Pyinstaller打包Scrapy项目的实现步骤
Sep 22 Python
浅析Python的命名空间与作用域
Nov 25 Python
如何利用opencv判断两张图片是否相同详解
Jul 07 Python
Python实现Linux监控的方法
May 16 #Python
计算机二级python学习教程(3) python语言基本数据类型
May 16 #Python
Python实现爬取亚马逊数据并打印出Excel文件操作示例
May 16 #Python
Python实现的文轩网爬虫完整示例
May 16 #Python
计算机二级python学习教程(2) python语言基本语法元素
May 16 #Python
计算机二级python学习教程(1) 教大家如何学习python
May 16 #Python
详解Python传入参数的几种方法
May 16 #Python
You might like
PHP执行linux系统命令的常用函数使用说明
2010/04/27 PHP
php写入数据到CSV文件的方法
2015/03/14 PHP
深入理解PHP中的Streams工具
2015/07/03 PHP
PHP排序算法之基数排序(Radix Sort)实例详解
2018/04/21 PHP
PHP 观察者模式深入理解与应用分析
2019/09/25 PHP
用JTrackBar实现的模拟苹果风格的滚动条
2007/08/06 Javascript
js获取class的所有元素
2013/03/28 Javascript
JSON序列化与解析原生JS方法且IE6和chrome测试通过
2013/09/05 Javascript
JavaScript简单下拉菜单实例代码
2015/09/07 Javascript
JS使用正则表达式除去字符串中重复字符的方法
2015/11/05 Javascript
Bootstrap的popover(弹出框)2秒后定时消失的实现代码
2017/02/27 Javascript
JS实现简单获取最近7天和最近3天日期的方法
2018/04/18 Javascript
详解如何在nuxt中添加proxyTable代理
2018/08/10 Javascript
JS实现换肤功能的方法实例详解
2019/01/30 Javascript
小程序实现列表多个批量倒计时
2021/01/29 Javascript
微信小程序保存多张图片的实现方法
2019/03/05 Javascript
解决layui的input独占一行的问题
2019/09/10 Javascript
jQuery实现轮播图效果demo
2020/01/11 jQuery
[13:40]TI3青蛙君全程回顾 DOTA2我们为梦想再战
2013/09/13 DOTA
python通过ftplib登录到ftp服务器的方法
2015/05/08 Python
django轻松使用富文本编辑器CKEditor的方法
2017/03/30 Python
Python开发的HTTP库requests详解
2017/08/29 Python
详解PANDAS 数据合并与重塑(join/merge篇)
2019/07/09 Python
python django 原生sql 获取数据的例子
2019/08/14 Python
Python 70行代码实现简单算式计算器解析
2019/08/30 Python
python pygame实现球球大作战
2019/11/25 Python
使用PyOpenGL绘制三维坐标系实例
2019/12/24 Python
python两种获取剪贴板内容的方法
2020/11/06 Python
python RSA加密的示例
2020/12/09 Python
canvas实现有递增动画的环形进度条的实现方法
2019/07/10 HTML / CSS
Ancheer官方户外和运动商店:销售电动自行车
2019/08/07 全球购物
车贷收入证明范本
2014/01/09 职场文书
消防工作实施方案
2014/06/09 职场文书
学习党的群众路线教育实践活动剖析材料
2014/10/13 职场文书
2014年销售员工作总结
2014/12/01 职场文书
2016党员干部廉政准则学习心得体会
2016/01/20 职场文书