We propose a weighted residual nonnegative matrix factorization (NMF) with spatial regularization to unmix hyperspectral data. NMF decomposes a matrix into the product of two nonnegative matrices. However, NMF is known to be generally sensitive to noise, which makes difficult to retrieve the global minimum of the underlying objective function. To overcome this limitation, we include a residual weighting mechanism in the conventional NMF formulation. This strategy treats each row of the residual based on the weighting factor. In this manner, residuals with large values are penalized less and residuals with small values are penalized more to make NMF based unmixing problem more robust. Furthermore, we include a weight term in the form of an $\ell_{1}$ norm regularizer to provide spatial information of the abundance matrix.
The model is described in the paper published in IEEE Geoscience Remote Sensing Letters:
The MATLAB code is available online: