Multi-wavelength solar images in the EUV are routinely used for analysing solar features such as coronal holes, filaments, and flares. However, images taken in different bands often look remarkably similar as each band receives contributions coming from regions with a range of different temperatures. This has motivated the search for empirical techniques that may unmix these contributions and concentrate salient morphological features of the corona in a smaller set of less redundant source images. Blind Source Separation (BSS) precisely does this. Here we show how this novel concept also provides new insight into the physics of the solar corona, using observations made by SDO/AIA. The source images are extracted using a Bayesian positive source separation technique. We show how observations made in six spectral bands, corresponding to optically thin emissions, can be reconstructed by linear combination of three sources. These sources have a narrower temperature response and allow for considerable data reduction since the pertinent information from all six bands can be condensed in only one single composite picture. In addition, they give access to empirical temperature maps of the corona. The limitations of the BSS technique and some applications are briefly discussed.

Bayesian Positive Source Separation (BPSS) is a useful unsupervised approach for hyperspectral data unmixing, where numerical non-negativity of spectra and abundances has to be ensured, such in remote sensing. Moreover, it is sensible to impose a sum-to-one (full additivity) constraint to the estimated source abundances in each pixel. Even though non-negativity and full additivity are two necessary properties to get physically interpretable results, the use of BPSS algorithms has been so far limited by high computation time and large memory requirements due to the Markov chain Monte Carlo calculations. An implementation strategy which allows one to apply these algorithms on a full hyperspectral image, as typical in Earth and Planetary Science, is introduced. Effects of pixel selection, the impact of such sampling on the relevance of the estimated component spectra and abundance maps, as well as on the computation times, are discussed. For that purpose, two different dataset have been used: a synthetic one and a real hyperspectral image from Mars.

Many studies assume that the solar irradiance in the EUV can be decomposed into different contributions, which makes the modelling of the spectral variability considerably easier. We consider a different approach, in which these contributions are not imposed a priori but are effectively and robustly inferred from spectral irradiance measurements. This is a source separation problem with a positivity constraint, for which we use a Bayesian solution. Using five years of daily EUV spectra recorded by the TIMED/SEE satellite, we show that the spectral irradiance can be decomposed into three elementary spectra. Our results suggest that they describe different layers of the solar atmosphere rather than specific regions. The temporal variability of these spectra is discussed.