code: def get_boundary_coords(x, y, spatial_pad=0):
      r
      west = np.min(x) - spatial_pad
      east = np.max(x) + spatial_pad
      north = np.max(y) + spatial_pad
      south = np.min(y) - spatial_pad
  
      return {'west': west, 'south': south, 'east': east, 'north': north}
docs: r"""Return bounding box based on given x and y coordinates assuming northern hemisphere.
  
      x: numeric
          x coordinates.
      y: numeric
          y coordinates.
      spatial_pad: numeric
          Number of meters to add to the x and y dimensions to reduce
          edge effects.
  
      Returns
      -------
      bbox: dictionary
          dictionary containing coordinates for corners of study area