docs: def get_bearing(origin_point, destination_point):
      """
      Calculate the bearing between two lat-long points. Each tuple should
      represent (lat, lng) as decimal degrees.
  
      Parameters
      ----------
      origin_point : tuple
      destination_point : tuple
  
      Returns
      -------
      bearing : float
          the compass bearing in decimal degrees from the origin point
          to the destination point
      """
  
      if not (isinstance(origin_point, tuple) and isinstance(destination_point, tuple)):
          raise TypeError('origin_point and destination_point must be (lat, lng) tuples')
  
      # get latitudes and the difference in longitude, as radians
      lat1 = math.radians(origin_point[0])
      lat2 = math.radians(destination_point[0])
      diff_lng = math.radians(destination_point[1] - origin_point[1])
  
      # calculate initial bearing from -180 degrees to +180 degrees
      x = math.sin(diff_lng) * math.cos(lat2)
      y = math.cos(lat1) * math.sin(lat2) - (math.sin(lat1) * math.cos(lat2) * math.cos(diff_lng))
      initial_bearing = math.atan2(x, y)
  
      # normalize initial bearing to 0 degrees to 360 degrees to get compass bearing
      initial_bearing = math.degrees(initial_bearing)
      bearing = (initial_bearing + 360) % 360
  
      return bearing