MatLab理想低通濾波及高通濾波實現(xiàn)高頻增強(qiáng)實驗報告

1、實驗作業(yè)5：1、 用理想低通濾波器在頻率域?qū)崿F(xiàn)低通濾波程序代碼：I=imread('d:/3.jpg');I=rgb2gray(I);figure(1),imshow(I);title('原圖像');s=fftshift(fft2(I);figure(2);imshow(abs(s),);title('圖像傅里葉變換所得頻譜');a,b=size(s);a0=round(a/2);b0=round(b/2);d=10;for i=1:a for j=1:b distance=sqrt(i-a0)2+(j-b0)2); if distance<

2、;=d h=1; else h=0; end; s(i,j)=h*s(i,j); end;end;s=uint8(real(ifft2(ifftshift(s);figure(3);imshow(s);title('低通濾波所得圖像'); 濾波d=102、用理想高通濾波器在頻率域?qū)崿F(xiàn)高頻增強(qiáng)I=imread('d:/3.jpg');I=rgb2gray(I);figure(1),imshow(I);title('原圖像');s=fftshift(fft2(I);figure(2);imshow(abs(s),);title('圖像傅里葉變

3、換所得頻譜');figure(3);imshow(log(abs(s),);title('圖像傅里葉變換取對數(shù)所得頻譜');a,b=size(s);a0=round(a/2);b0=round(b/2);d=10;p=0.2;q=0.5;for i=1:a for j=1:b distance=sqrt(i-a0)2+(j-b0)2); if distance<=d h=0; else h=1; end; s(i,j)=(p+q*h)*s(i,j); end;end;s=uint8(real(ifft2(ifftshift(s);figure(4);imshow(

4、s);title('高通濾波所得圖像');figure(5);imshow(s+I);title('高通濾波所得高頻增強(qiáng)圖像'); 心得體會：這次按照低通濾波器和高通濾波器的定義，按照低通濾波的過程，一步一步寫，先是進(jìn)行傅里葉變換，再對其頻譜進(jìn)行平移，使其中心位于中心，再對此時的頻譜進(jìn)行圓形濾波，剛開始糾結(jié)于公式，要怎么想出一個H（s）的濾波器的表達(dá)式，然后再進(jìn)行相乘，后來沒有想出來，就直接在每次循環(huán)里面直接進(jìn)行乘法運算，如：for i=1:a for j=1:b distance=sqrt(i-a0)2+(j-b0)2); if distance<=d h=1; else h=0; end; s(i,j)=h*s(i,j); end;end;后來低通濾波這么一寫，高通濾波也就挺簡單的了，沒怎么想，把上面那個h=1和h=0的位置調(diào)換一下就變成了高通濾波器，至于高頻增強(qiáng)，就改用了s(i,j)=(p+q*h)