function image = myIRadon(radonim,theta,filter)
% =========================================================================
% image = myiradon(radonim,theta,filter)
% 
% Funkce rekonstruuje obraz z radonova obrazu
% 
% PARAMETRY:
% * radonim - vstupni sedotonovy obrazek o velikosti [p theta]
% * theta - vektor ze ktereho jen obraz rekonstruovan
% * filter - filter pro filtraci p v FT spektru
% 
% VYSTUP:
% * image - rekonstrukce CT obrazu
% 
% PRIKLAD:
%     filter = designFilter(flt.type, flt.len, flt.d);
%     image = myIRadon(data.imgCT,10,filter);
%     figure; subplot(131); imshow(data.imgCT); title('CT obrazek');
%     subplot(132); imagesc(image); title('Muje rekonstrukce');
%     subplot(133); imagesc(data.phantom); title('Original');
%   NEBO
%     image = myIRadon(data.imgCT,10);
%     figure; subplot(131); imshow(data.imgCT); title('CT obrazek');
%     subplot(132); imagesc(image); title('Muje rekonstrukce');
%     subplot(133); imagesc(data.phantom); title('Original');
% -------------------------------------------------------------------------

if nargin<2,  error('too few parameters');  end;

    % urceni velikosti obrazku
    
    % x = ceil( p / sqrt(2) );

    % inicializace vysledneho obrazku

    % vytvoreni gridu pro praci s interp2

    % LOOP   

        % filtrace P vectoru
        % IF % pouziti filtru
        
        % ELSE % bez fitru 

        % prepocet souradnic    

        % interpolace

        % integral
