On approximating the non-negative rank: Applications to unsupervised image reduction