ROZ2 – cvičení 2
description
Transcript of ROZ2 – cvičení 2
ROZ2 – cvičení 2
Image restoration- výsledky
Maska gaussiánu
function h = gauss(N, sigma)
PI = 4*atan(1);
npul = (N-1)/2;
[x,y] = meshgrid(-1*npul:npul);
h = 1/(2*PI*sigma^2) * exp(-1*(x.^2 + y.^2)/(2*sigma^2));
h = h / sum(h(:));
gauss(11,5)
gauss(21,7)
Poškození obrázku
function g = poskod(f, h, SNR)
g1 = conv2(f, h);
MinI = min(g1(:));
MaxI = max(g1(:));
var_f = var(f(:));
var_n = var_f / 10^(SNR / 10);
g = g1 + var_n*randn(size(g1));
g(g<MinI) = MinI;
g(g>MaxI) = MaxI;
Výsledky
gauss(11,3) gauss(21,7)
BSNR
1/0
60
50
Inverzní filtrfunction f = inverz(g, h)
G = fft2(g);
H = fft2(h, size(g,1), size(g,2));
F = G ./ H;
f = ifft2(F);
SNR
200
150
120 130 140
Wienerův filtrfunction f = wiener(g, h, konst)
G = fft2(g);
H = fft2(h, size(g,1), size(g,2));
R = conj(H) ./ (abs(H).^2 + konst);
F = R .* G;
f = ifft2(F);
Gauss(11,3)SNR 60
konst 0,1 0,01 0,001 0,0001
Rozmazání pohybem
h = ones(1,10);
Rozmazání pohybemm1 = log(abs(fft2(f).^2));
m2 = real(fft2(m1));
mi = min(m2(:));
m3 = m2 < 0.9*mi;
m3 m1 fftshift+výřez
Rozmazání defokusací
h =kruh(20,10);
Rozmazání defokusacím1 = log(abs(fft2(f).^2));
m2 = real(fft2(m1));
mi = min(m2(:));
m3 = m2 < 0.9*mi;
fftshift(m3) fftshift(m1)