import pyproj
geodesic = pyproj.Geod(ellps='WGS84')
fwd_azimuth, back_azimuth, distance = geodesic.inv(1, 1, 2, 2) # lon, lat