I have several Sentinel-2 products connected to ODC. My purpose is to use a shapefile with parcels in order to extract the average value of each parcel (creation of object based dataset). In addition, there is a raster (SLC) which shows the probability of cloud on each pixel. Therefore I do want to mask and leave out all the pixels mentioned as clouds.
How can I make the aggregation of the pixels per parcel and at the same time the cloud masking based on SLC (Until now this process takes place with gdal and pure Python and I need to implement this via ODC)?