Diffusion bridges for stochastic Hamiltonian systems with applications to shape analysis

Stochastically evolving geometric systems are studied in geometric mechanics for modelling turbulence parts of multi-scale fluid flows and in shape analysis for stochastic evolutions of shapes of, e.g., human organs. Recently introduced models involve stochastic differential equations that govern the dynamics of a diffusion process $X$. In applications $X$ is only partially observed at times $0$ and $T>0$. Conditional on these observations, interest lies in inferring parameters in the dynamics of the diffusion and reconstructing the path $(X_t,\, t\in [0,T])$. The latter problem is known as bridge simulation. We develop a general scheme for bridge sampling in the case of finite dimensional systems of shape landmarks and singular solutions in fluid dynamics. This scheme allows for subsequent statistical inference of properties of the fluid flow or the evolution of observed shapes. It covers stochastic landmark models for which no suitable simulation method has been proposed in the literature, that removes restrictions of earlier approaches, improves the handling of the nonlinearity of the configuration space leading to more effective sampling schemes and allows to generalise the common inexact matching scheme to the stochastic setting.

[1]  E. Hairer,et al.  Geometric Numerical Integration: Structure Preserving Algorithms for Ordinary Differential Equations , 2004 .

[2]  Darryl D. Holm,et al.  Universitet A Geometric Framework for Stochastic Shape Analysis , 2018 .

[3]  Bo Markussen,et al.  Large deformation diffeomorphisms with application to optic flow , 2007, Comput. Vis. Image Underst..

[4]  Franccois-Xavier Vialard Extension to Infinite Dimensions of a Stochastic Second-Order Model associated with the Shape Splines , 2010, 1003.3957.

[5]  Bjarne Kjær Ersbøll,et al.  Extending and Applying Active Appearance Models for Automated, High Precision Segmentation* , 2001 .

[6]  Gareth O. Roberts,et al.  Importance sampling techniques for estimation of diffusion models , 2009 .

[7]  Miles Lubin,et al.  Forward-Mode Automatic Differentiation in Julia , 2016, ArXiv.

[8]  L. Younes Shapes and Diffeomorphisms , 2010 .

[9]  Hamiltonian Monte Carlo on Symmetric and Homogeneous Spaces via Symplectic Reduction , 2019, 1903.02699.

[10]  M. Yor,et al.  Grossissements de filtrations: exemples et applications , 1985 .

[11]  P. Fearnhead,et al.  Exact and computationally efficient likelihood‐based estimation for discretely observed diffusion processes (with discussion) , 2006 .

[12]  Akshay Pai,et al.  A Stochastic Large Deformation Model for Computational Anatomy , 2016, IPMI.

[13]  Darryl D. Holm,et al.  The Euler–Poincaré Equations and Semidirect Products with Applications to Continuum Theories , 1998, chao-dyn/9801015.

[14]  Franccois-Xavier Vialard,et al.  Shape Splines and Stochastic Shape Evolutions: A Second Order Point of View , 2010, 1003.3895.

[15]  Darryl D. Holm,et al.  An integrable shallow water equation with peaked solitons. , 1993, Physical review letters.

[16]  Mogens Bladt,et al.  Simulation of multivariate diffusion bridges , 2014 .

[17]  Darryl D. Holm,et al.  Stochastic discrete Hamiltonian variational integrators , 2016, BIT Numerical Mathematics.

[18]  G. Roberts,et al.  On inference for partially observed nonlinear diffusion models using the Metropolis–Hastings algorithm , 2001 .

[19]  Moritz Schauer,et al.  Continuous-discrete smoothing of diffusions , 2017, Electronic Journal of Statistics.

[20]  Stefan Sommer,et al.  Computational Anatomy in Theano , 2017, GRAIL/MFCA/MICGen@MICCAI.

[21]  L. Rogers,et al.  Diffusions, Markov processes, and martingales , 1979 .

[22]  Jean-Louis Marchand,et al.  Conditioning diffusions with respect to partial observations , 2011, 1105.1608.

[23]  F. Meulen,et al.  Bayesian estimation of discretely observed multi-dimensional diffusion processes using guided proposals , 2014, Electronic Journal of Statistics.

[24]  B. Delyon,et al.  Simulation of conditioned diffusion and application to parameter estimation , 2006 .

[25]  Michael I. Miller,et al.  Landmark matching via large deformation diffeomorphisms , 2000, IEEE Trans. Image Process..

[26]  Gareth O. Roberts,et al.  Examples of Adaptive MCMC , 2009 .

[27]  Darryl D. Holm,et al.  Momentum maps and measure-valued solutions (peakons, filaments, and sheets) for the EPDiff equation , 2003, nlin/0312048.

[28]  Harry van Zanten,et al.  Guided proposals for simulating multi-dimensional diffusion bridges , 2013, 1311.3606.

[29]  Fabrice Baudoin,et al.  Conditioned stochastic differential equations: theory, examples and application to finance , 2002 .

[30]  Stefan Sommer,et al.  Bridge Simulation and Metric Estimation on Landmark Manifolds , 2017, GRAIL/MFCA/MICGen@MICCAI.

[31]  Mark A. Girolami,et al.  Geometry and Dynamics for Markov Chain Monte Carlo , 2017, ArXiv.

[32]  Michael I. Miller,et al.  Evolutions equations in computational anatomy , 2009, NeuroImage.

[33]  M. Girolami,et al.  Riemann manifold Langevin and Hamiltonian Monte Carlo methods , 2011, Journal of the Royal Statistical Society: Series B (Statistical Methodology).

[34]  Jialin Hong,et al.  Explicit pseudo-symplectic methods for stochastic Hamiltonian systems , 2018 .

[35]  Darryl D. Holm Variational principles for stochastic fluid dynamics , 2014, Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[36]  Stephen R. Marsland,et al.  Langevin Equations for Landmark Image Registration with Uncertainty , 2016, SIAM J. Imaging Sci..

[37]  Darryl D. Holm,et al.  Soliton dynamics in computational anatomy , 2004, NeuroImage.

[38]  Evelyn Buckwar,et al.  Spectral density-based and measure-preserving ABC for partially observed diffusion processes. An illustration on Hamiltonian SDEs , 2019, Statistics and Computing.

[39]  G. N. Milstein,et al.  Numerical Methods for Stochastic Systems Preserving Symplectic Structure , 2002, SIAM J. Numer. Anal..

[40]  LINE KÜHNEL,et al.  Differential geometry and stochastic dynamics with deep learning numerics , 2017, Appl. Math. Comput..

[41]  Martin Bauer,et al.  Overview of the Geometries of Shape Spaces and Diffeomorphism Groups , 2013, Journal of Mathematical Imaging and Vision.

[42]  D. Kendall SHAPE MANIFOLDS, PROCRUSTEAN METRICS, AND COMPLEX PROJECTIVE SPACES , 1984 .