Page 195 -

โครงการหนังสืออิเล็กทรอนิกส์ เฉลิมพระเกียรติสมเด็จพระเทพรัตนราชสุดาฯ สยามบรมราชกุมารี

186

รูปที่ 6.22 ภาพที่ไดจากการผานตัวกรองแบบเมื่อ (a) D =20 (b) D =30
0
0

คําสั่งเปลี่ยนตัวกรองในโดเมนพื้นที่ใหเปนตัวกรองในโดเมนความถี่

เมื่อเราตองการใชตัวกรองที่ไดออกแบบไวแลวในโดเมนพื้นที่ เราสามารถใชคําสั่ง
freqz2(h, rows, cols) เพื่อใชแปลงตัวกรอง h ในโดเมนพื้นที่ ใหมาอยูในโดเมน

ความถี่ หลังการแปลงตัวกรองจะมีขนาดเทากับ rows×cols ตัวอยางเชน เราตองการ

กรองภาพดวยตัวกรองคาเฉลี่ยขนาด 3×3 เราสามารถทําการกรองในโดเมนความถี่ได
ดวยคําสั่ง MATLAB ตอไปนี้

>> im = imread (‘lfp_test.jpg’); %อานภาพอินพุต

>> h = fspecial(‘average’); %สรางตัวกรองคาเฉลี่ยขนาด3×3ในโดเมนพื้นที่
%แปลงตัวกรองใหอยูในโดเมนความถี่และใหมีขนาดเทากับขนาดภาพอินพุต

>> H = freqz2(h , size(im,2), size(im,1));
>> lpf_im = fftshift(fft2(im)).*H; % ทําการกรองในโดเมนความถี่

>> im_out = ifft2(lpf_im); %แปลงภาพที่ไดกลับใหมาอยูในโดเมนพื้นที่

>> mesh (H); % แสดงภาพตัวกรองแบบสามมิติ
>> title(‘Average lowpass filter’); axis tight % เขียนอธิบายรูป

>> figure, colormap(gray);imagesc(abs(im_out)); %แสดงผล