docs: def _tiles_from_cleaned_bounds(tp, bounds, zoom):
      """Return all tiles intersecting with bounds."""
      lb = _tile_from_xy(tp, bounds.left, bounds.bottom, zoom, on_edge_use="rt")
      rt = _tile_from_xy(tp, bounds.right, bounds.top, zoom, on_edge_use="lb")
      for tile_id in product([zoom], range(rt.row, lb.row + 1), range(lb.col, rt.col + 1)):
          yield tp.tile(*tile_id)