code: def bounds_window(bounds, affine):
      
      w, s, e, n = bounds
      row_start, col_start = rowcol(w, n, affine)
      row_stop, col_stop = rowcol(e, s, affine, op=math.ceil)
      return (row_start, row_stop), (col_start, col_stop)
docs: Create a full cover rasterio-style window