Uncertainty Quantification for LDDMM Using a Low-Rank Hessian Approximation
نویسندگان
چکیده
This paper presents an approach to estimate the uncertainty of registration parameters for the large displacement diffeomorphic metric mapping (LDDMM) registration framework. Assuming a local multivariate Gaussian distribution as an approximation for the registration energy at the optimal registration parameters, we propose a method to approximate the covariance matrix as the inverse of the Hessian of the registration energy to quantify registration uncertainty. In particular, we make use of a low-rank approximation to the Hessian to accurately and efficiently estimate the covariance matrix using few eigenvalues and eigenvectors. We evaluate the uncertainty of the LDDMM registration results for both synthetic and real imaging data.
منابع مشابه
Analysis of the Hessian for Inverse Scattering Problems. Part II: Inverse Medium Scattering of Acoustic Waves
We address the inverse problem for scattering of acoustic waves due to an inhomogeneous medium. We derive and analyze the Hessian in both Hölder and Sobolev spaces. Using an integral equation approach based on Newton potential theory and compact embeddings in Hölder and Sobolev spaces, we show that the Hessian can be decomposed into two components, both of which are shown to be compact operator...
متن کاملUncertainty quantification of reservoir performance predictions using a stochastic optimization algorithm
The simultaneous perturbation stochastic approximation (SPSA) algorithm is modified to obtain a stochastic Gaussian search direction (SGSD) algorithm for automatic history matching. The search direction in the SGSD algorithm is obtained by simultaneously perturbing the reservoir model with unconditional realizations from a Gaussian distribution. This search direction has two nice properties: (1...
متن کاملAnalysis of the Hessian for inverse scattering problems: I. Inverse shape scattering of acoustic waves
We derive expressions for the shape Hessian operator of the data misfit functional corresponding to the inverse problem of inferring the shape of a scatterer from reflected acoustic waves, using a Banach space setting and the Lagrangian approach. The shape Hessian is then analyzed in both Hölder and Sobolev spaces. Using an integral equation approach and compact embeddings in Hölder and Sobolev...
متن کاملAnalysis of the Hessian for Inverse Scattering Problems. Part I: Inverse Shape Scattering of Acoustic Waves
We derive expressions for the shape Hessian operator of the data misfit functional corresponding to the inverse problem of inferring the shape of a scatterer from reflected acoustic waves, using a Banach space setting and the Lagrangian approach. The shape Hessian is then analyzed in both Hölder and Sobolev spaces. Using an integral equation approach and compact embeddings in Hölder and Sobolev...
متن کاملAnalysis of the Hessian for Inverse Scattering Problems. Part Iii: Inverse Medium Scattering of Electromagnetic Waves in Three Dimensions
Continuing our previous work [6, Inverse Problems, 2012, 28, 055002] and [5, Inverse Problems, 2012, 28, 055001], we address the ill-posedness of the inverse scattering problem of electromagnetic waves due to an inhomogeneous medium by studying the Hessian of the data misfit. We derive and analyze the Hessian in both Hölder and Sobolev spaces. Using an integral equation approach based on Newton...
متن کامل