I don't think it would make sense to get a cloud score from an imageCollection that has been reduced to a single median image. Instead, by removing the cloud mask, median reducer, and sorting step when you assign your imageCollection, you can map the getCloudScores function to all images in the collection. Note that I include a different district because ...
You do not need to do any further atmospheric correction with Landsat OLI/TIRS Level-2 data products as they are already corrected to surface reflectance. These data will be sufficient for time-series analysis as well as work involving multiple scenes. However, make sure any older imagery (i.e. earlier than Landsat 8 OLI/TIRS) is also corrected to surface ...
For clarification, filterBounds does not clip images, it only removes all images not intersect with the bound. If you want to clip every image in your collection, just use map and clip function as suggested in the first snippet of code suggested by Mathieu Gravey.
For exporting all NDVI images (there should be 71 of them, I guess), Mathieu Gravey's answer ...
if you want to add the mean to each object of the collection you can do:
And if you want the global mean (mean of means) :
var globalMean= colWithMeans.aggregate_mean('requestedMean')
to clip your collection you can do something like that
but if it's to export all the images you can simply try something like that:
Where you specify your geometry for the ...