pyacs.lib.coordinates.xyz2geo

pyacs.lib.coordinates.xyz2geo.xyz2geo(x, y, z, A=6378137.0, E2=0.006694380022903, unit='radians')[source]

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

Parameters:

  • x (float or array_like) – X coordinate(s) in meters.

  • y (float or array_like) – Y coordinate(s) in meters.

  • z (float or array_like) – Z coordinate(s) in meters.

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

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

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

Returns:

  • long (float or ndarray) – Longitude in requested unit.

  • lat (float or ndarray) – Latitude in requested unit.

  • he (float or ndarray) – Height above ellipsoid in meters.

Notes

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