geod2geoc

Convert geodetic latitude to geocentric latitude

Syntax

geocentricLatitude = geod2geoc(geodeticLatitude,height)

[geocentricLatitude,radii] = geod2geoc(geodeticLatitude,height)

geocentricLatitude = geod2geoc(geodeticLatitude,height,model)

[geocentricLatitude,radii] = geod2geoc(geodeticLatitude,height,model)

geocentricLatitude = geod2geoc(geodeticLatitude,height,flattening,Re)

[geocentricLatitude,radii] = geod2geoc(geodeticLatitude,height,flattening,Re)

Description

WGS84 Ellipsoid Planet

geocentricLatitude = geod2geoc(geodeticLatitude,height) converts an array of geodetic latitudes, geodeticLatitude, and an array of heights above the planetary ellipsoid, height, into an array of geocentric latitudes, geocentricLatitude.

[geocentricLatitude,radii] = geod2geoc(geodeticLatitude,height) returns an array of radii radii from the center of the planet to the center of gravity.

Specific Ellipsoid Planet

geocentricLatitude = geod2geoc(geodeticLatitude,height,model) converts from geodetic to geocentric latitude for a specific ellipsoid planet.

[geocentricLatitude,radii] = geod2geoc(geodeticLatitude,height,model) returns the radius radii from the center of the planet to the center of gravity.

Custom Ellipsoid Planet

geocentricLatitude = geod2geoc(geodeticLatitude,height,flattening,Re) converts from geodetic to geocentric latitude for a custom ellipsoid planet defined by flattening, flattening, and the equatorial radius, equatorialRadius.

[geocentricLatitude,radii] = geod2geoc(geodeticLatitude,height,flattening,Re) returns the radius radii from the center of the planet to the center of gravity.

Examples

Determine Geocentric Latitude and Radius with Geodetic Latitude and Altitude

This example shows how to determine geocentric latitude given a geodetic latitudes and ellipsoidal altitudes.

[gc,r] = geod2geoc(45,1000)

gc = 
44.8076

r = 
6.3685e+06

Determine Geocentric Latitudes and Radii at Multiple Geodetic Latitudes and Altitudes, Specifying WGS84 Ellipsoid Model

This example shows how to determine geocentric latitude at multiple geodetic latitudes and ellipsoidal altitudes, specifying WGS84 ellipsoid model.

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

gc = 1×3

         0   44.8076   90.0000

r = 1×3
106 ×

    6.3791    6.3675    6.3588

Determine Geocentric Latitudes and Radii at Multiple Geodetic Latitudes, Given Altitude and Custom Ellipsoid Model

This example shows how to determine geocentric latitude at multiple geodetic latitudes, given multiple geodetic latitudes that specify a custom ellipsoid model.

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

gc = 1×3

         0   44.7084   90.0000

r = 1×3
106 ×

    3.3990    3.3904    3.3817

Input Arguments

geodeticLatitude — Geodetic latitude

array

Geodetic latitude, specified as an array in degrees.

Data Types: double

height — Height above the planetary ellipsoid

scalar

Height above the planetary ellipsoid, specified as a scalar in meters.

Data Types: double

model — Specific ellipsoid planet model

'WGS84'

Specific ellipsoid planet model, specified as 'WGS84'.

Data Types: double

flattening — Flattening

scalar

Flattening at each pole, specified as a scalar.

Data Types: double

Re — Equatorial radius

scalar

Equatorial radius, specified as a scalar in meters.

Data Types: double

Output Arguments

geocentricLatitude — Geocentric latitudes

array

Geocentric latitudes, returned as an array in degrees. Latitude values can be any value. However, values of +90 and -90 may return unexpected values because of singularity at the poles.

radii — Radii from center of planet

array

Radii from the center of the planet, returned as an array in meters.

Limitations

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

References

[1] Stevens, Brian L., Frank L. Lewis, Aircraft Control and Simulation, Wiley–Interscience, 1992.

Version History

Introduced in R2006b