傅里叶变化频谱图及频域滤波

close all; clc; pic = imread('Lena.jpg'); %-----------------------理想高通------------------------------ image20 = ideal_highpass(pic,10); image30 = ideal_highpass(pic,20); image40 = ideal_highpass(pic,6); image50 = ideal_highpass(pic,10); image60 = ideal_highpass(pic,16); % 显示图像 subplot(2,3,1), imshow(pic), title('原图像(Id)'); subplot(2,3,2), imshow(image20), title('截止频率D0 = 10'); subplot(2,3,3), imshow(image30), title('截止频率D0 = 30'); subplot(2,3,4), imshow(image40), title('截止频率D0 = 60'); subplot(2,3,5), imshow(image50), title('截止频率D0 = 100'); subplot(2,3,6), imshow(image60), title('截止频率D0 = 160'); %-----------------------巴特沃斯高通-------------------------- image21 = Bhighfilter(pic, 10, 2);%截止频率10,阶数2 image31 = Bhighfilter(pic, 20, 2);%截止频率30 image41 = Bhighfilter(pic, 60, 2);%截止频率30 image51 = Bhighfilter(pic, 100, 2);%截止频率30 image61 = Bhighfilter(pic, 160, 2);%截止频率30 % 显示图像 figure subplot(2,3,1), imshow(pic), title('原图像(Bw)'); subplot(2,3,2), imshow(image21), title('截止频率D0 = 10'); subplot(2,3,3), imshow(image31), title('截止频率D0 = 100'); subplot(2,3,4), imshow(image41), title('截止频率D0 = 60'); subplot(2,3,5), imshow(image51), title('截止频率D0 = 160'); subplot(2,3,6), imshow(image61), title('截止频率D0 = 360'); %------------------------高斯高通--------------------------- image22 = gauss_highpass(pic, 10);%截止频率10 image32 = gauss_highpass(pic, 20);%截止频率30 image42 = gauss_highpass(pic, 60);%截止频率30 image52 = gauss_highpass(pic, 100);%截止频率30 image62 = gauss_highpass(pic, 160);%截止频率30 % 显示图像 figure subplot(2,3,1), imshow(pic), title('原图像(Gau)'); subplot(2,3,2), imshow(image22), title('截止频率D0 = 10'); subplot(2,3,3), imshow(image32), title('截止频率D0 = 30'); subplot(2,3,4), imshow(image42), title('截止频率D0 = 60'); subplot(2,3,5), imshow(image52), title('截止频率D0 = 160'); subplot(2,3,6), imshow(image62), title('截止频率D0 = 360');