function[Lat,Lon]=DF17_poloha_surface(DF,latS,lonS)
%DF=binární 112 bitová zpráva
%latS= 50.10083333 %Referencni zemepisna sirka
%lonS= 14.26000000 %Referencni zemepisna delka 

CPR_format=DF(54);% Formát CPR odd/even
CPR_lat=DF(55:71);
CPR_lon=DF(72:88);

YZ=polyval(CPR_lat,2); % pøevede binární èíslo na decinální
XZ=polyval(CPR_lon,2);

NZ=15; % Pocet sektoru
Nb=length(CPR_lat);
%------------------------------Zemepisna sirka-----------------------------
Dlat=90/(4*NZ-CPR_format);

j1=latS/Dlat;
j2=0.5+(mod(latS,Dlat))/Dlat-YZ/(2^17);
j=floor(j1)+floor(j2);

Lat=Dlat*(j+YZ/(2^Nb))

%------------------------------Zemepisna delka-----------------------------
NL_Rlat1=1-(1-cos(pi/(2*NZ)))/((cos((pi/180)*abs(Lat)))^2);
NL_Rlat2=acos(NL_Rlat1);
NL_Rlat3=NL_Rlat2^(-1);
NL_Rlat4=NL_Rlat3*2*pi;
NL_Rlat=floor(NL_Rlat4);

NL_Rlat_Dlon=NL_Rlat-CPR_format;

if NL_Rlat_Dlon==0
    Dlon=90;
else
    Dlon=90/NL_Rlat_Dlon;
end
m1=lonS/Dlon;
m2=0.5+(mod(lonS,Dlon))/Dlon-XZ/(2^17);
m=floor(m1)+floor(m2);

Lon=Dlon*(m+XZ/(2^17))
end
