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)