docs: def _get_voronoi_poly_points(vert_index_list, voronoi_vertices,
                               voronoi_centroid):
      """
      This function returns the corner points for a
      polygon from scipy voronoi information
      """
      voronoi_poly_points = []
      if -1 not in vert_index_list and len(vert_index_list) > 3:
          voronoi_poly_points = voronoi_vertices[vert_index_list]
      elif vert_index_list.size > 0:
          # ASSUME RECTANGLE
          vert_index_list = vert_index_list[vert_index_list >= 0]
          voronoi_poly_points = voronoi_vertices[vert_index_list]
          # CASE 1: 2 valid voronoi vertices
          if vert_index_list.size == 2:
              center_lon = voronoi_centroid[0]
              center_lat = voronoi_centroid[1]
              corner_lon1 = voronoi_poly_points[0][0]
              corner_lat1 = voronoi_poly_points[0][1]
              corner_lon2 = voronoi_poly_points[1][0]
              corner_lat2 = voronoi_poly_points[1][1]
  
              # check if need to add points in lon or lat
              if abs(corner_lon1-corner_lon2) > abs(corner_lat1-corner_lat2):
                  dLat = center_lat - corner_lat1
                  # append the corners in order
                  voronoi_poly_points = np.array([
                      [corner_lon1, corner_lat1],
                      [corner_lon2, corner_lat2],
                      [corner_lon2, center_lat + dLat],
                      [corner_lon1, center_lat + dLat]
                  ])
              else:
                  dLon = center_lon - corner_lon1
                  # append the corners in order
                  voronoi_poly_points = np.array([
                      [corner_lon1, corner_lat1],
                      [corner_lon2, corner_lat2],
                      [center_lon + dLon, corner_lat2],
                      [center_lon + dLon, corner_lat1]
                  ])
          # CASE 2: 1 valid voronoi vertex
          elif vert_index_list.size == 1:
              center_lon = voronoi_centroid[0]
              center_lat = voronoi_centroid[1]
              corner_lon = voronoi_poly_points[0][0]
              corner_lat = voronoi_poly_points[0][1]
              dLat = center_lat - corner_lat
              dLon = center_lon - corner_lon
              # append the corners in order
              voronoi_poly_points = np.array([
                  [corner_lon, corner_lat],
                  [center_lon + dLon, corner_lat],
                  [center_lon + dLon, center_lat + dLat],
                  [corner_lon, center_lat + dLat]
              ])
  
      return voronoi_poly_points