Rasterio Hint

get-available-tiles.py

import mercantile
import rasterio
from rasterio.session import OSSSession
from rasterio.warp import transform_bounds
from rio_tiler import constants

with rasterio.open(prefix) as src:
        bounds = transform_bounds(
            src.crs, constants.WGS84_CRS, *src.bounds, densify_pts=21
        )
        mintile = mercantile.tile(bounds[0], bounds[3], 10)
        maxtile = mercantile.tile(bounds[2], bounds[1], 10)
        print(mintile,maxtile)