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

1、實(shí)驗(yàn)作業(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:afor j=1:bdistance=sqrt(i-a0)A2+(j-b0)A2);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('低通濾波所得圖像'); 工祝"5 "1： L1S1A- ” 工：工濾波d=10圖像詣里葉玄也班浮頒諳圣像博里葉變換班對(duì)數(shù),所諄翅語(yǔ)2、用理想高通濾波器在頻率域?qū)崿F(xiàn)高頻增強(qiáng)I=imread('d:/3.jpg');I=rgb2gray(I);figure(1),imshow(I);title('原圖像');s=fftshift(f

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

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