Create_an_AOI.py

# Request the administrative boundaries as a bytes dataset.
france_admin_url = "https://www.data.gouv.fr/en/datasets/r/17062524-991f-4e13-9bf0-b410cc2216fd"
france_admin = requests.get(france_admin_url, stream=True)
france_admin_bytes = bytes(france_admin.content)

# Import administrative boundaries into a Geodataframe
with fiona.BytesCollection(france_admin_bytes) as f:
    crs = f.crs
    france_df = gpd.GeoDataFrame.from_features(f, crs=crs)
    
# Clear variables
france_admin_bytes = None
france_admin = None

# Select the municipalities of interest based on their Insee code.
municipalities = france_df[france_df['insee'].isin(["17024", "17268", "17473"])]

aoi = municipalities.dissolve()