Index Page
cspice_sphlat
A  B  C  D  E  F  G  H  I  J  K  L  M  N  O  P  Q  R  S  T  U  V  W  X 

Abstract
I/O
Examples
Particulars
Required Reading
Version
Index_Entries

Abstract


   CSPICE_SPHLAT converts spherical coordinates to latitudinal
   coordinates.

I/O


   Given:

      r        the value(s) describing the distance of the position
               from the origin.

               [1,n] = size(r); double = class(r)

      colat    the value(s) describing the angle between the point and the
               positive z-axis, measured in radians (also referred to
               as the polar angle).

               [1,n] = size(colat); double = class(colat)

      lons     the value(s) describing the angle of the projection of the
               point to the XY plane from the positive X-axis, measured
               in radians, with range:

                   -pi < lons <= pi

               The positive Y-axis is at longitude PI/2 radians.

               [1,n] = size(lons); double = class(lons)

   the call:

      [radius, lon, lat] = cspice_sphlat(r, colat, lons)

   returns:

      radius   the value(s) describing the distance of the position
               from the origin.

               [1,n] = size(radius); double = class(radius)

      lon      the value(s) describing the angle of the position from
               the XZ plane measured in radians.

               [1,n] = size(lon); double = class(lon)

      lat      the value(s) describing the angle of the position from the
               XY plane measured in radians.

               [1,n] = size(lat); double = class(lat)

               The argument 'radius' returns in the same units associated
               with 'r'.

               'radius', 'lon', and 'lat' return with the same
                vectorization measure (N) as the 'r', 'colat',
                and 'lons'.

Examples


   Any numerical results shown for this example may differ between
   platforms as the results depend on the SPICE kernels used as input
   and the machine specific arithmetic implementation.

      %
      % Load an SPK, leapseconds, and PCK kernel set.
      %
      cspice_furnsh( 'standard.tm' )

      %
      % Create a vector of scalar times.
      %
      et = [0:2]*2.*cspice_spd;

      %
      % Retrieve the position of the moon seen from earth at 'et'
      % in the J2000 frame without aberration correction.
      %
      [pos, et] = cspice_spkpos( 'MOON', et, 'J2000', 'NONE', 'EARTH' );

      %
      % Convert the array of position vectors 'pos' to latitudinal
      % coordinates.
      %
      [r, colat, lons] = cspice_recsph(pos);

      %
      % Convert the latitudinal coords to spherical.
      %
      [ radius, lon, lat] = cspice_sphlat(r, colat, lons);

      %
      % Convert the spherical coords to rectangular.
      %
      [rectan] = cspice_latrec( radius, lon, lat);

      %
      % Calculate the relative error against the original position
      % vectors.
      %
      (rectan-pos) ./ pos

   MATLAB outputs:

      1.0e-14 *

                     0  -0.05552320600838   0.63783453323816
      0.02182376758148                  0  -0.01531271963894
      0.01912147275010  -0.02427608514229   0.02039513446643

   Example(2):

      %
      % Define six sets of spherical coordinates, 'lon' and 'colat'
      % expressed in degrees - converted to radians by use of cspice_rpd.
      %
      r     = [  1.,  1., 1.4142, 1.4142, 1.  , 0. ];
      colat = [ 90., 90., 45.   , 135.  , 0.  , 0. ] * cspice_rpd;
      lons  = [  0., 90., 180.  , 180.  , 180., 33.] * cspice_rpd;

      %
      % ...convert the latitudinal coordinates to spherical coordinates
      %
      [rad, lon, lat] = cspice_sphlat(r, colat, lons);

      %
      % ...convert angular measure to degrees.
      %
      colat = colat * cspice_dpr;
      lon   = lon   * cspice_dpr;
      lons  = lons  * cspice_dpr;
      lat   = lat   * cspice_dpr;

      %
      % Output banner.
      %
      disp('    r         colat       lons       rad        lon        lat   ')
      disp('  --------   --------   --------   --------   --------   --------')

      %
      % Create an array of values for output.
      %
      output = [ r; colat; lons; rad; lon; lat];
      txt    = sprintf( '%10.4f %10.4f %10.4f %10.4f %10.4f %10.4f\n', output);
      disp( txt )

      %
      % It's always good form to unload kernels after use,
      % particularly in MATLAB due to data persistence.
      %
      cspice_kclear

   MATLAB outputs:

       r         colat       lons       rad        lon        lat
     --------   --------   --------   --------   --------   --------
       1.0000    90.0000     0.0000     1.0000     0.0000     0.0000
       1.0000    90.0000    90.0000     1.0000    90.0000     0.0000
       1.4142    45.0000   180.0000     1.4142   180.0000    45.0000
       1.4142   135.0000   180.0000     1.4142   180.0000   -45.0000
       1.0000     0.0000   180.0000     1.0000   180.0000    90.0000
       0.0000     0.0000    33.0000     0.0000    33.0000    90.0000

Particulars


   None.

Required Reading


   For important details concerning this module's function, please refer to
   the CSPICE routine sphlat_c.

   MICE.REQ

Version


   -Mice Version 1.0.1, 01-DEC-2014, EDW (JPL)

       Edited I/O section to conform to NAIF standard for Mice documentation.

   -Mice Version 1.0.0, 12-DEC-2005, EDW (JPL)

Index_Entries


   spherical to latitudinal coordinates


Wed Apr  5 18:00:34 2017