def test_deprecation_warning_geodetic_to_geocentric(self):
"""Test if geodetic_to_geocentric in coords is deprecated"""
with warnings.catch_warnings(record=True) as war:
coords.geodetic_to_geocentric(45.0, lon_in=8.0)
assert len(war) >= 1
assert war[0].category == DeprecationWarning
def test_geodetic_to_geocentric_inverse(self):
"""Tests the reversibility of geodetic to geocentric conversions"""
lat1 = 37.5
lon1 = 117.3
lat2, lon2, rad_e = coords.geodetic_to_geocentric(lat1, lon_in=lon1,
inverse=False)
lat3, lon3, rad_e = coords.geodetic_to_geocentric(lat2, lon_in=lon2,
inverse=True)
assert abs(lon1-lon3) < 1.0e-6
assert abs(lat1-lat3) < 1.0e-6
def test_geodetic_to_geocentric_single(self):
"""Test conversion from geodetic to geocentric coordinates"""
lat, lon, rad = coords.geodetic_to_geocentric(45.0, lon_in=8.0)
assert abs(lat - 44.807576784018046) < 1.0e-6
assert abs(lon - 8.0) < 1.0e-6
assert abs(rad - 6367.489543863465) < 1.0e-6
def test_geocentric_to_geodetic_single(self):
"""Test conversion from geocentric to geodetic coordinates"""
lat, lon, rad = coords.geodetic_to_geocentric(45.0, lon_in=8.0,
inverse=True)
assert abs(lat - 45.192423215981954) < 1.0e-6
assert abs(lon - 8.0) < 1.0e-6
assert abs(rad - 6367.345908499981) < 1.0e-6
def test_geodetic_to_geocentric_mult(self):
"""Test array conversion from geodetic to geocentric coordinates"""
arr = np.ones(shape=(10,), dtype=float)
lat, lon, rad = coords.geodetic_to_geocentric(45.0*arr, lon_in=8.0*arr)
assert lat.shape == arr.shape
assert lon.shape == arr.shape
assert rad.shape == arr.shape
assert abs(lat - 44.807576784018046).max() < 1.0e-6
assert abs(lon - 8.0).max() < 1.0e-6
assert abs(rad - 6367.489543863465).max() < 1.0e-6
def test_geocentric_to_geodetic_mult(self):
"""Test array conversion from geocentric to geodetic coordinates"""
arr = np.ones(shape=(10,), dtype=float)
lat, lon, rad = coords.geodetic_to_geocentric(45.0*arr, lon_in=8.0*arr,
inverse=True)
assert lat.shape == arr.shape
assert lon.shape == arr.shape
assert rad.shape == arr.shape
assert abs(lat - 45.192423215981954).max() < 1.0e-6
assert abs(lon - 8.0).max() < 1.0e-6
assert abs(rad - 6367.345908499981).max() < 1.0e-6
