Abderrahim Halimi

Nonlinear Unmixing of Hyperspectral Images Using a Generalized Bilinear Model

Nonlinear models have recently shown interesting properties for spectral unmixing. This paper studies a generalized bilinear model and a hierarchical Bayesian algorithm for unmixing hyperspectral images. The proposed model is a generalization not only of the accepted linear mixing model but also of a bilinear model that has been recently introduced in the literature. Appropriate priors are chosen for its parameters to satisfy the positivity and sum-to-one constraints for the abundances. The joint posterior distribution of the unknown parameter vector is then derived. Unfortunately, this posterior is too complex to obtain analytical expressions of the standard Bayesian estimators. As a consequence, a Metropolis-within-Gibbs algorithm is proposed, which allows samples distributed according to this posterior to be generated and to estimate the unknown model parameters. The performance of the resulting unmixing strategy is evaluated via simulations conducted on synthetic and real data.

Supervised Nonlinear Spectral Unmixing Using a Postnonlinear Mixing Model for Hyperspectral Imagery

This paper presents a nonlinear mixing model for hyperspectral image unmixing. The proposed model assumes that the pixel reflectances are nonlinear functions of pure spectral components contaminated by an additive white Gaussian noise. These nonlinear functions are approximated using polynomial functions leading to a polynomial postnonlinear mixing model. A Bayesian algorithm and optimization methods are proposed to estimate the parameters involved in the model. The performance of the unmixing strategies is evaluated by simulations conducted on synthetic and real data.

Unmixing hyperspectral images using the generalized bilinear model

Nonlinear models have recently shown interesting properties for spectral unmixing. This paper considers a generalized bilinear model recently introduced for unmixing hyperspectral images. Different algorithms are studied to estimate the parameters of this bilinear model. The positivity and sum-to-one constraints for the abundances are ensured by the proposed algorithms. The performance of the resulting unmixing strategy is evaluated via simulations conducted on synthetic and real data.

Unsupervised Unmixing of Hyperspectral Images Accounting for Endmember Variability

This paper presents an unsupervised Bayesian algorithm for hyperspectral image unmixing, accounting for endmember variability. The pixels are modeled by a linear combination of endmembers weighted by their corresponding abundances. However, the endmembers are assumed random to consider their variability in the image. An additive noise is also considered in the proposed model, generalizing the normal compositional model. The proposed algorithm exploits the whole image to benefit from both spectral and spatial information. It estimates both the mean and the covariance matrix of each endmember in the image. This allows the behavior of each material to be analyzed and its variability to be quantified in the scene. A spatial segmentation is also obtained based on the estimated abundances. In order to estimate the parameters associated with the proposed Bayesian model, we propose to use a Hamiltonian Monte Carlo algorithm. The performance of the resulting unmixing strategy is evaluated through simulations conducted on both synthetic and real data.

Parameter Estimation for Peaky Altimetric Waveforms

Much attention has been recently devoted to the analysis of coastal altimetric waveforms. When approaching the coast, altimetric waveforms are sometimes corrupted by peaks caused by high reflective areas inside the illuminated land surfaces or by the modification of the sea state close to the shoreline. This paper introduces a new parametric model for these peaky altimetric waveforms. This model assumes that the received altimetric waveform is the sum of a Brown echo and an asymmetric Gaussian peak. The asymmetric Gaussian peak is parameterized by a location, an amplitude, a width, and an asymmetry coefficient. A maximum-likelihood estimator is studied to estimate the Brown plus peak model parameters. The Cramér-Rao lower bounds of the model parameters are then derived providing minimum variances for any unbiased estimator, i.e., a reference in terms of estimation error. The performance of the proposed model and the resulting estimation strategy are evaluated via many simulations conducted on synthetic and real data. Results obtained in this paper show that the proposed model can be used to retrack efficiently standard oceanic Brown echoes as well as coastal echoes corrupted by symmetric or asymmetric Gaussian peaks. Thus, the Brown with Gaussian peak model is useful for analyzing altimetric measurements closer to the coast.

A Semi-Analytical Model for Delay/Doppler Altimetry and Its Estimation Algorithm

The concept of delay/Doppler (DD) altimetry (DDA) has been under study since the mid-1990s, aiming at reducing the measurement noise and increasing the along-track resolution in comparison with the conventional pulse-limited altimetry. This paper introduces a new model for the mean backscattered power waveform acquired by a radar altimeter operating in synthetic aperture radar mode, as well as an associated least squares (LS) estimation algorithm. As in conventional altimetry (CA), the mean power can be expressed as the convolution of three terms: the flat surface impulse response (FSIR), the probability density function of the heights of the specular scatterers, and the time/frequency point target response of the radar. An important contribution of this paper is to derive an analytical formula for the FSIR associated with DDA. This analytical formula is obtained for a circular antenna pattern, no mispointing, no vertical speed effect, and a uniform scattering. The double convolution defining the mean echo power can then be computed numerically, resulting in a 2-D semi-analytical model called the DD map (DDM). This DDM depends on three altimetric parameters: the epoch, the sea surface wave height, and the amplitude. A multi-look model is obtained by summing all the reflected echoes from the same along-track surface location of interest after applying appropriate delay compensation (range migration) to align the DDM on the same reference. The second contribution of this paper concerns the estimation of the parameters associated with the multi-look semi-analytical model. An LS approach is investigated by means of the Levenberg-Marquardt algorithm. Simulations conducted on simulated altimetric waveforms allow the performance of the proposed estimation algorithm to be appreciated. The analysis of Cryosat-2 waveforms shows an improvement in parameter estimation when compared to the CA.

Hyperspectral Unmixing in Presence of Endmember Variability, Nonlinearity, or Mismodeling Effects

This paper presents three hyperspectral mixture models jointly with Bayesian algorithms for supervised hyperspectral unmixing. Based on the residual component analysis model, the proposed general formulation assumes the linear model to be corrupted by an additive term whose expression can be adapted to account for nonlinearities (NLs), endmember variability (EV), or mismodeling effects (MEs). The NL effect is introduced by considering a polynomial expression that is related to bilinear models. The proposed new formulation of EV accounts for shape and scale endmember changes while enforcing a smooth spectral/spatial variation. The ME formulation considers the effect of outliers and copes with some types of EV and NL. The known constraints on the parameter of each observation model are modeled via suitable priors. The posterior distribution associated with each Bayesian model is optimized using a coordinate descent algorithm, which allows the computation of the maximum a posteriori estimator of the unknown model parameters. The proposed mixture and Bayesian models and their estimation algorithms are validated on both synthetic and real images showing competitive results regarding the quality of the inferences and the computational complexity, when compared with the state-of-the-art algorithms.

Including Antenna Mispointing in a Semi-Analytical Model for Delay/Doppler Altimetry

Delay/Doppler altimetry (DDA) aims at reducing the measurement noise and increasing the along-track resolution in comparison with conventional pulse-limited altimetry. In a previous paper, we have proposed a semi-analytical model for DDA, which considers some simplifications as the absence of mispointing antenna. This paper first proposes a new analytical expression for the flat surface impulse response (FSIR), considering antenna mispointing angles, a circular antenna pattern, no vertical speed effect, and uniform scattering. The 2-D delay/Doppler map is then obtained by a numerical computation of the convolution between the proposed analytical function, the probability density function of the heights of the specular scatterers, and the time/frequency point target response of the radar. The approximations used to obtain the semi-analytical model are analyzed, and the associated errors are quantified by analytical bounds for these errors. The second contribution of this paper concerns the estimation of the parameters associated with the multilook semi-analytical model. Two estimation strategies based on the least squares procedure are proposed. The proposed model and algorithms are validated on both synthetic and real waveforms. The obtained results are very promising and show the accuracy of this generalized model with respect to the previous model assuming zero antenna mispointing.

Bayesian Estimation of Smooth Altimetric Parameters: Application to Conventional and Delay/Doppler Altimetry

This paper proposes a new Bayesian strategy for the smooth estimation of altimetric parameters. The altimetric signal is assumed to be corrupted by a thermal and speckle noise distributed according to an independent and non-identically Gaussian distribution. We introduce a prior enforcing a smooth temporal evolution of the altimetric parameters which improves their physical interpretation. The posterior distribution of the resulting model is optimized using a gradient descent algorithm which allows us to compute the maximum a posteriori estimator of the unknown model parameters. This algorithm has a low computational cost that is suitable for real-time applications. The proposed Bayesian strategy and the corresponding estimation algorithm are evaluated using both synthetic and real data associated with conventional and delay/Doppler altimetry. The analysis of real Jason-2 and CryoSat-2 waveforms shows an improvement in parameter estimation when compared to state-of-the-art estimation algorithms.

Supervised nonlinear spectral unmixing using a polynomial post nonlinear model for hyperspectral imagery

This paper studies a hierarchical Bayesian model for nonlinear hyperspectral image unmixing. The proposed model assumes that the pixel reflectances are polynomial functions of linear mixtures of pure spectral components contaminated by an additive white Gaussian noise. The parameters involved in this model satisfy constraints that are naturally expressed within a Bayesian framework. A Gibbs sampler allows one to sample the unknown abundances and nonlinearity parameters according to the joint posterior of interest. The performance of the resulting unmixing strategy is evaluated thanks to simulations conducted on synthetic and real data.