BZPEERReading, trust and collaboration for researchDiscover
Workspaces
Network
Opportunities
Account
Researcher profile
Salem Said contributes to research discovery and scholarly infrastructure.
Trust snapshot
Actions
Identity and collaboration
Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.
Log in to claimDirect collaboration
Claim this author entity first to unlock direct invitations.
Research graph
Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph 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.