Transform geodetic coordinates to geocentric Earth-centered Earth-fixed
[X,Y,Z] = geodetic2ecef(lat,lon,h,spheroid) is supported
but not recommended. Unlike the previous syntaxes, specify
lon in radians. Specify
spheroid as either a reference spheroid or an ellipsoid
vector of the form
[semimajor_axis, eccentricity]. Specify
h in the same units as the length unit of the
spheroid argument. Additionally, the outputs
in the same units as the length unit of the
Find the ECEF coordinates of Paris, France, using its geodetic coordinates.
First, specify the reference spheroid as WGS84 with length units measured in kilometers. For more information about WGS84, see Reference Spheroids. The units for the ellipsoidal height and ECEF coordinates must match the units specified by the
LengthUnit property of the reference spheroid.
wgs84 = wgs84Ellipsoid('kilometer');
Specify the geodetic coordinates of Paris. Specify
h as ellipsoidal height in kilometers.
lat = 48.8562; lon = 2.3508; h = 0.0674;
Then, calculate the ECEF coordinates of Paris. In this example,
y display in scientific notation.
[x,y,z] = geodetic2ecef(wgs84,lat,lon,h)
x = 4.2010e+03
y = 172.4603
z = 4.7801e+03
Reverse the transformation using the
[lat,lon,h] = ecef2geodetic(wgs84,x,y,z)
lat = 48.8562
lon = 2.3508
h = 0.0674
The geocentric Cartesian (ECEF) coordinate system is fixed with respect to the Earth, with its origin at the center of the spheroid and its positive x-, y-, and z-axes intersecting the surface at the following points:
|x-axis||0||0||Equator at the Prime Meridian|
|y-axis||0||90||Equator at 90-degrees East|