基于Matlab的各种图像滤波Filter算法（代码开源）_matlab 图像滤波_ct图像 滤波 matlab源码

1.2 图像噪声

在数字图像处理中，由于成像系统、传输介质和记录设备等的不完善，数字图像在其形成、传输和记录过程中往往会受到多种噪声的污染。**这些噪声在图像上常表现为一些孤立像素点或像素块，它们以无用的信息形式出现，扰乱了图像的可观测信息。**因此，图像滤波核心目的是消除这些噪声，以便更好地提取出图像的特征。

★常见的图像噪声分类：（1）椒盐噪声；（2）高斯噪声；（3）泊松噪声

1.2.1 椒盐噪声

椒盐噪声(salt-and-pepper noise)又称脉冲噪声，它随机改变一些像素值，在二值图像上表现为使一些像素点变白，一些像素点变黑。 是由图像传感器，传输信道，解码处理等产生的黑白相间的亮暗点噪声，也就是老人们比较熟悉的所谓“雪花”。作者补充：RGB图像上的椒盐噪声也可以为红绿蓝（RGB）的像素点。

matlab代码：

%%椒盐噪声
close all;
clear all;
clc;
I=imread('test.png');
I=im2double(I);
% J = imnoise(I,‘salt & pepper’,d),默认d=0.05
% 在添加类型为Salt & Pepper的噪声时,符号&的前面和后面必须有空格,否则系统会出错
J=imnoise(I,'salt & pepper',0.01);
K=imnoise(I,'salt & pepper',0.05);
% 图像中黑色的像素点为椒盐噪声,白色的像素点为盐噪声
subplot(131),imshow(I);
subplot(132),imshow(J);
subplot(133),imshow(K);
1
2
3
4
5
6
7
8
9
10
11
12
13
14

1.2.2 高斯噪声

高斯噪声是一种源于电子电路噪声和由低照明度或高温带来的传感器噪声。高斯噪声又称为正态噪声，是自然界中最常见的噪声。高斯噪声可以通过平滑滤波技术或图像复原技术来消除。

matlab代码：

%%高斯滤波
close all;
clear all;
clc;
I=imread('test1.jpg');
h=0:0.1:1; % h为在[0,1]之间的向量,表示图像的亮度值
v=0.01:-0.001:0; % v为一个长度和h相同,表示与h中亮度对应的高斯噪声的方差
% J=imnoise(I,'localvar',h,v)在图像的不同亮度值上叠加不同方差的高斯噪声,向量h中没有的亮度值将自动插值得到
J=imnoise(I,'localvar',h,v);
subplot(121),imshow(I);
subplot(122),imshow(J);
1
2
3
4
5
6
7
8
9
10
11

1.2.3 泊松噪声

泊松噪声则是满足泊松分布的噪音，你会觉得它和正态分布很相似，其实如果我们采集的数据越多，精度越密，其形态上它越发接近高斯分布函数，也就是正态分布，是常见的一种满足指数函数分布的离散模型。

作者补充：噪声的添加效果一定程度上依赖初始图像，泊松噪声和高斯噪声看似一模一样，其实还是有些区别的。椒盐噪声和高斯噪声（指数型噪声）存在明显的区别，特别有随机性的感觉。读者朋友可以使用自己的图像进行测试，好好熟悉一下这些噪声！

matlab代码：

%%泊松噪声
close all;
clear all;
clc;
I=imread('test3.png');
% J = imnoise(I,‘poisson’)添加泊松噪声给图像I
J=imnoise(I,'poisson');
subplot(121),imshow(I);
subplot(122),imshow(J);
1
2
3
4
5
6
7
8
9

1.3 图像滤波分类

常见的线性滤波器：
**（1）低通滤波器：**允许低频率通过；
**（2）高通滤波器：**允许高频率通过； 
**（3）带通滤波器：**允许一定区域的频率通过；
**（4）带阻滤波器：**阻止一定范围内的频率并且允许其他频率通过；
**（5）全通滤波器：**允许所有频率通过，仅仅改变相位；
**（6）陷波滤波器 Band stop filter：**阻止一个狭窄频率范围通过的特殊带阻滤波器。

★根据数学性质，图像滤波算法通常分为：线性滤波与非线性滤波

均值滤波–>        Blur函数 –>                        线性滤波
高斯滤波–>        GaussianBlur函数 –>         线性滤波
中值滤波–>        MedianBlur函数 –>            非线性滤波
双边滤波–>        BilateralFilter函数 –>         非线性滤波
引导滤波–>        GuideFilter函数 –>            非线性滤波

1.4 线性滤波原理

图像线性滤波原理：每个像素的输出值是输入像素的加权和。即利用给定像素周围的像素值的决定此像素的最终输出值的一种算子，本质上就是卷积核的使用！！！卷积就是图像处理中平滑（smoothing）/模糊（bluring）最大的幕后BOSS！

假设有 7x7 的图像像素点矩阵。卷积滤波过程：7x7 上面是个3x3的窗口，从左向右，从上向下移动，将图像中像素点与滤波窗口中的数据进行卷积计算得到数据作为新的像素点值保留。每次移动一个像素格**（滤波窗口与步伐可以自己设定）。**

二、均值滤波

2.1 算法概述

均值滤波是最简单的平滑滤波，通常图片中一个方块区域（一般为 3*3）内，中心点的像素为全部点像素值的平均值。均值滤波就是对于整张图片进行以上操作。

取卷积核（Kernel）区域下所有像素的平均值并替换中心元素，如下公式：

均值滤波算法案例图：

2.2 代码与分析

matlab代码：

%%均值滤波 
close all;
clear all;
clc;
I = imread('Noise Image.jpg');  % 读取图像 
%I = rgb2gray(I);  % 转换为灰度图像（如果需要）  
% 定义均值滤波器  
h = ones(3,3) / 9;  % 3x3的均值滤波器   
J = imfilter(I, h);  % 应用滤波器   
% 显示原始图像和滤波后的图像  
subplot(1,2,1), imshow(I), title('初始图像');  
subplot(1,2,2), imshow(J), title('均值滤波后的图像');
1
2
3
4
5
6
7
8
9
10
11
12

均值滤波图：

**优点：**均值滤波算法简单，计算速度快，算法容易实现。

**缺点：**由于均值滤波是通过对邻域像素的平均值来替代中心像素的值，因此容易导致图像细节的丢失；均值滤波会平滑图像中的高频部分，包括边缘和细节，使得图像边缘变得模糊。

三、中值滤波

3.1 算法概述

中值滤波是一种非线性信号处理技术，它将像素邻域内的所有像素值进行排序，并将中值作为输出。中值滤波对去除椒盐噪声特别有效。

3.2 代码与分析

matlab代码：

%%中值滤波
close all;
clear all;
clc;
I = imread('Noise Image.jpg');  % 读取图像  
I = rgb2gray(I);  %转换为灰度图像（如果需要）  
J = medfilt2(I);  % 使用2D中值滤波器  % 应用中值滤波    
% 显示原始图像和中值滤波后的图像
subplot(1,2,1), imshow(I), title('原始图像');  
subplot(1,2,2), imshow(J), title('中值滤波后的图像');
1
2
3
4
5
6
7
8
9
10

中值滤波图：

**优点：**中值滤波可以有效地去除图像中的椒盐噪声、加性高斯噪声等噪声类型，使图像更清晰；相比其他滤波方法，中值滤波能够在去噪的同时较好地保留图像的边缘信息；中值滤波的计算过程也是简单容易的。

**缺点：**对于图像中含有较多点、线、尖角细节的情况，中值滤波可能不适用，因为它可能会改变这些细节的像素值；中值滤波是一种非线性滤波方法，它通过将像素邻域内的所有像素值排序来选择中值，这可能导致图像的不连续性。

四、高斯滤波

4.1 算法概述

高斯滤波(Gaussian filter) 包含许多种，包括低通、带通和高通等，我们通常图像上说的高斯滤波，指的是高斯模糊(Gaussian Blur) ，是一种高斯低通滤波 ，其过滤调图像高频成分（图像细节部分），保留图像低频成分（图像平滑区域），所以对图像进行 ‘高斯模糊’ 后，图像会变得模糊。高斯模糊对于抑制高斯噪声 (服从正态分布的噪声) 非常有效。

★高斯滤波和均值滤波一样，都是利用一个滤波窗口和图像进行卷积求解，不同之处在于：

• 均值滤波器的滤波窗口系数都是相同的，为1。
• 而高斯滤波器的滤波窗口系数，随着距离模板中心距离的增大，系数减小（服从二维高斯分布）。

高斯函数分布：

公式：

分布：高斯分布（正态分布）

**正态分布（高斯分布）**是一种钟形曲线，越接近中心，取值越大，越远离中心，取值越小。计算平滑结果时，只需要将"中心点"作为原点，其他点按照其在正态曲线上的位置，分配权重，就可以得到一个加权平均值。

4.2 代码与分析

matlab代码：

%%高斯滤波
close all;
clear all;
clc; 
I = imread('Noise Image.jpg');  % 读取图像   
%I = rgb2gray(I);  % 转换为灰度图像（如果需要）    
J = imgaussfilt(I, 2);  % 使用标准差为2的高斯滤波器  % 应用高斯滤波器    
% 显示原始图像和高斯滤波后的图像  
subplot(1,2,1), imshow(I), title('原始图像');  
subplot(1,2,2), imshow(J), title('高斯滤波后的图像');
1
2
3
4
5
6
7
8
9
10

高斯滤波图：

**优点：**高斯滤波能够平滑图像，去除噪声，并保留图像的整体结构；高斯滤波可以通过调整高斯核的标准差来控制，灵活度高；高斯滤波在离散化窗口滑窗卷积实现时计算速度快，适用于实时处理。

**缺点：**高斯滤波器会平滑图像中的高频部分，包括边缘和细节，可能导致图像边缘模糊；高斯滤波需要进行卷积运算，计算复杂度相对较高，可能会消耗大量的计算资源和时间。

五、双边滤波

5.1 算法概述

双边滤波由 C. Tomasi 在1998 年提出，是一种经典的非线性空间滤波方法。在滤波器稀疏的制定上，双边滤波同时考虑到了输出像素与邻域内其它像素的欧氏距离和取值的差异，即：同时考虑到了空间域和值域间的差别。

（1）空间域核

由像素位置欧式距离决定的模板权值 。

为模板窗口的其他系数的坐标；

为模板窗口的中心坐标点；

**是空间域上的高斯核函数的标准差，用于控制像素位置的权值。

使用公式生成的滤波器模板和高斯滤波器使用的模板是没有区别的。权值称为定义域核，也称为空间系数、空间域。显示由的计算公式可知，它是计算临近点q到中心点p临近程度，因此定义域核是用于衡量空间临近的程度。

（2）值域核

由像素值的差值决定的模板权值 。

其中，为模板窗口的其他系数的坐标， 表示图像在点处的像素值；

为模板窗口的中心坐标点，对应的像素值为；

是像素值域上的高斯核函数的标准差，用于控制像素值的权值。

一般会将权值称为值域核，或像素值域，但不管是值域核还是空间域核，其大小都在[0 1]之间。

最后，将上述两个模板相乘就得到了双边滤波器的滑窗权值：

因此，双边滤波器的数据公式可以表示如下：

双边滤波器中的一边即为上述的高斯平滑滤波，两个像素物理距离越大则权值越小，反之则权值越大;双边滤波器中另一边由像素值值域决定，两个像素，灰度值相差越小，越不可能是边缘，那么越应该对其进行平滑处理，也就是应该提高其在滤波器中的权值,反之，像素值相差越大则越有可能是边缘，则应该尽力保留。

双边滤波的过程如下图，在突变的边缘上，使用了像素差权重，所以很好的保留了边缘：

5.2 代码与分析

matlab代码：

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%彩色图双边滤波%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function g = bfilt_rgb(f,r,a,b)
% f灰度图；r滤波半径；a全局方差；b局部方差
[x,y] = meshgrid(-r:r);
w1 = exp(-(x.^2+y.^2)/(2*a^2));
f = tofloat(f);%f = im2double(f);
 
h = waitbar(0,'Applying bilateral filter...');
set(h,'Name','Bilateral Filter Progress');
 
fr = f(:,:,1);
fg = f(:,:,2);
fb = f(:,:,3);
[m,n] = size(fr);
**自我介绍一下，小编13年上海交大毕业，曾经在小公司待过，也去过华为、OPPO等大厂，18年进入阿里一直到现在。**

**深知大多数嵌入式工程师，想要提升技能，往往是自己摸索成长或者是报班学习，但对于培训机构动则几千的学费，着实压力不小。自己不成体系的自学效果低效又漫长，而且极易碰到天花板技术停滞不前！**

**因此收集整理了一份《2024年嵌入式&物联网开发全套学习资料》，初衷也很简单，就是希望能够帮助到想自学提升又不知道该从何学起的朋友，同时减轻大家的负担。**

![img](https://img-blog.csdnimg.cn/img_convert/2853d19cc6f7d4de0d0ea571cdd52bcd.png)

![img](https://img-blog.csdnimg.cn/img_convert/26299cdd8450226449a877f480849237.jpeg)

![img](https://img-blog.csdnimg.cn/img_convert/6eafcf4a0cb038d5ff2f1b1db9f57a42.png)

 **既有适合小白学习的零基础资料，也有适合3年以上经验的小伙伴深入学习提升的进阶课程，基本涵盖了95%以上嵌入式&物联网开发知识点，真正体系化！**

![img](https://img-blog.csdnimg.cn/img_convert/a6935ca6fa8a6e732ad9baf0b9769291.png)

![img](https://img-blog.csdnimg.cn/img_convert/ee28c93ba40130cc438189b3d5ac6744.png)

 

**由于文件比较大，这里只是将部分目录大纲截图出来，每个节点里面都包含大厂面经、学习笔记、源码讲义、实战项目、讲解视频，并且后续会持续更新**

**如果你觉得这些内容对你有帮助，可以+V：Vip1104z获取！！！ （备注：嵌入式）**

<img src="https://img-community.csdnimg.cn/images/73bb5de17851459088c6af944156ee24.jpg" alt="img" style="zoom: 67%;" />



# 最后

**资料整理不易，觉得有帮助的朋友可以帮忙点赞分享支持一下小编~**

**你的支持，我的动力；祝各位前程似锦，offer不断，步步高升！！！**

伴深入学习提升的进阶课程，基本涵盖了95%以上嵌入式&物联网开发知识点，真正体系化！**

[外链图片转存中...(img-IS7BFwaR-1712382846745)]

[外链图片转存中...(img-6xHzmQct-1712382846745)]

 

**由于文件比较大，这里只是将部分目录大纲截图出来，每个节点里面都包含大厂面经、学习笔记、源码讲义、实战项目、讲解视频，并且后续会持续更新**

**如果你觉得这些内容对你有帮助，可以+V：Vip1104z获取！！！ （备注：嵌入式）**

<img src="https://img-community.csdnimg.cn/images/73bb5de17851459088c6af944156ee24.jpg" alt="img" style="zoom: 67%;" />



# 最后

**资料整理不易，觉得有帮助的朋友可以帮忙点赞分享支持一下小编~**

**你的支持，我的动力；祝各位前程似锦，offer不断，步步高升！！！**

**[更多资料点击此处获qu！！](https://bbs.csdn.net/topics/618376385)**
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71