Aerospace Toolbox

geod2geoc - Convert geodetic latitude to geocentric latitude

Syntax

gc = geod2geoc(gd, h) gc = geod2geoc(gd, h, model) gc = geod2geoc(gd, h, f, Re)

Description

gc = geod2geoc(gd, h) converts an array of m geodetic latitudes, gd, and an array of mean sea level altitudes, h, into an array of m geocentric latitudes, gc. Both gc and gd are in degrees. h is in meters.

gc = geod2geoc(gd, h, model) is an alternate method for converting from geodetic to geocentric latitude for a specific ellipsoid planet. Currently only 'WGS84' is supported for model.

gc = geod2geoc(gd, h, f, Re) is another alternate method for converting from geodetic to geocentric latitude for a custom ellipsoid planet defined by flattening, f, and the equatorial radius, Re, in meters.

Examples

Determine geocentric latitude given a geodetic latitude and altitude:

gc = geod2geoc(45, 1000)


gc =

   44.8076

Determine geocentric latitude at multiple geodetic latitudes and altitudes, specifying WGS84 ellipsoid model:

gc = geod2geoc([0 45 90], [1000 0 2000], 'WGS84')


gc =

         0
   44.8076
   90.0000

Determine geocentric latitude at multiple geodetic latitudes, given an altitude and specifying custom ellipsoid model:

f = 1/196.877360;
Re = 3397000;
gc = geod2geoc([0 45 90], 2000, f, Re)


gc =

         0
   44.7084
   90.0000

Assumptions and Limitations

This implementation generates a geocentric latitude that lies between ±90 degrees.

References