jistiak

import geopandas as gpd
import rioxarray as rxr
from shapely.geometry import Polygon

aoi = gpd.read_file("/PATH/TO/.shp")

image = rxr.open_rasterio("/PATH/TO/.tif")
xmin, ymin, xmax, ymax  = image.rio.bounds()
im_poly = Polygon([(xmin, ymin), (xmin, ymax), (xmax, ymax), (xmax, ymin)])
intersection = aoi.intersection(im_poly)