#vypocet sumu pozadia
load "pozadie.dat"
pozadie = pozadie(:,3);
p = mean(pozadie)

# bllac
load "bl.dat"

#priemer stvorca
velkost = sqrt(size(bl(:,1),1));

#matica s odpocitanym sumom
M = [bl(:,1) bl(:,2)  bl(:,3) - p];

#definovane premenne
r = velkost/2;
x0 = velkost/2;
y0 = velkost/2;
c = 0;
#prepocet toku v danom filtri R
d = 3e8*140e-9/((666e-9)**2);
#f0 pre R filter
f0 = 3.08e-26;

#suma hodnot v kruhu vyrezaneho obrazka
for i=1:velkost
for j=1:velkost
if (sqrt((i-x0)**2+(j-y0)**2)<r)
c = c + M((i+velkost*(j-1)),3) ;
endif
endfor
endfor

#vypocitana hodnota na meter stvorcovi za sekundu
C_bl = c/(0.3**2*pi*60)

#Pomocí Pogsonovy rovnice vypočteme hustoty toky pro známé objekty


# B
load "B.dat"

#priemer stvorca
velkost = sqrt(size(B(:,1),1));

#matica s odpocitanym sumom
M = [B(:,1) B(:,2)  B(:,3) - p];

#definovane premenne
r = velkost/2;
x0 = velkost/2;
y0 = velkost/2;
c = 0;

#suma hodnot v kruhu vyrezaneho obrazka
for i=1:velkost
for j=1:velkost
if (sqrt((i-x0)**2+(j-y0)**2)<r)
c = c + M((i+velkost*(j-1)),3) ;
endif
endfor
endfor

#vypocitana hodnota na meter stvorcovi za sekundu
C_B = c/(0.3**2*pi*60)

#Pomocí Pogsonovy rovnice vypočteme hustoty toky pro známé objekty
m = 11.93
f = f0*10**(-0.4*m)

#hustota pre filter
F = d*f

#priemer pomeru medzi inštrumentalnymi a kalibracnymi tokmi
r=F/C_B 

# C
load "C.dat"

#priemer stvorca
velkost = sqrt(size(C(:,1),1));

#matica s odpocitanym sumom
M = [C(:,1) C(:,2)  C(:,3) - p];

#definovane premenne
r = velkost/2;
x0 = velkost/2;
y0 = velkost/2;
c = 0;

#suma hodnot v kruhu vyrezaneho obrazka
for i=1:velkost
for j=1:velkost
if (sqrt((i-x0)**2+(j-y0)**2)<r)
c = c + M((i+velkost*(j-1)),3) ;
endif
endfor
endfor

#vypocitana hodnota na meter stvorcovi za sekundu
C_C = c/(0.3**2*pi*60)

#Pomocí Pogsonovy rovnice vypočteme hustoty toky pro známé objekty
m = 13.69
f = f0*10**(-0.4*m)

#hustota pre filter
F = d*f

#priemer pomeru medzi inštrumentalnymi a kalibracnymi tokmi
r=F/C_C

# K
load "K.dat"

#priemer stvorca
velkost = sqrt(size(K(:,1),1));

#matica s odpocitanym sumom
M = [K(:,1) K(:,2)  K(:,3) - p];

#definovane premenne
r = velkost/2;
x0 = velkost/2;
y0 = velkost/2;
c = 0;

#suma hodnot v kruhu vyrezaneho obrazka
for i=1:velkost
for j=1:velkost
if (sqrt((i-x0)**2+(j-y0)**2)<r)
c = c + M((i+velkost*(j-1)),3) ;
endif
endfor
endfor

#vypocitana hodnota na meter stvorcovi za sekundu
C_K = c/(0.3**2*pi*60)

#Pomocí Pogsonovy rovnice vypočteme hustoty toky pro známé objekty
m = 14.88
f = f0*10**(-0.4*m)

#hustota pre filter
F = d*f

#priemer pomeru medzi inštrumentalnymi a kalibracnymi tokmi
r=F/C_K

# H
load "H.dat"

#priemer stvorca
velkost = sqrt(size(H(:,1),1));

#matica s odpocitanym sumom
M = [H(:,1) H(:,2)  H(:,3) - p];

#definovane premenne
r = velkost/2;
x0 = velkost/2;
y0 = velkost/2;
c = 0;

#suma hodnot v kruhu vyrezaneho obrazka
for i=1:velkost
for j=1:velkost
if (sqrt((i-x0)**2+(j-y0)**2)<r)
c = c + M((i+velkost*(j-1)),3) ;
endif
endfor
endfor

#vypocitana hodnota na meter stvorcovi za sekundu
C_H = c/(0.3**2*pi*60)

#Pomocí Pogsonovy rovnice vypočteme hustoty toky pro známé objekty
m = 13.60
f = f0*10**(-0.4*m)

#hustota pre filter
F = d*f

#priemer pomeru medzi inštrumentalnymi a kalibracnymi tokmi
r=F/C_H

