pyacs.lib.coordinates.geo2xyz

pyacs.lib.coordinates.geo2xyz.geo2xyz(llambda, phi, he, unit='radians', A=6378137.0, E2=0.006694380022903)[source]

Convert geodetic coordinates (lon, lat, h) to geocentric XYZ.

Parameters:

  • llambda (float or array_like) – Longitude.

  • phi (float or array_like) – Latitude.

  • he (float or array_like) – Ellipsoidal height in meters.

  • unit ({'radians', 'dec_deg'}, optional) – Units for longitude and latitude. Default is ‘radians’.

  • A (float, optional) – Semi-major axis (equatorial radius), meters. Default 6378137.0 (GRS80).

  • E2 (float, optional) – Squared eccentricity. Default 0.006694380022903 (GRS80).

Returns:

x, y, z – Geocentric cartesian coordinates in meters.

Return type:

float or ndarray

Notes

Default ellipsoid is GRS80 (WGS84): A=6378137 m, E2=0.006694380022903, flattening F = 1 - sqrt(1-E2).