Code covered by the BSD License  

Highlights from
Geodesics on an ellipsoid of revolution

image thumbnail

Geodesics on an ellipsoid of revolution

by

 

20 Nov 2012 (Updated )

Solutions of the direct and inverse geodesic problems

defaultellipsoid
function ellipsoid = defaultellipsoid
%DEFAULTELLIPSOID  Return the WGS84 ellipsoid
%
%   ELLIPSOID = DEFAULTELLIPSOID
%
%   returns a vector of the equatorial radius and eccentricity for the
%   WGS84 ellipsoid.  Use ECC2FLAT and FLAT2ECC to convert between
%   the eccentricity and the flattening.
%
%   See also ECC2FLAT, FLAT2ECC.

  a = 6378137;
  f = 1/298.257223563;
  e = flat2ecc(f);
  ellipsoid = [a, e];
end

Contact us