geographiclib_wgs84_direct_inverse_geodesic_calculation.py

1from geographiclib.geodesic import Geodesic
2
3# The geodesic object for the WGS84 ellipsoid
4geod = Geodesic.WGS84
5
6# Solve the direct problem:
7# Compute the destination given a starting point, an azimuth, and a distance.
8# Input: lat1, lon1, azi1, s12
9d = geod.Direct(40.6, -73.8, 51, 5500e3)
10print("The distance is {:.3f} m.".format(d['s12']))
11print("The destination is ({:.5f}, {:.5f}).".format(d['lat2'], d['lon2']))
12
13# Solve the inverse problem:
14# Compute the distance and azimuths between two points.
15# Input: lat1, lon1, lat2, lon2
16g = geod.Inverse(40.6, -73.8, 51.6, -0.5)
17print("The distance is {:.3f} m.".format(g['s12']))
18print("The azimuths are {:.5f} and {:.5f} degrees.".format(g['azi1'], g['azi2']))