docs: def get_catalogue_bounding_polygon(catalogue):
      '''
      Returns a polygon containing the bounding box of the catalogue
      '''
      upper_lon = np.max(catalogue.data['longitude'])
      upper_lat = np.max(catalogue.data['latitude'])
      lower_lon = np.min(catalogue.data['longitude'])
      lower_lat = np.min(catalogue.data['latitude'])
  
      return Polygon([Point(lower_lon, upper_lat), Point(upper_lon, upper_lat),
                      Point(upper_lon, lower_lat), Point(lower_lon, lower_lat)])