%% Funkce prevodu snimku z kartezskych do polarnich
function [normalizovany,R_image,R_step, D_step] = car_pol(xs,ys,radius,I)
  % xs a ys  ...    souradnice stredu zornice (vnitri hranice duhovky)
  % radius   ...    vysledny polomer zornice
  % I       ...    puvodni sedotonovy snimek
  % normalizovan
%     pom = (size(im,1))-(ys + radius); % jaka je sirka duhovky od vnitri hranice z okraji snimku
    outer = 40;    % rucne zadana sirka duhovky
    degrees = 361; % pocet prevadenych stupnu
    R_image = I;  % puvodni snimek
    R_step = 0.65; % krok pro polomer
    D_step = 0.5;  % krok pro stupne
    rows = floor(outer/R_step);            % pocet radku normalizovaneho snimku
    columns = floor (degrees/D_step)-1;    % pocet sloupcu
    normalizovany = zeros(rows,columns); % preddefinovani normalizovaneho   
 
    %% zajmova oblast normalizovany = 0-360 stupnu
    index = 0;
    
    for int1 = 1:R_step:outer % pocet prochazenych polomeru
        index = index +1;
        index2 = 0;

        for int2 = 0:D_step:360  % stupne 
            index2 = index2 +1;

             x1 = xs + round(((radius + int1)*cos((int2*pi)/180))); % prevod na radiany (i*pi)180 
             y1 = ys + round(((radius + int1)*sin((int2*pi)/180)));

             normalizovany(index,index2) = im(y1, x1);

            R_image(y1,x1) = 255;

        end
    end
        
end