def test_geocentric2cart2geocentric(self):
"""Test conversion from geocentric to cartesian system and back."""
ref = (1, -13, 42)
cart = geodesy.geocentric2cart(*ref)
geo = geodesy.cart2geocentric(*cart)
assert np.allclose(ref, geo)
def test_geocentric2cart2geocentric(self):
"""Test conversion from geocentric to cartesian system and back."""
ref = (1, -13, 42)
cart = geodesy.geocentric2cart(*ref)
geo = geodesy.cart2geocentric(*cart)
assert np.allclose(ref, geo)
def test_cart2geocentric(self):
"""Test conversion from cartesian to geocentric system."""
cartesian = (np.array([1, 0, 0]), np.array([0, 1, 0]), np.array([0, 0, 1])) # x # y # z
reference = (np.array([1, 1, 1]), np.array([0, 0, 90]), np.array([0, 90, 0])) # r # lat # lon
conversion = geodesy.cart2geocentric(*cartesian)
assert np.allclose(conversion, reference)
def test_cart2geocentric(self):
"""Test conversion from cartesian to geocentric system."""
cartesian = (np.array([1, 0, 0]), # x
np.array([0, 1, 0]), # y
np.array([0, 0, 1]), # z
)
reference = (np.array([1, 1, 1]), # r
np.array([0, 0, 90]), # lat
np.array([0, 90, 0]), # lon
)
conversion = geodesy.cart2geocentric(*cartesian)
assert np.allclose(conversion, reference)
