Skip to content

Instantly share code, notes, and snippets.

@sandyUni
Forked from klucar/ECEF2LLA.matlab
Created June 13, 2017 06:38
Show Gist options
  • Save sandyUni/4e05196277518b3d778bc664be4f40f7 to your computer and use it in GitHub Desktop.
Save sandyUni/4e05196277518b3d778bc664be4f40f7 to your computer and use it in GitHub Desktop.
LLA2ECEF matlab
% LLA2ECEF - convert latitude, longitude, and altitude to
% earth-centered, earth-fixed (ECEF) cartesian
%
% USAGE:
% [x,y,z] = lla2ecef(lat,lon,alt)
%
% x = ECEF X-coordinate (m)
% y = ECEF Y-coordinate (m)
% z = ECEF Z-coordinate (m)
% lat = geodetic latitude (radians)
% lon = longitude (radians)
% alt = height above WGS84 ellipsoid (m)
%
% Notes: This function assumes the WGS84 model.
% Latitude is customary geodetic (not geocentric).
%
% Source: "Department of Defense World Geodetic System 1984"
% Page 4-4
% National Imagery and Mapping Agency
% Last updated June, 2004
% NIMA TR8350.2
%
% Michael Kleder, July 2005
function [x,y,z]=lla2ecef(lat,lon,alt)
% WGS84 ellipsoid constants:
a = 6378137;
e = 8.1819190842622e-2;
% intermediate calculation
% (prime vertical radius of curvature)
N = a ./ sqrt(1 - e^2 .* sin(lat).^2);
% results:
x = (N+alt) .* cos(lat) .* cos(lon);
y = (N+alt) .* cos(lat) .* sin(lon);
z = ((1-e^2) .* N + alt) .* sin(lat);
return
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment