geographiclib_geodesic_inverse_direct_distance_azimuth_calculation.py

1from geographiclib.geodesic import Geodesic
2
3# The geodesic object for the WGS84 ellipsoid
4geod = Geodesic.WGS84
5
6# Solve the inverse problem: distance and azimuth between two points
7# Wellington, New Zealand: 41.2833° S, 174.7667° E
8# Salamanca, Spain: 40.9667° N, 5.6667° W
9g = geod.Inverse(-41.2833, 174.7667, 40.9667, -5.6667)
10print("The distance is {:.3f} m.".format(g['s12']))
11
12# Solve the direct problem: destination point given start, azimuth, and distance
13# JFK Airport: 40.64° N, 73.78° W
14# Paris, France: azimuth 51°, distance 5849 km
15g = geod.Direct(40.64, -73.78, 51, 5849e3)
16print("The destination is ({:.4f}, {:.4f}).".format(g['lat2'], g['lon2']))