Pyproj

import pyproj

geodesic = pyproj.Geod(ellps='WGS84')

fwd_azimuth, back_azimuth, distance = geodesic.inv(1, 1, 2, 2) # lon, lat