def _test_geodetic2cart2geodetic(self, ellipsoid):
"""Test conversion from geodetic to cartesian system and back."""
ref = (1, -13, 42)
cart = geodesy.geodetic2cart(*ref, ellipsoid)
geod = geodesy.cart2geodetic(*cart, ellipsoid)
assert np.allclose(ref, geod)
def test_geodetic2cart2geodetic(self):
"""Test conversion from geodetic to cartesian system and back."""
ellipsoid = geodesy.ellipsoidmodels().get('WGS84')
ref = (1, -13, 42)
cart = geodesy.geodetic2cart(*ref, ellipsoid)
geod = geodesy.cart2geodetic(*cart, ellipsoid)
assert np.allclose(ref, geod)
def test_cart2geodetic(self):
"""Test conversion from cartesian to geodetic system."""
r = geodesy.ellipsoidmodels().get('WGS84')[0]
cartesian = (np.array([r, 0]), # x
np.array([0, r]), # y
np.array([0, 0]), # z
)
reference = (np.array([0, 0]), # r
np.array([0, 00]), # lat
np.array([0, 90]), # lon
)
conversion = geodesy.cart2geodetic(*cartesian)
assert np.allclose(conversion, reference)