Main Content

geod2geoc

Convert geodetic latitude to geocentric latitude

collapse all in page

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.

example

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.

example

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.

example

Examples

collapse all

Open Live Script

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
Open Live Script

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
Open Live Script

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

collapse all

Geodetic latitude, specified as an array in degrees.

Data Types: double

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

Data Types: double

Specific ellipsoid planet model, specified as 'WGS84'.

Data Types: double

Flattening at each pole, specified as a scalar.

Data Types: double

Equatorial radius, specified as a scalar in meters.

Data Types: double

Output Arguments

collapse all

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 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

See Also

| | | |