Adding in desaturation functionality

[wafels]

It is common for AIA images to show evidence of pixel saturation during flare events. Certain pixels receive too many photons and the charge 'bleeds over' into neighboring pixels, obscuring the details of the flaring regions.

It is possible to compensate for this, and compensation methods have been implemented in SSWIDL. For example, Schwartz et al (2015)

https://ui.adsabs.harvard.edu/abs/2015A%26C....13..117S/abstract

Citations to this paper reveal subsequent approaches:

https://ui.adsabs.harvard.edu/abs/2015A%26C....13..117S/citations

A simple desaturation method implemented in Python would allow aiapy users to at least get an idea of the flaring regions.