I have a raster stack (numpy 3d array) of daily MODIS LST quality bands. The quality bands have pixel values that correspond to the quality of the LST data in each pixel. For example, 2 = no data received due to clouds, 3 = possible cloud contamination, 17 = low quality for other reason, etc.
I would like to create a line plot (x axis = proportion of pixels, y axis = day of year) that depicts the proportion of pixels within the MODIS tile that have a given quality value on any day of the year. The line plot would have multiple lines, each representing a value in the quality band. I would like to do this in Python, if possible.
As an example, lets say I have a 2d array, which would represent a MODIS tile for one day, with the following information: [[2,3,0], [17,3,0], [2,2,0]]
For this array I would need to calculate the proportion of pixels that correspond to each value, and then have the output in some form (table or csv file) that I could then plot. I would then repeat this process for each day of the year.