图像处理实验之频域低通，高通，带通，带阻滤波器，同态滤波器

clear all；clc
f=imread(‘D:\桌面\图像处理实验\图像\fingerprint.tif’);subplot(3,2,1),figure(1),imshow(f); title(‘原始图’); F=fft2(f); %对图像进行傅里叶变换subplot(3,2,2);imshow(F),title(‘傅里叶图’);Fc=fftshift(F);[M,N]=size(Fc);%分别返回Fc的行数到M中，列数到N中D0=500; %初始化D0for i=1:M for j=1:N d=sqrt((i-M)^2+(j-N)2); %点（i,j）到傅里叶变换中心的距离if d<=D0 %点（i,j）在通带内的情况H(i,j)=1; %通带变换函数elseH(i,j)=0; %阻带变换函数 end end endGc=H.*Fc; %理想低通滤波器滤波后的频域表示G=ifftshift(Gc); %去中心化g=ifft2(G); %反傅里叶变换subplot(3,2,3);imshow(g),title(‘低通滤波’);D0=50; %初始化D0for i=1:Mfor j=1:Nd=sqrt((i-M)^2+(j-N)2); %点（i,j）到傅里叶变换中心的距离if d<=D0 %点（i,j）在通带内的情况H(i,j)=0; %通带变换函数elseH(i,j)=1; %阻带变换函数endendendGc=H.*Fc; %理想低通滤波器滤波后的频域表示G=ifftshift(Gc); %去中心化g=ifft2(G); %反傅里叶变换subplot(3,2,4);imshow(g),title(‘高通滤波’);D0=50; D1=500;%初始化D0for i=1:Mfor j=1:Nd=sqrt((i-M