code: def uniform_random_global_points(n=100):
      
      xs = np.random.uniform(-180, 180, n)
      ys = np.random.uniform(-90, 90, n)
      return [shapely.geometry.Point(x, y) for x, y in zip(xs, ys)]
docs: Returns an array of `n` uniformally distributed `shapely.geometry.Point` objects. Points are coordinates
      distributed equivalently across the Earth's surface.