1. Imports, setting up, data sources, data preparation
import rasterio
from rasterio.plot import show
from rasterio.mask import mask
import matplotlib.colors as mcolors
from shapely.geometry import box
import geopandas as gpd
import numpy as np
import matplotlib.pyplot as plt
import math
from matplotlib.colors import LinearSegmentedColormap
# File path t…