前言

卷积：函数空间域的卷积的傅里叶变换是函数傅里叶变换的乘积。对应地，频率域的卷积与空间域的乘积存在对应关系。

给定频率域滤波器，可对其进行傅里叶逆变换得到对应的空域滤波器；滤波在频域更为直观，但空域适合使用更小的滤波模板以提高滤波速度。因为相同尺寸下，频域滤波器效率高于空域滤波器，故空域滤波需要一个更小尺寸的模板近似得到需要的滤波结果。

一、空域滤波

将模板在图像中逐像素移动，将卷积核的每个元素分别和图像矩阵对应位置元素相乘并将结果累加，累加和作为

模板中心对应像素点的卷积结果。通俗的讲，卷积就是对整幅图像进行加权平均的过程，每一个像素点的值，都

由其本身和邻域内的其他像素值经过加权平均后得到。

二、频域滤波

频率域是由傅里叶变换和频率变量 (u,v)定义的空间，频域均匀滤波处理过程：先对图像进行傅里叶变换，转换至频率域，在频域使用滤波函数进行滤波，最后将结果反变换至空间域。即：
（1）计算原始图像f(x,y)的DFT，得到F(u,v)。
（2）创建均值滤波器，并使用fft2计算其DFT得到H(u,v)
（3）计算滤波器函数H(u,v)与F(u,v)的乘积G(u,v)。
（4）计算第（3）步计算结果的傅里叶反变换g(x,y)。
（5）取g(x,y)的实部作为最终滤波后的结果图像。
按照该步骤，在MATLAB中很容易编程实现频域滤波。滤波能否取得理想结果的关键取决于频域滤波函数H（u，v），常常称之为滤波器，或滤波器传递函数。因为它在滤波中抑制或滤除了频谱中某些频率的分量，而保留其他一些频率不受影响。

注意：空域的滤波器x转变的频域时用的函数是fft2(x,m,n),m*n是需要处理的图片大小。

三，matlab代码

t=imread('E:\matlab DMP\lenaG.bmp');
[m,n,z]=size(t);y1=0+10*randn(m,n);%二维高斯分布矩阵 0是均值 20是标准差
y2=0+20*randn(m,n);%二维高斯分布矩阵 0是均值 20是标准差
y1=uint8(y1);
y2=uint8(y2);
%加上噪声
t1=t+y1;
t2=t+y2;
figure('name','添加高斯噪声');;
subplot(2,2,1),imshow(t),title('原图');
subplot(2,2,3),imshow(t1),title('加入均值为0，标准差为10的高斯噪声后');
subplot(2,2,4),imshow(t2),title('加入均值为0，标准差为20的高斯噪声后');K1=conv2(fspecial('average',5),t2 )/255;
K2=conv2(fspecial('average',10),t2 )/255 ;
K3=conv2(fspecial('average',20),t2 )/255;figure('name','空域均值滤波');subplot(2,2,1),imshow(t2),title('原图');
subplot(2,2,2),imshow(K1),title('均值滤波,sizes  5x5');
subplot(2,2,3),imshow(K2),title('均值滤波,sizes  10x10');
subplot(2,2,4),imshow(K3),title('均值滤波,sizes  20x20');Fc = fft2(t2);h1 = fspecial('average',5);
h2 = fspecial('average',10);
h3 = fspecial('average',20);H1 = fft2(h1,512,512);
H2 = fft2(h2,512,512);
H3 = fft2(h3,512,512);G1 = H1.*Fc;
G2 = H2.*Fc;
G3 = H3.*Fc;g1 = ifft2(G1);
g1 = im2uint8(mat2gray(abs(g1)));g2 = ifft2(G2);
g2 = im2uint8(mat2gray(abs(g2)));g3 = ifft2(G3);
g3 = im2uint8(mat2gray(abs(g3)));figure('name','频域均值滤波');
subplot(2,2,1),imshow(t2),title('高斯噪声');
subplot(2,2,2),imshow(g1),title('高斯过滤,size 5*5');
subplot(2,2,3),imshow(g2),title('高斯过滤,size 10*10');
subplot(2,2,4),imshow(g3),title('高斯过滤,size 20*20');Kk1 = medfilt2(t2,[3,3]);
Kk2 = medfilt2(t2,[5,5]);
Kk3 = medfilt2(t2,[7,7]);
figure('name','频域中值滤波');
subplot(2,2,1),imshow(t2),title('高斯噪声');
subplot(2,2,2),imshow(Kk1),title('高斯过滤,size 3*3');
subplot(2,2,3),imshow(Kk2),title('高斯过滤,size 5*5');
subplot(2,2,4),imshow(Kk3),title('高斯过滤,size 7*7');

结果

matlab：lena图像在空域和频域的均匀滤波相关推荐

matlab图像算子变换,拉布拉斯算子锐化图像的空域和频域--matlab实现
本文将用matlab实现拉布拉斯算子提取边缘,锐化图像.分别采用空域和频域两种方法. 空域步骤: 1.读取图像/(标定) 2.生成空间滤波器(fspecial('laplacian',0)) 3.卷 ...
【图像去噪】基于空域和频域两种滤波实现图像去噪含Matlab源码
1 简介本文是基于MATLAB的GUI来对不同的(彩色或灰色)图像进行图像增强的相关处理,主要包括以下四点: 频域变换:傅里叶变换和离散余弦变换及其反变换,可将图像在空域上的特性转化到频域上: 叠加 ...
matlab实现频域率滤波,基于Matlab的图像的频域滤波实现及研究.doc
摘要:图像的频域滤波是图像增强的一种方法.图像增强是图像处理的方法之一,有频率域法和空间域法.频率域法把图像看成一种二维信号,对其进行二维傅里叶变换的信号增强,采用低通滤波法可以去掉图像的噪声:采用高 ...
MATLAB图像的频域低通滤波（灰度图像滤波+彩色图像滤波）
数字图像处理完整MATLAB代码在我的资源可以看到,为方便下载,下面是百度网盘资源: 链接:https://pan.baidu.com/s/17S7PZJwwvb3PFMFVxqEY5w 提取码:H ...
matlab怎么添加条纹噪声,基于频域的图像条纹噪声消除方法
基于频域的图像条纹噪声消除方法徐花艳孙瑞宁 (长安大学陕西西安 710064) 摘要: 主要以图像去噪为主的研究对象,利用图像的频谱图对图像中的条纹噪声的频率分布进行分析,并使用巴特沃斯带阻 ...
详解空域图像增强+频域图像增强
空域图像增强+频域图像增强文章目录空域图像增强+频域图像增强前言空域图像增强一.灰度映射:点处理 1.图像求反(将灰度值翻转) 2.对数变换(动态范围压缩) 3.幂次变换 4.阶梯量化 5. ...
如何保持空域与频域滤波结果的一致性
我们知道,给定一个空(时)间滤波器,既可以在空(时)域直接完成数字信号的滤波,也可以在频域完成.空域滤波的数学运算为卷积/相关,对应频域则为点乘/频域数据的共轭(G*)与滤波器(H)的乘积. 简单起见 ...
matlab空域图像增强,图像处理的MATLAB实现实验一空域图像增强.doc
图像处理的MATLAB实现实验一空域图像增强图像处理的MATLAB实现实验一空域图像增强一.实验目的 (1)掌握基本的空域图像增强方法,观察图像增强的效果,加深理解: (2)了解空域平滑模板的 ...
【老生谈算法】matlab实现图像滤波处理算法源码——图像滤波处理算法
matlab图像滤波处理算法详解 1.文档下载: 本算法已经整理成文档如下,有需要的朋友可以点击进行下载序号文档(点击下载) 本项目文档 [老生谈算法]matlab图像滤波处理.doc 2.算法详 ...
加权平均法融合图像matlab,基于MATLAB的图像融合算法
内容介绍原文档由会员小甜甜发布基于MATLAB的图像融合算法 1.9万字 34页摘要图像融合能够将不同类型传感器获取的同一对象的图像数据进行空间配准.并且采用一定的算法将各图像数据所含的信 ...

matlab：lena图像在空域和频域的均匀滤波

文章目录

前言

一、空域滤波

二、频域滤波

三，matlab代码

结果

matlab：lena图像在空域和频域的均匀滤波相关推荐

最新文章

热门文章