频域滤波是一种图像处理技术,可以通过在频域中增加或减弱某些频率分量,从而实现图像去噪、锐化、平滑等功能。常见的频域滤波包括频域低通滤波、频域高通滤波和频域同态滤波。 在使用这些滤波器进行频域处理时,通常需要选择合适的参数、预处理和后处理方法,以达到理想的处理效果。同时,不同的图像场景和需求也需要选择适合的滤波器和处理方法。
1.频域低通滤波
频域低通滤波:低通滤波器将保留图像中低频信息,同时削弱高频信息,从而实现图像模糊、平滑等效果。低通滤波器的频率响应通常为一个圆形区域,圆形区域内表示通过的频率分量,圆形区域外表示被滤掉的频率分量。低通滤波器适合处理噪声较多或需要平滑图像的场景。
以下是实现频域低通滤波中的理想低通滤波:
close all; %关闭当前所有图形窗口
clear all; %清空工作空间变量
clc; %清屏
I=imread('lena.bmp');
I=rgb2gray(I);
I=im2double(I);
[m,n]=size(I);
M=2*m;N=2*n;%滤波器的行数和列数
u=-M/2:(M/2-1);
v=-N/2:(N/2-1);
[U,V]=meshgrid(u,v);%meshgrid(u,v)的作用是分别产生以向量u为行,向量v为裂的两个大小相同的矩阵
D=sqrt(U.^2+V.^2);%设置频率点(U,V)与频域中心的距离为D(U,V)
D0=40; %截止频率
H=double(D<=D0); %理想低通滤波器
J=fftshift(fft2(I,size(H,1),size(H,2)));%通过傅里叶变换将时域图像转换到频域图像,并移动到中心位置
K=J.*H; %滤波处理
L=ifft2(ifftshift(K)); %傅里叶反转换
L=L(1:m,1:n); %设置图像的大小
figure;
subplot(131),imshow(J);
subplot(132),imshow(I);
subplot(133),imshow(L);
文章来源地址https://www.toymoban.com/news/detail-652012.html
2.频域高通滤波
频域高通滤波:高通滤波器则削弱低频信息,同时保留高频信息,从而实现图像增强、锐化等效果。高通滤波器的频率响应通常为一环形区域,环形区域内表示被滤掉的频率分量,环形区域外表示通过的频率分量。高通滤波器适合处理模糊或需要突出细节、边缘等特定信息的场景。
以下是实现频域高通滤波中的理想高通滤波:
close all; %关闭当前所有图形窗口
clear all; %清空工作空间变量
clc; %清屏
I=imread('lena.bmp');
I=rgb2gray(I);
I=im2double(I);
[m,n]=size(I);
M=2*m;N=2*n;%滤波器的行数和列数
u=-M/2:(M/2-1);
v=-N/2:(N/2-1);
[U,V]=meshgrid(u,v);%meshgrid(u,v)的作用是分别产生以向量u为行,向量v为裂的两个大小相同的矩阵
D=sqrt(U.^2+V.^2);%设置频率点(U,V)与频域中心的距离为D(U,V)
D0=40; %截止频率
H=double(D>D0); %理想高通滤波器
J=fftshift(fft2(I,size(H,1),size(H,2)));%通过傅里叶变换将时域图像转换到频域图像,并移动到中心位置
K=J.*H; %滤波处理
L=ifft2(ifftshift(K)); %傅里叶反转换
L=L(1:m,1:n); %设置图像的大小
figure;
subplot(131),imshow(J);
subplot(132),imshow(I);
subplot(133),imshow(L);
3.频域同态滤波
频域同态滤波:同态滤波器可以同时处理低频和高频信息,通过在频域中对图像进行对数变换,将乘积域转化为加法域,从而实现亮度和反差的调整。同态滤波器常用于处理光照不均和阴影等问题,同时也有一定的去噪和增强效果。
close all; %关闭当前所有图形窗口
clear all; %清空工作空间变量
clc; %清屏
I=imread('lena.bmp');
I=rgb2gray(I);
J=log(im2double(I)+1); %取对数
K=fft2(J); %傅里叶变换
n=5; %同态滤波参数
D0=0.1*pi;
rh=0.8;
rl=0.3;
[row,column]=size(J);
for i=1:row
for j=1:column
D1=sqrt(i^2+j^2);
H(i,j)=rl+(rh/(1+(D0/D1)^(2*n)));
end
end
L=K.*H; %滤波
M=ifft2(L); %傅里叶反变换
N=exp(M)-1; %取指数
figure;
subplot(121),imshow(I);
subplot(122),imshow(real(N));
文章来源:https://www.toymoban.com/news/detail-652012.html
到了这里,关于matlab实现图像频域处理(低通滤波、高通滤波、同态滤波)的文章就介绍完了。如果您还想了解更多内容,请在右上角搜索TOY模板网以前的文章或继续浏览下面的相关文章,希望大家以后多多支持TOY模板网!