pyart.core.cartesian_to_geographic_aeqd#

pyart.core.cartesian_to_geographic_aeqd(x, y, lon_0, lat_0, R=6370997.0)[source]#

Azimuthal equidistant Cartesian to geographic coordinate transform.

Transform a set of Cartesian/Cartographic coordinates (x, y) to geographic coordinate system (lat, lon) using a azimuthal equidistant map projection [1].

$$\begin{array}{r}\begin{array}{c}lat=\arcsin (\cos (c)\ast \sin (la{t}_0)+(y\ast \sin (c)\ast \cos (la{t}_0)/\rho ))\\ lon=lo{n}_0+\arctan 2(x\ast \sin (c),\rho \ast \cos (la{t}_0)\ast \cos (c)-y\ast \sin (la{t}_0)\ast \sin (c))\\ \rho =\sqrt{(}{x}^2+y^2)\\ c=\rho /R\end{array}\end{array}$$

Where x, y are the Cartesian position from the center of projection; lat, lon the corresponding latitude and longitude; lat_0, lon_0 are the latitude and longitude of the center of the projection; R is the radius of the earth (defaults to ~6371 km). lon is adjusted to be between -180 and 180.

Parameters:

• x, y (array-like) – Cartesian coordinates in the same units as R, typically meters.

• lon_0, lat_0 (float) – Longitude and latitude, in degrees, of the center of the projection.

• R (float, optional) – Earth radius in the same units as x and y. The default value is in units of meters.

Returns:

lon, lat (array) – Longitude and latitude of Cartesian coordinates in degrees.

References