Convert geocentric (ECEF) to local vertical coordinates
Note:

[xl,yl,zl] = ecef2lv(x,y,z,phi0,lambda0,h0,ellipsoid)
[xl,yl,zl] = ecef2lv(x,y,z,phi0,lambda0,h0,ellipsoid)
converts
geocentric point locations specified by the coordinate arrays x
, y
,
and z
to the local vertical coordinate system,
with its origin at geodetic latitude phi0
, longitude lambda0
,
and ellipsoidal height h0
. The arrays x
, y
,
and z
may be of any shape, as long as they all
match in size. phi0
, lambda0
,
and h0
must be scalars. ellipsoid
is
a referenceEllipsoid
(oblateSpheroid
) object, a referenceSphere
object, or a vector
of the form [semimajor axis, eccentricity]
). x
, y
, z
,
and h0
must have the same length units as the semimajor
axis. phi0
and lambda0
must
be in radians. The output coordinate arrays, xl
, yl
,
and zl
are the local vertical coordinates of the
input points. They have the same size as x
, y
,
and z
and have the same length units as the semimajor
axis.
In the local vertical Cartesian system defined by phi0
, lambda0
, h0
,
and ellipsoid
, the xl
axis is
parallel to the plane tangent to the ellipsoid at (phi0
,lambda0
)
and points due east. The yl
axis is parallel to
the same plane and points due north. The zl
axis
is normal to the ellipsoid at (phi0
,lambda0
)
and points outward into space. The local vertical system is sometimes
referred to as eastnorthup or ENU.
For a definition of the geocentric system,
also known as EarthCentered, EarthFixed (ECEF), see the help for geodetic2ecef
.