Page 224 -

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

215

การฟนฟูภาพโดยใชตัวกรองแบบปรับตัวชนิด LMS ขนาด 3×3 ตามคําสั่ง

MATLAB ดานลาง รูปที่ 7.5(ค) แสดงภาพที่ไดจากการฟนฟูโดยใชตัวกรองแบบ
ปรับตัวชนิด LMS ขนาด 7×7 ซึ่งจะเห็นไดวาภาพที่ไดจะเบลอมากขึ้นเมื่อขนาดของ

ตัวกรองใหญขึ้น รูปที่ 7.5(ง) แสดงภาพที่ไดจากการฟนฟูภาพโดยใชตัวกรองเฉลี่ย

(average filter) ขนาด 3×3 เมื่อเปรียบเทียบกับภาพที่ 7.5(ข) จะเห็นไดอยางชัดเจนวา

ตัวกรองแบบปรับตัวจะคงความคมชัดบริเวณขอบภาพไดดีกวาตัวกรองคาเฉลี่ยมาก

>> f = imread('lotus.jpg'); % อานภาพอินพุต

>> f = f(:,:,1); % พิจารณาเฉพาะภาพเพลนสีแดง

>> varN = 0.5;
% เพิ่มสัญญาณรบกวนเกาสเซีย ที่มีคาmean =จ คา var = 0.varN

>> g =imnoise(f, 'gaussian',0,var_noise); >> imshow(g);

>> g = double(g);
>> m = 3; n = 3; %ขนาดของหนากากที่ใช

>> mean = imfilter(g, ones(m,n), 'circular')/(m*n); %คํานวณคา local mean

>> fun = @(x) std2(x)^2; %กําหนดฟงกชันคํานวณคาความแปรปรวน
>> var = nlfilter(g, [m,n],fun); %คํานวณคา local variance

% ฟนฟูภาพ คําสั่ง eps จะใหคาเล็กๆอออกเพิ่มเขามาเพื่อหลีกเลี่ยงกรณีหารดวยศูนย

>> f_hat = g -((varN./(var+eps)).*[g-mean]);
>> figure,colormap(gray), imshow(uint8(f_hat));

7.4 การฟนฟูภาพเสื่อมจากสัญญาณรบกวนแบบคาบในโดเมนความถี่

(Periodic Noise Removal by Frequency Domain Filtering)

เนื่องจากสัญญาณรบกวนแบบคาบที่กลาวรายละเอียดไวในบทที่ 7.2 มีคุณสมบัติที่

เฉพาะในโดเมนความถี่ คือสัญญาณรบกวนแบบคาบจะปรากฎเปนกลุมจุดภาพใน
โดเมนความถี่เปนกลุมจุดที่มีขนาดของสัญญาณสูง ดังแสดงตัวอยางในภาพที่ 7.2