[论文解读] Estimation under group actions: recovering orbits from invariants

Afonso S. Bandeira, Ben Blum-Smith+4 more|arXiv (Cornell University)|Dec 29, 2017

Advanced Electron Microscopy Techniques and Applications被引用 10

一句话总结

本文通过不变量理论与实代数几何工具，建立了从独立同分布的噪声观测中恢复信号在群作用下的轨道所需样本复杂度的匹配上下界。研究证明，在冷冻电镜（cryo-EM）中，所需图像数量随噪声方差 $\sigma^2$ 的六次方增长，即 $\sigma^6$，并提出了一个计算高效的代数不变量认证程序。

ABSTRACT

We study a class of orbit recovery problems in which we observe independent copies of an unknown element of $\mathbb{R}^p$, each linearly acted upon by a random element of some group (such as $\mathbb{Z}/p$ or $\mathrm{SO}(3)$) and then corrupted by additive Gaussian noise. We prove matching upper and lower bounds on the number of samples required to approximately recover the group orbit of this unknown element with high probability. These bounds, based on quantitative techniques in invariant theory, give a precise correspondence between the statistical difficulty of the estimation problem and algebraic properties of the group. Furthermore, we give computer-assisted procedures to certify these properties that are computationally efficient in many cases of interest. The model is motivated by geometric problems in signal processing, computer vision, and structural biology, and applies to the reconstruction problem in cryo-electron microscopy (cryo-EM), a problem of significant practical interest. Our results allow us to verify (for a given problem size) that if cryo-EM images are corrupted by noise with variance $σ^2$, the number of images required to recover the molecule structure scales as $σ^6$. We match this bound with a novel (albeit computationally expensive) algorithm for ab initio reconstruction in cryo-EM, based on invariant features of degree at most 3. We further discuss how to recover multiple molecular structures from mixed (or heterogeneous) cryo-EM samples.

研究动机与目标

• 确定在独立同分布的噪声观测下，恢复 $\mathbb{R}^p$ 中未知信号在随机群作用下的群轨道所需的最少观测次数。
• 建立轨道恢复的统计难度与群 $G$ 的代数结构之间的精确对应关系。
• 开发计算高效的代数不变量认证程序，以验证对样本复杂度边界至关重要的代数不变量。
• 将该框架应用于冷冻电镜，验证在从头重建中图像数量遵循 $\sigma^6$ 的标度律。
• 将理论扩展至异质性冷冻电镜，其中单一样品中存在多种分子结构。

提出的方法

• 将轨道恢复问题形式化为：从独立同分布的观测 $y = g \cdot \theta + \xi$ 中估计信号 $\theta \in \mathbb{R}^p$ 的 $G$-轨道，其中 $g \sim \text{Haar}(G)$ 且 $\xi \sim \mathcal{N}(0, \sigma^2 I)$。
• 利用傅里叶分析与特征函数证明，不同轨道产生不同的分布，从而确保可识别性。
• 应用不变量理论中的技术——特别是实零点定理（real Nullstellensatz）与实根理想成员判定——以刻画不变量分离轨道的条件。
• 采用半定规划与符号计算（基于 Reid、Wang 与 Wolkowicz 的算法）来认证有限个不变量是否能分离轨道。
• 提出一种新颖的从头重建算法，基于次数不超过 3 的不变量特征，其样本复杂度边界与 $\sigma^6$ 完全匹配。
• 利用分离曲面（separating variety）与希尔伯特零点集（Hilbert nullcone）将轨道分离问题转化为多项式环中的理想成员判定问题。

实验结果

研究问题

• RQ1 在给定独立同分布的噪声观测下，恢复紧致群作用下信号轨道的最优样本复杂度是多少？
• RQ2 样本复杂度如何依赖于群结构与噪声方差 $\sigma^2$？
• RQ3 群作用的代数不变量能否用于认证有限个不变量可分离轨道？
• RQ4 在所提出的理论框架下，冷冻电镜图像数量的 $\sigma^6$ 标度律是否成立？
• RQ5 该理论能否扩展至异质性冷冻电镜，即样品中存在多种分子结构混合的情况？

主要发现

• 在高噪声区域（$\sigma \to \infty$），轨道恢复的样本复杂度与 $\sigma^6$ 成正比，且上下界完全匹配。
• 作者证明，由于哈尔测度的右不变性，信号 $\theta \in \mathbb{R}^p$ 的轨道可从 $y = g \cdot \theta + \xi$ 的分布中唯一识别。
• 通过半定规划实现的实根理想成员判定，开发了一种计算高效的认证程序，用于验证有限个不变量是否能分离轨道。
• 本文提出一种基于次数不超过 3 的不变量特征的新型从头重建算法，其样本复杂度边界与 $\sigma^6$ 完全匹配。
• 该框架适用于异质性冷冻电镜，可从单一混合样品中恢复多种分子结构。
• 理论结果基于实零点定理与分离曲面的几何性质，将统计估计与代数不变量紧密联系。