geographiclib_geodesic_distance_azimuth_inverse_direct_calculation.py

1from geographiclib.geodesic import Geodesic
2
3# The WGS84 ellipsoid
4geod = Geodesic.WGS84
5
6# Solve the inverse problem: compute the distance and azimuth from 
7# Wellington, New Zealand, to Salamanca, Spain
8# Wellington: 41.32S, 174.81E; Salamanca: 40.96N, 5.66W
9g = geod.Inverse(-41.32, 174.81, 40.96, -5.66)
10print("The distance is {:.3f} m.".format(g['s12']))
11print("The azimuth is {:.3f} degrees.".format(g['azi1']))
12
13# Solve the direct problem: compute the location of a point 10000 km
14# northeast of JFK Airport (40.64N, 73.78W)
15# JFK: 40.64N, 73.78W; Azimuth: 45 degrees; Distance: 10,000,000 m
16g = geod.Direct(40.64, -73.78, 45, 10e6)
17print("The position is ({:.8f}, {:.8f}).".format(g['lat2'], g['lon2']))
18print("The azimuth is {:.3f} degrees.".format(g['azi2']))