function FNorm = normalize(odata)
% entfernt das Rauschen und ersetzt die Intensi�ten mit Photonenanzahl

total_rows = uint32(size(odata,1))
z1 = input('Rauschschwelle des Akzeptor Kanals: ');
a1 = input('Gernzwert einzelnes Akzeptor Photonen: ');
b1 = input('Grenzwert von h�chstens 2 Akzeptor Photonen: ');
z2 = input('Rauschschwelle des Donator Kanals: ');
a2 = input('Gernzwert einzelnes Donator Photonen: ');
b2 = input('Grenzwert von h�chstens 2 Donator Photonen: ');
% alles dar�ber sind erstmal 3 Photonen
a=a1;
b=b1;
z=z1;
tic;
for m = 1:2
    for n = 1:total_rows
        if  odata(n,m) < z
            FNorm(n,m) = z;
         elseif odata(n,m) < a
                FNorm(n,m) = 1;      
         elseif odata(n,m) < b
                FNorm(n,m) = 2;
         else
            FNorm(n,m) = 3;
        end
    end
    a=a2
    b=b2
    z=z2
end;
toc
t=toc;

%schneller mit odata(odata<z1=0); odata(odata>a1=1...) 

