class test_GeodaticDatum(unittest.TestCase):
def setUp(self):
self.datum = GeodeticDatum('WGS84')
def test_toECEF(self):
x, y, z = self.datum.toECEF(0, 90)
nptest.assert_almost_equal(np.array([0, 0, self.datum.geod.b]),
np.array([x, y, z]),
decimal=5)
x, y, z = self.datum.toECEF(0, 0)
nptest.assert_almost_equal(np.array([self.datum.geod.a, 0, 0]),
np.array([x, y, z]),
decimal=5)
def test_ParallelRadi(self):
r = self.datum.ParallelRadi(0.)
nptest.assert_almost_equal(r, self.datum.geod.a, decimal=5)
r = self.datum.ParallelRadi(90.)
nptest.assert_almost_equal(r, 0., decimal=5)
def test_GeocentricDistance(self):
r = self.datum.GeocentricDistance(0., 0.)
nptest.assert_almost_equal(r, self.datum.geod.a, decimal=5)
r = self.datum.GeocentricDistance(0., 90.)
nptest.assert_almost_equal(r, self.datum.geod.b, decimal=5)
def test_ParallelArcDist(self):
dist = self.datum.ParallelArcDist(0., 0., 360.)
nptest.assert_almost_equal(dist, self.datum.geod.a * 2 * np.pi)
class test_GeodaticDatum(unittest.TestCase):
def setUp(self):
self.datum = GeodeticDatum('WGS84')
def test_toECEF(self):
x, y, z = self.datum.toECEF(0, 90)
nptest.assert_almost_equal(np.array([0, 0, self.datum.geod.b]),
np.array([x, y, z]),
decimal=5)
x, y, z = self.datum.toECEF(0, 0)
nptest.assert_almost_equal(np.array([self.datum.geod.a, 0, 0]),
np.array([x, y, z]),
decimal=5)
def test_ParallelRadi(self):
r = self.datum.ParallelRadi(0.)
nptest.assert_almost_equal(r, self.datum.geod.a, decimal=5)
r = self.datum.ParallelRadi(90.)
nptest.assert_almost_equal(r, 0., decimal=5)
def test_GeocentricDistance(self):
r = self.datum.GeocentricDistance(0., 0.)
nptest.assert_almost_equal(r, self.datum.geod.a, decimal=5)
r = self.datum.GeocentricDistance(0., 90.)
nptest.assert_almost_equal(r, self.datum.geod.b, decimal=5)
def test_ParallelArcDist(self):
dist = self.datum.ParallelArcDist(0., 0., 360.)
nptest.assert_almost_equal(dist, self.datum.geod.a * 2 * np.pi)
lat, lon1, lon2 = 45., -5., 5.
__, __, great_circle_dist = self.datum.geod.inv(lon1, lat, lon2, lat)
parallel_dist = self.datum.ParallelArcDist(lat, lon1, lon2)
assert great_circle_dist < parallel_dist, \
(great_circle_dist, parallel_dist)
