Region of Interest (ROI) is Luxembourg

We use the communes adminitrative division which is standardized across Europe by EUROSTAT at:
https://ec.europa.eu/eurostat/web/gisco/geodata/reference-data/administrative-units-statistical-units
This is roughly equivalent to the notion municipalities in most countries.

From the link above, communes definition are taken from COMM_RG_01M_2016_4326.shp and country borders
are taken from NUTS_RG_01M_2021_3035.shp.

images: Sentinel2 rgb from 2020-01-01 to 2020-31-12
        filtered out pixels with >5% clouds acoording to MSK_CLDPRB and then we take the
        median at each pixel for the remaining of acquisitions during the year.
        see https://developers.google.com/earth-engine/datasets/catalog/COPERNICUS_S2_SR_HARMONIZED

        Google Earth Engine Python code
        
            def maskS2clouds(image):
                qa = image.select('MSK_CLDPRB')
                mask = qa.lt(5)
                return image.updateMask(mask)

            image_s2 = ee.ImageCollection('COPERNICUS/S2_SR_HARMONIZED')
                              .filterDate('2020-01-01', '2020-12-31')
                              .map(maskS2clouds)
                              .select('B4', 'B3', 'B2')
                              .median()
                              .visualize(min=0, max=4000)        


labels: ESA WorldCover 10m V100
        labels mapped to the interval [1,11] according to the following map
        { 0:0, 10: 1, 20:2, 30:3, 40:4, 50:5, 60:6, 70:7, 80:8, 90:9, 95:10, 100:11 }
        pixel value zero is reserved for invalid data (no cloudless pixels, errors
        in data, etc. )
        see https://developers.google.com/earth-engine/datasets/catalog/ESA_WorldCover_v100
        
        Google Earth Engine Python code
        
            image_esawc = ee.ImageCollection("ESA/WorldCover/v100").first()   
