实验五 频域变换.pptVIP

第1页，共17页，星期日，2025年，2月5日1.傅里叶变换(1)简单的人工二值图像clearall,closeallf=zeros(50,50);f(15:35,23:28)=1;figure(1),imshow(f,notruesize)F=fft2(f,128,128);F1=fftshift(F);figure(2),imshow(log(abs(F1)),[-15]);colormap(gray);colorbarfigure(3),mesh(1:128,1:128,abs(F1));colormap(gray);colorbarF2=fft2(imrotate(f,90),128,128);F3=fftshift(F2);figure(4),imshow(imrotate(f,90),notruesize)figure(5),imshow(log(abs(F3)),[-15]);colormap(gray);colorbarfigure(6),mesh(1:128,1:128,abs(F3));colormap(gray);colorbar第2页，共17页，星期日，2025年，2月5日（2）实际图像傅里叶变换clearall,closeallA=imread(rice.tif);imshow(A);FT=fft2(A);AbsFT=abs(FT);F2=log(AbsFT);figure,imshow(F2,[-110],notruesize);FinvT=ifft2(FT);AbsFinvT=abs(FinvT);F3=log(AbsFinvT);figure,imshow(mat2gray(FinvT));第3页，共17页，星期日，2025年，2月5日(3)频域滤波closeall,clearallbw=imread(text.tif);imshow(bw);C=real(ifft2(fft2(bw).*fft2(fspecial(average),256,256)));figure,imshow(C,[]);第4页，共17页，星期日，2025年，2月5日2.离散余弦变换(1)对cameraman.tif进行变换和逆变换重建clearall,closeallf=imread(cameraman.tif);imshow(f);F=dct2(f);AbsFT=abs(F);figure,imshow(log(AbsFT));FinvT=idct2(F);figure,imshow(mat2gray(FinvT));第5页，共17页，星期日，2025年，2月5日%保留左上角50*50个数据，反变换后看效果F1=F;[m,n]=size(F1);F1(50:m,50:n)=0;AbsFT=abs(F1);figure,imshow(log(AbsFT));FinvT=idct2(F1);figure,imshow(mat2gray(FinvT));第6页，共17页，星期日，2025年，2月5日I1=imread(saturn.tif);figure,imshow(I1)I2=imnoise(I1,saltpepper);figure,imshow(I2)f=double(I2);g=fft2(f);g=fftshift(g);[N1,N2]=size(g);n=2;d0=50;n1=fix(N1/2);n2=fix(N2/2);fori=1:N1forj=1:N2d=sqrt((i-n1)^2+(j-n2)^2);h=1/(1+0.414*(d/d0)^(2*n));%计算Butterworth低通转换函数result(i,j)=h*g(i,j);endendresult=ifftshift(result);X2=ifft2(result);X3=uint8(real(X2));figure,imshow(X3)Butterworth低通滤波器第7页，共17页，星期日，2025年，2月5日I1=imread(blood1.tif);figure(1),imshow(I1)f=