+ All Categories
Home > Documents > ROZ2 – cvičení 2

ROZ2 – cvičení 2

Date post: 12-Jan-2016
Category:
Upload: naiya
View: 50 times
Download: 0 times
Share this document with a friend
Description:
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. - PowerPoint PPT Presentation
10
ROZ2 – cvičení 2 Image restoration - výsledky
Transcript
Page 1: ROZ2 –  cvičení 2

ROZ2 – cvičení 2

Image restoration- výsledky

Page 2: ROZ2 –  cvičení 2

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)

Page 3: ROZ2 –  cvičení 2

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;

Page 4: ROZ2 –  cvičení 2

Výsledky

gauss(11,3) gauss(21,7)

BSNR

1/0

60

50

Page 5: ROZ2 –  cvičení 2

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

Page 6: ROZ2 –  cvičení 2

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

Page 7: ROZ2 –  cvičení 2

Rozmazání pohybem

h = ones(1,10);

Page 8: ROZ2 –  cvičení 2

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

Page 9: ROZ2 –  cvičení 2

Rozmazání defokusací

h =kruh(20,10);

Page 10: ROZ2 –  cvičení 2

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)


Recommended