In hyperspectral image unmixing, a collection of pure spectra, the so-called endmembers, is identified and their abundance fractions are estimated at each pixel. While endmembers are often extracted using a geometric approach, the abundances are usually estimated using a least-squares approach by solving an inverse problem. In this paper, we tackle the problem of abundance estimation by using a...