Provable and Robust Wavefront Sensing via Self-Reference Interferometry¶

会议: ECCV 2026
论文: ECCV 官方页面
领域: 医学图像
关键词: 波前传感、自参考干涉测量、相位恢复、互质位移、图传播算法

一句话总结¶

本文提出了一种基于空间平移自参考干涉测量的非迭代波前传感框架，通过证明连续互质位移构建具备紧凑跃点直径的连通图，并结合等跃点路径平均与频域最小二乘闭式精炼，在仅需 8 至 16 幅干涉测量下实现了兼具理论收敛保证与高抗噪性的全局相位恢复。

研究背景与动机¶

波前传感旨在恢复光场的复振幅分布（包含强度与相位信息），是自适应光学、定量相位显微、浑浊介质成像以及全息成像等领域的核心使能技术。然而，常规光电探测器仅能对光波的时间平均光强响应，导致波前相位信息完全丢失，使得非凸的相位恢复问题成为波前传感的核心瓶颈。现有计算相位恢复方案主要依赖交替投影（如 Gerchberg-Saxton 算法）或非凸梯度下降优化，这类方法高度依赖初始猜测、计算耗时漫长且容易陷入局部极小；而深度学习方法与即插即用去噪先验（PnP）虽提升了恢复质量，但缺乏可解释的理论界定与测量物理可溯源性。

传统相移干涉测量（PSI）虽能通过引入已知参考光场推导解析恢复公式，但外部参考光路极其脆弱，在实际复杂环境、活体显微与非受控场景中极易受到机械振动与空气湍流的扰动，使得干涉条纹对比度与相位步进严重失真。点衍射干涉（PDI）等共光路设计虽利用小孔空间滤波入射光自身产生参考光以抵抗环境干扰，但极小的针孔透光导致参考光能量极度衰减，引发探测信噪比（SNR）急剧下降。横向剪切干涉测量（LSI）通过波前自相干测量剪切差分，然而经典差分积分算法对噪声极为敏感，长距离累积误差会导致相位边缘严重畸变，且缺乏对剪切量选取的系统性理论指导。

本文的切入角度是：将入射波前与其空间位移副本之间的自相干干涉建模为离散像素图上的相量差分积分，通过数论与图论工具重新设计测量采样的位移向量。核心 idea：通过引入连续互质位移构建强连通像素图，严格证明其达成最小可能跃点距离以界定最坏情况下的误差累积，并提出单次非迭代的等跃点路径相量传播与闭式频域最小二乘精炼算法，实现兼具可证明稳健性与高吞吐的波前重建。

方法详解¶

整体框架¶

本文方法的输入为待测未知复光场 $x = |x| \odot e^{j\phi}$ 在空间光调制器（SLM）调控下与若干空间位移副本叠加后的四步正交相移干涉强度图 $y_{k,q}$，输出为精确重建的二维全局相位图 $\hat{\phi}_{LS}$。整个求解流程避免了传统非凸优化中的前向与伴随算子反复迭代，分为自参考相量差分提取、互质位移图构建、等跃点路径相量传播以及频域全局最小二乘精炼四个主要阶段。

%%{init: {'flowchart': {'rankSpacing': 24, 'nodeSpacing': 28, 'padding': 6, 'wrappingWidth': 400}}}%%
flowchart TD
    In["输入入射光场与移相干涉强度<br/>x, y_{k,q}"] --> D1["自参考相量差分提取<br/>四步相移正交抵消求取相对相量 p_k"]
    D1 --> D2["互质位移图构建<br/>连续互质位移 s, t 最小化图直径 h*"]
    D2 --> D3["等跃点路径相量传播<br/>BFS 遍历与同距离相量平均"]
    D3 --> D4["频域全局最小二乘精炼<br/>整数缠绕展开与 FFT 闭式去噪解"]
    D4 --> Out["输出高保真全波前相位<br/>phi_{LS}"]

关键设计¶

1. 自参考相量差分提取：规避外部参考光并解耦幅相调制 传统共光路点衍射方案通过极细针孔滤波生成单点参考光，导致参考通道能量微弱且信噪比受限。本文设计了全口径空间平移自干涉测量架构，使未知波前 $x$ 与经平移算子 $\mathcal{S}_{\Delta_k}$ 作用后的自身副本发生干涉。对于给定的二维位移向量 $\Delta_k$，系统引入四个正交步进相位 $\varphi_q = q\frac{\pi}{2}$ ($q \in \{0,1,2,3\}$)，捕获干涉强度： $$y_{k,q} = |x + e^{j\varphi_q}\mathcal{S}_{\Delta_k}(x)|^2 + \eta_{k,q} = |x|^2 + |\mathcal{S}_{\Delta_k}(x)|^2 + 2\text{Re}\left(|x| \odot |\mathcal{S}_{\Delta_k}(x)| e^{j(\phi - \mathcal{S}_{\Delta_k}(\phi) - \varphi_q)}\right) + \eta_{k,q}$$ 通过对四步正交测量施加相位加权求和，直流背景能量与自强分量被完全抵消，直接解耦出携带像素间相位差的复交叉项。进一步对提取出的复数场进行逐元素模长归一化，消除了空间非均匀光照幅度的调制影响，解析导出纯粹由两像素相位差决定的单位相量图 $p_k(i) \approx e^{j(\phi(i+\Delta_k) - \phi(i))}$。该过程完全不需要外部分束光路，从物理层面抵抗共模环境振动，同时保持全孔径采光效率。

2. 互质位移图构建：以数论性质保全图连通与跃点下界 在传统横向剪切或邻域一阶差分中，像素点仅通过单位位移 $\Delta=1$ 连接，图拓扑退化为简单链状图，远离参考点的边缘像素需要遍历多达 $O(N)$ 个中间节点，各边独立测量噪声方差随路径长度线性累加发散。本文将测量位移组建模为定义在离散网格 $V$ 上的图边集 $E = \{(i, i+s_k)\}$。针对一维 $N$ 节点（以及张量积扩展至二维）的连通性问题，论文从数论上证明了：只要两个位移量 $s, t \in \mathbb{N}$ 互质（即 $\gcd(s,t)=1$）且满足 $s+t \le N$，通过模 $(s+t)$ 剩余类全覆盖定理与滑动窗口引理，图上任意节点必可经有限条边跳转至全局基准参考点。更为关键的是，为了最小化噪声沿路径传播的扩散，论文确立了任意双位移图最大跃点距离的普适理论下界： $$h^\star \ge \left\lceil \frac{-1 + \sqrt{2N - 1}}{2} \right\rceil$$ 本文严格证明了选取连续互质整数对 $s = \lfloor \sqrt{N/2} \rfloor$ 与 $t = \lfloor \sqrt{N/2} \rfloor + 1$ 可以渐近达到这一理论下界。以 $N=512$ 分辨率为例，仅需设置位移量 $s=16, t=17$，即可将从中心参考点到全图任意像素的最大跃点距离由单步剪切的 256 次骤降至 16 次，从根本上锁死最坏情形下的误差累积上界。

3. 等跃点路径相量传播：抑制累积方差的非迭代波前遍历 在获得互质位移差分相量后，全图相位积分需在复数相量域完成以规避三角函数反求带来的缠绕伪影。本文提出基于广度优先搜索（BFS）的波前推进策略。算法将参考节点固定为零相位（即基准相量 $\hat{p}(i_0) = 1$），维护一个从近及远的波前队列。当多个不同方向或位移组合的多条路径同时到达同一候选节点时，若盲目进行全路径平均，较长路径所携带的额外累积噪声反而会污染估计精度。为此，本文设计了等跃点约束平均准则：仅当新路径与已到达路径具有完全相同的跃点计数 $h(i') = h(i) + 1$ 时，才对估计相量执行矢量平均；一旦出现跃点数更大的长路径则直接丢弃。这一机制在严格保持无偏性的同时最大化降低了局部估计方差，并在 GPU 上通过张量化波前扩张实现全图并行推进。

4. 频域全局最小二乘精炼：消除缠绕环路积分误差的闭式求解 广度优先相量传播虽然速度极快且有效避免了长链发散，但因其树状遍历特性未能充分利用闭合环路上的冗余差分约束。本文将离散差分算子表示为 $D_{\Delta_k} = \mathcal{S}_{\Delta_k} - I$。利用图传播初步得到的全图相位 $\hat{\phi}$，首先精确计算各差分测量值的整周期缠绕跳变数： $$m_k = \text{round}\left( \frac{D_{\Delta_k}\hat{\phi} - \nabla_k\phi}{2\pi} \right)$$ 从而将原始非线性缠绕差分 $\nabla_k\phi = \angle p_k \in (-\pi, \pi]$ 线性化展开为平滑差分 $\nabla_k\phi^\text{unwrap} = \nabla_k\phi + 2\pi m_k$。随后，构建联合拟合所有测量位移的能量泛函 $\min_\phi \sum_k \|D_{\Delta_k}\phi - \nabla_k\phi^\text{unwrap}\|_2^2$。由于前向差分算子在空间域等价于循环卷积核 $d_{\Delta_k}$，该最小二乘问题可在二维傅里叶变换域完美对角化，推导出全局解析闭式解： $$\hat{\phi}_{LS} = \mathcal{F}^{-1}\left( \frac{\sum_k \overline{\mathcal{F}(d_{\Delta_k})} \odot \mathcal{F}(\nabla_k\phi^\text{unwrap})}{\sum_k |\mathcal{F}(d_{\Delta_k})|^2 + \lambda} \right)$$ 式中 $\mathcal{F}$ 为二维离散傅里叶变换，上划线表示复共轭，$\lambda$ 为微小正则化系数。该精炼步骤仅需若干次二维快速傅里叶变换（FFT），便能在极短时间内将全图微弱残余噪声降至理论极小。

实验关键数据¶

主实验¶

在合成数据集上，研究使用 DIV2K 验证集的 100 幅图像作为真实光强分布，结合二次相位（Quadratic）、随机相位（Random）及平滑峰值相位（Smooth peaks）生成复振幅场，并施加泊松-高斯混合噪声模拟探测器在 22 dB SNR 条件下的退化。对比基线包括随机初值梯度下降（GD-Rand）、谱初始化梯度下降（GD-Spec）、结合预训练 DRUNet 的即插即用近端算法（PnP-FISTA）、基于空间相移改进的交替投影传感算法（WISH）以及未训练深度图像先验（DIP）。

算法	二次相位 (16测)	二次相位 (32测)	随机相位 (16测)	随机相位 (32测)	平滑相位 (16测)	平滑相位 (32测)
GD-Rand	0.695 ± 0.365	0.470 ± 0.434	0.138 ± 0.148	0.079 ± 0.152	0.791 ± 0.361	0.570 ± 0.404
GD-Spec	0.737 ± 0.400	0.400 ± 0.432	0.691 ± 0.367	0.449 ± 0.440	0.686 ± 0.397	0.371 ± 0.374
WISH	0.852 ± 0.299	0.423 ± 0.366	0.854 ± 0.320	0.453 ± 0.377	0.898 ± 0.330	0.425 ± 0.367
PnP-FISTA	0.669 ± 0.341	0.456 ± 0.304	0.735 ± 0.329	0.425 ± 0.264	0.663 ± 0.341	0.479 ± 0.338
DIP	1.486 ± 0.031	1.491 ± 0.051	1.495 ± 0.037	1.447 ± 0.065	1.442 ± 0.063	1.423 ± 0.081
本文方法 (Ours)	0.158 ± 0.104	0.099 ± 0.083	0.156 ± 0.099	0.094 ± 0.075	0.155 ± 0.096	0.097 ± 0.071
本文精炼 (Ours+LS)	0.183 ± 0.082	0.126 ± 0.065	0.158 ± 0.075	0.121 ± 0.055	0.205 ± 0.063	0.159 ± 0.050

消融实验¶

为验证互质位移理论与跃点距离优化对相位重建精度的直接影响，下表对比了在 $N=512$ 空间网格不同位移配置下的图拓扑性质与重构性能表现，以及不同重建算法在单张 NVIDIA RTX 6000 Ada GPU 上的单帧重构耗时与加速倍率。

位移配置 / 算法	最大跃点数 $h$	图拓扑连通性	累积误差趋势	单帧耗时 (s)	相对加速比
单步剪切 ($\Delta_1=1$)	256	线性链状图连通	随几何跨度线性恶化发散	-	-
非互质位移 (如 16, 18)	$\infty$	不连通 (分裂为 2 个不连通子图)	存在未覆盖盲区，完全失效	-	-
远距次优互质 (如 2, 3)	86	稀疏长跨度连通	局部环路长，跳步噪声放大	-	-
最优互质位移 (16, 17)	16	全连通紧致图 (达到理论下界 $h^\star$)	跃点高度集中，方差最小	-	-
PnP-FISTA 迭代重构	-	-	依赖迭代去噪循环	10.7	1.0× (基准)
DIP 神经网络拟合	-	-	需测试时反向传播优化权重	8.7	1.2×
WISH 交替投影算法	-	-	多轮前后向傅里叶投影	3.6	3.0×
梯度下降 (GD-Rand/Spec)	-	-	依赖非凸能量函数迭代下降	1.5	7.1×
本文方法 (Ours / +LS)	-	-	解析波前推进 + 闭式 FFT 最小二乘	1.2	8.9×

关键发现¶

在仅有 16 幅干涉测量（2 个最优位移 $\times$ 水平/垂直 2 个维度 $\times$ 4 步正交相移）的极端欠定条件下，所有传统优化与深度先验基线在平滑与二次连续波前上的重构误差均高于 0.6 rad，而本文方法直接压缩至 0.155 rad 左右，成为唯一能精确刻画波前细节与条纹拓扑的算法。
模拟信噪比扫描实验（SNR 跨越 10 dB 至 30 dB）显示，频域全局最小二乘精炼（Ours with LS）在低信噪比恶劣工况下提升尤为显著，有效熨平了纯图遍历可能偶发的微弱路径条纹。
硬件光学台搭建验证了透镜相位测量、纯计算数值自动重聚焦（Fresnel 数值衍射传播在 $z=67\text{ mm}$ 处清晰呈现星标测试靶锐利边缘）以及透过散射介质（Thorlabs DG10-120-A 磨砂玻璃扩散片）透过成像，成功在散斑光场中校正波前相位畸变并恢复隐藏目标结构。

亮点与洞察¶

将光学干涉工程巧妙映射为图论中的度量遍历问题，跳出了传统非凸反问题“前向建模-梯度下降”的思维定势，给出了首个兼具可证明连通性与误差界的最优位移解析解。
建立在连续互质数对基础上的测量设计，不仅在纯数学上实现了最少跳步覆盖，在光学物理上亦与空间光调制器（SLM）的相位斜率调制特性完美契合，使得纯干涉自参考方案无需引入精密针孔即可获得极高能量利用率。
提出的等跃点平均准则极具启发性：在信息融合中并非“路径越多越好”，由于噪声随步数单调递增，截断劣质长路径、仅对相同置信度层级的相量进行矢量合成，是保证全局估计无偏且方差最小的优雅策略。

局限与展望¶

作者指出当前算法对相同跃点的路径赋予了相等的权重，然而在实际光学系统中，各像素点的局部干涉条纹对比度与光照强度存在差异；引入基于质量引导（Quality-guided）或置信度加权的相量平均有望进一步提升极端斑图下的稳健性。
在纯高频随机离散相位分布下，由于相位梯度在像素间跳变剧烈且缺乏空间连续性支撑，单纯图传播方法相较于无先验随机梯度下降的优势有所收敛。
目前搭建的硬件原型机依赖空间光调制器串行加载四步相移与位移相位光栅，暂无法捕获纳秒级或微秒级超快动态波前畸变；未来与偏振分光阵列或单次曝光多波长全息相结合，有望走向实时活体显微与动态自适应光学闭环。

评分¶

新颖性: ⭐⭐⭐⭐⭐ 创造性地融合数论互质理论与图论最短路模型解决自参考波前传感难题，理论完备且优雅。
实验充分度: ⭐⭐⭐⭐⭐ 从严谨的数学界定、全流程仿真、多算法耗时基准，到实物搭建 4f 光路、透镜重建、散射介质穿透与自动重聚焦，链条闭环极强。
写作质量: ⭐⭐⭐⭐⭐ 理论推导清晰透彻，定理与引理陈述严谨，实验图表排布极具说服力。
价值: ⭐⭐⭐⭐⭐ 为无参考光环境下的高分辨波前重构提供了高速解析新范式，在生物活体显微与天文学自适应光学中有广阔落地前景。

位移配置 / 算法	最大跃点数 \(h\)	图拓扑连通性	累积误差趋势	单帧耗时 (s)	相对加速比
单步剪切 (\(\Delta_1=1\))	256	线性链状图连通	随几何跨度线性恶化发散	-	-
非互质位移 (如 16, 18)	\(\infty\)	不连通 (分裂为 2 个不连通子图)	存在未覆盖盲区，完全失效	-	-
远距次优互质 (如 2, 3)	86	稀疏长跨度连通	局部环路长，跳步噪声放大	-	-
最优互质位移 (16, 17)	16	全连通紧致图 (达到理论下界 \(h^\star\))	跃点高度集中，方差最小	-	-
PnP-FISTA 迭代重构	-	-	依赖迭代去噪循环	10.7	1.0× (基准)
DIP 神经网络拟合	-	-	需测试时反向传播优化权重	8.7	1.2×
WISH 交替投影算法	-	-	多轮前后向傅里叶投影	3.6	3.0×
梯度下降 (GD-Rand/Spec)	-	-	依赖非凸能量函数迭代下降	1.5	7.1×
本文方法 (Ours / +LS)	-	-	解析波前推进 + 闭式 FFT 最小二乘	1.2	8.9×