import holoviews as hv xmin,ymin,xmax,ymax = box = (22,65,32,75) square = hv.Rectangles(data=box).opts(fill_alpha=0.1) crop = da.sel(x=slice(xmin,xmax),y=slice(ymin,ymax),channel = ['R','G','B']) cropview = crop.hvplot.rgb('x','y','value',bands = 'channel',flip_yaxis=True,data_aspect=1) show((imview*square+cropview).opts(shared_axes=False).cols(1))