Regularising orientation estimation in Cryo-EM 3D map refinement through measure-based lifting over Riemannian manifolds

Abstract

Motivated by the trade-off between noise-robustness and data-consistency for joint 3D map reconstruction and rotation estimation in single particle cryogenic-electron microscopy (Cryo-EM), we propose ellipsoidal support lifting (ESL), a measure-based lifting scheme for regularising and approximating the global minimiser of a smooth function over a Riemannian manifold. Under a uniqueness assumption on the minimiser we show several theoretical results, in particular well-posedness of the method and an error bound due to the induced bias with respect to the global minimiser. Additionally, we use the developed theory to integrate the measure-based lifting scheme into an alternating update method for joint homogeneous 3D map reconstruction and rotation estimation, where typically tens of thousands of manifold-valued minimisation problems have to be solved and where regularisation is necessary because of the high noise levels in the data. The joint recovery method is used to test both the theoretical predictions and algorithmic performance through numerical experiments with Cryo-EM data. In particular, the induced bias due to the regularising effect of ESL empirically estimates better rotations, i.e., rotations closer to the ground truth, than global optimisation would.

References in corpus (11)

• Single-particle cryo-electron microscopy: Mathematical theory, computational challenges, and opportunities
• A Second Order Non-Smooth Variational Model for Restoring Manifold-Valued Images
• Restoration of Manifold-Valued Images by Half-Quadratic Minimization
• Total variation regularization for manifold-valued data
• Consensus-Based Optimization on Hypersurfaces: Well-Posedness and Mean-Field Limit
• Method of moments for 3-D single particle ab initio modeling with non-uniform distribution of viewing angles
• A Parallel Douglas Rachford Algorithm for Minimizing ROF-like Functionals on Images with Values in Symmetric Hadamard Manifolds
• Total Generalized Variation for Manifold-valued Data
• Priors with Coupled First and Second Order Differences for Manifold-Valued Image Processing
• An Inexact Semi-smooth Newton Method on Riemannian Manifolds with Application to Duality-based Total Variation Denoising
• Formal uniqueness in Ewald sphere corrected single particle analysis