高通滤波:
function [ one ] = Untitled( )
img = imread('D:\图片\文件图片\Fig1(moon).tif');
[f1,f2]=freqspace(size(img),'meshgrid');
d = 0.3;
r = f1.^2+f2.^2;
for i =1:size(img,1)
for j = 1:size(img,2)
t = r(i,j)/(d*d);
hd(i,j) = 1-exp(-t);
end
end
y = fft2(double(img));
y = fftshift(y);
ya = y.*hd;
ya = fftshift(ya);
yb = real(ifft2(ya));
subplot(2,2,1),imshow(uint8(img));
title('原图图像');
subplot(2,2,2),imshow(uint8(yb));
title('高通滤波图像');
end
贷通滤波:
img = imread('D:\图片\文件图片\Fig1(moon).tif');
[f1,f2]=freqspace(size(img),'meshgrid');
d = 0.3;
r = f1.^2+f2.^2;
for i =1:size(img,1)
for j = 1:size(img,2)
t = r(i,j)/(d*d);