High_Freq.rar_high_滤波程序_高通滤波

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)