Амалий иш №9. Частота муҳитида фильтрларни қуриш
Амалий иш мақсади: Частота муҳитида фильтрларни қуришни ўрганиш ва уларни MATLAB дастурида қуриш.
Частота ва фазовий муҳитида чизиқли фильтрлашнинг асосини свертка теоремаси ташкил этади ва уни қуйидаги кўринишда ёзиш мумкин:
ва тескари ҳолатда ҳам
Б у ерда икки функциянинг свертка операциясини англатади, стрелка эса Фурье ўзгартиришда мос жуфтликларни аниқлайди.
Частота муҳитида фильтрлашнинг асосий ғоясини бериладиган фильтрни функцияларини танлаб олишдир, қайсики, F(u,v) функцияси билан ўзига хос шаклга келтиришдир. Частота муҳитида фильтрлаш қуйидаги кетма-кетлик асосидй амалга оширилади.
Мисол:
image = imread('E:\Магистратура\ТРИБ\TIBV\rasm1.tif');
[М, N] = size(image);
F = fft2(image);
sig = 10;
Н = lpfilter('gaussian', M, N, sig);
G = H.*F;
g = real(ifft(2(G));
imshow(g, [ ]);
Мисол:
F = fft2(f);
S = fftshift(log(l + abs(F)));
S = gscale(S) ;
imshow(S);
h = fspecial('sobel');
freqz2(h)
PQ = paddedsize(size(f));
Н = freqz2(h, PQA), PQB));
HI = ifftshift(H);
imshow(abs(H) , [ ]) ;
figure, imshow(abs(Hl) , [ ]) ;
gs = imfilter (double(f) , h) ;
gf = dftfilt(f, HI) ;
imshow(gs, [ ]) ;
figure, imshow(gf, [ ]) ;
figure, imshow(abs(gs) , [ ]) ;
figure, imshow (abs (gf), [ ]) ;
figure, imshow(abs(gs) > 0.2*abs(max(gs(:))));
figure, imshow(abs(gf) > 0.2*abs(max(gf(:))));
Do'stlaringiz bilan baham: |