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)