function [deltaN,detpair]=petMeasure(fantom,ti,tb,te,Nd,r,n,sigma2,taur)
% =========================================================================
% [deltaN,detpair]=petMeasure(fantom,ti,tb,te,Nd,r,sigma2,taur)
% 
% funkce ktera spocita dvojice tekovanych rozpadu spolu s poctem vsech
% rozpadu po dobu mereni
% 
% PARAMETRY:
% * fantom - realny fantom (obrazek) o velikosti [M N]
% * ti - cas ve kterem byla latka pacientovy podana v s
% * tb - cas kdy se zacal pacient snimat
% * te - cas kdy skoncilo snimani pacienta
% * Nd - pocet detektoru na celem obvodu kruznice
% * r - polomer krunice na ktere jsou rozmisteny detektory
% * n - mnozstvi podaneho radiofarmaka
% * sigma2 - rozptyl, jak daleho doleti castice pred svym rozpadem
% * taur - polocas rozpdu dane latky
% 
% VYSTUP:
% * deltaN - pocet rozpadlych castic po dobu mereni
% * detpair - dvojice detekovanych rozpadu zvolenymi tetektory
% 
% PRIKLAD:
%     load ph_def;
%     phantom=phantom(activity_def);
%     r=sqrt(sum((size(activity_phantom)./2).^2));
%     [deltaN,detpair]=petMeasure(phantom,1800,2100,2340,100,r,5e-12,3,6600);
%     figure; imagesc(detpair); 
% 
% -------------------------------------------------------------------------

    % pocet rozpadu

    % uprava fantomu podle poctu rozpadu

    % vektor pro vsechny hodnoty fantomu
    
    % nahodne smery vsech rozpadu

    % vygenrovani prislucnych souradnic k nahodnym smerum
    
    % urceni parametru primky - a,b,c

    % spocteni pruseciku s kruznici

    % zaokrouhleni (urceni ktery detektor detoval dany rozpad)

    % sumace vsech pruseciku pro danou kruznici
    % vytvoreni mapy vsech dvojic
