[論文レビュー] 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$ による軌道を、独立同分布のノイズあり観測から回復するのに必要な最小標本数を特定すること。
• 軌道回復の統計的難易度と群 $G$ の代数的構造との間の明確な対応関係を確立すること。
• 標本複雑度の境界に不可欠な重要な代数的不変量を検証する計算的に効率的な手続きを開発すること。
• 理論枠組みを低温電子顕微鏡（cryo-EM）に適用し、アビニシオ再構成における画像数の $\sigma^6$ スケーリング則を検証すること。
• 複数の分子構造が同一のサンプルに存在する不均一な低温電子顕微鏡（heterogeneous cryo-EM）へ理論を拡張すること。

提案手法

• 信号 $\theta \in \mathbb{R}^p$ の $G$-軌道を推定する問題を、$g \sim \text{Haar}(G)$ および $\xi \sim \mathcal{N}(0, \sigma^2 I)$ のもとで、独立同分布の観測 $y = g \cdot \theta + \xi$ から定式化する。
• フーリエ解析と特性関数を用いて、異なる軌道が異なる分布をもたらすことを証明し、同定可能性を保証する。
• 不変理論の技術、特に実代数的ゼロ点定理（real Nullstellensatz）と実根の属性テストを用いて、不変量が軌道を分離する条件を特徴づける。
• 半正定値計画法と記号計算（Reid, Wang, および Wolkowicz のアルゴリズムを用いて）を用い、有限個の不変量が軌道を分離することを検証する。
• 3次以下の不変特徴量に基づく、アビニシオ再構成のための新規アルゴリズムを構築し、$\sigma^6$ の標本複雑度境界と一致させる。
• 分離多様体（separating variety）とヒルベルトの零点多様体（Hilbert nullcone）を用いて、軌道分離問題を多項式環における理想属性テストに還元する。

実験結果

リサーチクエスチョン

• RQ1 ノイズのある独立同分布観測のもとで、コンpakト群作用下の信号の軌道を回復するための最適な標本複雑度は何か？
• RQ2 標本複雑度は、群の構造とノイズ分散 $\sigma^2$ にどのように依存するか？
• RQ3 群作用の代数的不変量を用いて、有限個の不変量が軌道を分離することを検証できるか？
• RQ4 提案された理論枠組みのもとで、低温電子顕微鏡（cryo-EM）における画像数の $\sigma^6$ スケーリング則は成立するか？
• RQ5 複数の分子構造が混合したサンプルにおいて、理論を不均一な低温電子顕微鏡（heterogeneous cryo-EM）に拡張できるか？

主な発見

• 高ノイズ領域（$\sigma \to \infty$）において、軌道回復の標本複雑度は $\sigma^6$ に比例し、上限と下限が一致する。
• 著者らは、ハア測度の右不変性のおかげで、$\theta \in \mathbb{R}^p$ の軌道が $y = g \cdot \theta + \xi$ の分布から一意に同定可能であることを証明した。
• 半正定値計画法による実根の属性テストを用いて、計算的に効率的な検証手順を開発し、有限個の不変量が軌道を分離することを確認した。
• 本稿は、3次以下の不変特徴量に基づく、低温電子顕微鏡（cryo-EM）におけるアビニシオ再構成のための新規アルゴリズムを提供し、$\sigma^6$ の標本複雑度境界を達成した。
• この枠組みは不均一な低温電子顕微鏡（heterogeneous cryo-EM）にも適用可能であり、1つの混合サンプルから複数の分子構造を回復可能である。
• 理論的結果は実代数的ゼロ点定理と分離多様体の幾何学に基づいており、統計的推定と代数的不変量を結びつける。