26 Downloads (Pure)

Abstract

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
Volume1
Issue number3
Early online date23 Sep 2015
DOIs
Publication statusPublished - Sep 2015

Keywords

  • stat.ME

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

Cite this