function [eqh, eqobr] = equal_hist(obrazek);
% funkce ekvalizuje histogram obrazku zadaneho jako parametr funkce a 
% transformuje obrazek podle ekvalizovaneho histogramu

%vypocet histogramu
histogram = hist(obrazek(:),0:255);

%provedeni ekvalizace podle vzorce ze skript (viz Hlavac Vaclav,
%Sedlacek Milos: Zaklady zpracovani signalu a obrazu, str. 145
pocet_radku = length(obrazek(:,1));
pocet_sloupcu = length(obrazek(1,:));

for i=1:pocet_radku
    for j=1:pocet_sloupcu
       eqobr(i,j) = uint8((255/(pocet_radku*pocet_sloupcu))*sum(histogram(1:obrazek(i,j))));
    end;
end; 

%vypocet ekvlizovaneho histogramu
eqh = hist(eqobr(:),0:255);
