Initial bearing from one geo position to another
GeoDirection[{lat1, lon1}, {lat2, lon2}]
GeoDirection[loc1, loc2]
GeoDirection[{37, -109}, {40.112981, -88.261227}]
→ Quantity[72.78588619288448, AngularDegrees]GeoDirection[{46.42, -21.02}, {-30.23, 40.51}]
→ Quantity[130.1976535041487, AngularDegrees]GeoDirection[GeoPosition[{30., 40.}], GeoPosition[{-50., -60.}]]
→ Quantity[-133.71145709858484, AngularDegrees]