clc
clear all
close all
[filename, pathname] = uigetfile('*.bmp;*.pgm;*.tif;*.jpg;*.*', 'Load input image');%读取文件
filename = [pathname filename];
I=imread(filename);
[height,width]=size(I);
width=width/3;
figure(10)
imshow(I)
I1=I(:,:,1);
I2=I(:,:,2);
I3=I(:,:,3);
J1=fft2(I1);
J1_1=fftshift(J1);
JM1=log(abs(J1_1));
JB1=mat2gray(JM1);
%figure(4)
%imshow(JB1)
J2=fft2(I2);
J2_1=fftshift(J2);
JM2=log(abs(J2_1));
JB2=mat2gray(JM2);
%figure(5)
%imshow(JB2)
J3=fft2(I3);
J3_1=fftshift(J3);
JM3=log(abs(J3_1));
JB3=mat2gray(JM3);
%figure(6)
%imshow(JB3)
y=-1:2/(height-1):1;
x=-1:2/(width-1):1;
[X,Y]=meshgrid(x,y);
Z=1-(X.^2+Y.^2)/2;
%figure(333)
%mesh(X,Y,Z)
JBD1=abs(ifft2((J1).*Z));%这里有没有加ffshift分别对应低频和高频
JBDF1=JBD1/max(max(JBD1));
%figure(7)
%imshow(JBDF1)
JBD2=abs(ifft2((J2).*Z));
JBDF2=JBD2/max(max(JBD2));
%figure(8)
%imshow(JBDF2)
JBD3=abs(ifft2((J3).*Z));
JBDF3=JBD3/max(max(JBD3));
%figure(9)
%imshow(JBDF3)
Final=zeros(height,width,3);
for i=1:height
for j=1:width
Final(i,j,1)=JBDF1(i,j);
end
end
for i=1:height
for j=1:width
Final(i,j,2)=JBDF2(i,j);
end
end
for i=1:height
for j=1:width
Final(i,j,3)=JBDF3(i,j);
end
end
%Final=abs(Final)/max(max(max(abs(Final))));
figure(10)
imshow(Final)
