OpenCV-Python 图像平滑处理1:卷积函数filter2D详解及用于均值滤波的案例

2023-11-09

一、图像平滑处理简介

图像平滑处理属于图像空间滤波的一种,用于模糊处理和降低噪声。模糊处理经常用于图像预处理任务中,例如在(大)目标提取之前去除图像中的一些琐碎细节,以及桥接直线或曲线的缝隙。模糊处理后的图像,可以通过阈值处理、形态处理等方式进行再加工,从而去除一些噪点。

平滑滤波器包括线性滤波器和非线性滤波器平滑线性空间滤波器的输出(响应)是包含在滤波器模板邻域内的像素的简单平均值。平滑线性空间滤波器有时也称为均值滤波器,它们属于低通滤波器

平滑线性滤波器的基本概念非常直观。它使用滤波器模板确定的邻域内像素的平均/加权平均灰度值代替图像中每个像素的值。所有系数都相等(非加权平均)的空间均值滤波器也称为盒状滤波器

非线性滤波器可能有多种,统计排序滤波器是常用的,如中值滤波、最小值滤波(如图像腐蚀)、最大值滤波(如图像膨胀)都属于统计排序滤波器。

更多关于图像平滑处理知识的介绍请参考《数字图像处理:线性和非线性滤波的平滑空间滤波器(Smoothing Spatial Filters)》的介绍。

二、filter2D介绍

2.1、简介

filter2D是OpenCV使用卷积核对图像进行卷积运算的函数,该函数能对图像进行任意的线性滤波处理,具体滤波方式由核矩阵确认。

该函数其实执行的是相关操作而不是卷积操作,计算公式如下:
在这里插入图片描述
关于相关和卷积的关系请参考《《数字图像处理》空间滤波学习感悟2:空间相关与卷积的概念、区别及联系》,不过对应系数相等的盒状滤波来说,由于核矩阵的对称性,卷积和相关的处理结果相同。

2.2、语法说明

语法
dst	=	cv.filter2D(	src, ddepth, kernel[, dst[, anchor[, delta[, borderType]]]]	)
参数说明
  • src:输入图像

  • ddepth:目标图像深度(请参考《图像表示的相关概念:图像深度、像素深度、位深的区别和关系》),如果目标图像深度和输入图像深度相同,则传值-1,老猿测试在Python中此时取值None、0效果也一样。针对输入图像对应的目标图像,该参数的可选传值对应关系如下:
    在这里插入图片描述

  • kernel:卷积核(convolution kernel ),如上概述所述,实际上是相关核(correlation kernel),为一个单通道的浮点数矩阵,如果针对图像不同通道需要使用不同核,则需要将图像进行split拆分成单通道并使用对应核逐个进行处理

  • dst:结果图像

  • anchor:核矩阵的锚点,用于定位核距中与当前处理像素点对齐的点,默认值(-1,-1)表示锚点位于内核中心,否则就是核矩阵锚点位置坐标,锚点位置对卷积处理的结果会有非常大的影响;

  • delta:在将卷积处理后的像素值存储到dst之前,向其添加的可选值,老猿测试验证当有值时,卷积后的像素结果值会与delta相加,得到的结果作为最终输出的像素值,注意这个加法是饱和运算,超过255的被置为255;

  • borderType:当要扩充输入图像矩阵边界时的像素取值方法,当核矩阵锚点与像素重合但核矩阵覆盖范围超出到图像外时,函数可以根据指定的边界模式进行插值运算。可选模式包括:
    在这里插入图片描述

注意

  • BORDER_WRAP在此不支持;
  • 经老猿测试,默认值为BORDER_DEFAULT ,与BORDER_REFLECT_101 、BORDER_REFLECT101相同

2.4、返回值

返回值为结果图像矩阵,因此输入参数中的dst参数无需输入。

三、使用案例

下面的案例脱胎于OpenCV帮助文档,代码对输入图像进行均值滤波:

import cv2
import numpy as np
import matplotlib.pyplot as plt

def smoothingByOpenCV():
    img = imread('f:\\pic\\opencvLogo.JPG')
    kernal = np.ones((5, 5), np.float32) / 25
    dst = cv2.filter2D(img, None, kernal)
    plt.subplot(121), plt.imshow(img), plt.title('Original')
    plt.xticks([]), plt.yticks([])
    plt.subplot(122), plt.imshow(dst), plt.title('Averaging')
    plt.xticks([]), plt.yticks([])
    plt.show()
    return dst
    
smoothingByOpenCV()

结果输出:
在这里插入图片描述

可以看到输出图像比输入图像变模糊了。

上面的代码kernal中各元素相加结果为1,没有改变图像整体的亮度,当将其改为:

kernal = np.ones((5, 5), np.float32) / 5

此时,输出结果如下:
在这里插入图片描述

可以看到图像的整体零度提升,并扩展了前景色范围,这是因为卷积过程中卷积核的元素值变大导致卷积结果值相比原值整体变大导致的。

如果不改变kernal,而改变delta参数,如:

	kernal = np.ones((5, 5), np.float32) / 25
    dst = cv2.filter2D(img, None, kernal,delta=250)

则输出图像为:
在这里插入图片描述

这是因为delta设置为250后,导致结果图像大部分像素值达到饱和导致的。

当然filter2D不只是用于均值滤波,所有线性滤波都可以实现,只需要将核矩阵根据滤波任务预置不同的元素即可。

四、小结

本文介绍了图像平滑处理及均值滤波等基础概念,并详细介绍了卷积函数filter2D的Python语法及参数,并用之进行了对图像的均值滤波处理,可以看到卷积核元素值以及相关参数如delta等对卷积处理结果的影响。

整体卷积过程是将核矩阵和处理图像从左到右、从上到下移动逐一计算像素的卷积结果过程,为了更直观的了解卷积处理,老猿用Python、numpy矩阵运算以及OpenCV-Python的图像基础操作模拟实现了一个卷积程序,其效果与filter2D基本功能完全等价。在这个过程中用到了一些小技巧,有兴趣的同好请参考《卷积处理过程模拟:用Python实现OpenCV函数filter2D等效的卷积功能》一文的介绍。

更多图像处理请参考专栏OpenCV-Python图形图像处理》、《图像处理基础知识》以及付费专栏《OpenCV-Python初学者疑难问题集》的介绍。

对于缺乏Python基础的同仁,可以通过老猿的免费专栏《 专栏:Python基础教程目录》从零开始学习Python。

如对文章内容存在疑问,可在博客评论区留言,或关注:老猿Python 微信公号发消息咨询:
老猿Python微信公号

也可以扫描博客左边栏目的微信公号二维码咨询。

写博不易,敬请支持:

如果阅读本文于您有所获,敬请点赞、评论、收藏,谢谢大家的支持!

关于老猿的付费专栏

  1. 付费专栏《https://blog.csdn.net/laoyuanpython/category_9607725.html 使用PyQt开发图形界面Python应用》专门介绍基于Python的PyQt图形界面开发基础教程,对应文章目录为《 https://blog.csdn.net/LaoYuanPython/article/details/107580932 使用PyQt开发图形界面Python应用专栏目录》;
  2. 付费专栏《https://blog.csdn.net/laoyuanpython/category_10232926.html moviepy音视频开发专栏 )详细介绍moviepy音视频剪辑合成处理的类相关方法及使用相关方法进行相关剪辑合成场景的处理,对应文章目录为《https://blog.csdn.net/LaoYuanPython/article/details/107574583 moviepy音视频开发专栏文章目录》;
  3. 付费专栏《https://blog.csdn.net/laoyuanpython/category_10581071.html OpenCV-Python初学者疑难问题集》为《https://blog.csdn.net/laoyuanpython/category_9979286.html OpenCV-Python图形图像处理 》的伴生专栏,是笔者对OpenCV-Python图形图像处理学习中遇到的一些问题个人感悟的整合,相关资料基本上都是老猿反复研究的成果,有助于OpenCV-Python初学者比较深入地理解OpenCV,对应文章目录为《https://blog.csdn.net/LaoYuanPython/article/details/109713407 OpenCV-Python初学者疑难问题集专栏目录
  4. 付费专栏《https://blog.csdn.net/laoyuanpython/category_10762553.html Python爬虫入门 》站在一个互联网前端开发小白的角度介绍爬虫开发应知应会内容,包括爬虫入门的基础知识,以及爬取CSDN文章信息、博主信息、给文章点赞、评论等实战内容。

前两个专栏都适合有一定Python基础但无相关知识的小白读者学习,第三个专栏请大家结合《https://blog.csdn.net/laoyuanpython/category_9979286.html OpenCV-Python图形图像处理 》的学习使用。

对于缺乏Python基础的同仁,可以通过老猿的免费专栏《https://blog.csdn.net/laoyuanpython/category_9831699.html 专栏:Python基础教程目录)从零开始学习Python。

如果有兴趣也愿意支持老猿的读者,欢迎购买付费专栏。

老猿Python,跟老猿学Python!

☞ ░ 前往老猿Python博文目录 https://blog.csdn.net/LaoYuanPython

本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系:hwhale#tublm.com(使用前将#替换为@)

OpenCV-Python 图像平滑处理1:卷积函数filter2D详解及用于均值滤波的案例 的相关文章

  • 从一个数组中删除另一个数组中的元素

    假设我有这些二维数组A and B 我怎样才能从中删除元素A那些在B 集合论中的补集 A B A np asarray 1 1 1 1 1 2 1 1 3 1 1 4 B np asarray 0 0 0 1 0 2 1 0 3 1 0 4
  • 在Linux中的端口80上运行flask[重复]

    这个问题在这里已经有答案了 也许以前有过这个问题的答案 所以请重定向我 如果是这样的话 我正在考虑在端口 80 上运行 Flask 所以我检查了是否有任何东西正在使用端口 80 因为事实证明端口 80 没有运行 所以当我输入以下内容时 if
  • matplotlib get_color 用于子图

    我正在按照这里的教程进行操作 https matplotlib org gallery ticks and spines multiple yaxis with spines html https matplotlib org galler
  • python 打开相对文件夹中所有以.txt结尾的文件

    我需要打开并解析文件夹中的所有文件 但我必须使用相对路径 类似于 input files 我知道在 JavaScript 中你可以使用 path 库来解决这个问题 我怎样才能在Python中做到这一点 这样您就可以获得路径中的文件列表作为列
  • 如何逐行替换(更新)文件中的文本

    我试图通过读取每一行 测试它 然后写入是否需要更新来替换文本文件中的文本 我不想保存为新文件 因为我的脚本已经先备份文件并对备份进行操作 这是我到目前为止所拥有的 我从 os walk 获取路径 并且保证 pathmatch var 正确返
  • 将元组列表转换为字符串 Python

    例如 我用 python 编写了一个返回列表的函数 1 1 2 2 3 3 但我希望输出为字符串 这样我就可以用另一个字符替换逗号 这样输出就是 1 1 2 2 3 3 有什么简单的方法可以解决这个问题吗 感谢您提前提供任何提示 这看起来像
  • PRAW 出现 SSLError?

    我正在尝试开始使用 PRAW 但在使用 login 时遇到问题 我有以下代码 import praw r praw Reddit This is a test bot r login myRedditUsername password 我收
  • 在 Pyinstaller、语音识别和 Pyttsx3 中使用“-w”时,PySimpleGUI 中出现“OSError:[WinError 6] 句柄无效”

    所以我用 PySimpleGUI 创建了一个程序 然后用 Pyinstaller 从它创建了 exe 文件 这是我的命令 pyinstaller hidden import pyttsx3 drivers hidden import pyt
  • Django 单元测试数据库没有被拆除?

    我编写了一些单元测试来测试我的 Django 应用程序 特别是一个测试套件中包含大量代码setUp 功能 所述代码的目的是为数据库创建测试数据 是的 我了解固定装置 并且选择在这种情况下不使用它们 当我运行单元测试套件时 运行的第一个测试通
  • 提取二值图像中的最中心区域

    我正在处理二进制图像 之前使用此代码来查找二进制图像中的最大区域 Use the hue value to convert to binary thresh 20 thresh thresh img cv2 threshold h thre
  • Python 日志记录 - 如何检查记录器是否为空

    我刚刚在我的应用程序中实现了日志记录 我想知道是否有一种方法可以检查记录器是否为空 我的想法是在我的脚本中设置两个处理程序 一个用于带水平仪的控制台WARNING 一个用于带级别的文件DEBUG 在脚本的最后 我需要检查是否CONSOLE记
  • Python:使用Excel CSV文件仅读取某些列和行

    虽然我可以读取 csv 文件而不是读取整个文件 但如何仅打印某些行和列 想象一下这是 Excel A B C D E State Heart Disease Rate Stroke Death Rate HIV Diagnosis Rate
  • 有什么理由不在Python中混合使用多处理和线程模块

    我正在考虑使用Python来实现一个需要大量多线程的程序 另一个要求是它将在桌面上运行 因此拥有许多进程将使应用程序显得混乱且难以杀死 在任务管理器中 因此 我正在考虑使用线程和多处理模块来减少进程数量 据我了解 GIL 仅适用于单个进程
  • 从 Cython 代码生成 SIMD 指令

    我需要概述在高性能数字代码中使用 Cython 可以获得的性能 我感兴趣的事情之一是找出优化的 C 编译器是否可以对 Cython 生成的代码进行矢量化 所以我决定写下面的小例子 import numpy as np cimport num
  • 在 CSV 文件的最上面一行写入

    我有这个sample csv 文件 a 1 apple b 2 banana c 3 cranberry d 4 durian e 5 eggplant 并有以下代码 samplefile open sample csv rb rows s
  • 无法从源 pylance 解析导入烧瓶

    我正在学习 Python 课程的一部分是使用 Flask 设置网络服务器 我按照 Flask 安装文档执行了步骤 由于某种原因 flask 模块带有下划线 如下所示 当我将鼠标悬停时 我会得到如下附加信息 无法从源 pylance 解析导入
  • 从 Python 访问 802.11 无线管理帧

    我想从 Linux 上的 Python 嗅探 802 11 管理 探测请求 帧 这可以从 Scapy 中实现 如下所示 coding utf 8 from scapy all import def proc p if p haslayer
  • 如何从集合中检索元素而不删除它?

    假设如下 gt gt gt s set 1 2 3 我如何获得一个值 任何值 s不做s pop 我想将该项目保留在集合中 直到我确定可以删除它 这只有在异步调用另一个主机之后才能确定 又快又脏 gt gt gt elem s pop gt
  • Pandas 将时间序列数据重新采样为 15 分钟和 45 分钟 - 使用多索引或列

    我有一些时间序列数据作为 Pandas 数据框 它从每小时过去 15 分钟和过去 45 分钟 时间间隔为 30 分钟 的观察开始 然后将频率更改为每分钟 我想对数据进行重新采样 以便整个数据帧的频率为每 30 分钟一次 15 点和 45 点
  • Django Python - LDAP 身份验证

    我目前正在研究 Django Python 我的目标是从 Ldap 目录对用户进行身份验证 我确实有 python 代码来访问 ldap 目录并检索信息 Code import ldap try l ldap open ldap forum

随机推荐

  • Storybook(一)

    TOC 欢迎使用Storybook 你好 这是你第一次使用 Storybook 这个工具 如果你想学习如何使用Storybook 可以仔细阅读这篇文章 了解一下Storybook的基本知识 基本用法 我们对Storybook进行了一些说明
  • 简单的个人介绍网页-主页面【附代码】

    主页面 代码1 style css nav height 41px border top 3px solid b4fffa border bottom 1px solid edeef0 background color fcfcfc lin
  • 微机原理中操作数位置总结

    立即数 立即数在CS里 CS在存储器里 但和存储器操作数有本质的区别 但是不属于存储器数的一种 本质上是立即数寻址而不是存储器寻址 寄存器数 在cup内部 AX BX CX DX属于EU部件 存储器操作数 1看他有没有段超越前缀 超越到了什
  • 在职的我竟然一次通过了注册测绘师考试(注册测绘师备考经验分享)

    先上一波成绩吧 首先说一下我是万万没想到案例能过的 本来做好二战准备的 认为案例答的不够全面 时间不够了 既然过了那我就回想下从几个方面简短的整理下经验吧 没过的话叫教训 希望帮助给后面的同学 时间准备 确切的开始时间记不清了 但应该是20
  • Unity SpriteAtlas 打包AssetBundle的一些问题总结

    1 新版sprite要打包成图集 需要手动创建SpriteAtlas资源 然后选择需要打包进图集资源的Sprite目录或资源文件 参考 https docs unity3d com Manual class SpriteAtlas html
  • 从入门到入土:[SEED-Lab]- Packet Sniffing and Spoofing Lab

    此博客仅用于记录个人学习进度 学识浅薄 若有错误观点欢迎评论区指出 欢迎各位前来交流 部分材料来源网络 若有侵权 立即删除 本人博客所有文章纯属学习之用 不涉及商业利益 不合适引用 自当删除 若被用于非法行为 与我本人无关 SEED Lab
  • UE4 后处理材质 镜头光晕(Lensflare)效果

    介绍 基于PostProcessMaterial实现镜头光晕 Lensflare 的效果 不使用UE自带的Lensflare 在ShaderToy上发现一个好看的 ShaderToy地址 ShaderToyLensflareSample 将
  • Java面向对象编程

    上网的时候 访问某个网页却突然出现了某个运营商的网页 如联通 电信 出现此问题可能的原因是 A DNS劫持 B DDoS攻击 C MAC地址欺骗 D 伪造DHCP服务器 答案 A TCP建立连接的三次握手中 第二次握手发送的包会包含的标记
  • Vue3 使用 jsx

    使用 jsx 格式文件 和 defineComponent 引用自定义组件 传递属性 基本使用 child jsx 组件的配置 demo jsx setup 函数 父组件引用 import demo from components demo
  • 课堂实验-【集成学习】

    实验名称 实验 集成学习 实验目的 1 了解决策树 随机森林理论基础 2 平台实现算法 3 编程实现决策树 随机森林算法 实验原理 决策树 Decision Tree 是在已知各种情况发生概率的基础上 通过构成决策树来求取净现值的期望值大于
  • 【pandas】TypeError: concat() got multiple values for argument 'axis'

    TypeError concat got multiple values for argument axis 原代码 var LotArea data pd concat df train SalePrice df train var ax
  • Hive--添加/更改表字段类型

    Hive添加表字段 新增字段表 alter table 表名 add columns 字段名 数据类型 alter table table name add columns c time string comment 当前时间 正确 添加在
  • Mock.js配合vuex进行登录获取token

    新项目 后端数据还没准备好 没办法直接从后端获取接口数据 于是我利用Mock简单实现一下这个功能 以作为测试用 首先 利用 cnpm install mockjs save dev 指令安装完mock后 在项目里创建文件夹mock 再创建i
  • 安装python软件找不到setuptools的解决方案

    from http blog csdn net newjueqi article details 47446965 在新装的centos系统中 用源码安装supervisor提示如下 plain view plain copy root 1
  • jeecgboot使用例子

    文章目录 java代码示例 方法请求例子 switch case 语句 判断包含多个实体的List中 是实体之间的userid重复 java8如何用Stream查List对象某属性是否有重复 JAVA判断字符串是否为null 是否为空 my
  • 恭喜广工荣获2023CUBAL冠军

    恭喜广工 作者 良知犹存 转载授权以及围观 欢迎添加微信公众号 羽林君 1 竞技比赛就有这样的刺激 也许到最后一刻才能决出胜负 一个赛季努力在最后几秒进行抉择 让人热血沸腾 尤其是人们也喜欢弱者的逆袭 再一次面对清华 一支来自由全国顶级篮球
  • Vue 对象数组修改数组中的对象的某个元素,元素更新了视图不更新的问题解决

    问题描述 当我修改数组对象中的某个元素的时候 元素更新了 但是视图没有更新 例如下数组对象 attendanceTime workDays 周一 shiftId 1 shiftTime workDays 周二 shiftId 2 shift
  • 使用uni-app开发叮咚买菜的一些笔记

    本项目纯开源免费 商品信息均来自网络 如有侵权 请联系作者 我会及时删除 前一段 看到一个大佬写的叮咚买菜 挺好 于是乎自己在家也琢磨了一下 自己搞了一版本 记录一些uni app的开发事项 纯属个人总结 不对之处 还请指教 功能不是很全面
  • 硬件工程师-MOS管

    目录 DS等效电阻 Rdson MOS管 GS电容 MOS管的损耗 MOS管体二极管 讨论MOSFET的GS电容问题 MOSFET的GS电容对管子开通特性的影响 MOS管GS下拉电阻 MOS管的等效模型 MOSFET是如何进行导通的 Rds
  • OpenCV-Python 图像平滑处理1:卷积函数filter2D详解及用于均值滤波的案例

    一 图像平滑处理简介 图像平滑处理属于图像空间滤波的一种 用于模糊处理和降低噪声 模糊处理经常用于图像预处理任务中 例如在 大 目标提取之前去除图像中的一些琐碎细节 以及桥接直线或曲线的缝隙 模糊处理后的图像 可以通过阈值处理 形态处理等方