基于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实现随机密码字典生成器示例
Apr 09 Python
pycharm 使用心得(七)一些实用功能介绍
Jun 06 Python
Python编程入门的一些基本知识
May 13 Python
使用Python发送各种形式的邮件的方法汇总
Nov 09 Python
利用python微信库itchat实现微信自动回复功能
May 18 Python
python递归打印某个目录的内容(实例讲解)
Aug 30 Python
python读取word 中指定位置的表格及表格数据
Oct 23 Python
python 字典套字典或列表的示例
Dec 16 Python
浅谈python的elementtree模块处理中文注意事项
Mar 06 Python
浅谈keras 模型用于预测时的注意事项
Jun 27 Python
Django中F函数的使用示例代码详解
Jul 06 Python
Python3以GitHub为例来实现模拟登录和爬取的实例讲解
Jul 30 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
echo, print, printf 和 sprintf 区别
2006/12/06 PHP
PHP实现手机号码中间四位用星号(*)隐藏的自定义函数分享
2014/09/27 PHP
php的ddos攻击解决方法
2015/01/08 PHP
php微信浏览器分享设置以及回调详解
2016/08/01 PHP
基本jquery的控制tabs打开的数量的代码
2010/10/17 Javascript
javascript密码强度校验代码(两种方法)
2015/08/10 Javascript
JS正则表达式比较常见用法
2016/01/26 Javascript
对js中回调函数的一些看法
2016/08/29 Javascript
JavaScript Uploadify文件上传实例
2017/02/28 Javascript
JS实现颜色动态淡化效果
2017/03/06 Javascript
JS实现图片预加载之无序预加载功能代码
2017/05/12 Javascript
基于Vuex无法观察到值变化的解决方法
2018/03/01 Javascript
Bootstrap的aria-label和aria-labelledby属性实例详解
2018/11/02 Javascript
原生JS检测CSS3动画是否结束的方法详解
2019/01/27 Javascript
vue中watch和computed为什么能监听到数据的改变以及不同之处
2019/12/27 Javascript
Node.js API详解之 timer模块用法实例分析
2020/05/07 Javascript
JS数据类型判断的几种常用方法
2020/07/07 Javascript
antd日期选择器禁止选择当天之前的时间操作
2020/10/29 Javascript
[03:00]《DAC最前线》之欧美新秀VS老将
2015/02/01 DOTA
python实现360皮肤按钮控件示例
2014/02/21 Python
利用Python脚本生成sitemap.xml的实现方法
2017/01/31 Python
python opencv之分水岭算法示例
2018/02/24 Python
python+opencv识别图片中的圆形
2020/03/25 Python
python中利用numpy.array()实现俩个数值列表的对应相加方法
2019/08/26 Python
Python with关键字,上下文管理器,@contextmanager文件操作示例
2019/10/17 Python
Python代码执行时间测量模块timeit用法解析
2020/07/01 Python
如何基于Python实现word文档重新排版
2020/09/29 Python
购买瑞典当代设计的腕表和太阳眼镜:TRIWA
2016/10/30 全球购物
Fanatics英国官网:美国体育电商
2018/11/06 全球购物
澳大利亚拥有最好的家具和家居用品在线目的地:Nestz
2019/02/23 全球购物
入学申请自荐信范文
2014/02/26 职场文书
企业金融服务方案
2014/06/03 职场文书
事业单位个人查摆问题及整改措施
2014/10/28 职场文书
乐山大佛导游词
2015/02/02 职场文书
法院个人总结
2015/03/03 职场文书
解决Pytorch修改预训练模型时遇到key不匹配的情况
2021/06/05 Python