function [lat,lon,h]=ltp2geo(x,y,z,lat0,lon0,h0)
% in ENU coordinates x-->Northing
%                    y-->Easting
h0=h0;
R=rotct2lg(lat0*pi/180,lon0*pi/180);
[x0,y0,z0]=ell2xyz(lat0*pi/180,lon0*pi/180,h0);
p0=[x0
    y0
    z0];
a(1,:)=x;
a(2,:)=y;
a(3,:)=z;
p=R'*a+p0*ones(length(x));
[lat,lon,h]=xyz2ell3(p(1,:),p(2,:),p(3,:));
lat=lat*180/pi;
lon=lon*180/pi;
return