Ptyhon: how to compute mean and percentiles over a list of arrays avoiding zeros

Question

I am working with satellite images that are geotiff files in form of arrays. Let say that I have multiple images with blank spaces (because of clouds or other eleme