主頁 > 後端開發 > Python 計算機視覺(九)—— OpenCV進行影像平滑

Python 計算機視覺(九)—— OpenCV進行影像平滑

2021-10-30 09:36:34 後端開發

參考的一些文章以及論文我都會給大家分享出來 —— 鏈接就貼在原文,論文我上傳到資源中去,大家可以免費下載學習,如果當天資源區找不到論文,那就等等,可能正在審核,審核完后就可以下載了,大家一起學習,一起進步!加油!!

目錄

1. 影像平滑概述

2. 均值濾波

(1)基本概念

(2)代碼實作

圖片添加噪聲

?進行均值濾波

3. 方框濾波

(1)基本概念

(2)代碼實作

4. 高斯濾波

(1)基本概念

(2)代碼實作

5. 中值濾波

(1)基本概念

(2)代碼實作

6.結束語

7.參考論文


1. 影像平滑概述

此處參考論文:《幾種影像平滑處理方法比較研究_關雪梅》

處理和傳輸數字影像的程序中可能會受到不同的噪聲的干擾,影像噪聲會引起影像質量降低、影像變得模糊,使影像的特征被淹沒、而通過除去噪聲來達到影像增強的目的的操作被稱為影像平滑處理;

一般來說,污染影像的噪聲的來源主要有以下三種方式:

(1)第一種是在光電、電磁轉換程序中產生的人為噪聲

(2)第二種是大氣層電(磁)暴、閃 電、電壓、浪涌等引起的強脈沖性沖激噪聲

(3)第三種是自然起伏性噪聲,由物理量的不連續性或粒子性所引起

2. 均值濾波

(1)基本概念

假設被噪聲污染的數字影像空間含有 N*N 個像素點,而這幅數字影像中每一個像素點值可以通過領域內(可以理解為使用矩形括起來)幾個像素點的平均值計算出來,采用平滑影像處理方法就可以得到一幅新的影像,這個程序就被稱為均值濾波

如果對于這種二維中的影像的平滑處理理解有點難度的話,可以參考我之前的一篇關于語音信號經過滑動平均濾除噪聲的文章,應該理解起來輕松一點:M-point moving-average(M點滑動平均)Matlab 實作

(2)代碼實作

圖片添加噪聲

首先肯定就是讀取影像資訊了:

"""
Author:XiaoMa
date:2021/10/28
"""
import cv2
import numpy as np

#讀取影像資訊
img0 = cv2.imread("E:\From Zhihu\For the desk\cvnine0.jpg")
img1 = cv2.cvtColor(img0, cv2.COLOR_BGR2GRAY)
h, w = img0.shape[:2]
print(h, w)
cv2.namedWindow("W0")
cv2.imshow("W0", img0)
cv2.waitKey(delay = 0)

得到的影像資訊如下:450 720

然后在原圖中加入噪聲:

使用 skimage.util.random_noise(image, mode=‘gaussian’, seed=None, clip=True, **kwargs)

image: 添加噪聲的影像

mode: 添加的噪聲型別(‘gaussian’ 高斯加性噪聲

‘localvar’ 高斯加性噪聲,每點具有特定的區域方差

‘poisson’ 泊松分布的噪聲

'salt’鹽噪聲,隨機用1替換像素,屬于高灰度噪聲

‘pepper’ 胡椒噪聲,隨機用0或-1替換像素,屬于低灰度噪聲

‘s&p’椒鹽噪聲,兩種噪聲同時出現,呈現出黑白雜點)

在這篇文章中我們就加入簡單的高斯白噪聲即可:

#給影像增加隨機噪聲
img2 = skimage.util.random_noise(img0, mode='gaussian')
cv2.namedWindow("W1")
cv2.imshow("W1", img2)
cv2.waitKey(delay = 0)

進行均值濾波

使用 cv2.blur() 函式

其中主要有兩個引數,需要處理的影像操作的核

#進行均值濾波
img3 = cv2.blur(img2, (6, 5))
cv2.namedWindow("W2")
cv2.imshow("W2", img3)
cv2.waitKey(delay = 0)

可以看出來濾波后的影像對于加了噪聲的影像來說確實好了一點,但沒辦法達到原圖那樣的效果;按照我個人的理解,加入噪聲就是使得一些像素點的像素值發生了改變,所以整個影像中的像素總值也發生了改變,即便求平均也會對原圖產生影響,而且還把那些之前沒有被改變的像素點也給影響了

3. 方框濾波

(1)基本概念

前面提到的均值濾波是對整個領域內的像素點求平均,再將平均值賦給目標像素點,但是方框濾波可以選擇是否對像素值進行歸一化,即可以將得到的像素點的總值賦給一個像素點

(2)代碼實作

使用 cv2.boxFilter(src,ddepth,ksize,normalize)

src: 原圖

ddepth: 目標影像深度

ksize: 核大小(和均值濾波類似)

normalize: 是否進行歸一化,使用布林值進行選擇

若我們進行歸一化處理,效果就會和均值濾波一樣,所以在這里我們只演示不進行歸一化操作:

#進行方框濾波
img4 = cv2.boxFilter(img2, ddepth = -1, ksize = (2, 2), normalize = False)
cv2.namedWindow("W3")
cv2.imshow("W3", img4)
cv2.waitKey(delay = 0)

在這里我們將核的大小設定為(2, 2),得出結果:
設定為(3, 3):
再次擴大為(5, 5):

整個影像變白了,這是為什么呢?

我們知道不進行歸一化就是將方框內的所有像素點值賦給一個像素點,而它的值是有最大值255的,如果超過這個值像素點就溢位了,變成了白色;前面我們將核值設為(x, y),那么就有 x*y個像素點的值賦給一個像素點,超過某個界限肯定會導致所有的點都溢位了,整個影像自然也就變白了

4. 高斯濾波

(1)基本概念

均值濾波中我們是將方框中的每一個像素點的值都加起來求平均然后賦給原來的像素點,這樣每一個像素點的權值都是1,但是在高斯濾波中,每一個像素點的權值以和該像素點的距離成高斯分布;即距離越遠的像素點的權值越小

(2)代碼實作

#進行高斯濾波
img5 = cv2.GaussianBlur(img2, (5, 5), 0)  #引數分別為:原圖、濾波器大小、標準差(標準差設定為0則會自動按照大小進行選擇)
cv2.namedWindow("W4")
cv2.imshow("W4", img5)
cv2.waitKey(delay = 0)

上圖中 W4 是高斯濾波的結果,W2 是均值濾波的結果(反正我本人看不出區別,哈哈哈哈)

5. 中值濾波

(1)基本概念

中值濾波將一個滑動視窗內所有像素點灰度值進行排序, 將滑動視窗中中心像素點的灰度值用中值來替代; 比如一個視窗內的像素點值為[90,89,76,89,94,93,23,54,87],那么將會使用中值89來代替原先的中點94

(2)代碼實作

使用函式 cv2.medianblur()

此處參考:opencv學習(二十一)之中值濾波medianBlur

一般引數有:

輸入影像:可規定影像深度,與濾波器大小有關

濾波器尺寸大小:必須為奇數

當濾波器大小為 3、5時,輸入影像的深度為 CV_8U、CV_16U、CV_32F中的一個,當濾波器尺寸過大時,輸入影像的深度只能是 CV_8U

CV_8U - 8-bit unsigned integers ( 0..255 )
CV_8S - 8-bit signed integers ( -128..127 )
CV_16U - 16-bit unsigned integers ( 0..65535 )
CV_16S - 16-bit signed integers ( -32768..32767 )
CV_32S - 32-bit signed integers ( -2147483648..2147483647 )
CV_32F - 32-bit floating-point numbers ( -FLT_MAX..FLT_MAX, INF, NAN )
CV_64F - 64-bit floating-point numbers ( -DBL_MAX..DBL_MAX, INF, NAN )

#進行中值濾波
img6 = cv2.imread("E:\From Zhihu\For the desk\cvnine0_0.png") #讀取的是添加了噪聲的影像
img7 = cv2.medianBlur(np.uint8(img6), 9)
cv2.namedWindow("W5")
cv2.imshow("W5", img7)
cv2.waitKey(delay = 0)

個人認為并不見得效果有多好,通過調整濾波器尺寸可以進行多次嘗試,來達到最優的一個影像平滑效果

6.結束語

這篇文章主要介紹了影像平滑中常用到的幾種濾波方法,比如均值濾波、方框濾波、高斯濾波以及中值濾波,影像平滑的主要目的是減小影像中的一些噪聲,對影像進行優化處理;當然影像平滑的代價就是失去一些影像中的細節變化,如果大家學過《信號與系統》就會知道,影像變換域中決定細節部分的是高頻部分,所以通過影像平滑,影像肯定丟失了高頻部分,這表明我們在本篇文章中構建的這幾個濾波器本質上都屬于低通濾波器,肯定會使影像變得模糊的,至于去噪和保留細節這兩個對立的面的平衡自己通過調整引數進行尋找即可,在本篇文章中參考的文章鏈接我都粘貼上去了,論文既上傳到資源部分也粘貼到了下面,需要的直接查看就可以了,加油!

7.參考論文

轉載請註明出處,本文鏈接:https://www.uj5u.com/houduan/341894.html

標籤:python

上一篇:Pythont打開 txt 格式的檔案

下一篇:7000字整理: 全網最詳細Pandas合并資料集操作總結

標籤雲
其他(157675) Python(38076) JavaScript(25376) Java(17977) C(15215) 區塊鏈(8255) C#(7972) AI(7469) 爪哇(7425) MySQL(7132) html(6777) 基礎類(6313) sql(6102) 熊猫(6058) PHP(5869) 数组(5741) R(5409) Linux(5327) 反应(5209) 腳本語言(PerlPython)(5129) 非技術區(4971) Android(4554) 数据框(4311) css(4259) 节点.js(4032) C語言(3288) json(3245) 列表(3129) 扑(3119) C++語言(3117) 安卓(2998) 打字稿(2995) VBA(2789) Java相關(2746) 疑難問題(2699) 细绳(2522) 單片機工控(2479) iOS(2429) ASP.NET(2402) MongoDB(2323) 麻木的(2285) 正则表达式(2254) 字典(2211) 循环(2198) 迅速(2185) 擅长(2169) 镖(2155) 功能(1967) .NET技术(1958) Web開發(1951) python-3.x(1918) HtmlCss(1915) 弹簧靴(1913) C++(1909) xml(1889) PostgreSQL(1872) .NETCore(1853) 谷歌表格(1846) Unity3D(1843) for循环(1842)

熱門瀏覽
  • 【C++】Microsoft C++、C 和匯編程式檔案

    ......

    uj5u.com 2020-09-10 00:57:23 more
  • 例外宣告

    相比于斷言適用于排除邏輯上不可能存在的狀態,例外通常是用于邏輯上可能發生的錯誤。 例外宣告 Item 1:當函式不可能拋出例外或不能接受拋出例外時,使用noexcept 理由 如果不打算拋出例外的話,程式就會認為無法處理這種錯誤,并且應當盡早終止,如此可以有效地阻止例外的傳播與擴散。 示例 //不可 ......

    uj5u.com 2020-09-10 00:57:27 more
  • Codeforces 1400E Clear the Multiset(貪心 + 分治)

    鏈接:https://codeforces.com/problemset/problem/1400/E 來源:Codeforces 思路:給你一個陣列,現在你可以進行兩種操作,操作1:將一段沒有 0 的區間進行減一的操作,操作2:將 i 位置上的元素歸零。最終問:將這個陣列的全部元素歸零后操作的最少 ......

    uj5u.com 2020-09-10 00:57:30 more
  • UVA11610 【Reverse Prime】

    本人看到此題沒有翻譯,就附帶了一個自己的翻譯版本 思考 這一題,它的第一個要求是找出所有 $7$ 位反向質數及其質因數的個數。 我們應該需要質數篩篩選1~$10^{7}$的所有數,這里就不慢慢介紹了。但是,重讀題,我們突然發現反向質數都是 $7$ 位,而將它反過來后的數字卻是 $6$ 位數,這就說明 ......

    uj5u.com 2020-09-10 00:57:36 more
  • 統計區間素數數量

    1 #pragma GCC optimize(2) 2 #include <bits/stdc++.h> 3 using namespace std; 4 bool isprime[1000000010]; 5 vector<int> prime; 6 inline int getlist(int ......

    uj5u.com 2020-09-10 00:57:47 more
  • C/C++編程筆記:C++中的 const 變數詳解,教你正確認識const用法

    1、C中的const 1、區域const變數存放在堆疊區中,會分配記憶體(也就是說可以通過地址間接修改變數的值)。測驗代碼如下: 運行結果: 2、全域const變數存放在只讀資料段(不能通過地址修改,會發生寫入錯誤), 默認為外部聯編,可以給其他源檔案使用(需要用extern關鍵字修飾) 運行結果: ......

    uj5u.com 2020-09-10 00:58:04 more
  • 【C++犯錯記錄】VS2019 MFC添加資源不懂如何修改資源宏ID

    1. 首先在資源視圖中,添加資源 2. 點擊新添加的資源,復制自動生成的ID 3. 在解決方案資源管理器中找到Resource.h檔案,編輯,使用整個專案搜索和替換的方式快速替換 宏宣告 4. Ctrl+Shift+F 全域搜索,點擊查找全部,然后逐個替換 5. 為什么使用搜索替換而不使用屬性視窗直 ......

    uj5u.com 2020-09-10 00:59:11 more
  • 【C++犯錯記錄】VS2019 MFC不懂的批量添加資源

    1. 打開資源頭檔案Resource.h,在其中預先定義好宏 ID(不清楚其實ID值應該設定多少,可以先新建一個相同的資源項,再在這個資源的ID值的基礎上遞增即可) 2. 在資源視圖中選中專案資源,按F7編輯資源檔案,按 ID 型別 相對路徑的形式添加 資源。(別忘了先把檔案拷貝到專案中的res檔案 ......

    uj5u.com 2020-09-10 01:00:19 more
  • C/C++編程筆記:關于C++的參考型別,專供新手入門使用

    今天要講的是C++中我最喜歡的一個用法——參考,也叫別名。 參考就是給一個變數名取一個變數名,方便我們間接地使用這個變數。我們可以給一個變數創建N個參考,這N + 1個變數共享了同一塊記憶體區域。(參考型別的變數會占用記憶體空間,占用的記憶體空間的大小和指標型別的大小是相同的。雖然參考是一個物件的別名,但 ......

    uj5u.com 2020-09-10 01:00:22 more
  • 【C/C++編程筆記】從頭開始學習C ++:初學者完整指南

    眾所周知,C ++的學習曲線陡峭,但是花時間學習這種語言將為您的職業帶來奇跡,并使您與其他開發人員區分開。您會更輕松地學習新語言,形成真正的解決問題的技能,并在編程的基礎上打下堅實的基礎。 C ++將幫助您養成良好的編程習慣(即清晰一致的編碼風格,在撰寫代碼時注釋代碼,并限制類內部的可見性),并且由 ......

    uj5u.com 2020-09-10 01:00:41 more
最新发布
  • Rust中的智能指標:Box<T> Rc<T> Arc<T> Cell<T> RefCell<T> Weak

    Rust中的智能指標是什么 智能指標(smart pointers)是一類資料結構,是擁有資料所有權和額外功能的指標。是指標的進一步發展 指標(pointer)是一個包含記憶體地址的變數的通用概念。這個地址參考,或 ” 指向”(points at)一些其 他資料 。參考以 & 符號為標志并借用了他們所 ......

    uj5u.com 2023-04-20 07:24:10 more
  • Java的值傳遞和參考傳遞

    值傳遞不會改變本身,參考傳遞(如果傳遞的值需要實體化到堆里)如果發生修改了會改變本身。 1.基本資料型別都是值傳遞 package com.example.basic; public class Test { public static void main(String[] args) { int ......

    uj5u.com 2023-04-20 07:24:04 more
  • [2]SpinalHDL教程——Scala簡單入門

    第一個 Scala 程式 shell里面輸入 $ scala scala> 1 + 1 res0: Int = 2 scala> println("Hello World!") Hello World! 檔案形式 object HelloWorld { /* 這是我的第一個 Scala 程式 * 以 ......

    uj5u.com 2023-04-20 07:23:58 more
  • 理解函式指標和回呼函式

    理解 函式指標 指向函式的指標。比如: 理解函式指標的偽代碼 void (*p)(int type, char *data); // 定義一個函式指標p void func(int type, char *data); // 宣告一個函式func p = func; // 將指標p指向函式func ......

    uj5u.com 2023-04-20 07:23:52 more
  • Django筆記二十五之資料庫函式之日期函式

    本文首發于公眾號:Hunter后端 原文鏈接:Django筆記二十五之資料庫函式之日期函式 日期函式主要介紹兩個大類,Extract() 和 Trunc() Extract() 函式作用是提取日期,比如我們可以提取一個日期欄位的年份,月份,日等資料 Trunc() 的作用則是截取,比如 2022-0 ......

    uj5u.com 2023-04-20 07:23:45 more
  • 一天吃透JVM面試八股文

    什么是JVM? JVM,全稱Java Virtual Machine(Java虛擬機),是通過在實際的計算機上仿真模擬各種計算機功能來實作的。由一套位元組碼指令集、一組暫存器、一個堆疊、一個垃圾回收堆和一個存盤方法域等組成。JVM屏蔽了與作業系統平臺相關的資訊,使得Java程式只需要生成在Java虛擬機 ......

    uj5u.com 2023-04-20 07:23:31 more
  • 使用Java接入小程式訂閱訊息!

    更新完微信服務號的模板訊息之后,我又趕緊把微信小程式的訂閱訊息給實作了!之前我一直以為微信小程式也是要企業才能申請,沒想到小程式個人就能申請。 訊息推送平臺🔥推送下發【郵件】【短信】【微信服務號】【微信小程式】【企業微信】【釘釘】等訊息型別。 https://gitee.com/zhongfuch ......

    uj5u.com 2023-04-20 07:22:59 more
  • java -- 緩沖流、轉換流、序列化流

    緩沖流 緩沖流, 也叫高效流, 按照資料型別分類: 位元組緩沖流:BufferedInputStream,BufferedOutputStream 字符緩沖流:BufferedReader,BufferedWriter 緩沖流的基本原理,是在創建流物件時,會創建一個內置的默認大小的緩沖區陣列,通過緩沖 ......

    uj5u.com 2023-04-20 07:22:49 more
  • Java-SpringBoot-Range請求頭設定實作視頻分段傳輸

    老實說,人太懶了,現在基本都不喜歡寫筆記了,但是網上有關Range請求頭的文章都太水了 下面是抄的一段StackOverflow的代碼...自己大修改過的,寫的注釋挺全的,應該直接看得懂,就不解釋了 寫的不好...只是希望能給視頻網站開發的新手一點點幫助吧. 業務場景:視頻分段傳輸、視頻多段傳輸(理 ......

    uj5u.com 2023-04-20 07:22:42 more
  • Windows 10開發教程_編程入門自學教程_菜鳥教程-免費教程分享

    教程簡介 Windows 10開發入門教程 - 從簡單的步驟了解Windows 10開發,從基本到高級概念,包括簡介,UWP,第一個應用程式,商店,XAML控制元件,資料系結,XAML性能,自適應設計,自適應UI,自適應代碼,檔案管理,SQLite資料庫,應用程式到應用程式通信,應用程式本地化,應用程式 ......

    uj5u.com 2023-04-20 07:22:35 more