Nâng cao chất lượng ảnh
185
M1 = 0.15, SIG1 = 0.05, M2 = 0.75, SIG2 = 0.05, A1 = 1, A2 = 0.07, K =
0.002
c1 = A1 * (1 / ((2 * pi) ^ 0.5) * sig1);
k1 = 2 * (sig1 ^ 2);
c2 = A2 * (1 / ((2 * pi) ^ 0.5) * sig2);
k2 = 2 * (sig2 ^ 2);
z = linspace(0, 1, 256);
p = k + c1 * exp(–((z – m1) .^ 2) ./ k1) + ...
c2 * exp(–((z – m2) .^ 2) ./ k2);
p = p ./ sum(p(:));
% p: dạng histogram mong muốn
g = histeq(f, p)
figure, imshow(g)
% Ảnh sau khi phối hợp histogram
figure, imhist(g)
% và histogram tương ứng
a)
c)
e)
b)
d)
f)
Hình 12.9.