Let’s see how to read a window instead of the whole bands of a given product and the potential time-consumption gain.

Warning: The durations shown hereunder may not be representative of your computer’s performances. Please take it as a hint about relative performances between constellations.

## Try with Landsat-8#

Let’s open a Landsat-8 OLI collection 2 tile.

# Imports
import os
import geopandas as gpd
import rasterio
from shapely.geometry import box

from rasterio.windows import Window, from_bounds