clc;
clear all;
close all;
J=imread('d:\111.PNG');
if size(J, 3)==3
J = rgb2gray(J);
end
subplot(221);imshow(uint8(J));
title('原图像');
J=double(J);
f=fft2(J); %采用傅里叶变换
g=fftshift(f); %数据矩阵平衡
[M,N]=size(f);
n1=floor(M/2);
n2=floor(N/2);
d0=20;
for i=1:M %进行理想高通滤波和理想高通加强滤波
for j=1:N
d=sqrt((i-n1)^2+(j-n2)^2);
if d>=d0
h1=1;
h2=1+0.5;
else
h1=0;
h2=0.5;
end
g1(i,j)=h1*g(i,j);
g2(i,j)=h2*g(i,j);
end
end
g1=ifftshift(g1);
g1=uint8(real(ifft2(g1)));
subplot(222);imshow(g1); %显示理想高通滤波结果
title('理想高通滤波');
g2=ifftshift(g2);
g2=uint8(real(ifft2(g2)));
subplot(223);imshow(g2); %显示理想高通加强滤波结果
title('理想高通加强滤波');
n=2;
d0=20;
评论0