dg2lg.m

Geodetic tools to be used in navigation analysis

function [dx,dy]=dg2lg(dlat,dlon,lat,h,a,e2)
% DG2LG  Converts 恖at,恖on,恏 to local geodetic coordinates.
%   Local origin at lat,lon.  If astronomic lat,h input,
%   then output is in local astronomic system.  Vectorized.
%   See also LG2DG.
% Version: 18 Jan 96
% Useage:  [dx,dy]=dg2lg(dlat,dlon,lat,h,a,e2)
% Input:   dlat - vector of latitude differences (rad)
%          dlon - vector of longitude differences (rad)
%          lat  - vector of lats of local system origins (rad)
%          h    - vector of hts of local system origins
%          a    - ref. ellipsoid major semi-axis
%          e2   - ref. ellipsoid eccentricity squared
% Output:  dx   - vector of x (N) coordinates in local system
%          dy   - vector of y (E) coordinates in local system
v=a./sqrt(1-e2.*sin(lat).^2);
r=v.*(1-e2)./(1-e2.*sin(lat).^2);
dx=dlat.*(r+h);
dy=dlon.*cos(lat).*(v+h);