cart2sphvec

Convert vector from Cartesian components to spherical representation

Syntax

vs = cart2sphvec(vr,az,el)

Description

vs = cart2sphvec(vr,az,el) converts the components of a vector or set of vectors, vr, from their representation in a local Cartesian coordinate system to a spherical basis representation contained in vs. A spherical basis representation is the set of components of a vector projected into a basis given by $$({\widehat{e}}_{az},{\widehat{e}}_{el},{\widehat{e}}_R)$$. The orientation of a spherical basis depends upon its location on the sphere as determined by azimuth, az, and elevation, el.

Examples

Spherical Representation of Unit Z-Vector

Start with a vector in Cartesian coordinates pointing along the z-direction and located at 45° azimuth, 45° elevation. Compute its components with respect to the spherical basis at that point.

vr = [0;0;1];
vs = cart2sphvec(vr,45,45)

vs = 3×1

         0
    0.7071
    0.7071

Input Arguments

vr — Vector in Cartesian basis representation

3-by-1 column vector | 3-by-N matrix

Vector in Cartesian basis representation specified as a 3-by-1 column vector or 3-by-N matrix. Each column of vr contains the three components of a vector in the right-handed Cartesian basis x,y,x.

Example: [4.0; -3.5; 6.3]

Data Types: double
Complex Number Support: Yes

az — Azimuth angle in degrees

scalar in range [–180, 180]

Azimuth angle in degrees, specified as a scalar in the closed range [–180, 180]. To define the azimuth angle of a point on a sphere, construct a vector from the origin to the point. The azimuth angle is the angle in the xy-plane from the positive x-axis to the vector's orthogonal projection into the xy-plane. As examples, zero azimuth angle and zero elevation angle specify a point on the x-axis while an azimuth angle of 90° and an elevation angle of zero specify a point on the y-axis.

Example: 45

Data Types: double

el — Elevation angle in degrees

scalar in range [–90, 90]

Elevation angle in degrees, specified as a scalar in the closed range [–90, 90]. To define the elevation of a point on the sphere, construct a vector from the origin to the point. The elevation angle is the angle from its orthogonal projection into the xy-plane to the vector itself. As examples, zero elevation angle defines the equator of the sphere and ±90° elevation define the north and south poles, respectively.

Example: 30

Data Types: double

Output Arguments

vs — Vector in spherical basis

3-by-1 column vector | 3-by-N matrix

Spherical representation of a vector returned as a 3-by-1 column vector or 3-by-N matrix having the same dimensions as vs. Each column of vs contains the three components of the vector in the right-handed $$({\widehat{e}}_{az},{\widehat{e}}_{el},{\widehat{e}}_R)$$ basis.

More About

Spherical basis representation of vectors

Spherical basis vectors are a local set of basis vectors which point along the radial and angular directions at any point in space.

The spherical basis is a set of three mutually orthogonal unit vectors $$({\widehat{e}}_{az},{\widehat{e}}_{el},{\widehat{e}}_R)$$ defined at a point on the sphere. The first unit vector points along lines of azimuth at constant radius and elevation. The second points along the lines of elevation at constant azimuth and radius. Both are tangent to the surface of the sphere. The third unit vector points radially outward.

The orientation of the basis changes from point to point on the sphere but is independent of R so as you move out along the radius, the basis orientation stays the same. The following figure illustrates the orientation of the spherical basis vectors as a function of azimuth and elevation:

For any point on the sphere specified by az and el, the basis vectors are given by:

$$\begin{array}{ll}{\widehat{e}}_{az}\hfill & =-\sin (az)\widehat{i}+\cos (az)\widehat{j}\hfill \\ {\widehat{e}}_{el}\hfill & =-\sin (el)\cos (az)\widehat{i}-\sin (el)\sin (az)\widehat{j}+\cos (el)\widehat{k}\hfill \\ {\widehat{e}}_R\hfill & =\cos (el)\cos (az)\widehat{i}+\cos (el)\sin (az)\widehat{j}+\sin (el)\widehat{k}\text{.}\hfill \end{array}$$

Any vector can be written in terms of components in this basis as $$v=v_{az}{\widehat{e}}_{az}+v_{el}{\widehat{e}}_{el}+v_R{\widehat{e}}_R$$. The transformations between spherical basis components and Cartesian components take the form

$$\left[\begin{array}{c}{v}_x\\ v_y\\ v_z\end{array}\right]=\left[\begin{array}{ccc}-\sin (az)& -\sin (el)\cos (az)& \cos (el)\cos (az)\\ \cos (az)& -\sin (el)\sin (az)& \cos (el)\sin (az)\\ 0& \cos (el)& \sin (el)\end{array}\right]\left[\begin{array}{c}{v}_{az}\\ v_{el}\\ v_R\end{array}\right]$$

.

and

$$\left[\begin{array}{c}{v}_{az}\\ v_{el}\\ v_R\end{array}\right]=\left[\begin{array}{ccc}-\sin (az)& \cos (az)& 0\\ -\sin (el)\cos (az)& -\sin (el)\sin (az)& \cos (el)\\ \cos (el)\cos (az)& \cos (el)\sin (az)& \sin (el)\end{array}\right]\left[\begin{array}{c}{v}_x\\ v_y\\ v_z\end{array}\right]$$.

Extended Capabilities

C/C++ Code Generation

Generate C and C++ code using MATLAB® Coder™.

Usage notes and limitations:

Does not support variable-size inputs.

Version History

Introduced in R2013a