A method performed by an electronic device, comprising: obtaining at least two images of at least two fluorophores marking different biomolecules, respectively; and unmixing the obtained images into images, each containing the signals of a single fluorophore, while iteratively minimizing mutual information shared by the obtained images; wherein the obtaining of at least two images of at least two fluorophores comprises obtaining two images of two fluorophores, and wherein the unmixing of the obtained images comprises obtaining two new images updated from the obtained images by applying a variable for minimizing the mutual information calculated with respect to the obtained images, whereby images, each containing the signals of a single fluorophore, are obtained; wherein the obtaining of two new images updated from the obtained images comprises processing the obtained images as low-resolution images having lower resolution than the obtained images, respectively, as in Equations (i) below; v 1 = pixelBinning u 1, k v 2 = pixelBinning u 2, k wherein the u 1 and the u 2 indicate the obtained images, respectively, the v 1 and the v 2 indicate the low-resolution images, respectively, and the k indicates a factor applied in order to process the obtained images as the low-resolution images, respectively, calculating values α̂ and β̂ of variables α and β respectively from the low-resolution images as in Equations (iii) below by using loss functions L 1 (α) and L 2 (β) defined as in Equations (ii) below; L 1 α = I v 1 − α × v 2; v 2 L 2 β = I v 1; v 2 − β × v 1. α ^ = arg min α L 1 α = arg min α I v 1 − α × v 2; v 2 β ^ = arg min β L 1 β = arg min β I v 1; v 2 − β × v 1 wherein α̂ and β̂ indicate the values of variables α and β which minimise the values of L 1 (α) and L 2 (β) respectively; whereby I(X;Y) represents the mutual information shared by two images X and Y; and obtaining updated new images from the obtained images by applying the values of variables α̂ and β̂ and an update ratio λ, having a value between 0 and 1, to the obtained images as in Equations (iv) below, u 1 ← u 1 − λ α ^ u 2 u 2 ← u 2 − λ β ^ u 1. obtaining at least two images of at least two fluorophores marking different biomolecules, respectively; and unmixing the obtained images into images, each containing the signals of a single fluorophore, while iteratively minimizing mutual information shared by the obtained images; wherein the obtaining of at least two images of at least two fluorophores comprises obtaining two images of two fluorophores, and wherein the unmixing of the obtained images comprises obtaining two new images updated from the obtained images by applying a variable for minimizing the mutual information calculated with respect to the obtained images, whereby images, each containing the signals of a single fluorophore, are obtained; wherein the obtaining of two new images updated from the obtained images comprises processing the obtained images as low-resolution images having lower resolution than the obtained images, respectively, as in Equations (i) below; v 1 = pixelBinning u 1, k v 2 = pixelBinning u 2, k wherein the u 1 and the u 2 indicate the obtained images, respectively, the v 1 and the v 2 indicate the low-resolution images, respectively, and the k indicates a factor applied in order to process the obtained images as the low-resolution images, respectively, calculating values α̂ and β̂ of variables α and β respectively from the low-resolution images as in Equations (iii) below by using loss functions L 1 (α) and L 2 (β) defined as in Equations (ii) below; L 1 α = I v 1 − α × v 2; v 2 L 2 β = I v 1; v 2 − β × v 1. α ^ = arg min α L 1 α = arg min α I v 1 − α × v 2; v 2 β ^ = arg min β L 1 β = arg min β I v 1; v 2 − β × v 1 wherein α̂ and β̂ indicate the values of variables α and β which minimise the values of L 1 (α) and L 2 (β) respectively; whereby I(X;Y) represents the mutual information shared by two images X and Y; and obtaining updated new images from the obtained images by applying the values of variables α̂ and β̂ and an update ratio λ, having a value between 0 and 1, to the obtained images as in Equations (iv) below, u 1 ← u 1 − λ α ^ u 2 u 2 ← u 2 − λ β ^ u 1.