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)