def test_GeodSolve0(self): point1 = GeoPoint(40.6, -73.8, frame=wgs84, degrees=True) point2 = GeoPoint(49.01666667, 2.55, frame=wgs84, degrees=True) s_ab, az_a, az_b = point1.distance_and_azimuth(point2, degrees=True) self.assertAlmostEqual(az_a, 53.47022, delta=0.5e-5) self.assertAlmostEqual(az_b, 111.59367, delta=0.5e-5) self.assertAlmostEqual(s_ab, 5853226, delta=0.5)
def test_GeodSolve0(self): point1 = GeoPoint(40.6, -73.8, frame=wgs84, degrees=True) point2 = GeoPoint(49.01666667, 2.55, frame=wgs84, degrees=True) s_ab, az_a, az_b = point1.distance_and_azimuth(point2, degrees=True) self.assertAlmostEqual(az_a, 53.47022, delta=0.5e-5) self.assertAlmostEqual(az_b, 111.59367, delta=0.5e-5) self.assertAlmostEqual(s_ab, 5853226, delta=0.5)
def test_wgs84_inverse(testcase): """Test inverse method on the WGS84 ellipsoid""" options = dict(frame=WGS84, degrees=True) (lat1, lon1, azi1, lat2, lon2, azi2, s12) = testcase[:7] point1 = GeoPoint(lat1, lon1, **options) point2 = GeoPoint(lat2, lon2, **options) s_ab, az_a, az_b = point1.distance_and_azimuth(point2, degrees=True) assert s_ab == approx(s12, rel=1e-14) assert az_a == approx(azi1, abs=2e-13) assert az_b == approx(azi2, abs=1e-13)
def test_inverse(self): options = dict(frame=wgs84, degrees=True) for l in GeodesicTest.testcases: (lat1, lon1, azi1, lat2, lon2, azi2, s12) = l[:7] point1 = GeoPoint(lat1, lon1, **options) point2 = GeoPoint(lat2, lon2, **options) s_ab, az_a, az_b = point1.distance_and_azimuth(point2, long_unroll=True, degrees=True) self.assertAlmostEqual(azi1, az_a, delta=1e-13) self.assertAlmostEqual(azi2, az_b, delta=1e-13) self.assertAlmostEqual(s12, s_ab, delta=1e-8)
def test_inverse(self): options = dict(frame=wgs84, degrees=True) for l in TESTCASES: (lat1, lon1, azi1, lat2, lon2, azi2, s12) = l[:7] point1 = GeoPoint(lat1, lon1, **options) point2 = GeoPoint(lat2, lon2, **options) s_ab, az_a, az_b = point1.distance_and_azimuth(point2, long_unroll=True, degrees=True) self.assertAlmostEqual(azi1, az_a, delta=1e-13) self.assertAlmostEqual(azi2, az_b, delta=1e-13) self.assertAlmostEqual(s12, s_ab, delta=1e-8)