BZPEERPapers, people, evidence and research pathsDiscover
Research tools
Network
Opportunities
Account
Source author record
Salem Said appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.
Catalog footprint
Research graph
Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
Hidden Markov models with observations in a Euclidean space play an important role in signal and image processing. Previous work extending to models where observations lie in Riemannian manifolds based on the Baum-Welch algorithm suffered from high memory usage and slow speed. Here we present an algorithm that is online, more accurate, and offers dramatic improvements in speed and efficiency.
Hidden Markov chain, or Markov field, models, with observations in a Euclidean space, play a major role across signal and image processing. The present work provides a statistical framework which can be used to extend these models, along with related, popular algorithms (such as the Baum-Welch algorithm), to the case where the observations lie in a Riemannian manifold. It is motivated by the potential use of hidden Markov chains and fields, with observations in Riemannian manifolds, as models for complex signals and images.
This paper studies fixed step-size stochastic approximation (SA) schemes, including stochastic gradient schemes, in a Riemannian framework. It is motivated by several applications, where geodesics can be computed explicitly, and their use accelerates crude Euclidean methods. A fixed step-size scheme defines a family of time-homogeneous Markov chains, parametrized by the step-size. Here, using this formulation, non-asymptotic performance bounds are derived, under Lyapunov conditions. Then, for any step-size, the corresponding Markov chain is proved to admit a unique stationary distribution, and to be geometrically ergodic. This result gives rise to a family of stationary distributions indexed by the step-size, which is further shown to converge to a Dirac measure, concentrated at the solution of the problem at hand, as the step-size goes to 0. Finally, the asymptotic rate of this convergence is established, through an asymptotic expansion of the bias, and a central limit theorem.
This entry contains the core material of my habilitation thesis, soon to be officially submitted. It provides a self-contained presentation of the original results in this thesis, in addition to their detailed proofs. The motivation of these results is the analysis of data which lie in Riemannian manifolds. Their aim is to bring about general, meaningful, and applicable tools, which can be used to model, and to learn from such "Riemannian data", as well as to analyse the various algorithms which may be required in this kind of pursuit (for sampling, optimisation, stochastic approximation, ...). The world of Riemannian data and algorithms can be quite different from its Euclidean counterpart, and this difference is the source of mathematical problems, addressed in my thesis. The first chapter provides some taylor-made geometric constructions, to be used in the thesis, while subsequent chapters (there are four more of them), address a series of issues, which arise from unresolved challenges, in the recent literature. A one-page guide, on how to read the thesis, is to be found right after the table of contents.
The Riemannian barycentre is one of the most widely used statistical descriptors for probability distributions on Riemannian manifolds. At present, existing algorithms are able to compute the Riemannian barycentre of a probability distribution, only if i.i.d. samples of this distribution are readily available. However, there are many cases where i.i.d. samples are quite difficult to obtain, and have to be replaced with non-independent samples, generated by a Markov chain Monte Carlo method. To overcome this difficulty, the present paper proposes a new Markov chain Monte Carlo algorithm for computing the Riemannian barycentre of a probability distribution on a Hadamard manifold (a simply connected, complete Riemannian manifold with non-positive curvature). This algorithm relies on two original propositions, proved in the paper. The first proposition states that the recursive barycentre of samples generated from a geometrically ergodic Markov chain converges in the mean-square to the Riemannian barycentre of the stationary distribution of this chain. The second proposition provides verifiable conditions which ensure a Metropolis-Hastings Markov chain, with its values in a symmetric Hadamard manifold, is geometrically ergodic. This latter result yields a partial solution, in the context of Riemannian manifolds, to the problem of geometric ergodicity of Metropolis-Hastings chains, which has previously attracted extensive attention when considered in Euclidean space. In addition to these two propositions, the new Markov chain Monte Carlo algorithm, proposed in this paper, is applied to a problem of Bayesian inference, arising from computer vision.
A new Riemannian geometry for the Compound Gaussian distribution is proposed. In particular, the Fisher information metric is obtained, along with corresponding geodesics and distance function. This new geometry is applied on a change detection problem on Multivariate Image Times Series: a recursive approach based on Riemannian optimization is developed. As shown on simulated data, it allows to reach optimal performance while being computationally more efficient.
Stochastic optimisation in Riemannian manifolds, especially the Riemannian stochastic gradient method, has attracted much recent attention. The present work applies stochastic optimisation to the task of recursive estimation of a statistical parameter which belongs to a Riemannian manifold. Roughly, this task amounts to stochastic minimisation of a statistical divergence function. The following problem is considered : how to obtain fast, asymptotically efficient, recursive estimates, using a Riemannian stochastic optimisation algorithm with decreasing step sizes? In solving this problem, several original results are introduced. First, without any convexity assumptions on the divergence function, it is proved that, with an adequate choice of step sizes, the algorithm computes recursive estimates which achieve a fast non-asymptotic rate of convergence. Second, the asymptotic normality of these recursive estimates is proved, by employing a novel linearisation technique. Third, it is proved that, when the Fisher information metric is used to guide the algorithm, these recursive estimates achieve an optimal asymptotic rate of convergence, in the sense that they become asymptotically efficient. These results, while relatively familiar in the Euclidean context, are here formulated and proved for the first time, in the Riemannian context. In addition, they are illustrated with a numerical application to the recursive estimation of elliptically contoured distributions.
This paper considers the problem of optimal filtering for partially observed signals taking values on the rotation group. More precisely, one or more components are considered not to be available in the measurement of the attitude of a 3D rigid body. In such cases, the observed signal takes its values on a Stiefel manifold. It is demonstrated how to filter the observed signal through the anti-development built from observations. A particle filter implementation is proposed to perform the estimation of the signal partially observed and corrupted by noise. The sampling issue is also addressed and interpolation methods are introduced. Illustration of the proposed technique on synthetic data demonstrates the ability of the approach to estimate the angular velocity of a partially observed 3D system partially observed.
Noncommutative harmonic analysis is used to solve a nonparametric estimation problem stated in terms of compound Poisson processes on compact Lie groups. This problem of decompounding is a generalization of a similar classical problem. The proposed solution is based on a char- acteristic function method. The treated problem is important to recent models of the physical inverse problem of multiple scattering.
We present a new model for the propagation of polarized light in a random birefringent medium. This model is based on a decomposition of the higher order statistics of the reduced Stokes parameters along the irreducible representations of the rotation group. We show how this model allows a detailed description of the propagation, giving analytical expressions for the probability densities of the Mueller matrix and the Stokes vector throughout the propagation. It also allows an exact description of the evolution of averaged quantities, such as the degree of polarization. We will also discuss how this model allows a generalization of the concepts of reduced Stokes parameters and degree of polarization to higher order statistics. We give some notes on how it can be extended to more general random media.