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

媒介：本文为手把手讲授 Matlab 平台下的各种图像滤波算法的教程，将编程代码与图像滤波知识相接洽，以实战为例！博客中图像滤波算法包含：均值滤波、中值滤波、高斯滤波、双边滤波、引导滤波。图像滤波算法是盘算机视觉领域CV必修课，被广泛运用于各行各业，尤其是科研领域！希望本篇博客能给读者朋侪的工程项目或科研生活给予些许帮助。（篇末代码开源！）
图像滤波算法总图：

一、图像滤波知识

1.1 根本概念

图像滤波是一种常见的图像处理技术，旨在抑制图像中的噪声，并在尽量保留图像细节特征的条件下改善图像质量。该技术是图像预处理中非常紧张的步调，将直接影响后续的图像数据挖掘与分析性能（例如：目标检测，目标分割与去雾去雨算法等）！实在如今大火的卷积神经网络 CNN 实在也是滤波的一种，都是用卷积核去提取图像的特征模式。不外传统的滤波，使用的卷积核是固定的参数，是由履历丰富的人去手动计划的，也称为手工特征。而卷积神经网络的卷积核参数是未知的，是根据差别的任务由数据驱动去学习得到的参数，更能顺应于差别的任务。作者补充：卷积是图像算法中很紧张的知识点，读者朋侪肯定须要熟练掌握!

1.2 图像噪声

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

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

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

matlab代码：

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

复制代码

1.2.2 高斯噪声

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

matlab代码：

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

复制代码

1.2.3 泊松噪声

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

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

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

复制代码

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代码：

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

复制代码

均值滤波图：

   优点：均值滤波算法简单，盘算速度快，算法容易实现。
  缺点：由于均值滤波是通过对邻域像素的平均值来替代中央像素的值，因此容易导致图像细节的丢失；均值滤波会平滑图像中的高频部分，包括边缘和细节，使得图像边缘变得含糊。
  三、中值滤波

3.1 算法概述

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

3.2 代码与分析

matlab代码：

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

复制代码

中值滤波图：

   优点：中值滤波可以有用地去除图像中的椒盐噪声、加性高斯噪声等噪声类型，使图像更清晰；相比其他滤波方法，中值滤波能够在去噪的同时较好地保留图像的边缘信息；中值滤波的盘算过程也是简单容易的。
  缺点：对于图像中含有较多点、线、尖角细节的情况，中值滤波大概不适用，因为它大概会改变这些细节的像素值；中值滤波是一种非线性滤波方法，它通过将像素邻域内的所有像素值排序来选择中值，这大概导致图像的不连续性。
  四、高斯滤波

4.1 算法概述

高斯滤波(Gaussian filter) 包含许多种，包括低通、带通和高通等，我们通常图像上说的高斯滤波，指的是高斯含糊(Gaussian Blur) ，是一种高斯低通滤波 ，其过滤调图像高频身分（图像细节部分），保留图像低频身分（图像平滑地区），以是对图像进行 ‘高斯含糊’ 后，图像会变得含糊。高斯含糊对于抑制高斯噪声 (服从正态分布的噪声) 非常有用。
   ★高斯滤波和均值滤波一样，都是使用一个滤波窗口和图像进行卷积求解，差别之处在于：
  

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

    高斯函数分布：
  公式：
  

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

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

matlab代码：

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

复制代码

高斯滤波图：

   优点：高斯滤波能够平滑图像，去除噪声，并保留图像的整体结构；高斯滤波可以通过调整高斯核的标准差来控制，灵活度高；高斯滤波在离散化窗口滑窗卷积实现时盘算速度快，适用于实时处理。
  缺点：高斯滤波器会平滑图像中的高频部分，包括边缘和细节，大概导致图像边缘含糊；高斯滤波须要进行卷积运算，盘算复杂度相对较高，大概会消耗大量的盘算资源和时间。
  五、双边滤波

5.1 算法概述

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

 。

   

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

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

是空间域上的高斯核函数的标准差，用于控制像素位置的权值。
  使用公式天生的滤波器模板和高斯滤波器使用的模板是没有区别的。权值称为定义域核，也称为空间系数、空间域。显示由的盘算公式可知，它是盘算邻近点q到中央点p邻近水平，因此定义域核是用于衡量空间邻近的水平。
  （2）值域核 
由像素值的差值决定的模板权值

 。

   其中，

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

 表现图像在点

处的像素值；
  

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

；
  

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

称为值域核，或像素值域，但不管是值域核

还是空间域核

，其巨细都在[0 1]之间。
  最后，将上述两个模板相乘就得到了双边滤波器的滑窗权值：

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

   双边滤波器中的一边即为上述的高斯平滑滤波，两个像素物理间隔越大则权值越小，反之则权值越大;双边滤波器中另一边由像素值值域决定，两个像素，灰度值相差越小，越不大概是边缘，那么越应该对其进行平滑处理，也就是应该进步其在滤波器中的权值,反之，像素值相差越大则越有大概是边缘，则应该尽力保留。
  双边滤波的过程如下图，在突变的边缘上，使用了像素差权重，以是很好的保留了边缘：

5.2 代码与分析

matlab代码：

1. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%彩色图双边滤波%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
2. function g = bfilt_rgb(f,r,a,b)
3. % f灰度图；r滤波半径；a全局方差；b局部方差
4. [x,y] = meshgrid(-r:r);
5. w1 = exp(-(x.^2+y.^2)/(2*a^2));
6. f = tofloat(f);%f = im2double(f);
8. h = waitbar(0,'Applying bilateral filter...');
9. set(h,'Name','Bilateral Filter Progress');
11. fr = f(:,:,1);
12. fg = f(:,:,2);
13. fb = f(:,:,3);
14. [m,n] = size(fr);
15. fr_temp = padarray(fr,[r r],'symmetric');
16. fg_temp = padarray(fg,[r r],'symmetric');
17. fb_temp = padarray(fb,[r r],'symmetric');
18. [gr,gg,gb] = deal(zeros(size(fr)));
21. for i = r+1:m+r
22.     for j = r+1:n+r
23.         temp1 = fr_temp(i-r:i+r,j-r:j+r);
24.         temp2 = fg_temp(i-r:i+r,j-r:j+r);
25.         temp3 = fb_temp(i-r:i+r,j-r:j+r);
26.         dr = temp1 - fr_temp(i,j);
27.         dg = temp2 - fg_temp(i,j);
28.         db = temp3 - fb_temp(i,j);
29.         w2 = exp(-(dr.^2+dg.^2+db.^2)/(2*b^2));
30.         w = w1.*w2;
31.         gr(i-r,j-r) = sum(sum(temp1.*w))/sum(w(:));
32.         gg(i-r,j-r) = sum(sum(temp2.*w))/sum(w(:));
33.         gb(i-r,j-r) = sum(sum(temp3.*w))/sum(w(:));
34.     end
35.     waitbar((i-r)/n);
36. end
37. g = cat(3,gr,gg,gb);
38. % g = revertclass(g);
39. end
41. %%%%%%%%%%%%%%%%%%%%%%%%%%%%可以用到的函数%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
42. function [out,revertclass]=tofloat(in)
43. identity=@(x) x;
44. tosingle=@im2single;
46. table={'uint8',tosingle,@im2uint8
47. 'uint16',tosingle,@im2uint16
48. 'int16',tosingle,@im2int16
49. 'logical',tosingle,@logical
50. 'double',identity,identity
51. 'single',identity,identity};
53. classIndex=find(strcmp(class(in),table(:,1)));
55. if isempty(classIndex)
56.     error('unsupported input immage class.');
57. end
59. out=table{classIndex,2}(in);
60. revertclass=table{classIndex,3};
61. end

复制代码

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

复制代码

双边滤波图：

   优点：双边滤波能够同时考虑空间邻域和灰度值差异，因此在平滑图像的同时较好地保留边缘信息；双边滤波能够有用地去除图像中的噪声，特别是椒盐噪声；双边滤波的参数可调，可以根据差别的应用场景调整参数以获得最佳效果。
  缺点：双边滤波须要进行卷积运算和排序操纵，盘算复杂度，会消耗非常多的盘算本钱；双边滤波的参数选择对滤波效果有很大影响，须要仔细调整以获得最佳效果。
  六、引导滤波

6.1 算法概述

引导滤波的过程如下图所示，其在保边性能与盘算复杂度上都到达了极优： 

6.2 代码与分析

matlab代码：

1. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
2. function q = guidedfilter(I, p, r, eps)
3. %   GUIDEDFILTER   O(1) time implementation of guided filter.
4. %
5. %   - guidance image: I (should be a gray-scale/single channel image)
6. %   - filtering input image: p (should be a gray-scale/single channel image)
7. %   - local window radius: r
8. %   - regularization parameter: eps
10. [hei, wid] = size(I);
11. N = boxfilter(ones(hei, wid), r); % the size of each local patch; N=(2r+1)^2 except for boundary pixels.
13. mean_I = boxfilter(I, r) ./ N;
14. mean_p = boxfilter(p, r) ./ N;
15. mean_Ip = boxfilter(I.*p, r) ./ N;
16. cov_Ip = mean_Ip - mean_I .* mean_p; % this is the covariance of (I, p) in each local patch.
18. mean_II = boxfilter(I.*I, r) ./ N;
19. var_I = mean_II - mean_I .* mean_I;
21. a = cov_Ip ./ (var_I + eps); % Eqn. (5) in the paper;
22. b = mean_p - a .* mean_I; % Eqn. (6) in the paper;
24. mean_a = boxfilter(a, r) ./ N;
25. mean_b = boxfilter(b, r) ./ N;
27. q = mean_a .* I + mean_b; % Eqn. (8) in the paper;
28. end
29. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
30. function imDst = boxfilter(imSrc, r)
32. %   BOXFILTER   O(1) time box filtering using cumulative sum
33. %
34. %   - Definition imDst(x, y)=sum(sum(imSrc(x-r:x+r,y-r:y+r)));
35. %   - Running time independent of r;
36. %   - Equivalent to the function: colfilt(imSrc, [2*r+1, 2*r+1], 'sliding', @sum);
37. %   - But much faster.
39. [hei, wid] = size(imSrc);
40. imDst = zeros(size(imSrc));
42. %cumulative sum over Y axis
43. imCum = cumsum(imSrc, 1);
44. %difference over Y axis
45. imDst(1:r+1, :) = imCum(1+r:2*r+1, :);
46. imDst(r+2:hei-r, :) = imCum(2*r+2:hei, :) - imCum(1:hei-2*r-1, :);
47. imDst(hei-r+1:hei, :) = repmat(imCum(hei, :), [r, 1]) - imCum(hei-2*r:hei-r-1, :);
49. %cumulative sum over X axis
50. imCum = cumsum(imDst, 2);
51. %difference over Y axis
52. imDst(:, 1:r+1) = imCum(:, 1+r:2*r+1);
53. imDst(:, r+2:wid-r) = imCum(:, 2*r+2:wid) - imCum(:, 1:wid-2*r-1);
54. imDst(:, wid-r+1:wid) = repmat(imCum(:, wid), [1, r]) - imCum(:, wid-2*r:wid-r-1);
55. end

复制代码

1. %%灰度图引导滤波
2. close all;
3. clear all;
4. clc;
5. I = imread('Noise Image.jpg');  % 读取图像
6. if length(size(I))>2
7.     I=rgb2gray(I);
8. end
9. I = double(I) / 255;
10. p = I;
11. r = 9;
12. eps = 0.4^2;      %0.021^2  0.04^2
13. J = guidedfilter(I,p,r,eps);           %引导滤波
14. % 显示原始图像和双边滤波后的图像  
15. subplot(1,2,1), imshow(I);  
16. subplot(1,2,2), imshow(J);

复制代码

引导滤波图：

   优点：引导滤波能够使用引导图像中的边缘信息，使输出图像的边缘与引导图像相似，因此在平滑图像的同时较好地保留边缘信息；由于引导滤波使用了引导图像的信息，因此对噪声具有肯定的鲁棒性。
  缺点：引导滤波须要进行卷积运算和迭代优化，大概会消耗很多的盘算资源和时间（比双边滤波消耗时间短）；如果输入图像中存在大范围的噪声，引导滤波的效果大概会受到影响。
  七、算法总结

图像滤波算法是盘算机视觉领域必须掌握的技能，从事算法科研亦或是嵌入式AI行业的工程职员都应该熟练掌握它！图像滤波算法随着时代科学的进步，如今出现了许多新型的图像滤波算法，例如：反盲卷积滤波，侧窗盒式滤波与维纳滤波等。希望读者朋侪在现实的工程项目中，不要盲目的追求各种牛逼的图像滤波算法，实在适合该工程的滤波就是好滤波！
八、代码开源

代码地址： 基于Matlab的各种图像滤波算法代码资源-CSDN文库
如果积分不够的朋侪，点波关注，评论区留下邮箱，作者无偿提供源码和后续问题解答。求求啦关注一波吧 ！！！

免责声明：如果侵犯了您的权益，请联系站长，我们会及时删除侵权内容，谢谢合作！更多信息从访问主页：qidao123.com:ToB企服之家，中国第一个企服评测及商务社交产业平台。