Helper function for cleaning Columbus P-1 datasets. Given lat and long coords in degree decimal, convert to radians and compute bearing.
calc_bearing(lat1, lon1, lat2, lon2)
latitude of starting point
longitude of starting point
latitude of ending point
longitude of ending point
bearing computed from given coordinates