Bayesian nonlinear hyperspectral unmixing with spatial residual component analysis

Research output: Contribution to journalArticlepeer-review

28 Citations (Scopus)
71 Downloads (Pure)


This paper presents a new Bayesian model and algorithm for nonlinear unmixing of hyperspectral images. The model proposed represents the pixel reflectances as linear combinations of the endmembers, corrupted by nonlinear (with respect to the endmembers) terms and additive Gaussian noise. Prior knowledge about the problem is embedded in a hierarchical model that describes the dependence structure between the model parameters and their constraints. In particular, a gamma Markov random field is used to model the joint distribution of the nonlinear terms, which are expected to exhibit significant spatial correlations. An adaptive Markov chain Monte Carlo algorithm is then proposed to compute the Bayesian estimates of interest and perform Bayesian inference. This algorithm is equipped with a stochastic optimisation adaptation mechanism that automatically adjusts the parameters of the gamma Markov random field by maximum marginal likelihood estimation. Finally, the proposed methodology is demonstrated through a series of experiments with comparisons using synthetic and real data and with competing state-of-the-art approaches.
Original languageEnglish
Pages (from-to)174 - 185
JournalIEEE Transactions on Computational Imaging
Issue number3
Early online date23 Sept 2015
Publication statusPublished - Sept 2015


  • stat.ME


Dive into the research topics of 'Bayesian nonlinear hyperspectral unmixing with spatial residual component analysis'. Together they form a unique fingerprint.

Cite this