Matlab频域高/低通滤波

建议参考书籍：数字图像处理_第三版冈萨雷斯

写在前面：
对于给定的低通滤波器的函数表达式，可以得到高通滤波器的函数表达式：

理想高/低通滤波器

理想高通

一个二维理想高通滤波器（IHPF）定义为：
其中，D0为截至频率。D的表达式为：（后同）
是频域率中点（u, v）与频率矩形中心的距离，其中P，Q为频率矩形的长宽。

具体实现代码：

clc;
clear;
close all;img=imread('test.jpg');
subplot(221),imshow(img),title('原图');
F=fft2(img);
[M,N]=size(F);
u=0:(M-1);v=0:(N-1);
idx = find(u > M/2);u(idx) = u(idx) - M;
idy = find(v > N/2);v(idy) = v(idy) - N;
[U,V]=meshgrid(u,v);
D=sqrt(U.^2+V.^2);
D0=M;% 可以自己调节截止频率
H=ones(M,N);
for i=1:Mfor j=1:Nif D(i,j)^2<=D0H(i,j)=0;endend
endshiftH=fftshift(H);
subplot(222),surf(shiftH(1:10:600,1:10:600)),title('理想高通滤波器透视图');
subplot(223),imshow(fftshift(H),[]),title('理想高通滤波器图像表示');
g=real(ifft2(F.*H));
subplot(224),imshow(g,[]),title('高通滤波');

代码效果：

理想低通

在以原点为圆心、以D0为半径的圆内，无衰减地通过所有频率，而在圆外‘切断‘所有频率的二维滤波器，称为理想低通滤波器（ILPF）；它由以下函数确定：

具体代码实现：

clc;
clear;
close all;img=imread('test.jpg');
subplot(221),imshow(img),title('原图');
F=fft2(img);
[M,N]=size(F);
u=0:(M-1);v=0:(N-1);
idx = find(u > M/2);u(idx) = u(idx) - M;
idy = find(v > N/2);v(idy) = v(idy) - N;
[U,V]=meshgrid(u,v);
D=sqrt(U.^2+V.^2);
D0=M;% 可以自己调节截止频率
H=zeros(M,N);
for i=1:Mfor j=1:Nif D(i,j)^2<=D0H(i,j)=1;endend
endshiftH=fftshift(H);
subplot(222),surf(shiftH(1:10:600,1:10:600)),title('理想低通滤波器透视图');
subplot(223),imshow(fftshift(H),[]),title('理想低通滤波器图像表示');
g=real(ifft2(F.*H));
subplot(224),imshow(g,[]),title('低通滤波');

代码效果：

高斯高/低通滤波

高斯高通

高斯高通滤波器函数表达式：

实现代码：

clc;
clear;
close all;img=imread('test.jpg');
F=fft2(img);
[M,N]=size(F);
subplot(221),imshow(img),title('原图');
u=0:(M-1);v=0:(N-1);
idx = find(u > M/2);u(idx) = u(idx) - M;
idy = find(v > N/2);v(idy) = v(idy) - N;
[U,V]=meshgrid(u,v);
D=sqrt(U.^2+V.^2);
D0=0.1*M;
H=1-exp(-(D.^2)/(2*(D0)^2));
shiftH=fftshift(H);
subplot(222),mesh(shiftH(1:10:600,1:10:600)),title('高斯高通滤波透视图');
subplot(223),imshow(fftshift(H),[]),title('高斯高通滤波器图像表示');
g=real(ifft2(F.*H));
subplot(224),imshow(g,[]),title('高斯高通滤波');

代码效果：

高斯低通

高斯低通滤波器函数表达式：

具体代码实现：

clc;
clear;
close all;img=imread('test.jpg');
F=fft2(img);
[M,N]=size(F);
subplot(221),imshow(img),title('原图');
u=0:(M-1);v=0:(N-1);
idx = find(u > M/2);u(idx) = u(idx) - M;
idy = find(v > N/2);v(idy) = v(idy) - N;
[U,V]=meshgrid(u,v);
D=sqrt(U.^2+V.^2);
D0=0.1*M;
H=exp(-(D.^2)/(2*(D0)^2));
shiftH=fftshift(H);
subplot(222),mesh(shiftH(1:10:600,1:10:600)),title('高斯低通滤波透视图');
subplot(223),imshow(fftshift(H),[]),title('高斯低通滤波器图像表示');
g=real(ifft2(F.*H));
subplot(224),imshow(g,[]),title('高斯低通滤波');

代码效果：

巴特沃斯高/低通滤波器

巴特沃斯高通

n阶巴斯沃特高通滤波器函数表达式为：

实现代码：

clc;
clear;
close all;img=imread('test.jpg');
F=fft2(img);
[M,N]=size(F);
n=2;
subplot(221),imshow(img),title('原图');
u=0:(M-1);v=0:(N-1);
idx = find(u > M/2);u(idx) = u(idx) - M;
idy = find(v > N/2);v(idy) = v(idy) - N;
[U,V]=meshgrid(u,v);
D=sqrt(U.^2+V.^2);
D0=0.05*M;
H=1./(1+(D0./D).^(2*n));
shiftH=fftshift(H);
subplot(222),mesh(shiftH(1:10:600,1:10:600)),title('巴特沃斯2阶高通滤波透视图');
subplot(223),imshow(fftshift(H),[]),title('巴特沃斯2阶高通滤波器图像表示');
g=real(ifft2(F.*H));
subplot(224),imshow(g,[]),title('巴特沃斯2阶高通滤波');

代码效果：

巴特沃斯低通

截止频率位于距原点D0处的n阶巴特沃斯低通滤波器（BLPF）的传递函数为：

具体实现代码：

clc;
clear;
close all;img=imread('test.jpg');
F=fft2(img);
[M,N]=size(F);
subplot(221),imshow(img),title('原图');
u=0:(M-1);v=0:(N-1);
idx = find(u > M/2);u(idx) = u(idx) - M;
idy = find(v > N/2);v(idy) = v(idy) - N;
[U,V]=meshgrid(u,v);
D=sqrt(U.^2+V.^2);
D0=0.1*M;
H=1-exp(-(D.^2)/(2*(D0)^2));
shiftH=fftshift(H);
subplot(222),mesh(shiftH(1:10:600,1:10:600)),title('高斯高通滤波透视图');
subplot(223),imshow(fftshift(H),[]),title('高斯高通滤波器图像表示');
g=real(ifft2(F.*H));
subplot(224),imshow(g,[]),title('高斯高通滤波');

代码效果展示：

Matlab频域高/低通滤波相关推荐

基于matlab的音频信号低通滤波实验
基于matlab的音频信号低通滤波实验一.实验任务录制语音信号,查看其时域频谱图. 加入高频噪声信号,并查看其时域频谱图. 设计低通滤波器,将带有高频噪声的语音信号经过低通滤波器,将噪声滤除,查看 ...
matlab 对图像进行低通滤波,matlab图像处理实现低通滤波
matlab图像处理实现低通滤波 matlab图像处理实现低通滤波 clc; clear all; img=imread('lena.jpg'); img_noise = imnoise(img, ' ...
VTK：频域处理——低通滤波（理想+巴特沃斯）
1.低通滤波器低通滤波是将频域图像中的高频部分滤除而通过低频部分.图像的边缘和噪声对应于频域图像中的高频部分,而低频滤波的作用即是减弱这部分作用,从而达到图像平滑去噪的目的. 2.理想低通滤波器最 ...
MATLAB图像处理--高斯低通滤波、高斯高通滤波（代码及示例）
本文目录常用的滤波器高斯低通滤波器 matlab代码高斯高通滤波器 Butterworth 低通滤波器常用的滤波器常用的滤波器有以下几种: 频域滤波器频域平滑滤波器理想低通滤波器巴特沃 ...
MATLAB巴特沃斯低通滤波图像
clc,clear,close all % 清理命令区.清理工作区.关闭显示图形 warning off % 消除警告 feature jit off % 加速代码运行 D0 = 20; % 阻止的频 ...
MATLAB高通滤波与低通滤波GUI实现代码
1.需搭建的GUI界面两个axes,一个用于导入原始图像,一个用于输出处理后的头像,并且记好他们的Tag,在代码里会用到 2.导入图片功能代码实现给一个String为导入图片的button添加ca ...
Matlab对灰度图像的频域进行高通滤波和低通滤波
1. 要求对灰度图像进行离散傅里叶变换(Discrete Fourier Transfom, DFT)变换,在频域上分别使用理想的高通和低通滤波器进行滤波,显示滤波后的频域图像,以及逆离散傅里叶变换 ...
数字图像处理-频域滤波-高通/低通滤波
频域滤波频域滤波是在频率域对图像做处理的一种方法.步骤如下: 滤波器大小和频谱大小相同,相乘即可得到新的频谱. 滤波后结果显示,低通滤波去掉了高频信息,即细节信息,留下的低频信息代表了概貌.常用的例 ...
图像理想低通滤波_图像处理之滤波（下）
[toc]目录一.常规滤波低通高通带通带阻二.非局部均值滤波三.维纳滤波四.卡尔曼滤波前言所谓滤波,其实就是从混合在一起的诸多信号中提取出所需要的信号. 信号的分类: 确定型信号, ...

Matlab频域高/低通滤波

理想高/低通滤波器

理想高通

理想低通

高斯高/低通滤波

高斯高通

高斯低通

巴特沃斯高/低通滤波器

巴特沃斯高通

巴特沃斯低通

Matlab频域高/低通滤波相关推荐

最新文章

热门文章