Phép l c thơng cao Gauss 2 chi u có d ng:
, = − −�
2 ,
�2 . −
Ta dùng phép l c này để l c nh qu bí ngơ sau
Hình 1.5 – 31. nh qu bí ngơ
Kh i t o b l c trong Matlab
Thuật toán 1.5 – 6. L c thông cao Gauss
Input: ảnh cần lọc
Output: ảnh sau khi dùng bộ lọc thông cao Gauss
1 P = size(f, 1); %Tạo bộ lọc có kích thước bằng với ảnh 2 Q = size(f, 2); 3 h = zeros(P, Q); 4 5 sig = 40; %Gán � 6 b = 2.*sig.*sig; 7 8 for i = 1:P
35 9 for j = 1:Q 10 D = (i - P./2).^2 + (j - Q./2).^2; 11 h(i, j) = 1 - exp(-D./b); 12 end 13 end
14 H = fftshift(h); %Chuyển bộ lọc sang miền tần số
15 P = size(f, 1); %Tạo bộ lọc có kích thước bằng với ảnh 16 Q = size(f, 2); 17 h = zeros(P, Q); 18 19 sig = 40; %Gán � 20 b = 2.*sig.*sig; 21 22 for i = 1:P 23 for j = 1:Q 24 D = (i - P./2).^2 + (j - Q./2).^2; 25 h(i, j) = 1 - exp(-D./b); 26 end 27 end
28 H = fftshift(h); %Chuyển bộ lọc sang miền tần số
Ta đ ợc nh b l c trong mi n không gian và mi n t n s nh sau
Hình 1.5 – 32. nh b l c thông cao Gauss trong mi n khơng gian (bên trái, có d u ch m
đen nh gi a hình) sang mi n t n s (bên ph i)
Sử dụng thuật toán 1.5 – 3, ta tính tích chập c a nh g c trong mi n t n s và b l c trong mi n t n s , sau đó l y chuyển đổi ng ợc.
Hình 1.5 – 33. nh qu bí ngơ sau khi l c bằng b l c thông cao Gauss Ta quan sát phổ c a nh g c và nh sau khi l c (trên – ph i xu ng d i – trái)
36
Hình 1.5 – 34. So sánh t n s c a nh tr c và sau khi l c. nh trên – trái: nh ban đ u. nh trên – Ph i: nh sau khi l c. nh d i – trái: t n s nh ban đ u. nh d i – ph i:
t n s nh sau khi l c.
Nhìn nh, ta th y rằng nh khơng có hi u ng chng, phổ nh cho th y các giá tr t n s th p dao đ ng gi a t n s cao (biên qu bí ngơ).