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)
