A Bayesian approach for polarimetric data reduction: the Mueller imaging case.

In this paper, we extend to the Mueller imaging framework a formerly introduced Bayesian approach dealing with polarimetric data reduction and robust clustering of polarization encoded images in the piecewise constant case. The extension was made possible thanks to a suitable writing of the observation model in the Mueller context that relies on the system's coherency matrix and Cholesky decomposition such that the admissibility constraints are easily captured. This generalization comes at the cost of nonlinearity with respect to the parameters that have to be estimated. This estimation-clustering problem is tackled in a Bayesian framework where a hierarchical stochastic model based on a Markov random field proposed by Potts is used. This fully unsupervised approach is extensively tested over synthetic data as well as real Mueller images.