Techniques for removing image autoflourescence from fluorescently stained biological images are provided herein. The techniques utilize non-negative matrix factorization that may constrain mixing coefficients to be non-negative. The probability of convergence to local minima is reduced by using smoothness constraints. The non-negative matrix factorization algorithm provides the advantage of removing both dark current and autofluorescence.