2025-11-13T19:07:10.620379

A max filtering local stability theorem with application to weighted phase retrieval and cryo-EM

Qaddura

Given an inner product space $V$ and a group $G$ of linear isometries, max filtering offers a rich class of convex $G$-invariant maps. In this paper, we identify sufficient conditions under which these maps are locally bilipschitz on $R(G)$, the set of orbits with maximal dimension, with respect to the quotient metric on the orbit space $V/G$. Central to our proof is a desingularization theorem, which applies to open, dense neighborhoods around each orbit in $R(G)/G$ and may be of independent interest. As an application, we provide guarantees for stable weighted phase retrieval. That is, we construct componentwise convex bilipschitz embeddings of weighted complex (resp.\ quaternionic) projective spaces. These spaces arise as quotients of direct sums of nontrivial unitary irreducible complex (resp.\ quaternionic) representations of the group of unit complex numbers $S^1\cong \operatorname{SO}(2)$ (resp.\ unit quaternions $S^3\cong \operatorname{SU}(2)$). We also discuss the relevance of such embeddings to a nearest-neighbor problem in single-particle cryogenic electron microscopy (cryo-EM), a leading technique for resolving the spatial structure of biological molecules.

academic

A max filtering local stability theorem with application to weighted phase retrieval and cryo-EM

基本信息

论文ID: 2403.14042
标题: A max filtering local stability theorem with application to weighted phase retrieval and cryo-EM
作者: Yousef Qaddura (The Ohio State University)
分类: math.FA cs.IT math.IT
发表时间: 2024年3月 (arXiv预印本，v3版本于2025年10月13日更新)
论文链接: https://arxiv.org/abs/2403.14042

摘要

本文研究了在内积空间 $V$ 和线性等距群 $G$ 的框架下，max filtering映射的局部双Lipschitz性质。作者识别了使这些凸 $G$ -不变映射在正则点集合 $R(G)$ （具有最大维数的轨道集合）上关于商空间 $V/G$ 的商度量局部双Lipschitz的充分条件。证明的核心是一个去奇异化定理，该定理适用于 $R(G)/G$ 中每个轨道周围的开稠密邻域。作为应用，论文为稳定加权相位恢复提供了保证，构造了加权复（四元数）射影空间的分量凸双Lipschitz嵌入，并讨论了这些嵌入在单粒子冷冻电镜（cryo-EM）最近邻问题中的相关性。

研究背景与动机

核心问题

现代机器学习算法通常设计用于欧几里得数据，但许多实际数据表示存在由正交对称群 $G \leq O(V)$ 引起的模糊性。例如：

冷冻电镜数据可能存在于有限维复向量空间 $\mathbb{C}^d$ 中，受对角圆作用 $S^1 \to \mathbb{C}^{d \times d}$ 诱导的模糊性影响
相位恢复问题中的复数等价关系 $x \sim e^{i\theta}x$

研究意义

为了利用基于欧几里得的机器学习方法，需要将轨道空间 $V/G$ 以双Lipschitz方式嵌入到欧几里得空间中。这种嵌入确保 $V/G$ 中的距离得到忠实保持，使得欧几里得算法能够稳健地转移到轨道空间。

现有方法局限性

对于有限群 $G$ ，已知每个单射max filter bank都是双Lipschitz的
对于无限群，只有三个例外情况被解决：复相位恢复、极坐标作用
一般无限群的双Lipschitz性质仍是开放问题

研究动机

本文旨在研究在给定足够多通用模板的条件下，max filter banks何时是双Lipschitz的，特别是在所有非零轨道具有常数维数的群作用情况下。

核心贡献

建立了max filter banks的局部双Lipschitz条件：在正则点集合 $R(G)$ 上，当模板数量超过 $2 \cdot \chi(G) \cdot (c-1)$ 时，通用max filter banks是局部双Lipschitz的
提出了去奇异化定理：适用于 $R(G)/G$ 中每个轨道周围的开稠密邻域，可能具有独立的数学价值
构造了稳定加权相位恢复的双Lipschitz嵌入：为加权复/四元数射影空间提供了分量凸双Lipschitz嵌入
发展了Voronoi单元分解理论：提供了主要点和正则点的几何刻画，建立了详细的Voronoi分解理论
应用于冷冻电镜：为cryo-EM中的最近邻问题提供了理论保证，改进了现有的双谱嵌入方法

方法详解

任务定义

给定内积空间 $V$ 和紧致群 $G \leq O(V)$ ，寻找模板 $z_1, \ldots, z_n \in V$ 使得max filter bank $\Phi([x]) := \{\langle\langle[x], [z_i]\rangle\rangle\}_{i=1}^n$ 是双Lipschitz映射，其中max filtering映射定义为： $\langle\langle[x], [z]\rangle\rangle := \sup_{p \in [x], q \in [z]} \langle p, q \rangle$

核心概念

正则Voronoi复杂度

对于紧致群 $G \leq O(d)$ ，定义：

正则点集合： $R(G) := \{x \in \mathbb{R}^d : \dim([x]) = \max_{y \in \mathbb{R}^d} \dim([y])\}$
正则Voronoi复杂度： $\chi(G) := \max_{x,p \in R(G)} \{|G_x/G_p| : G_p \leq G_x\}$

其中 $G_y$ 表示 $y$ 在 $G$ 中的稳定化子。

Voronoi单元分解

对于 $x \in \mathbb{R}^d$ ，定义：

Voronoi单元： $U_x := \{z \in \mathbb{R}^d : \{x\} = \arg\max_{p \in [x]} \langle p, z \rangle\}$
开Voronoi单元： $V_x := \text{relint}(U_x)$
开Voronoi图： $Q_x := \bigsqcup_{p \in [x]} V_p$

主要定理

定理4（局部双Lipschitz性）

设 $G \leq O(d)$ 是紧致群， $c := d - \max_{x \in \mathbb{R}^d} \dim([x])$ 。对于通用的 $z_1, \ldots, z_n \in \mathbb{R}^d$ ，当 $n > 2 \cdot \chi(G) \cdot (c-1)$ 时，max filter bank $\Phi$ 在每个 $x \in R(G)$ 处局部双Lipschitz。

定理5（全局双Lipschitz性）

设 $G \leq O(d)$ 是紧致群且 $\mathbb{R}^d - \{0\} \subseteq R(G)$ ， $c := d - \max_{x \in \mathbb{R}^d} \dim([x])$ 。对于通用的 $z_1, \ldots, z_n \in \mathbb{R}^d$ ，当 $n > 2 \cdot \chi(G) \cdot (c-1)$ 时，max filter bank $\Phi$ 是双Lipschitz的。

技术创新点

几何刻画方法：通过Voronoi分解提供了主要点和正则点的几何刻画
去奇ularization技术：为非流形轨道空间构造了局部流形结构
半代数几何分析：利用半代数集合的维数保持性质进行复杂度分析
Riemannian几何工具：结合测地线和割轨理论分析轨道空间的几何性质

实验设置

理论验证

论文主要是理论工作，通过以下方式验证结果：

具体例子分析：
- 三维旋转反射群的Voronoi分解
- 圆群在复空间上的酉表示
- 加权相位恢复的特殊情况
维数计算：
- 对于复相位恢复： $\chi(G) = 1$ , $c = 2d-1$
- 对于加权情况： $\chi(G) \leq k_{\max}$ , $c \leq p$

应用验证

冷冻电镜应用

问题规模： $L \times L$ 像素图像， $k_{\max} = O(L)$ ， $p = O(L^2)$
模板需求： $O(L^3)$ 个通用模板（相比双谱嵌入的 $O(L^5)$ 有显著改进）
理论保证：提供了双Lipschitz常数的显式界限

实验结果

主要结果

维数界限的精确性：
- 证明了"坏"模板集合的维数上界
- 建立了半代数集合的维数估计
Voronoi分解的完整性：
- 证明了 $U_x = V_x$ 当且仅当特定条件成立
- 提供了开Voronoi单元的完整刻画
应用效果：
- 冷冻电镜：从 $O(L^5)$ 降低到 $O(L^3)$ 的复杂度
- 加权相位恢复：提供了稳定性保证

理论发现

几何互换性：
- 主要点： $z \in V_x \Leftrightarrow x \in V_z$
- 正则点： $z \in V_x \Leftrightarrow x \in V_z^{\text{loc}}$
维数关系：
- 正则Voronoi复杂度与群结构的深层联系
- 半代数维数的保持性质